Source author record

Ashish A. Mahabal

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

27 published item(s)

preprint2022arXiv

A 62-minute orbital period black widow binary in a wide hierarchical triple

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution.

preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

preprint2022arXiv

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

preprint2022arXiv

Removing Atmospheric Fringes from Zwicky Transient Facility i-Band Images using Principal Component Analysis

The Zwicky Transient Facility is a time-domain optical survey that has substantially increased our ability to observe and construct massive catalogs of astronomical objects by use of its 47 square degree camera that can observe in multiple filters. However the telescope's i-band filter suffers from significant atmospheric fringes that reduce photometric precision, especially for faint sources and in multi-epoch co-additions. Here we present a method for constructing models of these atmospheric fringes using Principal Component Analysis that can be used to identify and remove these artifacts from contaminated images. In addition, we present the Uniform Background Indicator as a quantitative measurement of the reduced correlated background noise and photometric error present after removing fringes. We conclude by evaluating the effect of our method on measuring faint sources through the injection and recovery of artificial stars in both single-image epochs and co-additions. Our method for constructing atmospheric fringe models and applying those models to produce cleaned images is available for public download in the open source python package \href{https://github.com/MichaelMedford/fringez}{fringez}.

preprint2022arXiv

SNIascore: Deep Learning Classification of Low-Resolution Supernova Spectra

We present SNIascore, a deep-learning based method for spectroscopic classification of thermonuclear supernovae (SNe Ia) based on very low-resolution (R $\sim100$) data. The goal of SNIascore is fully automated classification of SNe Ia with a very low false-positive rate (FPR) so that human intervention can be greatly reduced in large-scale SN classification efforts, such as that undertaken by the public Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS). We utilize a recurrent neural network (RNN) architecture with a combination of bidirectional long short-term memory and gated recurrent unit layers. SNIascore achieves a $<0.6\%$ FPR while classifying up to $90\%$ of the low-resolution SN Ia spectra obtained by the BTS. SNIascore simultaneously performs binary classification and predicts the redshifts of secure SNe Ia via regression (with a typical uncertainty of $<0.005$ in the range from $z = 0.01$ to $z = 0.12$). For the magnitude-limited ZTF BTS survey ($\approx70\%$ SNe Ia), deploying SNIascore reduces the amount of spectra in need of human classification or confirmation by $\approx60\%$. Furthermore, SNIascore allows SN Ia classifications to be automatically announced in real-time to the public immediately following a finished observation during the night.

preprint2022arXiv

The Zwicky Transient Facility phase I sample of hydrogen-rich superluminous supernovae without strong narrow emission lines

We present a sample of 14 hydrogen-rich superluminous supernovae (SLSNe II) from the Zwicky Transient Facility (ZTF) between 2018 and 2020. We include all classified SLSNe with peaks $M_{g}<-20$ mag and with observed \emph{broad} but not narrow Balmer emission, corresponding to roughly 20 per cent of all hydrogen-rich SLSNe in ZTF phase I. We examine the light curves and spectra of SLSNe II and attempt to constrain their power source using light-curve models. The brightest events are photometrically and spectroscopically similar to the prototypical SN 2008es, while others are found spectroscopically more reminiscent of non-superluminous SNe II, especially SNe II-L. $^{56}$Ni decay as the primary power source is ruled out. Light-curve models generally cannot distinguish between circumstellar interaction (CSI) and a magnetar central engine, but an excess of ultraviolet (UV) emission signifying CSI is seen in most of the SNe with UV data, at a wide range of photometric properties. Simultaneously, the broad H$α$ profiles of the brightest SLSNe II can be explained through electron scattering in a symmetric circumstellar medium (CSM). In other SLSNe II without narrow lines, the CSM may be confined and wholly overrun by the ejecta. CSI, possibly involving mass lost in recent eruptions, is implied to be the dominant power source in most SLSNe II, and the diversity in properties is likely the result of different mass loss histories. Based on their radiated energy, an additional power source may be required for the brightest SLSNe II, however -- possibly a central engine combined with CSI.

preprint2022arXiv

Tracing Milky Way substructure with an RR Lyrae hierarchical clustering forest

RR Lyrae variable stars have long been reliable standard candles used to discern structure in the Local Group. With this in mind, we present a routine to identify groupings containing a statistically significant number of RR Lyrae variables in the Milky Way environment. RR Lyrae variable groupings, or substructures, with potential Galactic archaeology applications are found using a forest of agglomerative, hierarchical clustering trees, whose leaves are Milky Way RR Lyrae variables. Each grouping is validated by ensuring that the internal RR Lyrae variable proper motions are sufficiently correlated. Photometric information was collected from the Gaia second data release and proper motions from the (early) third data release. After applying this routine to the catalogue of 91234 variables, we are able to report sixteen unique RR Lyrae substructures with physical sizes of less than 1 kpc. Five of these substructures are in close proximity to Milky Way globular clusters with previously known tidal tails and/or a potential connection to Galactic merger events. One candidate substructure is in the neighbourhood of the Large Magellanic Cloud but is more distant (and older) than known satellites of the dwarf galaxy. Our study ends with a discussion of ways in which future surveys could be applied to the discovery of Milky Way stellar streams.

preprint2022arXiv

Zwicky Transient Facility and Globular Clusters: The RR Lyrae gri-Band Period-Luminosity-Metallicity and Period-Wesenheit-Metallicity Relations

Based on time-series observations collected from Zwicky Transient Facility (ZTF), we derived period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations for RR Lyrae located in globular clusters. We have applied various selection criteria to exclude RR Lyrae with problematic or spurious light curves. These selection criteria utilized information on the number of data points per light curve, amplitudes, colors, and residuals on the period-luminosity and/or period-Wesenheit relations. Due to blending, a number of RR Lyrae in globular clusters were found to be anomalously bright and have small amplitudes of their ZTF light curves. We used our final sample of ~750 RR Lyrae in 46 globular clusters covering a wide metallicity range (-2.36 < [Fe/H] < -0.54 dex) to derive PLZ and PWZ relations in gri bands. In addition, we have also derived the period-color-metallicity (PCZ) and for the first time, the PQZ relations where the Q-index is extinction-free by construction. We have compared our various relations to empirical and theoretical relations available in literature, and found a good agreement with most studies. Finally, we applied our derived PLZ relation to a dwarf galaxy, Crater II, and found its true distance modulus should be larger than the most recent determination.

preprint2021arXiv

Cataclysmic Variables in the Second Year of the Zwicky Transient Facility

Using a filter in the GROWTH Marshal based on color and the amplitude and the timescale of variability, we have identified 372 objects as known or candidate cataclysmic variables (CVs) during the second year of operation of the Zwicky Transient Facility (ZTF). From the available difference imaging data, we found that 93 are previously confirmed CVs, and 279 are strong candidates. Spectra of four of the candidates confirm them as CVs by the presence of Balmer emission lines, while one of the four has prominent HeII lines indicative of containing a magnetic white dwarf. Gaia EDR3 parallaxes are available for 154 of these systems, resulting in distances from 108-2096 pc and absolute magnitudes in the range of 7.5-15.0, with the largest number of candidates between 10.5-12.5. The total numbers are 21% higher than from the previous year of the survey with a greater number of distances available but a smaller percentage of systems close to the Galactic plane. Comparison of these findings with a machine learning method of searching all the light curves reveals large differences in each dataset related to the parameters involved in the search process.

preprint2021arXiv

Gaia GraL: Gaia DR2 Gravitational Lens Systems. VII. XMM-Newton Observations of Lensed Quasars

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA.

preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

preprint2020arXiv

TDE Hosts are Green and Centrally Concentrated: Signatures of a Post-Merger System

We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact "green valley" galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as measured by Sersic index, we find this green valley over representation is even larger. That is, concentrated red sequence galaxies are not producing TDEs at elevated levels. We present host galaxy spectra which show that E+A galaxies are overrepresented in the ZTF sample by a factor of $\approx$22, which is lower than previous TDE host galaxy studies have found. We find that this overrepresentation can be fully accounted for when taking into account the masses, colors, and Sérsic indices of the ZTF TDE hosts. The combination of both green colors and high Sérsic index of the typical TDE host galaxy could be explained if the TDE rate is temporarily enhanced following a merger that leads to a higher central concentration of stars.

preprint2020arXiv

Zwicky Transient Facility constraints on the optical emission from the nearby repeating FRB 180916.J0158+65

The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3σ$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $η\lesssim 200$, corresponding to an optical energy $E_{\rm opt} \lesssim 3 \times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $η\lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2\pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.

preprint2015arXiv

A possible close supermassive black-hole binary in a quasar with optical periodicity

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model. The recent availability of large collections of astronomical time series of flux measurements (light curves) offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 $\pm$ 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of $\sim9$ years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.

preprint2015arXiv

A systematic search for close supermassive black hole binaries in the Catalina Real-Time Transient Survey

Hierarchical assembly models predict a population of supermassive black hole (SMBH) binaries. These are not resolvable by direct imaging but may be detectable via periodic variability (or nanohertz frequency gravitational waves). Following our detection of a 5.2 year periodic signal in the quasar PG 1302-102 (Graham et al. 2015), we present a novel analysis of the optical variability of 243,500 known spectroscopically confirmed quasars using data from the Catalina Real-time Transient Survey (CRTS) to look for close (< 0.1 pc) SMBH systems. Looking for a strong Keplerian periodic signal with at least 1.5 cycles over a baseline of nine years, we find a sample of 111 candidate objects. This is in conservative agreement with theoretical predictions from models of binary SMBH populations. Simulated data sets, assuming stochastic variability, also produce no equivalent candidates implying a low likelihood of spurious detections. The periodicity seen is likely attributable to either jet precession, warped accretion disks or periodic accretion associated with a close SMBH binary system. We also consider how other SMBH binary candidates in the literature appear in CRTS data and show that none of these are equivalent to the identified objects. Finally, the distribution of objects found is consistent with that expected from a gravitational wave-driven population. This implies that circumbinary gas is present at small orbital radii and is being perturbed by the black holes. None of the sources is expected to merge within at least the next century. This study opens a new unique window to study a population of close SMBH binaries that must exist according to our current understanding of galaxy and SMBH evolution.

preprint2015arXiv

AstroStat - A VO Tool for Statistical Analysis

AstroStat is an easy-to-use tool for performing statistical analysis on data. It has been designed to be compatible with Virtual Observatory (VO) standards thus enabling it to become an integral part of the currently available collection of VO tools. A user can load data in a variety of formats into AstroStat and perform various statistical tests using a menu driven interface. Behind the scenes, all analysis is done using the public domain statistical software - R and the output returned is presented in a neatly formatted form to the user. The analyses performable include exploratory tests, visualizations, distribution fitting, correlation & causation, hypothesis testing, multivariate analysis and clustering. The tool is available in two versions with identical interface and features - as a web service that can be run using any standard browser and as an offline application. AstroStat will provide an easy-to-use interface which can allow for both fetching data and performing power statistical analysis on them.

preprint2015arXiv

Infrared Time Lags for the Periodic Quasar PG 1302-102

The optical light curve of the quasar PG 1302-102 at $z = 0.278$ shows a strong, smooth 5.2 yr periodic signal, detectable over a period of $\sim 20$ yr. Although the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim 10^{5}$ yr due to gravitational wave emission alone. Here we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from {\it WISE} and {\it Akari}, we confirm that the periodic behavior reported in the optical light curve from Graham et al. (2015) is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6 μ$m time lags of $(2219 \pm 153, 2408 \pm 148)$ days for a near face-on orientation of the torus, or $(4103\pm 153, 4292 \pm 148)$ days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the further out dusty regions. The implied distance of this dusty, assumed toroidal region is $\sim$ 1.5 pc for a near face-on geometry, or $\sim$1.1 pc for the relativistic Doppler boosted case.

preprint2015arXiv

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

preprint2014arXiv

DAMEWARE: A web cyberinfrastructure for astrophysical data mining

Astronomy is undergoing through a methodological revolution triggered by an unprecedented wealth of complex and accurate data. The new panchromatic, synoptic sky surveys require advanced tools for discovering patterns and trends hidden behind data which are both complex and of high dimensionality. We present DAMEWARE (DAta Mining & Exploration Web Application REsource): a general purpose, web-based, distributed data mining environment developed for the exploration of large datasets, and finely tuned for astronomical applications. By means of graphical user interfaces, it allows the user to perform classification, regression or clustering tasks with machine learning methods. Salient features of DAMEWARE include its capability to work on large datasets with minimal human intervention, and to deal with a wide variety of real problems such as the classification of globular clusters in the galaxy NGC1399, the evaluation of photometric redshifts and, finally, the identification of candidate Active Galactic Nuclei in multiband photometric surveys. In all these applications, DAMEWARE allowed to achieve better results than those attained with more traditional methods. With the aim of providing potential users with all needed information, in this paper we briefly describe the technological background of DAMEWARE, give a short introduction to some relevant aspects of data mining, followed by a summary of some science cases and, finally, we provide a detailed description of a template use case.

preprint2013arXiv

A comparison of period finding algorithms

This paper presents a comparison of popular period finding algorithms applied to the light curves of variable stars from the Catalina Real-time Transient Survey (CRTS), MACHO and ASAS data sets. We analyze the accuracy of the methods against magnitude, sampling rates, quoted period, quality measures (signal-to-noise and number of observations), variability, and object classes. We find that measure of dispersion-based techniques - analysis-of-variance with harmonics and conditional entropy - consistently give the best results but there are clear dependencies on object class and light curve quality. Period aliasing and identifying a period harmonic also remain significant issues. We consider the performance of the algorithms and show that a new conditional entropy-based algorithm is the most optimal in terms of completeness and speed. We also consider a simple ensemble approach and find that it performs no better than individual algorithms.

preprint2013arXiv

A novel variability-based method for quasar selection: evidence for a rest frame ~54 day characteristic timescale

We compare quasar selection techniques based on their optical variability using data from the Catalina Real-time Transient Survey (CRTS). We introduce a new technique based on Slepian wavelet variance (SWV) that shows comparable or better performance to structure functions and damped random walk models but with fewer assumptions. Combining these methods with WISE mid-IR colors produces a highly efficient quasar selection technique which we have validated spectroscopically. The SWV technique also identifies characteristic timescales in a time series and we find a characteristic rest frame timescale of ~54 days, confirmed in the light curves of ~18000 quasars from CRTS, SDSS and MACHO data, and anticorrelated with absolute magnitude. This indicates a transition between a damped random walk and $P(f) \propto f^{-1/3}$ behaviours and is the first strong indication that a damped random walk model may be too simplistic to describe optical quasar variability.

preprint2013arXiv

A plausible (overlooked) super-luminous supernova in the SDSS Stripe 82 data

We present the discovery of a plausible super-luminous supernova (SLSN), found in the archival data of Sloan Digital Sky Survey (SDSS) Stripe 82, called PSN 000123+000504. The supernova peaked at M_g<-21.3 mag in the second half of September 2005, but was missed by the real-time supernova hunt. The observed part of the light curve (17 epochs) showed that the rise to the maximum took over 30 days, while the decline time lasted at least 70 days (observed frame), closely resembling other SLSNe of SN2007bi type. Spectrum of the host galaxy reveals a redshift of z=0.281 and the distance modulus of μ=40.77 mag. Combining this information with the SDSS photometry, we found the host galaxy to be an LMC-like irregular dwarf galaxy with the absolute magnitude of M_B=-18.2+/-0.2 mag and the oxygen abundance of 12+log[O/H]=8.3+/-0.2. Our SLSN follows the relation for the most energetic/super-luminous SNe exploding in low-metallicity environments, but we found no clear evidence for SLSNe to explode in low-luminosity (dwarf) galaxies only. The available information on the PSN 000123+000504 light curve suggests the magnetar-powered model as a likely scenario of this event. This SLSN is a new addition to a quickly growing family of super-luminous SNe.

preprint2013arXiv

Feature Selection Strategies for Classifying High Dimensional Astronomical Data Sets

The amount of collected data in many scientific fields is increasing, all of them requiring a common task: extract knowledge from massive, multi parametric data sets, as rapidly and efficiently possible. This is especially true in astronomy where synoptic sky surveys are enabling new research frontiers in the time domain astronomy and posing several new object classification challenges in multi dimensional spaces; given the high number of parameters available for each object, feature selection is quickly becoming a crucial task in analyzing astronomical data sets. Using data sets extracted from the ongoing Catalina Real-Time Transient Surveys (CRTS) and the Kepler Mission we illustrate a variety of feature selection strategies used to identify the subsets that give the most information and the results achieved applying these techniques to three major astronomical problems.

preprint2013arXiv

Machine-assisted discovery of relationships in astronomy

High-volume feature-rich data sets are becoming the bread-and-butter of 21st century astronomy but present significant challenges to scientific discovery. In particular, identifying scientifically significant relationships between sets of parameters is non-trivial. Similar problems in biological and geosciences have led to the development of systems which can explore large parameter spaces and identify potentially interesting sets of associations. In this paper, we describe the application of automated discovery systems of relationships to astronomical data sets, focussing on an evolutionary programming technique and an information-theory technique. We demonstrate their use with classical astronomical relationships - the Hertzsprung-Russell diagram and the fundamental plane of elliptical galaxies. We also show how they work with the issue of binary classification which is relevant to the next generation of large synoptic sky surveys, such as LSST. We find that comparable results to more familiar techniques, such as decision trees, are achievable. Finally, we consider the reality of the relationships discovered and how this can be used for feature selection and extraction.

preprint2012arXiv

Connecting the time domain community with the Virtual Astronomical Observatory

The time domain has been identified as one of the most important areas of astronomical research for the next decade. The Virtual Observatory is in the vanguard with dedicated tools and services that enable and facilitate the discovery, dissemination and analysis of time domain data. These range in scope from rapid notifications of time-critical astronomical transients to annotating long-term variables with the latest modeling results. In this paper, we will review the prior art in these areas and focus on the capabilities that the VAO is bringing to bear in support of time domain science. In particular, we will focus on the issues involved with the heterogeneous collections of (ancillary) data associated with astronomical transients, and the time series characterization and classification tools required by the next generation of sky surveys, such as LSST and SKA.

preprint2011arXiv

Real Time Classification of Transient Events in Synoptic Sky Surveys

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This problem will grow by orders of magnitude with the next generation of surveys. We are exploring a variety of novel automated classification techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. We describe briefly some of the methods used.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.