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

Kelley, L. A.

Normalized to: Kelley, L.

38 article(s) in total. 353 co-authors, from 1 to 15 common article(s). Median position in authors list is 17,0.

[1]  [pdf] - 2026596
Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays
Comments: Fixed typo in author list. Code that supports all calculations and figures is available at
Submitted: 2020-01-02, last modified: 2020-01-06
The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.
[2]  [pdf] - 2042263
The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory
Comments: 10 pages, 6 figures, submitted to ApJ
Submitted: 2019-11-19, last modified: 2019-12-06
The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.
[3]  [pdf] - 2025702
Probing Massive Black Hole Binary Populations with LISA
Comments: 18 pages, 9 figures, 3 tables
Submitted: 2019-08-15, last modified: 2019-11-01
ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.
[4]  [pdf] - 1981263
Forward Modeling of Double Neutron Stars: Insights from Highly-Offset Short Gamma-Ray Bursts
Comments: 15 pages, 7 figures, 1 table
Submitted: 2019-10-08, last modified: 2019-10-16
We perform a detailed analysis of two well-localized, highly-offset short gamma-ray bursts --- GRB070809 and GRB090515 --- investigating the prior kinematic evolution of their progenitors from compact object formation until merger. Calibrated to observations of their most probable galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of $\gtrsim\,250~\mathrm{km\,s}^{-1}$. For massive and quiescent hosts, we see $\sim\,20\%$ of mergers occur at $\gtrsim\,5$ effective radii from their hosts. Though the offset of GRB070809 is explainable with small supernova natal kicks, our modeling suggests that the second-born neutron star in the GRB090515 progenitor system received a natal kick of $\gtrsim\,200~\mathrm{km\,s}^{-1}$ at the $98\%$ credible level. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.
[5]  [pdf] - 2054030
The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics
Comments: 14 pages, 13 figures, fixed typo in abstract of earlier version
Submitted: 2019-09-18, last modified: 2019-09-20
An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.
[6]  [pdf] - 1969541
The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
Comments: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland (
Submitted: 2018-12-30, last modified: 2019-05-21
Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.
[7]  [pdf] - 1905635
The Astrophysics of Nanohertz Gravitational Waves
Comments: Accepted to The Astronomy and Astrophysics Review (A&ARv)
Submitted: 2018-11-15, last modified: 2019-05-17
Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.
[8]  [pdf] - 1862298
The IllustrisTNG Simulations: Public Data Release
Comments: The data is made available at (comments welcome)
Submitted: 2018-12-13, last modified: 2019-04-05
We present the full public release of all data from the TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~750 TB, including 1200 full volume snapshots and ~80,000 high time-resolution subbox snapshots. This will increase to ~1.1 PB with the future release of TNG50. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.
[9]  [pdf] - 1853139
Supermassive Black-hole Demographics & Environments With Pulsar Timing Arrays
Comments: Submitted to the Astro2020 Decadal Survey. One of 5 core white-papers authored by members of the NANOGrav Collaboration. 9 pages, 2 figures
Submitted: 2019-03-19
Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.
[10]  [pdf] - 1851819
Multi-Messenger Astrophysics with Pulsar Timing Arrays
Comments: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al
Submitted: 2019-03-18
Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.
[11]  [pdf] - 1846803
Massive BH Binaries as Periodically-Variable AGN
Comments: 18 pages, 13 figures, comments appreciated
Submitted: 2018-09-06, last modified: 2018-12-06
Massive black-hole (MBH) binaries, which are expected to form following the merger of their parent galaxies, produce gravitational waves which will be detectable by Pulsar Timing Arrays at nanohertz frequencies (year periods). While no confirmed, compact MBH binary systems have been seen in electromagnetic observations, a large number of candidates have recently been identified in optical surveys of AGN variability. Using a combination of cosmological, hydrodynamic simulations; comprehensive, semi-analytic binary merger models; and analytic AGN spectra and variability prescriptions; we calculate the expected electromagnetic detection rates of MBH binaries as periodically variable AGN. In particular, we consider two independent variability models: (i) Doppler boosting due to large orbital velocities, and (ii) hydrodynamic variability in which the fueling of MBH accretion disks is periodically modulated by the companion. Our models predict that numerous MBH binaries should be present and distinguishable in the existing data. In particular, our fiducial models produce an expectation value of $0.2$ (Doppler) and $5$ (hydrodynamic) binaries to be identifiable in CRTS, while $20$ and $100$ are expected after five years of LSST observations. The brightness variations in most systems are too small to be distinguishable, but almost $1\%$ of AGN at redshifts $z \lesssim 0.6$ could be in massive binaries. We analyze the predicted binary parameters of observable systems and their selection biases, and include an extensive discussion of our model parameters and uncertainties.
[12]  [pdf] - 1659501
Testing the binary hypothesis: pulsar timing constraints on supermassive black hole binary candidates
Comments: 14 pages, 11 figures, 3 tables. Resubmitted to the Astrophysical Journal after some major revision of the results including a proper estimate of the intrinsic mass of the binary candidates
Submitted: 2017-03-30, last modified: 2017-12-16
The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.
[13]  [pdf] - 1656164
Single Sources in the Low-Frequency Gravitational Wave Sky: properties and time to detection by pulsar timing arrays
Comments: 10 pages, 11 figures, comments welcome
Submitted: 2017-10-31
We calculate the properties, occurrence rates and detection prospects of individually resolvable 'single sources' in the low frequency gravitational wave (GW) spectrum. Our simulations use the population of galaxies and massive black hole binaries from the Illustris cosmological hydrodynamic simulations, coupled to comprehensive semi-analytic models of the binary merger process. Using mock pulsar timing arrays (PTA) with, for the first time, varying red-noise models, we calculate plausible detection prospects for GW single sources and the stochastic GW background (GWB). Contrary to previous results, we find that single sources are at least as detectable as the GW background. Using mock PTA, we find that these 'foreground' sources (also 'deterministic'/'continuous') are likely to be detected with $\sim 20 \, \textrm{yr}$ total observing baselines. Detection prospects, and indeed the overall properties of single sources, are only moderately sensitive to binary evolution parameters---namely eccentricity & environmental coupling, which can lead to differences of $\sim 5 \,\textrm{yr}$ in times to detection. Red noise has a stronger effect, roughly doubling the time to detection of the foreground between a white-noise only model ($\sim 10$ - $15 \, \textrm{yr}$) and severe red noise ($\sim 20$ - $30 \, \textrm{yr}$). The effect of red noise on the GWB is even stronger, suggesting that single source detections may be more robust. We find that typical signal-to-noise ratios for the foreground peak near $f = 0.1 \, \textrm{yr}^{-1}$, and are much less sensitive to the continued addition of new pulsars to PTA.
[14]  [pdf] - 1581420
The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays
Comments: 21 pages, 1 table, 18 figures; comments welcomed. v2: updated figures and calculations (times-to-detections changed), added author S.R.Taylor
Submitted: 2017-02-07, last modified: 2017-05-31
Pulsar Timing Arrays (PTA) around the world are using the incredible consistency of millisecond pulsars to measure low frequency gravitational waves from (super)Massive Black Hole (MBH) binaries. We use comprehensive MBH merger models based on cosmological hydrodynamic simulations to predict the spectrum of the stochastic Gravitational-Wave Background (GWB). We use real Time-of-Arrival (TOA) specifications from the European, NANOGrav, Parkes, and International PTA (IPTA) to calculate realistic times to detection of the GWB across a wide range of model parameters. In addition to exploring the parameter space of environmental hardening processes (in particular: stellar scattering efficiencies), we have expanded our models to include eccentric binary evolution which can have a strong effect on the GWB spectrum. Our models show that strong stellar scattering and high characteristic eccentricities enhance the GWB strain amplitude near the PTA sensitive "sweet-spot" (near the frequency $f = 1 \, \mathrm{yr}^{-1}$), slightly improving detection prospects in these cases. While the GWB $amplitude$ is degenerate between cosmological and environmental parameters, the location of a spectral turnover at low frequencies ($f \lesssim 0.1 \, \mathrm{yr}^{-1}$) is strongly indicative of environmental coupling. At high frequencies ($f\gtrsim 1 \, \mathrm{yr}^{-1}$), the GWB spectral index can be used to infer the number density of sources and possibly their eccentricity distribution. Even with merger models that use pessimistic environmental and eccentricity parameters, if the current rate of PTA expansion continues, we find that the International PTA is highly likely to make a detection within about 10 years.
[15]  [pdf] - 1530658
An Open Catalog for Supernova Data
Comments: 26 pages, 6 figures, 1 table. Updated in response to referee comments, added Luke Zoltan Kelley to author list. Catalog available at
Submitted: 2016-05-03, last modified: 2016-11-07
We present the Open Supernova Catalog, an online collection of observations and metadata for presently 36,000+ supernovae and related candidates. The catalog is freely available on the web (, with its main interface having been designed to be a user-friendly, rapidly-searchable table accessible on desktop and mobile devices. In addition to the primary catalog table containing supernova metadata, an individual page is generated for each supernova which displays its available metadata, light curves, and spectra spanning X-ray to radio frequencies. The data presented in the catalog is automatically rebuilt on a daily basis and is constructed by parsing several dozen sources, including the data presented in the supernova literature and from secondary sources such as other web-based catalogs. Individual supernova data is stored in the hierarchical, human- and machine-readable JSON format, with the entirety of each supernova's data being contained within a single JSON file bearing its name. The setup we present here, which is based upon open source software maintained via git repositories hosted on github, enables anyone to download the entirety of the supernova dataset to their home computer in minutes, and to make contributions of their own data back to the catalog via git. As the supernova dataset continues to grow, especially in the upcoming era of all-sky synoptic telescopes which will increase the total number of events by orders of magnitude, we hope that the catalog we have designed will be a valuable tool for the community to analyze both historical and contemporary supernovae.
[16]  [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
Submitted: 2016-10-10
The Astropy Project ( 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.
[17]  [pdf] - 1492245
Massive Black Hole Binary Mergers in Dynamical Galactic Environments
Comments: 28 pages, 28 figures, 3 tables, Monthly Notices of the Royal Astronomical Society 2016
Submitted: 2016-06-06, last modified: 2016-10-05
Gravitational Waves (GW) have now been detected from stellar-mass black hole binaries, and the first observations of GW from Massive Black Hole (MBH) Binaries are expected within the next decade. Pulsar Timing Arrays (PTA), which can measure the years long periods of GW from MBHB, have excluded many standard predictions for the amplitude of a stochastic GW Background (GWB). We use coevolved populations of MBH and galaxies from hydrodynamic, cosmological simulations ('Illustris') to calculate a predicted GWB. The most advanced predictions so far have included binary hardening mechanisms from individual environmental processes. We present the first calculation including all of the environmental mechanisms expected to be involved: dynamical friction, stellar 'loss-cone' scattering, and viscous drag from a circumbinary disk. We find that MBH binary lifetimes are generally multiple gigayears, and only a fraction coalesce by redshift zero. For a variety of parameters, we find all GWB amplitudes to be below the most stringent PTA upper limit of $A_{\textrm{yr}^{-1}} \approx 10^{-15}$. Our fairly conservative fiducial model predicts an amplitude of $A_{\textrm{yr}^{-1}} \approx 0.4\times 10^{-15}$---less than a factor of three below the current limit. At lower frequencies, we find $A_{0.1\,\textrm{yr}^{-1}} \approx 1.5\times 10^{-15}$ with spectral indices between $-0.4$ and $-0.6$---significantly flatter than the canonical value of $-2/3$ due to purely GW-driven evolution. Typical MBHB driving the GWB signal come from redshifts around $0.3$, with total masses of a few times $10^9\,M_\odot$, and in host galaxies with very large stellar masses. Even without GWB detections, our results can be connected to observations of dual AGN to constrain binary evolution.
[18]  [pdf] - 1338938
Recoiling black holes: prospects for detection and implications of spin alignment
Comments: 30 pages, 19 figures. Accepted to MNRAS after minor revisions
Submitted: 2015-08-06, last modified: 2015-12-10
Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.
[19]  [pdf] - 1222352
Tidal Disruption and Magnetic Flux Capture: Powering a Jet from a Quiescent Black Hole
Comments: MNRAS Accepted, September 30th, 2014. 20 pages, 14 figures. Animations online:
Submitted: 2014-10-01
The transient Swift J1644+57 is believed to have been produced by an unlucky star wandering too close to a supermassive black hole (BH) leading to a tidal disruption event. This unusual flare displayed highly super-Eddington X-ray emission which likely originated in a relativistic, collimated jet. This presents challenges to modern accretion and jet theory as upper limits of prior BH activity, which we obtain from the radio afterglow of this event, imply that both the pre-disruption BH and stellar magnetic fluxes fall many orders of magnitude short of what is required to power the observed X-ray luminosity. We argue that a pre-existing, "fossil" accretion disc can contain a sufficient reservoir of magnetic flux and that the stellar debris stream is capable of dragging this flux into the BH. To demonstrate this, we perform local, 3D magnetohydrodynamic simulations of the disc--stream interaction and demonstrate that the interface between the two is unstable to mixing. This mixing entrains a sufficient amount of fossil disc magnetic flux into the infalling stellar debris to power the jet. We argue that the interaction with the fossil disc can have a pronounced effect on the structure and dynamics of mass fallback and likely the resulting transient. Finally, we describe possible ramifications of these interactions on unresolved problems in tidal disruption dynamics, in particular, the efficiency of debris circularization, and effects of the disruption on the preexisting black hole system. Animations online:
[20]  [pdf] - 685325
Electromagnetic transients as triggers in searches for gravitational waves from compact binary mergers
Comments: 17 pages, 2 figures, 9 tables. Matches content of published PRD article
Submitted: 2012-09-13, last modified: 2013-06-26
The detection of an electromagnetic transient which may originate from a binary neutron star merger can increase the probability that a given segment of data from the LIGO-Virgo ground-based gravitational-wave detector network contains a signal from a binary coalescence. Additional information contained in the electromagnetic signal, such as the sky location or distance to the source, can help rule out false alarms, and thus lower the necessary threshold for a detection. Here, we develop a framework for determining how much sensitivity is added to a gravitational-wave search by triggering on an electromagnetic transient. We apply this framework to a variety of relevant electromagnetic transients, from short GRBs to signatures of r-process heating to optical and radio orphan afterglows. We compute the expected rates of multi-messenger observations in the Advanced detector era, and find that searches triggered on short GRBs --- with current high-energy instruments, such as Fermi --- and nucleosynthetic `kilonovae' --- with future optical surveys, like LSST --- can boost the number of multi-messenger detections by 15% and 40%, respectively, for a binary neutron star progenitor model. Short GRB triggers offer precise merger timing, but suffer from detection rates decreased by beaming and the high a priori probability that the source is outside the LIGO-Virgo sensitive volume. Isotropic kilonovae, on the other hand, could be commonly observed within the LIGO-Virgo sensitive volume with an instrument roughly an order of magnitude more sensitive than current optical surveys. We propose that the most productive strategy for making multi-messenger gravitational-wave observations is using triggers from future deep, optical all-sky surveys, with characteristics comparable to LSST, which could make as many as ten such coincident observations a year.
[21]  [pdf] - 1158931
Swift J1644+57 gone MAD: the case for dynamically-important magnetic flux threading the black hole in a jetted tidal disruption event
Comments: Submitted to MNRAS. 15 pages, 9 figures, uses mn2e.cls. Adjusted margins for clarity when printed
Submitted: 2013-01-09, last modified: 2013-01-10
The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]
[22]  [pdf] - 1091306
Towards improving the prospects for coordinated gravitational-wave and electromagnetic observations
Comments: Full text for proceedings of New Horizons in Time-domain Astronomy, IAU Symposium 285 [published proceedings will be abbreviated]
Submitted: 2011-10-31
We discuss two approaches to searches for gravitational-wave (GW) and electromagnetic (EM) counterparts of binary neutron star mergers. The first approach relies on triggering archival searches of GW detector data based on detections of EM transients. We introduce a quantitative approach to evaluate the improvement to GW detector reach due to the extra information gained from the EM transient and the increased confidence in the presence of a signal from a binary merger. We also advocate utilizing other transients in addition to short gamma ray bursts. The second approach involves following up GW candidates with targeted EM observations. We argue for the use of slower but optimal parameter-estimation techniques to localize the source on the sky, and for a more sophisticated use of astrophysical prior information, including galaxy catalogs, to find preferred followup locations.
[23]  [pdf] - 1041760
The Distribution of Coalescing Compact Binaries in the Local Universe: Prospects for Gravitational-Wave Observations
Comments: 9 pages, 3 figures; Accepted for publication in The Astrophysical Journal Letters
Submitted: 2010-11-04
Merging compact binaries are the most viable and best studied candidates for gravitational wave (GW) detection by the fully operational network of ground-based observatories. In anticipation of the first detections, the expected distribution of GW sources in the local universe is of considerable interest. Here we investigate the full phase space distribution of coalescing compact binaries at $z = 0$ using dark matter simulations of structure formation. The fact that these binary systems acquire large barycentric velocities at birth ("kicks") results in merger site distributions that are more diffusely distributed with respect to their putative hosts, with mergers occurring out to distances of a few Mpc from the host halo. Redshift estimates based solely on the nearest galaxy in projection can, as a result, be inaccurate. On the other hand, large offsets from the host galaxy could aid the detection of faint optical counterparts and should be considered when designing strategies for follow-up observations. The degree of isotropy in the projected sky distributions of GW sources is found to be augmented with increasing kick velocity and to be severely enhanced if progenitor systems possess large kicks as inferred from the known population of pulsars and double compact binaries. Even in the absence of observed electromagnetic counterparts, the differences in sky distributions of binaries produced by disparate kick-velocity models could be discerned by GW observatories, within the expected accuracies and detection rates of advanced LIGO--in particular with the addition of more interferometers.
[24]  [pdf] - 87077
Discovery of TeV Gamma-Ray Emission from the Cygnus Region of the Galaxy
Comments: Submitted to Astrophysical Journal Letters
Submitted: 2006-11-21
The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.
[25]  [pdf] - 71652
Constraints on Very High Energy gamma-ray emission from Gamma-Ray Bursts
Comments: submitted to ApJ
Submitted: 2005-03-11
The Milagro gamma-ray observatory employs a water Cherenkov detector to observe extensive air showers produced by high energy particles interacting in the Earth's atmosphere. Milagro has a wide field of view and high duty cycle, monitoring the northern sky almost continuously in the 100 GeV to 100 TeV energy range. Milagro is, thus, uniquely capable of searching for very high-energy emission from gamma-ray bursts (GRBs) during the prompt emission phase. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Twenty-five satellite-triggered GRBs occurred within the field of view of Milagro between January 2000 and December 2001. We have searched for counterparts to these GRBs and found no significant emission from any of the burst positions. Due to the absorption of high-energy gamma rays by the extragalactic background light, detections are only expected to be possible for redshifts less than ~0.5. Three of the GRBs studied have measured redshifts. GRB 010921 has a redshift low enough (0.45) to allow an upper limit on the fluence to place an observational constraint on potential GRB models.
[26]  [pdf] - 60969
Limits on Very High Energy Emission from Gamma-Ray Bursts with the Milagro Observatory
Comments: Accepted for publication in ApJL, minor revisions to the text, retypeset in emulateapj style
Submitted: 2003-11-17, last modified: 2004-02-05
The Milagro telescope monitors the northern sky for 100 GeV to 100 TeV transient emission through continuous very high energy wide-field observations. The large effective area and ~100 GeV energy threshold of Milagro allow it to detect very high energy (VHE) gamma-ray burst emission with much higher sensitivity than previous instruments and a fluence sensitivity at VHE energies comparable to that of dedicated gamma-ray burst satellites at keV to MeV energies. Even in the absence of a positive detection, VHE observations can place important constraints on gamma-ray burst (GRB) progenitor and emission models. We present limits on the VHE flux of 40 s -- 3 h duration transients nearby to earth, as well as sensitivity distributions which have been corrected for gamma-ray absorption by extragalactic background light and cosmological effects. The sensitivity distributions suggest that the typical intrinsic VHE fluence of GRBs is similar or weaker than the keV -- MeV emission, and we demonstrate how these sensitivity distributions may be used to place observational constraints on the absolute VHE luminosity of gamma-ray bursts for any GRB emission and progenitor model.
[27]  [pdf] - 50292
The high-energy gamma-ray fluence and energy spectrum of GRB 970417a from observations with Milagrito
Comments: Submitted to the Astrophysical Journal
Submitted: 2002-07-06
Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.
[28]  [pdf] - 33908
Evidence for TeV Emission from GRB 970417a
Comments: 10 pages, 3 figures
Submitted: 2000-01-07
Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.
[29]  [pdf] - 109103
First Results of a Study of TeV Emission from GRBs in Milagrito
Comments: 6 pages, 5 figures, to appear in "GeV-TeV Astrophysics: Toward a Major Atmospheric Cherenkov Telescope VI," Snowbird, Utah (August, 1999)
Submitted: 1999-10-29
Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.
[30]  [pdf] - 107439
Instrumented Water Tanks can Improve Air Shower Detector Sensitivity
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-07-15
Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.
[31]  [pdf] - 107138
Calibration of the Milagro Cosmic Ray Telescope
Comments: 4 pages, submitted to the XXVI International Cosmic Ray Conference
Submitted: 1999-06-25
The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.
[32]  [pdf] - 107105
Milagrito Detection of TeV Emission from Mrk 501
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.
[33]  [pdf] - 107108
Detection of 6 November 1997 Ground Level Event by Milagrito
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.
[34]  [pdf] - 107107
Search for Short Duration Bursts of TeV Gamma Rays with the Milagrito Telescope
Comments: submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.
[35]  [pdf] - 107106
Search for a TeV Component of GRBs using the Milagrito Detector
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.
[36]  [pdf] - 107109
Study of the Shadows of the Moon and the Sun with VHE Cosmic Rays
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.
[37]  [pdf] - 107104
Status of the Milagro Gamma Ray Observatory
Comments: 4 pages, submitted to XXVI International Cosmic Ray Conference, Salt Lake City
Submitted: 1999-06-24
The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.
[38]  [pdf] - 107090
An All-Sky Search for Steady VHE Gamma-Ray Sources
Comments: 4 pages
Submitted: 1999-06-23
The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.