-
The Solar System Notification Alert Processing System (SNAPS): Asteroid Population Outlier Detection
Authors:
Michael Gowanlock,
David E. Trilling,
Daniel Kramer,
Maria Chernyavskaya,
Andrew McNeill
Abstract:
The Solar System Notification Alert Processing System (SNAPS) is a ZTF and Rubin Observatory alert broker that will send alerts to the community regarding interesting events in the Solar System. SNAPS is actively monitoring Solar System objects and one of its functions is to compare objects (primarily main belt asteroids) to one another to find those that are outliers relative to the population. I…
▽ More
The Solar System Notification Alert Processing System (SNAPS) is a ZTF and Rubin Observatory alert broker that will send alerts to the community regarding interesting events in the Solar System. SNAPS is actively monitoring Solar System objects and one of its functions is to compare objects (primarily main belt asteroids) to one another to find those that are outliers relative to the population. In this paper, we use the SNAPShot1 dataset which contains 31,693 objects from ZTF and derive outlier scores for each of these objects. SNAPS employs an unsupervised approach; consequently, to derive outlier rankings for each object, we propose four different outlier metrics such that we can explore variants of outlier scores and add confidence to outlier rankings. We also provide outlier scores for each object in each permutation of 15 feature spaces, between 2 and 15 features, which yields 32,752 total feature spaces. We show that we can derive population outlier rankings each month at Rubin Observatory scale using four Nvidia A100 GPUs, and present several avenues of scientific investigation that can be explored using population outlier detection.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Removing Aliases in Time-Series Photometry
Authors:
Daniel Kramer,
Michael Gowanlock,
David Trilling,
Andrew McNeill,
Nicolas Erasmus
Abstract:
Ground-based, all-sky astronomical surveys are imposed with an inevitable day-night cadence that can introduce aliases in period-finding methods. We examined four different methods -- three from the literature and a new one that we developed -- that remove aliases to improve the accuracy of period-finding algorithms. We investigate the effectiveness of these methods in decreasing the fraction of a…
▽ More
Ground-based, all-sky astronomical surveys are imposed with an inevitable day-night cadence that can introduce aliases in period-finding methods. We examined four different methods -- three from the literature and a new one that we developed -- that remove aliases to improve the accuracy of period-finding algorithms. We investigate the effectiveness of these methods in decreasing the fraction of aliased period solutions by applying them to the Zwicky Transient Facility (ZTF) and the LSST Solar System Products Data Base (SSPDB) asteroid datasets. We find that the VanderPlas method had the worst accuracy for each survey. The mask and our newly proposed window method yields the highest accuracy when averaged across both datasets. However, the Monte Carlo method had the highest accuracy for the ZTF dataset, while for SSPDB, it had lower accuracy than the baseline where none of these methods are applied. Where possible, detailed de-aliasing studies should be carried out for every survey with a unique cadence.
△ Less
Submitted 26 April, 2023;
originally announced April 2023.
-
The Solar System Notification Alert Processing System (SNAPS): Design, Architecture, and First Data Release (SNAPShot1)
Authors:
David E. Trilling,
Michael Gowanlock,
Daniel Kramer,
Andrew McNeill,
Brian Donnelly,
Nat Butler,
John Kececioglu
Abstract:
We present here the design, architecture, and first data release for the Solar System Notification Alert Processing System (SNAPS). SNAPS is a Solar System broker that ingests alert data from all-sky surveys. At present, we ingest data from the Zwicky Transient Facility (ZTF) public survey, and we will ingest data from the forthcoming Legacy Survey of Space and Time (LSST) when it comes online. SN…
▽ More
We present here the design, architecture, and first data release for the Solar System Notification Alert Processing System (SNAPS). SNAPS is a Solar System broker that ingests alert data from all-sky surveys. At present, we ingest data from the Zwicky Transient Facility (ZTF) public survey, and we will ingest data from the forthcoming Legacy Survey of Space and Time (LSST) when it comes online. SNAPS is an official LSST downstream broker. In this paper we present the SNAPS design goals and requirements. We describe the details of our automatic pipeline processing in which physical properties of asteroids are derived. We present SNAPShot1, our first data release, which contains 5,458,459 observations of 31,693 asteroids observed by ZTF from July, 2018, through May, 2020. By comparing a number of derived properties for this ensemble to previously published results for overlapping objects we show that our automatic processing is highly reliable. We present a short list of science results, among many that will be enabled by our SNAPS catalog: (1) we demonstrate that there are no known asteroids with very short periods and high amplitudes, which clearly indicates that in general asteroids in the size range 0.3--20 km are strengthless; (2) we find no difference in the period distributions of Jupiter Trojan asteroids, implying that the L4 and L5 cloud have different shape distributions; and (3) we highlight several individual asteroids of interest. Finally, we describe future work for SNAPS and our ability to operate at LSST scale.
△ Less
Submitted 2 February, 2023;
originally announced February 2023.
-
GPU-Enabled Searches for Periodic Signals of Unknown Shape
Authors:
Michael Gowanlock,
Nathaniel R. Butler,
David E. Trilling,
Andrew McNeill
Abstract:
Recent and future generation observatories will enable the study of variable astronomical phenomena through their time-domain capabilities. High temporal fidelity will allow for unprecedented investigations into the nature of variable objects -- those objects that vary in brightness over time. A major bottleneck in data processing pipelines is constructing light curve solutions for catalogs of var…
▽ More
Recent and future generation observatories will enable the study of variable astronomical phenomena through their time-domain capabilities. High temporal fidelity will allow for unprecedented investigations into the nature of variable objects -- those objects that vary in brightness over time. A major bottleneck in data processing pipelines is constructing light curve solutions for catalogs of variable objects, as it is well-known that period finding algorithms are computationally expensive. Furthermore, there are many period finding algorithms that are often suited for specific science cases. In this paper, we present the first GPU-accelerated Super Smoother algorithm. Super Smoother is general purpose and uses cross-validation to fit line segments to a time series, and as such, is more computationally expensive than other algorithms, such as Lomb-Scargle. Because the algorithm requires making several scans over the input time series for a tested frequency, we also propose a novel generalized-validation variant of Super Smoother that only requires a single scan over the data. We compare the performance of our algorithms to analogous parallel multi-core CPU implementations on three catalogs of data, and show that it is generally advantageous to use the GPU algorithm over the CPU counterparts. Furthermore, we demonstrate that our single-pass variant of Super Smoother is roughly equally as accurate at finding correct period solutions as the original algorithm. Our software supports several features, such as batching the computation to eliminate the possibility of exceeding global memory on the GPU, processing a single object or batches of objects, and we allow for scaling the algorithm across multiple GPUs.
△ Less
Submitted 14 November, 2021;
originally announced November 2021.
-
Fast Period Searches Using the Lomb-Scargle Algorithm on Graphics Processing Units for Large Datasets and Real-Time Applications
Authors:
Michael Gowanlock,
Daniel Kramer,
David E. Trilling,
Nathaniel R. Butler,
Brian Donnelly
Abstract:
Computing the periods of variable objects is well-known to be computationally expensive. Modern astronomical catalogs contain a significant number of observed objects. Therefore, even if the period ranges for particular classes of objects are well-constrained due to expected physical properties, periods must be derived for a tremendous number of objects. In this paper, we propose a GPU-accelerated…
▽ More
Computing the periods of variable objects is well-known to be computationally expensive. Modern astronomical catalogs contain a significant number of observed objects. Therefore, even if the period ranges for particular classes of objects are well-constrained due to expected physical properties, periods must be derived for a tremendous number of objects. In this paper, we propose a GPU-accelerated Lomb-Scargle period finding algorithm that computes periods for single objects or for batches of objects as is necessary in many data processing pipelines. We demonstrate the performance of several optimizations, including comparing the use of shared and global memory GPU kernels and using multiple CUDA streams to copy periodogram data from the GPU to the host. Also, we quantify the difference between 32-bit and 64-bit floating point precision on two classes of GPUs, and show that the performance degradation of using 64-bit over 32-bit is greater on the CPU than a GPU designed for scientific computing. We find that the GPU algorithm achieves superior performance over the baseline parallel CPU implementation, achieving a speedup of up to 174.53$\times$. The Vera C. Rubin Observatory will carry out the Legacy Survey of Space and Time (LSST). We perform an analysis that shows we can derive the rotation periods of batches of Solar System objects at LSST scale in near real-time, which will be employed in a future LSST event broker. All source code has been made publicly available.
△ Less
Submitted 9 May, 2021;
originally announced May 2021.
-
The Breakthrough Listen Search For Intelligent Life Near the Galactic Center I
Authors:
Vishal Gajjar,
Karen I. Perez,
Andrew P. V. Siemion,
Griffin Foster,
Bryan Brzycki,
Shami Chatterjee,
Yuhong Chen,
James M. Cordes,
Steve Croft,
Daniel Czech,
David DeBoer,
Julia DeMarines,
Jamie Drew,
Michael Gowanlock,
Howard Isaacson,
Brian C. Lacki,
Matt Lebofsky,
David H. E. MacMahon,
Ian S. Morrison,
Cherry Ng,
Imke de Pater,
Danny C. Price,
Sofia Z. Sheikh,
Akshay Suresh,
Claire Webb
, et al. (1 additional authors not shown)
Abstract:
A line-of-sight towards the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys towards the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hours of deep observat…
▽ More
A line-of-sight towards the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys towards the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hours of deep observations across 0.7--93 GHz. We report preliminary results from our survey for ETI beacons across 1--8 GHz with 7.0 and 11.2 hours of observations with Parkes and GBT, respectively. With our narrowband drifting signal search, we were able to place meaningful constraints on ETI transmitters across 1--4 GHz and 3.9--8 GHz with EIRP limits of $\geq$4$\times$10$^{18}$ W among 60 million stars and $\geq$5$\times$10$^{17}$ W among half a million stars, respectively. For the first time, we were able to constrain the existence of artificially dispersed transient signals across 3.9--8 GHz with EIRP $\geq$1$\times$10$^{14}$ W/Hz with a repetition period $\leq$4.3 hours. We also searched our 11.2 hours of deep observations of the GC and its surrounding region for Fast Radio Burst-like magnetars with the DM up to 5000 pc cm$^{-3}$ with maximum pulse widths up to 90 ms at 6 GHz. We detected several hundred transient bursts from SGR J1745$-$2900, but did not detect any new transient burst with the peak luminosity limit across our observed band of $\geq$10$^{31}$ erg s$^{-1}$ and burst-rate of $\geq$0.23 burst-hr$^{-1}$. These limits are comparable to bright transient emission seen from other Galactic radio-loud magnetars, constraining their presence at the GC.
△ Less
Submitted 29 April, 2021;
originally announced April 2021.
-
The Habitability of our Evolving Galaxy
Authors:
Michael G. Gowanlock,
Ian S. Morrison
Abstract:
The notion of a Galactic Habitable Zone (GHZ), or regions of the Milky Way galaxy that preferentially maintain the conditions to sustain complex life, has recently gained attention due to the detection of numerous exoplanets and advances made in understanding habitability on the Earth and other environments. We discuss what a habitable environment means on large spatial and temporal scales, which…
▽ More
The notion of a Galactic Habitable Zone (GHZ), or regions of the Milky Way galaxy that preferentially maintain the conditions to sustain complex life, has recently gained attention due to the detection of numerous exoplanets and advances made in understanding habitability on the Earth and other environments. We discuss what a habitable environment means on large spatial and temporal scales, which necessarily requires an approximated definition of habitability to make an assessment of the astrophysical conditions that may sustain complex life. We discuss a few key exoplanet findings that directly relate to estimating the distribution of Earth-size planets in the Milky Way. With a broad notion of habitability defined and major observable properties of the Milky Way described, we compare selected literature on the GHZ and postulate why the models yield differing predictions of the most habitable regions at the present day, which include: (1) the majority of the galactic disk; (2) an annular ring between 7-9 kpc, and (3) the galactic outskirts. We briefly discuss the habitability of other galaxies as influenced by these studies. We note that the dangers to biospheres in the Galaxy taken into account in these studies may be incomplete and we discuss the possible role of Gamma-Ray Bursts and other dangers to life in the Milky Way. We speculate how changing astrophysical properties may affect the GHZ over time, including before the Earth formed, and describe how new observations and other related research may fit into the bigger picture of the habitability of the Galaxy.
△ Less
Submitted 20 February, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
Authors:
Kazunori Akiyama,
Shiro Ikeda,
Mollie Pleau,
Vincent L. Fish,
Fumie Tazaki,
Kazuki Kuramochi,
Avery Broderick,
Jason Dexter,
Monika Mościbrodzka,
Michael Gowanlock,
Mareki Honma,
Sheperd S. Doeleman
Abstract:
We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distributio…
▽ More
We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $λ/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.
△ Less
Submitted 3 March, 2017; v1 submitted 1 February, 2017;
originally announced February 2017.
-
Astrobiological Effects of Gamma-Ray Bursts in the Milky Way Galaxy
Authors:
Michael G. Gowanlock
Abstract:
A planet having protective ozone within the collimated beam of a Gamma Ray Burst (GRB) may suffer ozone depletion, potentially causing a mass extinction event to existing life on a planet's surface and oceans. We model the dangers of long GRBs to planets in the Milky Way and utilize a static statistical model of the Galaxy that matches major observable properties, such as the inside-out star forma…
▽ More
A planet having protective ozone within the collimated beam of a Gamma Ray Burst (GRB) may suffer ozone depletion, potentially causing a mass extinction event to existing life on a planet's surface and oceans. We model the dangers of long GRBs to planets in the Milky Way and utilize a static statistical model of the Galaxy that matches major observable properties, such as the inside-out star formation history, metallicity evolution, and 3-dimensional stellar number density distribution. The GRB formation rate is a function of both the star formation history and metallicity; however, the extent to which chemical evolution reduces the GRB rate over time in the Milky Way is still an open question. Therefore, we compare the damaging effects of GRBs to biospheres in the Milky Way using two models. One model generates GRBs as a function of the inside-out star formation history. The other model follows the star formation history, but generates GRB progenitors as a function of metallicity, thereby favoring metal-poor host regions of the Galaxy over time. If the GRB rate only follows the star formation history, the majority of the GRBs occur in the inner Galaxy. However, if GRB progenitors are constrained to low metallicity environments, then GRBs only form in the metal-poor outskirts at recent epochs. Interestingly, over the past 1 Gyr, the surface density of stars (and their corresponding planets) that survive a GRB is still greatest in the inner galaxy in both models. The present day danger of long GRBs to life at the solar radius ($R_\odot=8$ kpc) is low. We find that at least $\sim$65% of stars survive a GRB over the past 1 Gyr. Furthermore, when the GRB rate was expected to have been enhanced at higher redshifts, such as $z\gtrsim0.5$, our results suggest that a large fraction of planets would have survived these lethal GRB events.
△ Less
Submitted 29 September, 2016;
originally announced September 2016.
-
A Model of Habitability Within the Milky Way Galaxy
Authors:
Michael G. Gowanlock,
David R. Patton,
Sabine M. McConnell
Abstract:
We present a model of the Galactic Habitable Zone (GHZ), described in terms of the spatial and temporal dimensions of the Galaxy that may favour the development of complex life. The Milky Way galaxy is modelled using a computational approach by populating stars and their planetary systems on an individual basis using Monte-Carlo methods. We begin with well-established properties of the disk of the…
▽ More
We present a model of the Galactic Habitable Zone (GHZ), described in terms of the spatial and temporal dimensions of the Galaxy that may favour the development of complex life. The Milky Way galaxy is modelled using a computational approach by populating stars and their planetary systems on an individual basis using Monte-Carlo methods. We begin with well-established properties of the disk of the Milky Way, such as the stellar number density distribution, the initial mass function, the star formation history, and the metallicity gradient as a function of radial position and time. We vary some of these properties, creating four models to test the sensitivity of our assumptions. To assess habitability on the Galactic scale, we model supernova rates, planet formation, and the time required for complex life to evolve. Our study improves on other literature on the GHZ by populating stars on an individual basis and by modelling SNII and SNIa sterilizations by selecting their progenitors from within this preexisting stellar population. Furthermore, we consider habitability on tidally locked and non-tidally locked planets separately, and study habitability as a function of height above and below the Galactic midplane. In the model that most accurately reproduces the properties of the Galaxy, the results indicate that an individual SNIa is ~5.6 \times more lethal than an individual SNII on average. In addition, we predict that ~1.2% of all stars host a planet that may have been capable of supporting complex life at some point in the history of the Galaxy. Of those stars with a habitable planet, ~75% of planets are predicted to be in a tidally locked configuration with their host star. The majority of these planets that may support complex life are found towards the inner Galaxy, distributed within, and significantly above and below, the Galactic midplane.
△ Less
Submitted 6 July, 2011;
originally announced July 2011.