Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,\sigma$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$.

Variable sources probe a wide range of astrophysical phenomena. We present a catalog of over ten million variable source candidates found in Data Release 1 (DR1) of the Zwicky Transient Facility (ZTF). We perform a periodicity search up to a frequency of 160 day^-1, and we classify the light curves into erratic and smooth variables. We also present variability indicators and the results of a periodicity search, up to a frequency of 5 day^-1, for about 1 billion sources in the ZTF-DR1 light curve database. We present several new short-period (<90 min) candidates, and about 60 new dwarf nova candidates, including two candidate eclipsing systems. Both the 10 million variables catalog and ~1 billion source catalog are available online in catsHTM format.

We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

We present a sample of supernovae Type IIn (SNe IIn) from the untargeted, magnitude-limited surveys of the Palomar Transient Factory (PTF) and its successor, the intermediate PTF (iPTF). The SNe IIn found and followed by the PTF/iPTF were used to select a sample of 42 events with useful constraints on the rise times as well as with available post-peak photometry. The sample SNe were discovered in 2009-2016 and have at least one low-resolution classification spectrum, as well as photometry from the P48 and P60 telescopes at Palomar Observatory. We study the light-curve properties of these SNe IIn using spline fits (for the peak and the declining portion) and template matching (for the rising portion). We find that the typical rise times are divided into fast and slow risers at $20\pm6$ d and $50\pm11$ d, respectively. The decline rates are possibly divided into two clusters, but this division has weak statistical significance. We find no significant correlation between the peak luminosity of SNe IIn and their rise times, but the more luminous SNe IIn are generally found to be more long-lasting. Slowly rising SNe IIn are generally found to decline slowly. The SNe in our sample were hosted by galaxies of absolute magnitude $-22 \lesssim M_g \lesssim -13$ mag. The K-corrections at light-curve peak of the SNe IIn in our sample are found to be within 0.2 mag for the observer's frame $r$-band, for SNe at redshifts $z < 0.25$. By applying K-corrections and also including ostensibly "superluminous" SNe IIn, we find that the peak magnitudes are $M_{\rm peak}^{r} = -19.18\pm1.32$ mag. We conclude that the occurrence of conspicuous light-curve bumps in SNe IIn, such as in iPTF13z, are limited to $1.4^{+14.6}_{-1.0} \%$ of the SNe IIn. We also investigate a possible sub-type of SNe IIn with a fast rise to a $\gtrsim 50$ d plateau followed by a slow, linear decline.

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

The Zwicky Transient Facility (ZTF) is currently surveying the entire northern sky, including dense Galactic plane fields. Here, we present preliminary results of the search for gravitational microlensing events in the ZTF data collected from the beginning of the survey (March 20, 2018) through June 30, 2019.

In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust.

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{\gamma,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $\nu L_\nu \lesssim 10^{37}\,$erg sec$^{-1}$. Taken together, we find that the early ($<10\,$days) data are best explained by shock breakout in a massive shell of dense circumstellar material ($0.02\,M_\odot$) at large radii ($3 \times 10^{14}\,$cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time ($>10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply-imaged Type Ia supernova iPTF16geu, at $z = 0.409$, using $Hubble\,Space\,Telescope$ ($HST$) observations in combination with supporting ground-based data. Multi-band photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using $HST$ and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of $\mu = 67.8^{+2.6}_{-2.9}$, accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: -0.23 $\pm$ 0.99, -1.43 $\pm$ 0.74 and 1.36 $\pm$ 1.07 days. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: $\Delta m_1 = -3.88^{+0.07}_{-0.06}$, $\Delta m_2 = -2.99^{+0.09}_{-0.08}$, $\Delta m_3 = -2.19^{+0.14}_{-0.15}$ and $\Delta m_4 = -2.40^{+0.14}_{-0.12}$ mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro- or millilensing, in addition to the large scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

(Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$\sigma$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Weak lensing shear estimation typically results in per galaxy statistical errors significantly larger than the sought after gravitational signal of only a few percent. These statistical errors are mostly a result of shape-noise -- an estimation error due to the diverse (and a-priori unknown) morphology of individual background galaxies. These errors are inversely proportional to the limiting angular resolution at which localized objects, such as galaxy clusters, can be probed with weak lensing shear. In this work we report on our initial attempt to reduce statistical errors in weak lensing shear estimation using a machine learning approach -- training a multi-layered convolutional neural network to directly estimate the shear given an observed background galaxy image. We train, calibrate and evaluate the performance and stability of our estimator using simulated galaxy images designed to mimic the distribution of HST observations of lensed background sources in the CLASH galaxy cluster survey. Using the trained estimator, we produce weak lensing shear maps of the cores of 20 galaxy clusters in the CLASH survey, demonstrating an RMS scatter reduced by approximately 26% when compared to maps produced with a commonly used shape estimator. This is equivalent to a survey speed enhancement of approximately 60%. However, given the non-transparent nature of the machine learning approach, this result requires further testing and validation. We provide python code to train and test this estimator on both simulated and real galaxy cluster observations. We also provide updated weak lensing catalogues for the 20 CLASH galaxy clusters studied.

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN2018fif was surrounded by relatively small amounts of circumstellar material (CSM) compared to all previous cases. This particularity, coupled with the high cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman (2017). Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=1174_{-81}^{+208} solar radii, and an ejected mass of M=5.6_{-1.0}^{+9.1} solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The large radius differs from previously modeled objects, and the difference could be either intrinsic or due to the relatively small amount of CSM around SN2018fif, perhaps making it a "cleaner" candidate for applying shock-cooling analytical models.

The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$\sigma$ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline.

We present the Census of the Local Universe (CLU) narrow-band survey to search for emission-line (\ha) galaxies. CLU-\ha~has imaged $\approx$3$\pi$ of the sky (26,470~deg$^2$) with 4 narrow-band filters that probe a distance out to 200~Mpc. We have obtained spectroscopic follow-up for galaxy candidates in 14 preliminary fields (101.6~deg$^2$) to characterize the limits and completeness of the survey. In these preliminary fields, CLU can identify emission lines down to an \ha~flux limit of $10^{-14}$~$\rm{erg~s^{-1}~cm^{-2}}$ at 90\% completeness, and recovers 83\% (67\%) of the \ha~flux from catalogued galaxies in our search volume at the $\Sigma$=2.5 ($\Sigma$=5) color excess levels. The contamination from galaxies with no emission lines is 61\% (12\%) for $\Sigma$=2.5 ($\Sigma$=5). Also, in the regions of overlap between our preliminary fields and previous emission-line surveys, we recover the majority of the galaxies found in previous surveys and identify an additional $\approx$300 galaxies. In total, we find 90 galaxies with no previous distance information, several of which are interesting objects: 7 blue compact dwarfs, 1 green pea, and a Seyfert galaxy; we also identified a known planetary nebula. These objects show that the CLU-\ha~survey can be a discovery machine for objects in our own Galaxy and extreme galaxies out to intermediate redshifts. However, the majority of the CLU-\ha~galaxies identified in this work show properties consistent with normal star-forming galaxies. CLU-\ha~galaxies with new redshifts will be added to existing galaxy catalogs to focus the search for the electromagnetic counterpart to gravitational wave events.

We present a novel algorithm for scheduling the observations of time-domain imaging surveys. Our Integer Linear Programming approach optimizes an observing plan for an entire night by assigning targets to temporal blocks, enabling strict control of the number of exposures obtained per field and minimizing filter changes. A subsequent optimization step minimizes slew times between each observation. Our optimization metric self-consistently weights contributions from time-varying airmass, seeing, and sky brightness to maximize the transient discovery rate. We describe the implementation of this algorithm on the surveys of the Zwicky Transient Facility and present its on-sky performance.

We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at $z\le0.09$ in a well-controlled sample over three years of operation (2010-2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realisations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_v=2.43\pm0.29\,\text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, $\propto\mathrm{t}^{-\beta}$, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher-redshifts. We find $\beta{\sim}1$, consistent with a progenitor channel governed by the gravitational in-spiral of binary white dwarfs.

Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade.

I present a software tool for solving the astrometry of astronomical images. The code puts emphasis on robustness against failures for correctly matching the sources in the image to a reference catalog, and on the stability of the solutions over the field of view (e.g., using orthogonal polynomials for the fitted transformation). The code was tested on over 50,000 images from various sources, including the Palomar Transient Factory (PTF) and the Zwicky Transient Facility (ZTF). The tested images equally represent low and high Galactic latitude fields and exhibit failure/bad-solution rate of <2x10^-5. Running on PTF 60-s integration images, and using the GAIA-DR2 as a reference catalog, the typical two-axes-combined astrometric root-mean square (RMS) is 14 mas at the bright end, presumably due to astrometric scintillation noise and systematic errors. I discuss the effects of seeing, airmass and the order of the transformation on the astrometric accuracy. The software, available online, is developed in MATLAB as part of an astronomical image processing environment and it can be run also as a stand-alone code.

SN2010jl was a luminous Type IIn supernova (SN), detected in radio, optical, X-ray and hard X-rays. Here we report on its six year R- and g-band light curves obtained using the Palomar Transient Factory. The light curve was generated using a pipeline based on the proper image subtraction method and we discuss the algorithm performances. As noted before, the R-band light curve, up to about 300 days after maximum light is well described by a power-law decline with a power-law index of about -0.5. Between day 300 and day 2300 after maximum light, it is consistent with a power-law decline, with a power-law index of about -3.4. The longevity of the light curve suggests that the massive circum-stellar material around the progenitor was ejected on time scales of at least tens of years prior to the progenitor explosion.

We discuss the late time (tens of days) emission from the radioactive ejecta of mergers involving neutron stars, when the ionization energy loss time of beta-decay electrons and positrons exceeds the expansion time. We show that if the e$^\pm$ are confined to the plasma (by magnetic fields), then the time dependence of the plasma heating rate, $\dot{\varepsilon}_d$, and hence of the bolometric luminosity $L=\dot{\varepsilon}_d$, are given by $d\log L/d\log t\simeq-2.8$, nearly independent of the composition and of the instantaneous radioactive energy release rate, $\dot{\varepsilon}$. This universality of the late time behavior is due to the weak dependence of the ionization loss rate on composition and on e$^\pm$ energy. The late time IR and optical measurements of GW 170817 are consistent with this expected behavior provided that the ionization loss time exceeds the expansion time at $t>t_\varepsilon\approx 7$~d, as predicted based on the early (few day) electromagnetic emission.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5um luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few of the heaviest isotopes in the third abundance peak with half-life around 14 days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.

High proper motion stars probe several extreme stellar populations, including nearby objects, Galactic halo stars, and hyper-velocity stars. Extending the search for high proper motion stars, to faint limits can increase their numbers and help to identify interesting targets. We conduct a search for faint (r>19.5 mag) high proper motion stars (>~200 mas/yr) by comparing the Sloan Digital Sky Survey (SDSS) - Data Release (DR) 10 catalog to the Pan-STARRS1-DR1 stacked image catalog. Our main selection criterion is stars that moved >1.5 arcsec and up to 7 arcsec between the SDSS and PS1 epochs. We identify 2923 high proper motion stars, of which 826 do not have proper motion in the GAIA-DR2 catalog and 565 are not listed in the GAIA-DR2 catalog. Our SDSS-PS1 proper motions are consistent with the GAIA-measured proper motions with a robust rms of about 10 mas/yr.

We present the best 265 sampled R-band light curves of spectroscopically identified Type Ia supernovae (SNe) from the Palomar Transient Factory (PTF; 2009-2012) survey and the intermediate Palomar Transient Factory (iPTF; 2013-2017). A model-independent light curve template is built from our data-set with the purpose to investigate average properties and diversity in our sample. We searched for multiple populations in the light curve properties using machine learning tools. We also utilised the long history of our light curves, up to 4000 days, to exclude any significant pre- or post- supernova flares. From the shapes of light curves we found the average rise time in the R band to be $16.8^{+0.5}_{-0.6}$ days. Although PTF/iPTF were single-band surveys, by modelling the residuals of the SNe in the Hubble-Lema\^{i}tre diagram, we estimate the average colour excess of our sample to be $<$E$($B$-$V$)> \approx 0.05(2)$ mag and thus the mean corrected peak brightness to be $M_R = -19.02\pm0.02$ $+5 \log( {\rm H}_0 [{\rm km} \cdot{\rm s}^{-1} {\rm Mpc}^{-1}]/70)$ mag with only weakly dependent on light curve shape. The intrinsic scatter is found to be $\sigma_R = 0.186 \pm 0.033$ mag for the redshift range $0.05<z<0.1$, without colour corrections of individual SNe. Our analysis shows that Malmquist bias becomes very significant at z=0.13. A similar limitation is expected for the ongoing Zwicky Transient Facility (ZTF) survey using the same telescope, but new camera expressly designed for ZTF.

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a Type Ic supernova with a fast evolving light curve indicating an extremely low ejecta mass ($\approx 0.2$ solar masses) and low kinetic energy ($\approx 2 \times 10^{50}$ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended He-rich envelope, likely ejected in an intense pre-explosion mass loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems.

Identification of linear features (streaks) in astronomical images is important for several reasons, including: detecting fast-moving near-Earth asteroids; detecting or flagging faint satellites streaks; and flagging or removing diffraction spikes, pixel bleeding, line-like cosmic rays and bad-pixel features. Here we discuss an efficient and optimal algorithm for the detection of such streaks. The optimal method to detect streaks in astronomical images is by cross-correlating the image with a template of a line broadened by the point spread function of the system. To do so efficiently, the cross-correlation of the streak position and angle is performed using the Radon transform, which is the integral of pixel values along all possible lines through an image. A fast version of the Radon transform exists, which we here extend to efficiently detect arbitrarily short lines. While the brute force Radon transform requires an order of N^3 operations for an N by N image, the fast Radon transform has a complexity of the order of N^2 log(N). We apply this method to simulated images, recovering the theoretical signal-to-noise ratio, and to real images, finding long streaks of low-Earth-orbit satellites and shorter streaks of Global Positioning System satellites. We detect streaks that are barely visible to the eye, out of hundreds of images, without a-priori knowledge of the streaks' positions or angles. We provide implementation of this algorithm in Python and MATLAB.

Most telescope designs have a circular-shape aperture. We demonstrate that telescopes with an elongated pupil have better contrast, at lower separations, between a bright central star and a faint companion. We simulate images for an elongated-pupil telescope and for a circular-pupil telescope of equal aperture area and integration time, investigating specifically what is the maximal contrast for finding faint companions around bright stars as a function of angular separation. We show that this design gives better contrast at lower separation from a bright star. This is shown for diffraction-limited (for perfect and imperfect optics) and seeing-limited speckle images, assuming equal aperture area and observing time. We also show the results are robust to errors in measurement of the point spread function. To compensate for the wider point spread function of the short axis, images should be taken at different rotation angles, either by rotating the telescope around the optical axis or by allowing a stationary mirror array to scan different parallactic angles with time. Images taken at different rotation angles are added using the proper image coaddition algorithms developed by Zackay & Ofek. The final image has the same contrast in all angles, rather than in specific areas of diffraction nulls. We obtained speckle observations with a small, ground based elongated-aperture telescope and show the results are consistent with simulations.

We present a simple analytic model, that captures the key features of the emission of radiation from material ejected by the merger of neutron stars (NS), and construct the multi-band and bolometric luminosity light curves of the transient associated with GW170817, AT\,2017gfo, using all available data. The UV to IR emission is shown to be consistent with a single $\approx0.05$\,M$_\odot$ component ejecta, with a power-law velocity distribution between $\approx 0.1\,c$ and $>0.3\,c$, a low opacity, {$\kappa<1$\,cm$^2$\,g$^{-1}$}, and a radioactive energy release rate consistent with an initial $Y_{\rm e}<0.4$. The late time spectra require an opacity of $\kappa_\nu\approx0.1$\,cm$^2$\,g$^{-1}$ at 1 to $2\mu$m. If this opacity is provided entirely by Lanthanides, their implied mass fraction is $X_{\rm Ln}\approx10^{-3}$, approximately 30 times below the value required to account for the solar abundance. The inferred value of $X_{\rm Ln}$ is uncertain due to uncertainties in the estimates of IR opacities of heavy elements, which also do not allow the exclusion of a significant contribution to the opacity by other elements (the existence of a slower ejecta rich in Lanthanides, that does not contribute significantly to the luminosity, can also not be ruled out). The existence of a relatively massive, $\approx 0.05$\,M$_\odot$, ejecta with high velocity and low opacity is in tension with the results of numerical simulations of NS mergers.

We present the discovery of a slowly-evolving, extragalactic radio transient, FIRST J141918.9+394036, identified by comparing a catalog of radio sources in nearby galaxies against new observations from the Very Large Array Sky Survey. Analysis of other archival data shows that FIRST J141918.9+394036 faded by a factor of ~50 over 23 years, from a flux of ~26 mJy at 1.4 GHz in 1993 to an upper limit of 0.4 mJy at 3 GHz in 2017. FIRST J141918.9+394036 is likely associated with the small star-forming galaxy SDSS J141918.81+394035.8 at a redshift z=0.01957 (d=87 Mpc), which implies a peak luminosity $\nu L_\nu \gtrsim 3\times10^{38}$ erg s$^{-1}$. If interpreted as an isotropic synchrotron blast wave, the source requires an explosion of kinetic energy ~10^{51} erg some time prior to our first detection in late 1993. This explosion could plausibly be associated with a long gamma-ray burst (GRB) or the merger of two neutron stars. Alternatively, FIRST J141918.9+394036 could be the nebula of a newly-born magnetar. The radio discovery of any of these phenomena would be unprecedented. Joint consideration of the event light curve, host galaxy, lack of a counterpart gamma-ray burst, and volumetric rate suggests that FIRST J141918.9+394036 is the afterglow of an off-axis (`orphan') long GRB. The long time baseline of this event offers the best available constraint in afterglow evolution as the bulk of shock-accelerated electrons become non-relativistic. The proximity, age, and precise localization of FIRST J141918.9+394036 make it a key object for understanding the aftermath of rare classes of stellar explosion.

We present visible-light and ultraviolet (UV) observations of the supernova PTF12glz. The SN was discovered and monitored in near-UV and R bands as part of a joint GALEX and Palomar Transient Factory campaign. It is among the most energetic Type IIn supernovae observed to date (~10^51erg). If the radiated energy mainly came from the thermalization of the shock kinetic energy, we show that PTF12glz was surrounded by ~1 solar mass of circumstellar material (CSM) prior to its explosive death. PTF12glz shows a puzzling peculiarity: at early times, while the freely expanding ejecta are presumably masked by the optically thick CSM, the radius of the blackbody that best fits the observations grows at ~8000km/s. Such a velocity is characteristic of fast moving ejecta rather than optically thick CSM. This phase of radial expansion takes place before any spectroscopic signature of expanding ejecta appears in the spectrum and while both the spectroscopic data and the bolometric luminosity seem to indicate that the CSM is optically thick. We propose a geometrical solution to this puzzle, involving an aspherical structure of the CSM around PTF12glz. By modeling radiative diffusion through a slab of CSM, we show that an aspherical geometry of the CSM can result in a growing effective radius. This simple model also allows us to recover the decreasing blackbody temperature of PTF12glz. SLAB-Diffusion, the code we wrote to model the radiative diffusion of photons through a slab of CSM and evaluate the observed radius and temperature, is made available on-line.

Calcium rich gap transients represent an intriguing new class of faint and fast evolving supernovae that exhibit strong [Ca II] emission in their nebular phase spectra. In this paper, we present the discovery and follow-up observations of iPTF 16hgs -- an intermediate luminosity and fast evolving transient that exhibited a double peaked light curve. Exhibiting a typical Type Ib spectrum in the photospheric phase and an early transition to a [Ca II] dominated nebular phase, we show that iPTF 16hgs shows properties consistent with the class of Ca-rich gap transients, with two interesting exceptions. First, while the second peak of the light curve is similar to other Ca-rich gap transients (suggesting $M_{ej}$ of 0.4 M$_\odot$ and peak luminosity of $3 \times 10^{41}$ ergs s$^{-1}$), we show that the first blue and fast declining (over $2$ days) peak is unique to this source. Second, with Integral Field Unit observations of the host galaxy, we find that iPTF 16hgs occurred in the outskirts (projected offset of $6$ kpc $ = 1.9 R_{eff}$) of a low metallicity (0.4 Z$_\odot$), star forming, dwarf spiral galaxy. Using deep late-time VLA and uGMRT observations, we place stringent limits on the local environment of the source, ruling out a large parameter space of circumstellar densities and mass loss environments of the progenitor. If iPTF 16hgs shares explosion physics with the class of Ca-rich gap transients, the presence of the first peak can be explained by enhanced mixing of 0.01 M$_\odot$ of $^{56}$Ni into the outer layers the ejecta, reminiscent of some models of He-shell detonations on WDs. On the other hand, if iPTF 16hgs is physically unrelated to the class, the first peak is consistent with shock cooling emission (of an envelope with a mass of 0.08 M$_\odot$ and radius of 13 R$_\odot$) associated with a core-collapse explosion of a highly stripped massive star in a close binary system.

We are building a next-generation laser adaptive optics system, Robo-AO-2, for the UH 2.2-m telescope that will deliver robotic, diffraction-limited observations at visible and near-infrared wavelengths in unprecedented numbers. The superior Maunakea observing site, expanded spectral range and rapid response to high-priority events represent a significant advance over the prototype. Robo-AO-2 will include a new reconfigurable natural guide star sensor for exquisite wavefront correction on bright targets and the demonstration of potentially transformative hybrid AO techniques that promise to extend the faintness limit on current and future exoplanet adaptive optics systems.

The mass of single neutron stars (NSs) can be measured using astrometric microlensing events. In such events, the center-of-light motion of a star lensed by a NS will deviate from the expected non-lensed motion and this deviation can be used to measure the mass of the NS. I search for future conjunctions between pulsars, with measured proper motion, and stars in the GAIA-DR2 catalog. I identified one candidate event of a star that will possibly be lensed by a pulsar during the next ten years in which the expected light deflection of the background star will deviate from the non-lensed motion by more than 0.05 mas. Given the position and proper motion of PSR J0846-3533, it will possibly pass ~0.2" from a 19.0 G magnitude background star in 2022.9. Assuming a 1.4 solar mass NS, the expected maximum deviation of the background star images from the uniform-rate plus parallax motion will be 0.091 mas. This pulsar position has a relatively large uncertainty and therefore additional observations are required in order to verify this event. I briefly discuss the opposite case, in which a pulsar is being lensed by a star. Such events can be used to measure the stellar mass via pulsar timing measurements. I do not find good candidates for such events with predicted variations in the pulsar period derivative, divided by 1 s, exceeding 10^-20 s^-1. Since only about 10% of the known pulsars have measured proper motions, there is potential for an increase in the number of predicted pulsar lensing events.

We investigate the light-curve properties of a sample of 26 spectroscopically confirmed hydrogen-poor superluminous supernovae (SLSNe-I) in the Palomar Transient Factory (PTF) survey. These events are brighter than SNe Ib/c and SNe Ic-BL, on average, by about 4 and 2~mag, respectively. The peak absolute magnitudes of SLSNe-I in rest-frame $g$ band span $-22\lesssim M_g \lesssim-20$~mag, and these peaks are not powered by radioactive $^{56}$Ni, unless strong asymmetries are at play. The rise timescales are longer for SLSNe than for normal SNe Ib/c, by roughly 10 days, for events with similar decay times. Thus, SLSNe-I can be considered as a separate population based on photometric properties. After peak, SLSNe-I decay with a wide range of slopes, with no obvious gap between rapidly declining and slowly declining events. The latter events show more irregularities (bumps) in the light curves at all times. At late times, the SLSN-I light curves slow down and cluster around the $^{56}$Co radioactive decay rate. Powering the late-time light curves with radioactive decay would require between 1 and 10${\rm M}_\odot$ of Ni masses. Alternatively, a simple magnetar model can reasonably fit the majority of SLSNe-I light curves, with four exceptions, and can mimic the radioactive decay of $^{56}$Co, up to $\sim400$ days from explosion. The resulting spin values do not correlate with the host-galaxy metallicities. Finally, the analysis of our sample cannot strengthen the case for using SLSNe-I for cosmology.

Fast access to large catalogs is required for some astronomical applications. Here we introduce the catsHTM tool, consisting of several large catalogs reformatted into HDF5-based file format, which can be downloaded and used locally. To allow fast access, the catalogs are partitioned into hierarchical triangular meshes and stored in HDF5 files. Several tools are provided to perform efficient cone searches at resolutions spanning from a few arc seconds to degrees, within a few milliseconds time. The first released version includes the following catalogs (by alphabetical order): 2MASS, 2MASS extended sources, AKARI, APASS, Cosmos, DECaLS/DR5, FIRST, GAIA/DR1, GAIA/DR2, GALEX/DR6Plus7, HSC/v2, IPHAS/DR2, NED redshifts, NVSS, Pan-STARRS1/DR1, PTF photometric catalog, ROSAT faint source, SDSS sources, SDSS/DR14 spectroscopy, Spitzer/SAGE, Spitzer/IRAC galactic center, UCAC4, UKIDSS/DR10, VST/ATLAS/DR3, VST/KiDS/DR3, WISE and XMM. We provide Python code that allows to perform cone searches, as well as MATLAB code for performing cone searches, catalog cross-matching, general searches, as well as load and create these catalogs.

There has been speculation of a class of relativistic explosions with an initial Lorentz factor smaller than that of classical Gamma-Ray Bursts (GRBs). These "dirty fireballs" would lack prompt GRB emission but could be pursued via their optical afterglow, appearing as transients that fade overnight. Here we report a search for such transients (transients that fade by 5-$\sigma$ in magnitude overnight) in four years of archival photometric data from the intermediate Palomar Transient Factory (iPTF). Our search criteria yielded 45 candidates. Of these, two were afterglows to GRBs that had been found in dedicated follow-up observations to triggers from the Fermi GRB Monitor (GBM). Another (iPTF14yb; Cenko et al. 2015) was a GRB afterglow discovered serendipitously. Two were spurious artifacts of reference image subtraction and one was an asteroid. The remaining 37 candidates have red stellar counterparts in external catalogs. The photometric and spectroscopic properties of the counterparts identify these transients as strong flares from M dwarfs of spectral type M3-M7 at distances of d ~ 0.15-2.1 kpc; two counterparts were already spectroscopically classified as late-type M stars. With iPTF14yb as the only confirmed relativistic outflow discovered independently of a high-energy trigger, we constrain the all-sky rate of transients that peak at m = 18 and fade by $\Delta$2 mag in $\Delta$3 hr to be 680 per year with a 68% confidence interval of 119-2236 per year. This implies that the rate of visible dirty fireballs is at most comparable to that of the known population of long-duration GRBs.

We present an optical and infrared (IR) study of IC 10 X-2, a high-mass X-ray binary in the galaxy IC 10. Previous optical and X-ray studies suggest X-2 is a Supergiant Fast X-ray Transient: a large-amplitude (factor of $\sim$ 100), short-duration (hours to weeks) X-ray outburst on 2010 May 21. We analyze R- and g-band light curves of X-2 from the intermediate Palomar Transient Factory taken between 2013 July 15 and 2017 Feb 14 show high-amplitude ($\gtrsim$ 1 mag), short-duration ($\lesssim8$ d) flares and dips ($\gtrsim$ 0.5 mag). Near-IR spectroscopy of X-2 from Palomar/TripleSpec show He I, Paschen-$\gamma$, and Paschen-$\beta$ emission lines with similar shapes and amplitudes as those of luminous blue variables (LBVs) and LBV candidates (LBVc). Mid-IR colors and magnitudes from Spitzer/IRAC photometry of X-2 resemble those of known LBV/LBVcs. We suggest that the stellar companion in X-2 is an LBV/LBVc and discuss possible origins of the optical flares. Dips in the optical light curve are indicative of eclipses from optically thick clumps formed in the winds of the stellar counterpart. Given the constraints on the flare duration ($0.02 - 0.8$ d) and the time between flares ($15.1\pm7.8$ d), we estimate the clump volume filling factor in the stellar winds, $f_V$, to be $0.01 < f_V < 0.71$, which overlaps with values measured from massive star winds. In X-2, we interpret the origin of the optical flares as the accretion of clumps formed in the winds of an LBV/LBVc onto the compact object.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

Current time domain facilities are finding several hundreds of transient astronomical events a year. The discovery rate is expected to increase in the future as soon as new surveys such as the Zwicky Transient Facility (ZTF) and the Large Synoptic Sky Survey (LSST) come on line. At the present time, the rate at which transients are classified is approximately one order or magnitude lower than the discovery rate, leading to an increasing "follow-up drought". Existing telescopes with moderate aperture can help address this deficit when equipped with spectrographs optimized for spectral classification. Here, we provide an overview of the design, operations and first results of the Spectral Energy Distribution Machine (SEDM), operating on the Palomar 60-inch telescope (P60). The instrument is optimized for classification and high observing efficiency. It combines a low-resolution (R$\sim$100) integral field unit (IFU) spectrograph with "Rainbow Camera" (RC), a multi-band field acquisition camera which also serves as multi-band (ugri) photometer. The SEDM was commissioned during the operation of the intermediate Palomar Transient Factory (iPTF) and has already proved lived up to its promise. The success of the SEDM demonstrates the value of spectrographs optimized to spectral classification. Introduction of similar spectrographs on existing telescopes will help alleviate the follow-up drought and thereby accelerate the rate of discoveries.

A Type IIn supernova (SN) is dominated by the interaction of SN ejecta with the circumstellar medium (CSM). Some SNe IIn (e.g., SN 2006jd) have episodes of re-brightening ("bumps") in their light curves. We present iPTF13z, a SN IIn discovered by the intermediate Palomar Transient Factory (iPTF) and characterised by several bumps in its light curve. We analyse this peculiar behaviour trying to infer the properties of the CSM and of the SN explosion, as well as the nature of its progenitor star. We obtained multi-band optical photometry for over 1000 days after discovery with the P48 and P60 telescopes at Palomar Observatory. We obtained low-resolution optical spectra in the same period. We did an archival search for progenitor outbursts. We analyse our photometry and spectra, and compare iPTF13z to other SNe IIn. A simple analytical model is used to estimate properties of the CSM. iPTF13z was a SN IIn showing a light curve with five bumps during its decline phase. The bumps had amplitudes between 0.4 and 0.9 mag and durations between 20 and 120 days. The most prominent bumps appeared in all our different optical bands. The spectra showed typical SN IIn characteristics, with emission lines of H$\alpha$ (with broad component FWHM ~$10^{3}-10^{4} ~{\rm ~km ~s^{-1}}$ and narrow component FWHM ~$10^2 \rm ~km ~s^{-1}$) and He I, but also with Fe II, Ca II, Na I D and H$\beta$ P-Cygni profiles (with velocities of ~$10^{3}$ ${\rm ~km ~s^{-1}}$). A pre-explosion outburst was identified lasting $\gtrsim 50$ days, with $M_r \approx -15$ mag around 210 days before discovery. Large, variable progenitor mass-loss rates (~> 0.01 $M_{\odot} \rm ~yr^{-1}$) and CSM densities (~> 10$^{-16}$ g cm$^{-3}$) are derived. We suggest that the light curve bumps of iPTF13z arose from SN ejecta interacting with denser regions in the CSM, possibly produced by the eruptions of a luminous blue variable star.

Detection of templates (e.g., sources) embedded in low-number count Poisson noise is a common problem in astrophysics. Examples include source detection in X-ray images, gamma-rays, UV, neutrinos, and search for clusters of galaxies and stellar streams. However, the solutions in the X-ray-related literature are sub-optimal -- in some cases by considerable factors. Using the lemma of Neyman-Pearson we derive the optimal statistics for template detection in the presence of Poisson noise. We demonstrate that this method provides higher completeness, for a fixed false-alarm probability value, compared with filtering the image with the point-spread function (PSF). In turn, we find that filtering by the PSF is better than filtering the image using the Mexican-hat wavelet (used by wavedetect). For some background levels, our method improves the sensitivity of source detection by more than a factor of two over the popular Mexican-hat wavelet filtering. This filtering technique can also be used also for fast PSF photometry and flare detection, and it is efficient, as well as straight forward to implement. We provide an implementation in MATLAB.

We present observations of two new hydrogen-poor superluminous supernovae (SLSN-I), iPTF15esb and iPTF16bad, showing late-time H-alpha emission with line luminosities of (1-3)e+41 erg/s and velocity widths of (4000-6000) km/s. Including the previously published iPTF13ehe, this makes up a total of three such events to date. iPTF13ehe is one of the most luminous and the slowest evolving SLSNe-I, whereas the other two are less luminous and fast decliners. We interpret this as a result of the ejecta running into a neutral H-shell located at a radius of ~ 1.0e+16cm. This implies that violent mass loss must have occurred several decades before the supernova explosion. Such a short time interval suggests that eruptive mass loss could be common shortly prior to the death of a massive star as a SLSN. And more importantly, helium is unlikely to be completely stripped off the progenitor stars and could be present in the ejecta. It is a mystery why helium features are not detected, even though non-thermal energy sources, capable of ionizing He atoms, may exist as suggested by the O II absorption series in the early time spectra. At late times (+240d), our spectra appear to have intrinsically lower [O I]6300A luminosities than that of SN2015bn and SN2007bi, possibly an indication of smaller oxygen masses (<10-30Msun). The blue-shifted H-alpha emission relative to the hosts for all three events may be in tension with the binary star model proposed for iPTF13ehe. Finally, iPTF15esb has a peculiar light curve with three peaks separated from one another by ~ 22 days. The LC undulation is higher in bluer bands. One possible explanation is eject-CSM interaction.

The localization of the repeating fast radio burst (FRB), FRB 121102, suggests that it is associated with a persistent radio-luminous compact source in the FRB host galaxy. Using the FIRST radio catalog, I present a search for luminous persistent sources in nearby galaxies, with radio luminosities >10% of the FRB 121102 persistent source luminosity. The galaxy sample contains about 30% of the total galaxy g-band luminosity within <108 Mpc, in a footprint of 10,600 deg^2. After rejecting sources likely due to active galactic nuclei activity or background sources, I remain with 11 candidates that are presumably associated with galactic disks or star formation regions. At least some of these candidates are likely to be due to chance alignment. In addition, I find 85 sources within 1" of galactic nuclei. Assuming the radio persistent sources are not related to galactic nuclei and that they follow the galaxy g-band light, the 11 sources imply a 95% confidence upper limit on the space density of luminous persistent sources of <5x10^-5 Mpc^-3, and that at any given time only a small fraction of galaxies host a radio luminous persistent source (<10^-3 L_*^-1). Assuming persistent sources life time of 100 yr, this implies a birth rate of <5x10^-7 yr^-1 Mpc^-3. Given the FRB volumetric rate, and assuming that all FRBs repeat and are associated with persistent radio sources, this sets a lower limit on the rate of FRB events per persistent source of >0.8 yr^-1. I argue that these 11 candidates are good targets for FRB searches and I estimate the FRB event rate from these candidates.

Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, while the second one (Event B) occurred over one month later and reached M(r) ~ -18 mag. By inspecting archival data, a faint source at the position of SN 2016bdu is detectable several times in the past few years. We interpret these detections as signatures of a phase of erratic variability, similar to that experienced by SN 2009ip between 2008 and mid-2012, and resembling the currently observed variability of the luminous blue variable SN 2000ch in NGC 3432. Spectroscopic monitoring of SN 2016bdu during the second peak initially shows features typical of a SN IIn. One month after the Event B maximum, the spectra develop broad Balmer lines with P Cygni profiles and broad metal features. At these late phases, the spectra resemble those of a typical Type II SN. All members of this SN 2009ip-like group are remarkably similar to the Type IIn SN 2005gl. For this object, the claim of a terminal SN explosion is supported by the disappearance of the progenitor star. The similarity with SN 2005gl suggests that all members of this family may finally explode as genuine SNe, although the unequivocal detection of nucleosynthesised elements in their nebular spectra is still missing.

We present results derived from four stellar occultations by the plutino object (208996) 2003~AZ$_{84}$, detected at January 8, 2011 (single-chord event), February 3, 2012 (multi-chord), December 2, 2013 (single-chord) and November 15, 2014 (multi-chord). Our observations rule out an oblate spheroid solution for 2003~AZ$_{84}$'s shape. Instead, assuming hydrostatic equilibrium, we find that a Jacobi triaxial solution with semi axes $(470 \pm 20) \times (383 \pm 10) \times (245 \pm 8)$~km % axis ratios $b/a= 0.82 \pm 0.05$ and $c/a= 0.52 \pm 0.02$, can better account for all our occultation observations. Combining these dimensions with the rotation period of the body (6.75~h) and the amplitude of its rotation light curve, we derive a density $\rho=0.87 \pm 0.01$~g~cm$^{-3}$ a geometric albedo $p_V= 0.097 \pm 0.009$. A grazing chord observed during the 2014 occultation reveals a topographic feature along 2003~AZ$_{84}$'s limb, that can be interpreted as an abrupt chasm of width $\sim 23$~km and depth $> 8$~km or a smooth depression of width $\sim 80$~km and depth $\sim 13$~km (or an intermediate feature between those two extremes).

We report the discovery of a multiply-imaged gravitationally lensed Type Ia supernova, iPTF16geu (SN 2016geu), at redshift $z=0.409$. This phenomenon could be identified because the light from the stellar explosion was magnified more than fifty times by the curvature of space around matter in an intervening galaxy. We used high spatial resolution observations to resolve four images of the lensed supernova, approximately 0.3" from the center of the foreground galaxy. The observations probe a physical scale of $\sim$1 kiloparsec, smaller than what is typical in other studies of extragalactic gravitational lensing. The large magnification and symmetric image configuration implies close alignment between the line-of-sight to the supernova and the lens. The relative magnifications of the four images provide evidence for sub-structures in the lensing galaxy.

We report the first maximum-light far-Ultraviolet to near-infrared spectra (1000A - 1.62um, rest) of a H-poor superluminous supernova, Gaia16apd. At z=0.1018, it is one of the closest and the UV brightest such events, with 17.4 (AB) magnitude in Swift UV band (1928A) at -11days pre-maximum. Assuming an exponential form, we derived the rise time of 33days and the peak bolometric luminosity of 3x10^{44}ergs^-1. At maximum light, the estimated photospheric temperature and velocity are 17,000K and 14,000kms^-1 respectively. The inferred radiative and kinetic energy are roughly 1x10^{51} and 2x10^{52}erg. Gaia16apd is extremely UV luminous, emitting 50% of its total luminosity at 1000 - 2500A. Compared to the UV spectra (normalized at 3100A) of well studied SN1992A (Ia), SN2011fe(Ia), SN1999em (IIP) and SN1993J (IIb), it has orders of magnitude more far-UV emission. This excess is interpreted primarily as a result of weaker metal line blanketing due to much lower abundance of iron-group elements in the outer ejecta. Because these elements originate either from the natal metallicity of the star, or have been newly produced, our observation provides direct evidence that little of these freshly synthesized material, including 56Ni, was mixed into the outer ejecta, and the progenitor metallicity is likely sub-solar. This disfavors Pair-Instability Supernova (PISN) models with Helium core masses >=90Msun, where substantial 56Ni material is produced. Higher photospheric temperature of Gaia16apd than that of normal SNe may also contribute to the observed far-UV excess. We find some indication that UV luminous SLSNe-I like Gaia16apd could be common. Using the UV spectra, we show that WFIRST could detect SLSNe-I out to redshift of 8.

Recent studies of the stellar population in the solar neighborhood (<20 pc) suggest that there are undetected white dwarfs (WDs) in multiple systems with main sequence companions. Detecting these hidden stars and obtaining a more complete census of nearby WDs is important for our understanding of binary and galactic evolution, as well as the study of explosive phenomena. In an attempt to uncover these hidden WDs, we present intermediate resolution spectroscopy over the wavelength range 3000-25000 \AA\ of 101 nearby M dwarfs (dMs), observed with the Very Large Telescope X-Shooter spectrograph. For each star we search for a hot component superimposed on the dM spectrum. X-Shooter has excellent blue sensitivity and thus can reveal a faint hot WD despite the brightness of its red companion. Visual examination shows no clear evidence of a WD in any of the spectra. We place upper limits on the effective temperatures of WDs that may still be hiding by fitting dM templates to the spectra, and modeling WD spectra. On average our survey is sensitive to WDs hotter than about 5300 K. This suggests that the frequency of WD companions of T<5300 K with separation of order <50 AU among the local dM population is <3% at the 95% confidence level. The reduced spectra are made available on via WISeREP repository.

We present late-time optical $R$-band imaging data from the Palomar Transient Factory (PTF) for the nearby type Ia supernova SN 2011fe. The stacked PTF light curve provides densely sampled coverage down to $R\simeq22$ mag over 200 to 620 days past explosion. Combining with literature data, we estimate the pseudo-bolometric light curve for this event from 200 to 1600 days after explosion, and constrain the likely near-infrared contribution. This light curve shows a smooth decline consistent with radioactive decay, except over ~450 to ~600 days where the light curve appears to decrease faster than expected based on the radioactive isotopes presumed to be present, before flattening at around 600 days. We model the 200-1600d pseudo-bolometric light curve with the luminosity generated by the radioactive decay chains of $^{56}$Ni, $^{57}$Ni and $^{55}$Co, and find it is not consistent with models that have full positron trapping and no infrared catastrophe (IRC); some additional energy escape other than optical/near-IR photons is required. However, the light curve is consistent with models that allow for positron escape (reaching 75% by day 500) and/or an IRC (with 85% of the flux emerging in non-optical wavelengths by day 600). The presence of the $^{57}$Ni decay chain is robustly detected, but the $^{55}$Co decay chain is not formally required, with an upper mass limit estimated at 0.014 M$_{\odot}$. The measurement of the $^{57}$Ni/$^{56}$Ni mass ratio is subject to significant systematic uncertainties, but all of our fits require a high ratio >0.031 (>1.3 in solar abundances).

With the advent of new wide-field, high-cadence optical transient surveys, our understanding of the diversity of core-collapse supernovae has grown tremendously in the last decade. However, the pre-supernova evolution of massive stars, that sets the physical backdrop to these violent events, is theoretically not well understood and difficult to probe observationally. Here we report the discovery of the supernova iPTF13dqy = SN 2013fs, a mere ~3 hr after explosion. Our rapid follow-up observations, which include multiwavelength photometry and extremely early (beginning at ~6 hr post-explosion) spectra, map the distribution of material in the immediate environment (<~ 10^15 cm) of the exploding star and establish that it was surrounded by circumstellar material (CSM) that was ejected during the final ~1 yr prior to explosion at a high rate, around 10^-3 solar masses per year. The complete disappearance of flash-ionised emission lines within the first several days requires that the dense CSM be confined to within <~ 10^15 cm, consistent with radio non-detections at 70--100 days. The observations indicate that iPTF13dqy was a regular Type II SN; thus, the finding that the probable red supergiant (RSG) progenitor of this common explosion ejected material at a highly elevated rate just prior to its demise suggests that pre-supernova instabilities may be common among exploding massive stars.

We present the Palomar Transient Factory discoveries and the photometric and spectroscopic observations of PTF11kmb and PTF12bho. We show that both transients have properties consistent with the class of calcium-rich gap transients, specifically lower peak luminosities and rapid evolution compared to ordinary supernovae, and a nebular spectrum dominated by [Ca II] emission. A striking feature of both transients is their host environments: PTF12bho is an intra-cluster transient in the Coma Cluster, while PTF11kmb is located in a loose galaxy group, at a physical offset ~150 kpc from the most likely host galaxy. Deep Subaru imaging of PTF12bho rules out an underlying host system to a limit of $M_R > -8.0$ mag, while Hubble Space Telescope imaging of PTF11kmb reveals a marginal counterpart that, if real, could be either a background galaxy or a globular cluster. We show that the offset distribution of Ca-rich gap transients is significantly more extreme than that seen for Type Ia supernovae or even short-hard gamma-ray bursts (sGRBs). Thus, if the offsets are caused by a kick, they require larger kick velocities and/or longer merger times than sGRBs. We also show that almost all Ca-rich gap transients found to date are in group and cluster environments with elliptical host galaxies, indicating a very old progenitor population; the remote locations could partially be explained by these environments having the largest fraction of stars in the intra-group/intra-cluster light following galaxy-galaxy interactions.

Type Ibn supernovae (SNe) are a small yet intriguing class of explosions whose spectra are characterized by low-velocity helium emission lines with little to no evidence for hydrogen. The prevailing theory has been that these are the core-collapse explosions of very massive stars embedded in helium-rich circumstellar material (CSM). We report optical observations of six new SNe Ibn: PTF11rfh, PTF12ldy, iPTF14aki, iPTF15ul, SN 2015G, and iPTF15akq. This brings the sample size of such objects in the literature to 22. We also report new data, including a near-infrared spectrum, on the Type Ibn SN 2015U. In order to characterize the class as a whole, we analyze the photometric and spectroscopic properties of the full Type Ibn sample. We find that, despite the expectation that CSM interaction would generate a heterogeneous set of light curves, as seen in SNe IIn, most Type Ibn light curves are quite similar in shape, declining at rates around 0.1 mag/day during the first month after maximum light, with a few significant exceptions. Early spectra of SNe Ibn come in at least two varieties, one that shows narrow P Cygni lines and another dominated by broader emission lines, both around maximum light, which may be an indication of differences in the state of the progenitor system at the time of explosion. Alternatively, the spectral diversity could arise from viewing-angle effects or merely from a lack of early spectroscopic coverage. Together, the relative light curve homogeneity and narrow spectral features suggest that the CSM consists of a spatially confined shell of helium surrounded by a less dense extended wind.

We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.

We present the light curves of the hydrogen-poor superluminous supernovae (SLSNe-I) PTF12dam and iPTF13dcc, discovered by the (intermediate) Palomar Transient Factory. Both show excess emission at early times and a slowly declining light curve at late times. The early bump in PTF12dam is very similar in duration (~10 days) and brightness relative to the main peak (2-3 mag fainter) compared to those observed in other SLSNe-I. In contrast, the long-duration (>30 days) early excess emission in iPTF13dcc, whose brightness competes with that of the main peak, appears to be of a different nature. We construct bolometric light curves for both targets, and fit a variety of light-curve models to both the early bump and main peak in an attempt to understand the nature of these explosions. Even though the slope of the late-time light-curve decline in both SLSNe is suggestively close to that expected from the radioactive decay of $^{56}$Ni and $^{56}$Co, the amount of nickel required to power the full light curves is too large considering the estimated ejecta mass. The magnetar model including an increasing escape fraction provides a reasonable description of the PTF12dam observations. However, neither the basic nor the double-peaked magnetar model is capable of reproducing the iPTF13dcc light curve. A model combining a shock breakout in an extended envelope with late-time magnetar energy injection provides a reasonable fit to the iPTF13dcc observations. Finally, we find that the light curves of both PTF12dam and iPTF13dcc can be adequately fit with the circumstellar medium (CSM) interaction model.

[abridged] In large-scale time-domain surveys, the processing of data, from procurement up to the detection of sources, is generally automated. One of the main challenges is contamination by artifacts, especially in regions of strong unresolved emission. We present a novel method for identifying candidates for variables and transients from the outputs of such surveys' data pipelines. We use the method to systematically search for novae in iPTF observations of the bulge of M31. We demonstrate that most artifacts produced by the iPTF pipeline form a locally uniform background of false detections approximately obeying Poissonian statistics, whereas genuine variables and transients as well as artifacts associated with bright stars result in clusters of detections, whose spread is determined by the source localization accuracy. This makes the problem analogous to source detection on images produced by X-ray telescopes, enabling one to utilize tools developed in X-ray astronomy. In particular, we use a wavelet-based source detection algorithm from the Chandra data analysis package CIAO. Starting from ~2.5x10^5 raw detections made by the iPTF data pipeline, we obtain ~4000 unique source candidates. Cross-matching these candidates with the source-catalog of a deep reference image, we find counterparts for ~90% of them. These are either artifacts due to imperfect PSF matching or genuine variable sources. The remaining ~400 detections are transient sources. We identify novae among these candidates by applying selection cuts based on the expected properties of nova lightcurves. Thus, we recovered all 12 known novae registered during the time span of the survey and discovered three nova candidates. Our method is generic and can be applied for mining any target out of the artifacts in optical time-domain data. As it is fully automated, its incompleteness can be accurately computed and corrected for.

We present the analysis of the first circumbinary planet microlensing event, OGLE-2007-BLG-349. This event has a strong planetary signal that is best fit with a mass ratio of $q \approx 3.4\times10^{-4}$, but there is an additional signal due to an additional lens mass, either another planet or another star. We find acceptable light curve fits with two classes of models: 2-planet models (with a single host star) and circumbinary planet models. The light curve also reveals a significant microlensing parallax effect, which constrains the mass of the lens system to be $M_L \approx 0.7 M_\odot$. Hubble Space Telescope images resolve the lens and source stars from their neighbors and indicate excess flux due to the star(s) in the lens system. This is consistent with the predicted flux from the circumbinary models, where the lens mass is shared between two stars, but there is not enough flux to be consistent with the 2-planet, 1-star models. So, only the circumbinary models are consistent with the HST data. They indicate a planet of mass $m_c = 80\pm 13\,M_\oplus$, orbiting a pair of M-dwarfs with masses of $M_A = 0.41\pm 0.07 M_\odot$ and $M_B = 0.30\pm 0.07 M_\oplus$, which makes this the lowest mass circumbinary planet system known. The ratio of the separation between the planet and the center-of-mass to the separations of the two stars is $\sim 40$, so unlike most of the circumbinary planets found by Kepler, the planet does not orbit near the stability limit.

We describe the near real-time transient-source discovery engine for the intermediate Palomar Transient Factory (iPTF), currently in operations at the Infrared Processing and Analysis Center (IPAC), Caltech. We coin this system the IPAC/iPTF Discovery Engine (or IDE). We review the algorithms used for PSF-matching, image subtraction, detection, photometry, and machine-learned (ML) vetting of extracted transient candidates. We also review the performance of our ML classifier. For a limiting signal-to-noise ratio of 4 in relatively unconfused regions, "bogus" candidates from processing artifacts and imperfect image subtractions outnumber real transients by ~ 10:1. This can be considerably higher for image data with inaccurate astrometric and/or PSF-matching solutions. Despite this occasionally high contamination rate, the ML classifier is able to identify real transients with an efficiency (or completeness) of ~ 97% for a maximum tolerable false-positive rate of 1% when classifying raw candidates. All subtraction-image metrics, source features, ML probability-based real-bogus scores, contextual metadata from other surveys, and possible associations with known Solar System objects are stored in a relational database for retrieval by the various science working groups. We review our efforts in mitigating false-positives and our experience in optimizing the overall system in response to the multitude of science projects underway with iPTF.

While optical and radio transient surveys have enjoyed a renaissance over the past decade, the dynamic infrared sky remains virtually unexplored. The infrared is a powerful tool for probing transient events in dusty regions that have high optical extinction, and for detecting the coolest of stars that are bright only at these wavelengths. The fundamental roadblocks in studying the infrared time-domain have been the overwhelmingly bright sky background (250 times brighter than optical) and the narrow field-of-view of infrared cameras (largest is 0.6 sq deg). To begin to address these challenges and open a new observational window in the infrared, we present Palomar Gattini-IR: a 25 sq degree, 300mm aperture, infrared telescope at Palomar Observatory that surveys the entire accessible sky (20,000 sq deg) to a depth of 16.4 AB mag (J band, 1.25um) every night. Palomar Gattini-IR is wider in area than every existing infrared camera by more than a factor of 40 and is able to survey large areas of sky multiple times. We anticipate the potential for otherwise infeasible discoveries, including, for example, the elusive electromagnetic counterparts to gravitational wave detections. With dedicated hardware in hand, and a F/1.44 telescope available commercially and cost-effectively, Palomar Gattini-IR will be on-sky in early 2017 and will survey the entire accessible sky every night for two years. Palomar Gattini-IR will pave the way for a dual hemisphere, infrared-optimized, ultra-wide field high cadence machine called Turbo Gattini-IR. To take advantage of the low sky background at 2.5 um, two identical systems will be located at the polar sites of the South Pole, Antarctica and near Eureka on Ellesmere Island, Canada. Turbo Gattini-IR will survey 15,000 sq. degrees to a depth of 20AB, the same depth of the VISTA VHS survey, every 2 hours with a survey efficiency of 97%.

While dozens of stellar mass black holes have been discovered in binary systems, isolated black holes have eluded detection. Their presence can be inferred when they lens light from a background star. We attempt to detect the astrometric lensing signatures of three photometrically identified microlensing events, OGLE-2011-BLG-0022, OGLE-2011-BLG-0125, and OGLE-2012-BLG-0169 (OB110022, OB110125, and OB120169), located toward the Galactic Bulge. These events were selected because of their long durations, which statistically favors more massive lenses. Astrometric measurements were made over 1-2 years using laser-guided adaptive optics observations from the W. M. Keck Observatory. Lens model parameters were first constrained by the photometric light curves. The OB120169 light curve is well-fit by a single-lens model, while both OB110022 and OB110125 light curves favor binary-lens models. Using the photometric fits as prior information, no significant astrometric lensing signal was detected and all targets were consistent with linear motion. The significant lack of astrometric signal constrains the lens mass of OB110022 to 0.05-1.79 Msun in a 99.7% confidence interval, which disfavors a black hole lens. Fits to OB110125 yielded a reduced Einstein crossing time and insufficient observations during the peak, so no mass limits were obtained. Two degenerate solutions exist for OB120169, which have a lens mass between 0.2-38.8 Msun and 0.4-39.8 Msun for a 99.7% confidence interval. Follow-up observations of OB120169 will further constrain the lens mass. Based on our experience, we use simulations to design optimal astrometric observing strategies and show that, with more typical observing conditions, detection of black holes is feasible.

The progenitors of some supernovae (SNe) exhibit outbursts with super-Eddington luminosities prior to their final explosions. This behavior is common among Type IIn SNe, but the driving mechanisms of these precursors are not yet well understood. SNHunt 275 was announced as a possible new SN during May 2015. Here we report on pre-explosion observations of the location of this event by the Palomar Transient Factory (PTF) and report the detection of a precursor about 500 days prior to the 2015 May activity (PTF 13efv). The observed velocities in the 2015 transient and its 2013 precursor absorption spectra are low (1000-2000 km/s), so it is not clear yet if the recent activity indeed marks the final disruption of the progenitor. Regardless of the nature of this event, we use the PTF photometric and spectral observations, as well as Swift-UVOT observations, to constrain the efficiency of the radiated energy relative to the total kinetic energy of the precursor. We find that, using an order-of-magnitude estimate and under the assumption of spherical symmetry, the ratio of the radiated energy to the kinetic energy is in the range of 0.04 to 3400.

We report the discovery of 2 new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven H-alpha emitters were discovered using the H-alpha imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3-m telescope at Lick observatory. Based on their spectral types, three H-alpha emitters were confirmed as Be stars with H-alpha equivalent widths > -10 Angstrom. Two objects were also observed by the new spectrograph SED-Machine on the Palomar 60 inch Telescope. The SED-Machine results show strong H-alpha emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-Machine can provide rapid observations for Be stars in a comprehensive survey in the future.

Transient detection and flux measurement via image subtraction stand at the base of time domain astronomy. Due to the varying seeing conditions, the image subtraction process is non-trivial, and existing solutions suffer from a variety of problems. Starting from basic statistical principles, we develop the optimal statistic for transient detection, flux measurement and any image-difference hypothesis testing. We derive a closed-form statistic that: (i) Is mathematically proven to be the optimal transient detection statistic in the limit of background-dominated noise; (ii) Is numerically stable; (iii) For accurately registered, adequately sampled images, does not leave subtraction or deconvolution artifacts; (iv) Allows automatic transient detection to the theoretical sensitivity limit by providing credible detection significance; (v) Has uncorrelated white noise; (vi) Is a sufficient statistic for any further statistical test on the difference image, and in particular, allows to distinguish particle hits and other image artifacts from real transients; (vii) Is symmetric to the exchange of the new and reference images; (viii) Is at least an order of magnitude faster to compute than some popular methods; and (ix) Is straightforward to implement. Furthermore, we present extensions of this method that make it resilient to registration errors, color-refraction errors, and any noise source that can be modelled. In addition, we show that the optimal way to prepare a reference image is the proper image coaddition presented in Zackay \& Ofek (2015b). We demonstrate this method on simulated data and real observations from the Palomar Transient Factory data release 2. We provide an implementation of this algorithm in MATLAB and Python.

The recently discovered fast radio bursts (FRBs), presumably of extra-galactic origin, have the potential to become a powerful probe of the intergalactic medium (IGM). We point out a few such potential applications. We provide expressions for the dispersion measure and rotation measure as a function of redshift, and we discuss the sensitivity of these measures to the HeII reionization and the IGM magnetic field. Finally we calculate the microlensing effect from an isolate, extragalctic stellar-mass compact object on the FRB spectrum. The time delays between the two lensing images will induce constructive and destructive interference, leaving a specific imprint on the spectra of FRBs. With a high all-sky rate, a large statistical sample of FRBs is expected to make these applications feasible.

Stacks of digital astronomical images are combined in order to increase image depth. The variable seeing conditions, sky background and transparency of ground-based observations make the coaddition process non-trivial. We present image coaddition methods optimized for source detection and flux measurement, that maximize the signal-to-noise ratio (S/N). We show that for these purposes the best way to combine images is to apply a matched filter to each image using its own point spread function (PSF) and only then to sum the images with the appropriate weights. Methods that either match filter after coaddition, or perform PSF homogenization prior to coaddition will result in loss of sensitivity. We argue that our method provides an increase of between a few and 25 percent in the survey speed of deep ground-based imaging surveys compared with weighted coaddition techniques. We demonstrate this claim using simulated data as well as data from the Palomar Transient Factory data release 2. We present a variant of this coaddition method which is optimal for PSF or aperture photometry. We also provide an analytic formula for calculating the S/N for PSF photometry on single or multiple observations. In the next paper in this series we present a method for image coaddition in the limit of background-dominated noise which is optimal for any statistical test or measurement on the constant-in-time image (e.g., source detection, shape or flux measurement or star-galaxy separation), making the original data redundant. We provide an implementation of this algorithm in MATLAB.

Image coaddition is one of the most basic operations that astronomers perform. In Paper~I, we presented the optimal ways to coadd images in order to detect faint sources and to perfrom flux measurements under the assumption that the noise is approximately Gaussian. Here, we build on these results and derive from first principles a coaddition technique which is optimal for any hypothesis testing and measurement (e.g., source detection, flux or shape measurements and star/galaxy separation), in the background-noise-dominated case. This method has several important properties. The pixels of the resulting coadd image are uncorrelated. This image preserves all the information (from the original individual images) on all spatial frequencies. Any hypothesis testing or measurement that can be done on all the individual images simultaneously, can be done on the coadded image without any loss of information. The PSF of this image is typically as narrow, or narrower than the PSF of the best image in the ensemble. Moreover, this image is practically indistinguishable from a regular single image, meaning that any code that measures any property on a regular astronomical image can be applied to it unchanged. In particular, the optimal source detection statistic derived in paper~I is reproduced by matched filtering this image with its own PSF. This coaddition process, which we call proper coaddition, can be understood as a the maximum signal-to-noise ratio measurement of the Fourier transform of the image, weighted in such a way that the noise in the entire Fourier domain is of equal variance. This method has important implications for multi-epoch seeing-limited deep surveys, weak lensing galaxy shape measurements, and diffraction-limited imaging via speckle observations. The last topic will be covered in depth in future papers. We provide an implementation of this algorithm in MATLAB.

Supernovae (SNe) embedded in dense circumstellar material (CSM) may show prominent emission lines in their early-time spectra ($\leq 10$ days after the explosion), owing to recombination of the CSM ionized by the shock-breakout flash. From such spectra ("flash spectroscopy"), we can measure various physical properties of the CSM, as well as the mass-loss rate of the progenitor during the year prior to its explosion. Searching through the Palomar Transient Factory (PTF and iPTF) SN spectroscopy databases from 2009 through 2014, we found 12 Type II SNe showing flash-ionized (FI) signatures in their first spectra. All are younger than 10 days. These events constitute 14\% of all 84 SNe in our sample having a spectrum within 10 days from explosion, and 18\% of SNe~II observed at ages $<5$ days, thereby setting lower limits on the fraction of FI events. We classified as "blue/featureless" (BF) those events having a first spectrum which is similar to that of a black body, without any emission or absorption signatures. It is possible that some BF events had FI signatures at an earlier phase than observed, or that they lack dense CSM around the progenitor. Within 2 days after explosion, 8 out of 11 SNe in our sample are either BF events or show FI signatures. Interestingly, we found that 19 out of 21 SNe brighter than an absolute magnitude $M_R=-18.2$ belong to the FI or BF groups, and that all FI events peaked above $M_R=-17.6$ mag, significantly brighter than average SNe~II.

During the first few days after explosion, Type II supernovae (SNe) are dominated by relatively simple physics. Theoretical predictions regarding early-time SN light curves in the ultraviolet (UV) and optical bands are thus quite robust. We present, for the first time, a sample of $57$ $R$-band Type II SN light curves that are well monitored during their rise, having $>5$ detections during the first 10 days after discovery, and a well-constrained time of explosion to within $1-3$ days. We show that the energy per unit mass ($E/M$) can be deduced to roughly a factor of five by comparing early-time optical data to the model of Rabinak & Waxman (2011), while the progenitor radius cannot be determined based on $R$-band data alone. We find that Type II SN explosion energies span a range of $E/M=(0.2-20)\times 10^{51} \; \rm{erg/(10 M}_\odot$), and have a mean energy per unit mass of $\left\langle E/M \right\rangle = 0.85\times 10^{51} \; \rm{erg/(10 M}_\odot$), corrected for Malmquist bias. Assuming a small spread in progenitor masses, this indicates a large intrinsic diversity in explosion energy. Moreover, $E/M$ is positively correlated with the amount of $^{56}\rm{Ni}$ produced in the explosion, as predicted by some recent models of core-collapse SNe. We further present several empirical correlations. The peak magnitude is correlated with the decline rate ($\Delta m_{15}$), the decline rate is weakly correlated with the rise time, and the rise time is not significantly correlated with the peak magnitude. Faster declining SNe are more luminous and have longer rise times. This limits the possible power sources for such events.

The nature of fast radio bursts (FRB) has been extensively debated. Here we investigate FRB121102, detected at Arecibo telescope and remarkable for its unusually large spectral index. After extensive study we conclude that the spectral index is caused by a nebula with free-free absorption. We find that putative nebula must lie beyond the Milky Way. We conclude that FRBs are of extra-galactic origin and that they arise in dense star-forming regions. The challenge with extra-galactic models is the the high volumetric rate of FRBs. This high rate allows us to eliminate all models of catastrophic stellar deaths. Hyper-giant flares from young magnetars emerge as the most likely progenitors. Some of the consequences are: (i) Intergalactic FRB models can be safely ignored. (ii) The rich ISM environment of young magnetars can result in significant contribution to DM, Rotation Measure (RM) and in some cases to significant free-free optical depth. (iii) The star-forming regions in the host galaxies can contribute significantly to the DM. Including this contribution reduces the inferred distances to FRBs and correspondingly increases the volumetric rate of FRBs (and, in turn, may require that giant flares can also produce FRBs). (iv) FRBs are likely to be suppressed at lower frequencies. Conversely, searching for FRBs at higher frequencies (2-5 GHz) would be attractive. (v) The blast wave which produces the radio emission can undergo rapid deceleration if the circum-burst medium is dense (as maybe the case for FRB121102), leading to X-ray, radio and possibly gamma-ray emission. (vi) Galaxies with high star formation rate host will have a higher FRB rate. However, such FRBs will have differing DMs owing to differing local contributions. (vi) The DM and RM of FRBs will prove to be noisy probes of the intergalactic medium (density, magnetic field) and cosmography.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.

iPTF13ehe is a hydrogen-poor superluminous supernova (SLSN) at z=0.3434, with a slow-evolving light curve and spectral features similar to SN2007bi. It rises within (83-148)days (rest-frame) to reach a peak bolometric luminosity of 1.3x$10^{44}$erg/s, then decays very slowly at 0.015mag. per day. The measured ejecta velocity is 13000km/s. The inferred explosion characteristics, such as the ejecta mass (67-220$M_\odot$), the total radiative and kinetic energy ($10^{51}$ & 2x$10^{53}$erg respectively), is typical of a slow-evolving H-poor SLSN event. However, the late-time spectrum taken at +251days reveals a Balmer Halpha emission feature with broad and narrow components, which has never been detected before among other H-poor SLSNe. The broad component has a velocity width of ~4500km/s and has a ~300km/s blue-ward shift relative to the narrow component. We interpret this broad H$\alpha$ emission with luminosity of $\sim$2$\times10^{41}$\,erg\,s$^{-1}$ as resulting from the interaction between the supernova ejecta and a discrete H-rich shell, located at a distance of $\sim4\times10^{16}$\,cm from the explosion site. This ejecta-CSM interaction causes the rest-frame r-band LC to brighten at late times. The fact that the late-time spectra are not completely absorbed by the shock ionized CSM shell implies that its Thomson scattering optical depth is likely <1, thus setting upper limits on the CSM mass <30$M_\odot$ and the volume number density <4x$10^8cm^{-3}$. Of the existing models, a Pulsational Pair Instability Supernova model can naturally explain the observed 30$M_\odot$ H-shell, ejected from a progenitor star with an initial mass of (95-150)$M_\odot$ about 40 years ago. We estimate that at least $\sim$15\%\ of all SLSNe-I may have late-time Balmer emission lines.

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.

The progenitor stars of several Type IIb supernovae (SNe) show indications for extended hydrogen envelopes. These envelopes might be the outcome of luminous energetic pre-explosion events, so-called precursor eruptions. We use the Palomar Transient Factory (PTF) pre-explosion observations of a sample of 27 nearby Type IIb SNe to look for such precursors during the final years prior to the SN explosion. No precursors are found when combining the observations in 15-day bins, and we calculate the absolute-magnitude-dependent upper limit on the precursor rate. At the 90% confidence level, Type IIb SNe have on average $<0.86$ precursors as bright as absolute $R$-band magnitude $-14$ in the final 3.5 years before the explosion and $<0.56$ events over the final year. In contrast, precursors among SNe IIn have a $\gtrsim 5$ times higher rate. The kinetic energy required to unbind a low-mass stellar envelope is comparable to the radiated energy of a few-weeks-long precursor which would be detectable for the closest SNe in our sample. Therefore, mass ejections, if they are common in such SNe, are radiatively inefficient or have durations longer than months. Indeed, when using 60-day bins a faint precursor candidate is detected prior to SN 2012cs ($\sim2$% false-alarm probability). We also report the detection of the progenitor of SN 2011dh which does not show detectable variability over the final two years before the explosion. The suggested progenitor of SN 2012P is still present, and hence is likely a compact star cluster, or an unrelated object.

The IceCube neutrino observatory pursues a follow-up program selecting interesting neutrino events in real-time and issuing alerts for electromagnetic follow-up observations. In March 2012, the most significant neutrino alert during the first three years of operation was issued by IceCube. In the follow-up observations performed by the Palomar Transient Factory (PTF), a Type IIn supernova (SN) PTF12csy was found $0.2^\circ$ away from the neutrino alert direction, with an error radius of $0.54^\circ$. It has a redshift of $z=0.0684$, corresponding to a luminosity distance of about $300 \, \mathrm{Mpc}$ and the Pan-STARRS1 survey shows that its explosion time was at least 158 days (in host galaxy rest frame) before the neutrino alert, so that a causal connection is unlikely. The a posteriori significance of the chance detection of both the neutrinos and the SN at any epoch is $2.2 \, \sigma$ within IceCube's 2011/12 data acquisition season. Also, a complementary neutrino analysis reveals no long-term signal over the course of one year. Therefore, we consider the SN detection coincidental and the neutrinos uncorrelated to the SN. However, the SN is unusual and interesting by itself: It is luminous and energetic, bearing strong resemblance to the SN IIn 2010jl, and shows signs of interaction of the SN ejecta with a dense circumstellar medium. High-energy neutrino emission is expected in models of diffusive shock acceleration, but at a low, non-detectable level for this specific SN. In this paper, we describe the SN PTF12csy and present both the neutrino and electromagnetic data, as well as their analysis.

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are predicted to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. Here we summarize the results of the first systematic search for gamma-ray emission in Fermi-LAT data in the energy range from 100 MeV to 300 GeV from a large sample of SNe exploding in dense CSM. We search for a gamma-ray excess at the position of 147 SNe Type IIn in a one year time window after the optical peak time. In addition we combine the closest and optically brightest sources of our sample in a joint likelihood analysis in three different time windows (3, 6 and 12 months). No excess gamma-ray emission is found and limits on the gamma-ray luminosity and the ratio of gamma-ray to optical luminosity are presented.

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).

We present ultraviolet photometry (NUV band, 180--280 nm) of 405 asteroids observed serendipitously by the Galaxy Evolution Explorer (GALEX) from 2003--2012. All asteroids in this sample were detected by GALEX at least twice. Unambiguous visible-color-based taxonomic labels (C type versus S type) exist for 315 of these asteroids; of these, thermal-infrared-based diameters are available for 245. We derive NUV-V color using two independent models to predict the visual magnitude V at each NUV-detection epoch. Both V models produce NUV-V distributions in which the S types are redder than C types with more than 8-sigma confidence. This confirms that the S types' redder spectral slopes in the visible remain redder than the C types' into the NUV, this redness being consistent with absorption by silica-containing rocks. The GALEX asteroid data confirm earlier results from the International Ultraviolet Explorer, which two decades ago produced the only other sizeable set of UV asteroid photometry. The GALEX-derived NUV-V data also agree with previously published Hubble Space Telescope (HST) UV observations of asteroids 21 Lutetia and 1 Ceres. Both the HST and GALEX data indicate that NUV band is less useful than u band for distinguishing subgroups within the greater population of visible-color-defined C types (notably, M types and G types).

We fit 54,296 sparsely-sampled asteroid lightcurves in the Palomar Transient Factory to a combined rotation plus phase-function model. Each lightcurve consists of 20+ observations acquired in a single opposition. Using 805 asteroids in our sample that have reference periods in the literature, we find the reliability of our fitted periods is a complicated function of the period, amplitude, apparent magnitude and other attributes. Using the 805-asteroid ground-truth sample, we train an automated classifier to estimate (along with manual inspection) the validity of the remaining 53,000 fitted periods. By this method we find 9,033 of our lightcurves (of 8,300 unique asteroids) have reliable periods. Subsequent consideration of asteroids with multiple lightcurve fits indicate 4% contamination in these reliable periods. For 3,902 lightcurves with sufficient phase-angle coverage and either a reliably-fit period or low amplitude, we examine the distribution of several phase-function parameters, none of which are bimodal though all correlate with the bond albedo and with visible-band colors. Comparing the theoretical maximal spin rate of a fluid body with our amplitude versus spin-rate distribution suggests that, if held together only by self-gravity, most asteroids are in general less dense than 2 g/cm$^3$, while C types have a lower limit of between 1 and 2 g/cm$^3$, in agreement with previous density estimates. For 5-20km diameters, S types rotate faster and have lower amplitudes than C types. If both populations share the same angular momentum, this may indicate the two types' differing ability to deform under rotational stress. Lastly, we compare our absolute magnitudes and apparent-magnitude residuals to those of the Minor Planet Center's nominal $G=0.15$, rotation-neglecting model; our phase-function plus Fourier-series fitting reduces asteroid photometric RMS scatter by a factor of 3.

We present near-infrared (NIR) time-series spectroscopy, as well as complementary ultraviolet (UV), optical, and NIR data, of the Type Ia supernova (SN Ia) iPTF13ebh, which was discovered within two days from the estimated time of explosion. The first NIR spectrum was taken merely 2.3 days after explosion and may be the earliest NIR spectrum yet obtained of a SN Ia. The most striking features in the spectrum are several NIR C I lines, and the C I {\lambda}1.0693 {\mu}m line is the strongest ever observed in a SN Ia. Interestingly, no strong optical C II counterparts were found, even though the optical spectroscopic time series began early and is densely-cadenced. Except at the very early epochs, within a few days from the time of explosion, we show that the strong NIR C I compared to the weaker optical C II appears to be general in SNe Ia. iPTF13ebh is a fast decliner with {\Delta}m15(B) = 1.79 $\pm$ 0.01, and its absolute magnitude obeys the linear part of the width-luminosity relation. It is therefore categorized as a "transitional" event, on the fast-declining end of normal SNe Ia as opposed to subluminous/91bg-like objects. iPTF13ebh shows NIR spectroscopic properties that are distinct from both the normal and subluminous/91bg-like classes, bridging the observed characteristics of the two classes. These NIR observations suggest composition and density of the inner core similar to that of 91bg-like events, and a deep reaching carbon burning layer not observed in slower declining SNe Ia. There is also a substantial difference between the explosion times inferred from the early-time light curve and the velocity evolution of the Si II {\lambda}0.6355 {\mu}m line, implying a long dark phase of ~ 4 days.

We reanalyze microlensing events in the published list of anomalous events that were observed from the OGLE lensing survey conducted during 2004-2008 period. In order to check the existence of possible degenerate solutions and extract extra information, we conduct analyses based on combined data from other survey and follow-up observation and consider higher-order effects. Among the analyzed events, we present analyses of 8 events for which either new solutions are identified or additional information is obtained. We find that the previous binary-source interpretations of 5 events are better interpreted by binary-lens models. These events include OGLE-2006-BLG-238, OGLE-2007-BLG-159, OGLE-2007-BLG-491, OGLE-2008-BLG-143, and OGLE-2008-BLG-210. With additional data covering caustic crossings, we detect finite-source effects for 6 events including OGLE-2006-BLG-215, OGLE-2006-BLG-238, OGLE-2006-BLG-450, OGLE-2008-BLG-143, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. Among them, we are able to measure the Einstein radii of 3 events for which multi-band data are available. These events are OGLE-2006-BLG-238, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. For OGLE-2008-BLG-143, we detect higher-order effect induced by the changes of the observer's position caused by the orbital motion of the Earth around the Sun. In addition, we present degenerate solutions resulting from the known close/wide or ecliptic degeneracy. Finally, we note that the masses of the binary companions of the lenses of OGLE-2006-BLG-450 and OGLE-2008-BLG-210 are in the brown-dwarf regime.

The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 solar, explosion energies of 10^51 erg, and ejecta masses of 10 solar masses. Exploding blue supergiants and Wolf-Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1d after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission, is expected to find >100 SNe per year (~0.5 SN per deg^2), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life.

Impulsive radio bursts that are detectable across cosmological distances constitute extremely powerful probes of the ionized Inter-Galactic Medium (IGM), intergalactic magnetic fields, and the properties of space-time itself. Their dispersion measures (DMs) will enable us to detect the "missing" baryons in the low-redshift Universe and make the first measurements of the mean galaxy halo profile, a key parameter in models of galaxy formation and feedback. Impulsive bursts can be used as cosmic rulers at redshifts exceeding 2, and constrain the dark energy equation-of-state parameter, $w(z)$ at redshifts beyond those readily accessible by Type Ia SNe. Both of these goals are realisable with a sample of $\sim 10^4$ fast radio bursts (FRBs) whose positions are localized to within one arcsecond, sufficient to obtain host galaxy redshifts via optical follow-up. It is also hypothesised that gravitational wave events may emit coherent emission at frequencies probed by SKA1-LOW, and the localization of such events at cosmological distances would enable their use as cosmological standard sirens. To perform this science, such bursts must be localized to their specific host galaxies so that their redshifts may be obtained and compared against their dispersion measures, rotation measures, and scattering properties. The SKA can achieve this with a design that has a wide field-of-view, a substantial fraction of its collecting area in a compact configuration (80\% within a 3\,km radius), and a capacity to attach high-time-resolution instrumentation to its signal path.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

Astronomical radio bursts disperse while traveling through the interstellar medium. To optimally detect a short-duration signal within a frequency band, we have to precisely compensate for the pulse dispersion, which is a computationally demanding task. We present the Fast Dispersion Measure Transform (FDMT) algorithm for optimal detection of such signals. Our algorithm has a low theoretical complexity of 2N_f N_t+ N_t N_d log_2(N_f) where N_f, N_t and N_d are the numbers of frequency bins, time bins, and dispersion measure bins, respectively. Unlike previously suggested fast algorithms our algorithm conserves the sensitivity of brute force dedispersion. Our tests indicate that this algorithm, running on a standard desktop computer, and implemented in a high-level programming language, is already faster than the state of the art dedispersion codes running on graphical processing units (GPUs). We also present a variant of the algorithm that can be efficiently implemented on GPUs. The latter algorithm's computation and data transport requirements are similar to those of two-dimensional FFT, indicating that incoherent dedispersion can now be considered a non-issue while planning future surveys. We further present a fast algorithm for sensitive dedispersion of pulses shorter than normally allowed by incoherent dedispersion. In typical cases this algorithm is orders of magnitude faster than coherent dedispersion by convolution. We analyze the computational complexity of pulsed signal searches by radio interferometers. We conclude that, using our suggested algorithms, maximally sensitive blind searches for such pulses is feasible using existing facilities. We provide an implementation of these algorithms in Python and MATLAB.

We present an analysis of the early, rising light curves of 18 Type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ). We fit these early data flux using a simple power-law $(f(t) = {\alpha\times t^n})$ to determine the time of first light $({t_0})$, and hence the rise-time $({t_{rise}})$ from first light to peak luminosity, and the exponent of the power-law rise ($n$). We find a mean uncorrected rise time of $18.98 {\pm} 0.54$ days, with individual SN rise-times ranging from $15.98$ to $24.7$ days. The exponent n shows significant departures from the simple 'fireball model' of $n = 2$ (or ${f(t) \propto t^2}$) usually assumed in the literature. With a mean value of $n = 2.44 {\pm} 0.13$, our data also show significant diversity from event to event. This deviation has implications for the distribution of 56Ni throughout the SN ejecta, with a higher index suggesting a lesser degree of 56Ni mixing. The range of n found also confirms that the 56Ni distribution is not standard throughout the population of SNe Ia, in agreement with earlier work measuring such abundances through spectral modelling. We also show that the duration of the very early light curve, before the luminosity has reached half of its maximal value, does not correlate with the light curve shape or stretch used to standardise SNe Ia in cosmological applications. This has implications for the cosmological fitting of SN Ia light curves.

The AM CVn systems are a class of He-rich, post-period minimum, semi-detached, ultra-compact binaries. Their long-term light curves have been poorly understood due to the few systems known and the long (hundreds of days) recurrence times between outbursts. We present combined photometric light curves from the LINEAR, CRTS, and PTF synoptic surveys to study the photometric variability of these systems over an almost 10 yr period. These light curves provide a much clearer picture of the outburst phenomena that these systems undergo. We characterize the photometric behavior of most known outbursting AM CVn systems and establish a relation between their outburst properties and the systems' orbital periods. We also explore why some systems have only shown a single outburst so far and expand the previously accepted phenomenological states of AM CVn systems. We conclude that the outbursts of these systems show evolution with respect to the orbital period, which can likely be attributed to the decreasing mass transfer rate with increasing period. Finally, we consider the number of AM CVn systems that should be present in modeled synoptic surveys.

Millisecond duration bright radio pulses at 1.4-GHz with high dispersion measures (DM) were reported by Lorimer et al., Keane et al., and Thornton et al. Their all-sky rate is $\approx 10^4$/day above $\sim$1 Jy. Related events are "Perytons" -- similar pulsed, dispersed sources, but most certainly local. Suggested models of fast radio bursts (FRBs) can originate in the Earth's atmosphere, in stellar coronae, in other galaxies, and even at cosmological distances. Using physically motivated assumptions combined with observed properties, we explore these models. In our analysis, we focus on the Lorimer event: a 30 Jy, 5-ms duration burst with DM$=$ 375 cm$^{-3}$ pc, exhibiting a steep frequency-dependent pulse width (the {\it Sparker}). To be complete, we drop the assumption that high DMs are produced by plasma propagation and assume that the source produces pulses with frequency-dependent arrival time ("chirped signals"). Within this framework we explore a scenario in which Perytons, the {\it Sparker}, and the FRBs are all atmospheric phenomenon occurring at different heights. This model is {\it ad hoc} in that we cannot explain why Perytons at higher altitudes show greater DMs or exhibit narrower pulses. Nonetheless, we argue the {\it Sparker} may be a Peryton. We end with two remarks. First, the detection of a single FRB by an interferometer with a kilometer (or longer) baseline will prove that FRBs are of extra-terrestrial origin. Second, we urge astronomers to pursue observations and understanding of Perytons since they form (at least) a formidable foreground for the FRBs.

We present imaging and spectroscopy of a hydrogen-poor superluminous supernova (SLSN) discovered by the intermediate Palomar Transient Factory: iPTF13ajg. At a redshift of z=0.7403, derived from narrow absorption lines, iPTF13ajg peaked at an absolute magnitude M(u,AB)=-22.5, one of the most luminous supernovae to date. The uBgRiz light curves, obtained with the P48, P60, NOT, DCT, and Keck telescopes, and the nine-epoch spectral sequence secured with the Keck and the VLT (covering 3 rest-frame months), are tied together photometrically to provide an estimate of the flux evolution as a function of time and wavelength. The observed bolometric peak luminosity of iPTF13ajg is 3.2x10^44 erg/s, while the estimated total radiated energy is 1.3x10^51 erg. We detect narrow absorption lines of Mg I, Mg II, and Fe II, associated with the cold interstellar medium in the host galaxy, at two different epochs with X-shooter at the VLT. From Voigt-profile fitting, we derive the column densities log N(Mg I)=11.94+-0.06, log N(Mg II)=14.7+-0.3, and log N(Fe II)=14.25+-0.10. These column densities, as well as the Mg I and Mg II equivalent widths of a sample of hydrogen-poor SLSNe taken from the literature, are at the low end of those derived for gamma-ray bursts (GRBs), whose progenitors are also thought to be massive stars. This suggests that the environments of SLSNe and GRBs are different. From the nondetection of Fe II fine-structure absorption lines, we derive a strict lower limit on the distance between the supernova and the narrow-line absorbing gas of 50 pc. No host-galaxy emission lines are detected, leading to an upper limit on the unobscured star-formation rate of SFR([OII])<0.07 Msun/yr. Late-time imaging shows the host galaxy of iPTF13ajg to be faint, with g(AB)~27.0 and R(AB)>=26.0 mag, which roughly corresponds to M(B,Vega) >~ -17.7 mag. [abridged]

We present the results of a Palomar Transient Factory (PTF) archival search for blue transients which lie in the magnitude range between "normal" core-collapse and superluminous supernovae (i.e. with $-21\,{\leq}M_{R\,(peak)}\,{\leq}-19$). Of the six events found after excluding all interacting Type~IIn and Ia-CSM supernovae, three (PTF09ge, 09axc and 09djl) are coincident with the centers of their hosts, one (10iam) is offset from the center, and for two (10nuj and 11glr) a precise offset can not be determined. All the central events have similar rise times to the He-rich tidal disruption candidate PS1-10jh, and the event with the best-sampled light curve also has similar colors and power-law decay. Spectroscopically, PTF09ge is He-rich, while PTF09axc and 09djl display broad hydrogen features around peak magnitude. All three central events are in low star-formation hosts, two of which are E+A galaxies. Our spectrum of the host of PS1-10jh displays similar properties. PTF10iam, the one offset event, is different photometrically and spectroscopically from the central events and its host displays a higher star formation rate. Finding no obvious evidence for ongoing galactic nuclei activity or recent star formation, we conclude that the three central transients likely arise from the tidal disruption of a star by a super-massive black hole. We compare the spectra of these events to tidal disruption candidates from the literature and find that all of these objects can be unified on a continuous scale of spectral properties. The accumulated evidence of this expanded sample strongly supports a tidal disruption origin for this class of nuclear transients.

Almost every known low-luminosity Milky Way dwarf spheroidal (dSph) satellite galaxy contains at least one RR Lyrae star. Assuming that a fraction of distant (60 < d_{helio} < 100 kpc) Galactic halo RR Lyrae stars are members of yet to be discovered low-luminosity dSph galaxies, we perform a {\em guided} search for these low-luminosity dSph galaxies. In order to detect the presence of dSph galaxies, we combine stars selected from more than 123 sightlines centered on RR Lyrae stars identified by the Palomar Transient Factory. We find that this method is sensitive enough to detect the presence of Segue 1-like galaxies (M_V= -1.5^{+0.6}_{-0.8}, r_h=30 pc) even if only ~20 sightlines were occupied by such dSph galaxies. Yet, when our method is applied to the SDSS DR10 imaging catalog, no signal is detected. An application of our method to sightlines occupied by pairs of close (<200 pc) horizontal branch stars, also did not yield a detection. Thus, we place upper limits on the number of low-luminosity dSph galaxies with half-light radii from 30 pc to 120 pc, and in the probed volume of the halo. Stronger constraints on the luminosity function may be obtained by applying our method to sightlines centered on RR Lyrae stars selected from the Pan-STARRS1 survey, and eventually, from LSST. In the Appendix, we present spectroscopic observations of an RRab star in the Bo\"{o}tes 3 dSph and a light curve of an RRab star near the Bo\"{o}tes 2 dSph.

The explosive fate of massive stripped Wolf-Rayet (W-R) stars is a key open question in stellar physics. An appealing option is that hydrogen-deficient W-R stars are the progenitors of some H-poor supernova (SN) explosions of Types IIb, Ib, and Ic. A blue object, having luminosity and colors consistent with those of some W-R stars, has been recently identified at the location of a SN~Ib in pre-explosion images but has not yet been conclusively determined to have been the progenitor. Similar previous works have so far only resulted in nondetections. Comparison of early photometric observations of Type Ic supernovae with theoretical models suggests that the progenitor stars had radii <10^12 cm, as expected for some W-R stars. However, the hallmark signature of W-R stars, their emission-line spectra, cannot be probed by such studies. Here, we report the detection of strong emission lines in an early-time spectrum of SN 2013cu (iPTF13ast; Type IIb) obtained ~15.5 hr after explosion ("flash spectroscopy"). We identify W-R-like wind signatures suggesting a progenitor of the WN(h) subclass. The extent of this dense wind may indicate increased mass loss from the progenitor shortly prior to its explosion, consistent with recent theoretical predictions.

We present optical photometric and spectroscopic observations of the Type Ibn (SN 2006jc-like) supernova iPTF13beo. Detected by the intermediate Palomar Transient Factory ~3 hours after the estimated first light, iPTF13beo is the youngest and the most distant (~430 Mpc) Type Ibn event ever observed. The iPTF13beo light curve is consistent with light curves of other Type Ibn SNe and with light curves of fast Type Ic events, but with a slightly faster rise-time of two days. In addition, the iPTF13beo R-band light curve exhibits a double-peak structure separated by ~9 days, not observed before in any Type Ibn SN. A low-resolution spectrum taken during the iPTF13beo rising stage is featureless, while a late-time spectrum obtained during the declining stage exhibits narrow and intermediate-width He I and Si II features with FWHM ~ 2000-5000 km/s and is remarkably similar to the prototypical SN Ibn 2006jc spectrum. We suggest that our observations support a model of a massive star exploding in a dense He-rich circumstellar medium (CSM). A shock breakout in a CSM model requires an eruption releasing a total mass of ~0.1 Msun over a time scale of couple of weeks prior to the SN explosion.

Since the discovery of the unusual prototype SN 2002cx, the eponymous class of low-velocity, hydrogen-poor supernovae has grown to include at most another two dozen members identified from several heterogeneous surveys, in some cases ambiguously. Here we present the results of a systematic study of 1077 hydrogen-poor supernovae discovered by the Palomar Transient Factory, leading to nine new members of this peculiar class. Moreover we find there are two distinct subclasses based on their spectroscopic, photometric, and host galaxy properties: The "SN 2002cx-like" supernovae tend to be in later-type or more irregular hosts, have more varied and generally dimmer luminosities, have longer rise times, and lack a Ti II trough when compared to the "SN 2002es-like" supernovae. None of our objects show helium, and we counter a previous claim of two such events. We also find that these transients comprise 5.6+17-3.7% (90% confidence) of all SNe Ia, lower compared to earlier estimates. Combining our objects with the literature sample, we propose that these subclasses have two distinct physical origins.

The Palomar Transient Factory (PTF) is a multi-epochal robotic survey of the northern sky that acquires data for the scientific study of transient and variable astrophysical phenomena. The camera and telescope provide for wide-field imaging in optical bands. In the five years of operation since first light on December 13, 2008, images taken with Mould-R and SDSS-g' camera filters have been routinely acquired on a nightly basis (weather permitting), and two different H-alpha filters were installed in May 2011 (656 nm and 663 nm). The PTF image-processing and data-archival program at the Infrared Processing and Analysis Center (IPAC) is tailored to receive and reduce the data, and, from it, generate and preserve astrometrically and photometrically calibrated images, extracted source catalogs, and coadded reference images. Relational databases have been deployed to track these products in operations and the data archive. The fully automated system has benefited by lessons learned from past IPAC projects and comprises advantageous features that are potentially incorporable into other ground-based observatories. Both off-the-shelf and in-house software have been utilized for economy and rapid development. The PTF data archive is curated by the NASA/IPAC Infrared Science Archive (IRSA). A state-of-the-art custom web interface has been deployed for downloading the raw images, processed images, and source catalogs from IRSA. Access to PTF data products is currently limited to an initial public data release (M81, M44, M42, SDSS Stripe 82, and the Kepler Survey Field). It is the intent of the PTF collaboration to release the full PTF data archive when sufficient funding becomes available.

We present optical photometry and spectroscopy of the broad-lined Type Ic supernova (SNIc-BL) PTF10qts, which was discovered as part of the Palomar Transient Factory. The supernova was located in a dwarf galaxy of magnitude $r=21.1$ at a redshift $z=0.0907$. We find that the $R$-band light curve is a poor proxy for bolometric data and use photometric and spectroscopic data to construct and constrain the bolometric light curve. The derived bolometric magnitude at maximum light is $M_{\rm bol} = -18.51\pm0.2$ mag, comparable to that of SN 1998bw ($M_{\rm bol} = -18.7$ mag) which was associated with a gamma-ray burst (GRB). PTF10qts is one of the most luminous SNIc-BL observed without an accompanying GRB. We estimate the physical parameters of the explosion using data from our programme of follow-up observations, finding that it produced a larger mass of radioactive nickel compared to other SNeIc-BL with similar inferred ejecta masses and kinetic energies. The progenitor of the event was likely a $\sim20$M$_{\odot}$ star.

Interaction of supernova (SN) ejecta with the optically thick circumstellar medium (CSM) of a progenitor star can result in a bright, long-lived shock breakout event. Candidates for such SNe include Type IIn and superluminous SNe. If some of these SNe are powered by interaction, then there should be a relation between their peak luminosity, bolometric light-curve rise time, and shock-breakout velocity. Given that the shock velocity during shock breakout is not measured, we expect a correlation, with a significant spread, between the rise time and the peak luminosity of these SNe. Here, we present a sample of 15 SNe IIn for which we have good constraints on their rise time and peak luminosity from observations obtained using the Palomar Transient Factory. We report on a possible correlation between the R-band rise time and peak luminosity of these SNe, with a false-alarm probability of 3%. Assuming that these SNe are powered by interaction, combining these observables and theory allows us to deduce lower limits on the shock-breakout velocity. The lower limits on the shock velocity we find are consistent with what is expected for SNe (i.e., ~10^4 km/s). This supports the suggestion that the early-time light curves of SNe IIn are caused by shock breakout in a dense CSM. We note that such a correlation can arise from other physical mechanisms. Performing such a test on other classes of SNe (e.g., superluminous SNe) can be used to rule out the interaction model for a class of events.

We report on the discovery of SN2014J in the nearby galaxy M82. Given its proximity, it offers the best opportunity to date to study a thermonuclear supernova over a wide range of the electromagnetic spectrum. The first set of optical, near-IR and mid-IR observations of SN2014J, orchestrated by the intermediate Palomar Transient Factory (iPTF), show that SN2014J is a spectroscopically normal Type Ia supernova, albeit exhibiting high-velocity features in its spectrum and heavily reddened by dust in the host galaxy. Our earliest detections start just hours after the fitted time of explosion. We use high-resolution optical spectroscopy to analyze the dense intervening material and do not detect any evolution in the resolved absorption features during the lightcurve rise. Similarly to other highly reddened Type Ia supernovae, a low value of total-to-selective extinction, Rv < 2, provides the best match to our observations. We also study pre-explosion optical and near-IR images from HST with special emphasis on the sources nearest to the SN location.

We report on PTF1 J191905.19+481506.2, a newly discovered, partially eclipsing, outbursting AM CVn system found in the Palomar Transient Factory synoptic survey. This is only the second known eclipsing AM CVn system. We use high-speed photometric observations and phase-resolved spectroscopy to establish an orbital period of 22.4559(3) min. We also present a long-term light curve and report on the normal and super-outbursts regularly seen in this system, including a super-outburst recurrence time of 36.8(4) d. We use the presence of the eclipse to place upper and lower limits on the inclination of the system and discuss the number of known eclipsing AM CVn systems versus what would be expected.

We have identified spectral features in the late-time X-ray afterglow of the unusually long, slow-decaying GRB 130925A using NuSTAR, Swift-XRT, and Chandra. A spectral component in addition to an absorbed power-law is required at $>4\sigma$ significance, and its spectral shape varies between two observation epochs at $2\times10^5$ and $10^6$ seconds after the burst. Several models can fit this additional component, each with very different physical implications. A broad, resolved Gaussian absorption feature of several keV width improves the fit, but it is poorly constrained in the second epoch. An additive black body or second power-law component provide better fits. Both are challenging to interpret: the blackbody radius is near the scale of a compact remnant ($10^8$ cm), while the second powerlaw component requires an unobserved high-energy cutoff in order to be consistent with the non-detection by Fermi-LAT.

The optical and near-IR emission from some classes of supernovae (SNe), including Type IIn and possibly some super-luminous SNe, is likely powered by a collision between the SN ejecta and dense circumstellar material (CSM). We argue that for a range of CSM masses and their radii, a collisionless shock can form, allowing for efficient cosmic-ray (CR) acceleration. We show that pp collisions between these newly accelerated CRs and the CSM leads to not only gamma rays but also secondary electrons and positrons that radiate synchrotron photons in the high-frequency radio bands. Our estimates imply that various facilities including the Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA) may observe such SNe at Gpc distances by followup observations in months-to-years, although the detectability strongly depends on the CSM density as well as observed frequency. Detecting this signal would give us a unique probe of CR acceleration at early times, and even non-detections can put interesting limits on the possibility of CR ion acceleration. Following our previous work, we also show that GeV gamma rays can escape from the system without severe attenuation, encouraging point-source and stacking analyses with Fermi. We provide recipes for diagnosing interaction-powered SN scenario with multi-messenger (neutrino and gamma-ray) observations.

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova discovered by the Palomar Transient Factory. PTF11qcj is located at a distance of dL ~ 124 Mpc. Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated supernova 1998bw (L_{5GHz} ~ 10^{29} erg/s/Hz). PTF11qcj is also detected in X-rays with the Chandra observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the supernova interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of ~10^{-4} Msun/yr x (v_w/1000 km/s), and a velocity of ~(0.3-0.5)c for the fastest moving ejecta (at ~10d after explosion). However, these estimates are derived assuming the simplest model of supernova ejecta interacting with a smooth circumstellar material characterized by radial power-law density profile, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio light curve shows deviations from such a simple model, as well as a re-brightening at late times. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). An IR light echo from pre-existing dust is in agreement with our infrared data. Our analysis of pre-explosion data from the Palomar Transient Factory suggests that a precursor eruption of absolute magnitude M_r ~ -13 mag may have occurred ~ 2.5 yr prior to the supernova explosion. Based on our panchromatic follow-up campaign, we conclude that PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

There is a growing number of supernovae (SNe), mainly of Type IIn, which present an outburst prior to their presumably final explosion. These precursors may affect the SN display, and are likely related to some poorly charted phenomena in the final stages of stellar evolution. Here we present a sample of 16 SNe IIn for which we have Palomar Transient Factory (PTF) observations obtained prior to the SN explosion. By coadding these images taken prior to the explosion in time bins, we search for precursor events. We find five Type IIn SNe that likely have at least one possible precursor event, three of which are reported here for the first time. For each SN we calculate the control time. Based on this analysis we find that precursor events among SNe IIn are common: at the one-sided 99% confidence level, more than 50% of SNe IIn have at least one pre-explosion outburst that is brighter than absolute magnitude -14, taking place up to 1/3 yr prior to the SN explosion. The average rate of such precursor events during the year prior to the SN explosion is likely larger than one per year, and fainter precursors are possibly even more common. We also find possible correlations between the integrated luminosity of the precursor, and the SN total radiated energy, peak luminosity, and rise time. These correlations are expected if the precursors are mass-ejection events, and the early-time light curve of these SNe is powered by interaction of the SN shock and ejecta with optically thick circumstellar material.

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

Many photometric time-domain surveys are driven by specific goals, such as searches for supernovae or transiting exoplanets, which set the cadence with which fields are re-imaged. In the case of the Palomar Transient Factory (PTF), several sub-surveys are conducted in parallel, leading to non-uniform sampling over its $\sim$$20,000 \mathrm{deg}^2$ footprint. While the median $7.26 \mathrm{deg}^2$ PTF field has been imaged $\sim$40 times in \textit{R}-band, $\sim$$2300 \mathrm{deg}^2$ have been observed $>$100 times. We use PTF data to study the trade-off between searching for microlensing events in a survey whose footprint is much larger than that of typical microlensing searches, but with far-from-optimal time sampling. To examine the probability that microlensing events can be recovered in these data, we test statistics used on uniformly sampled data to identify variables and transients. We find that the von Neumann ratio performs best for identifying simulated microlensing events in our data. We develop a selection method using this statistic and apply it to data from fields with $>$10 $R$-band observations, $1.1\times10^9$ light curves, uncovering three candidate microlensing events. We lack simultaneous, multi-color photometry to confirm these as microlensing events. However, their number is consistent with predictions for the event rate in the PTF footprint over the survey's three years of operations, as estimated from near-field microlensing models. This work can help constrain all-sky event rate predictions and tests microlensing signal recovery in large data sets, which will be useful to future time-domain surveys, such as that planned with the Large Synoptic Survey Telescope.

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.

We analyze the planetary microlensing event MOA-2010-BLG-328. The best fit yields host and planetary masses of Mh = 0.11+/-0.01 M_{sun} and Mp = 9.2+/-2.2M_Earth, corresponding to a very late M dwarf and sub-Neptune-mass planet, respectively. The system lies at DL = 0.81 +/- 0.10 kpc with projected separation r = 0.92 +/- 0.16 AU. Because of the host's a-priori-unlikely close distance, as well as the unusual nature of the system, we consider the possibility that the microlens parallax signal, which determines the host mass and distance, is actually due to xallarap (source orbital motion) that is being misinterpreted as parallax. We show a result that favors the parallax solution, even given its close host distance. We show that future high-resolution astrometric measurements could decisively resolve the remaining ambiguity of these solutions.

We report the discovery of the optical afterglow of the gamma-ray burst (GRB) 130702A, identified upon searching 71 square degrees surrounding the Fermi Gamma-ray Burst Monitor (GBM) localization. Discovered and characterized by the intermediate Palomar Transient Factory (iPTF), iPTF13bxl is the first afterglow discovered solely based on a GBM localization. Real-time image subtraction, machine learning, human vetting, and rapid response multi-wavelength follow-up enabled us to quickly narrow a list of 27,004 optical transient candidates to a single afterglow-like source. Detection of a new, fading X-ray source by Swift and a radio counterpart by CARMA and the VLA confirmed the association between iPTF13bxl and GRB 130702A. Spectroscopy with the Magellan and Palomar 200-inch telescopes showed the afterglow to be at a redshift of z=0.145, placing GRB 130702A among the lowest redshift GRBs detected to date. The prompt gamma-ray energy release and afterglow luminosity are intermediate between typical cosmological GRBs and nearby sub-luminous events such as GRB 980425 and GRB 060218. The bright afterglow and emerging supernova offer an opportunity for extensive panchromatic follow-up. Our discovery of iPTF13bxl demonstrates the first observational proof-of-principle for ~10 Fermi-iPTF localizations annually. Furthermore, it represents an important step towards overcoming the challenges inherent in uncovering faint optical counterparts to comparably localized gravitational wave events in the Advanced LIGO and Virgo era.

We present our observations of SN 2010mb, a Type Ic SN lacking spectroscopic signatures of H and He. SN 2010mb has a slowly-declining light curve ($\sim600\,$days) that cannot be powered by $^{56}$Ni/$^{56}$Co radioactivity, the common energy source for Type Ic SNe. We detect signatures of interaction with hydrogen-free CSM including a blue quasi-continuum and, uniquely, narrow oxygen emission lines that require high densities ($\sim10^9$cm$^{-3}$). From the observed spectra and light curve we estimate that the amount of material involved in the interaction was $\sim3$M$_{\odot}$. Our observations are in agreement with models of pulsational pair-instability SNe described in the literature.

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our finely sampled light curves show a rise of 0.8mag within 2.5hr. Power-law fits [f(t)\sim(t-t_0)^n] to these data constrain the explosion date to within one day. We cannot rule out the expected quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic, and its high ejecta velocities (~30,000km/s four days after explosion) are closer to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) than to the observed velocities in normal Type Ic SNe. The high velocities are sustained through the SN early evolution, and are similar to those of GRB-SNe when the SN reach peak magnitude. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25-35 solar mass and a large (5-10E51 erg) kinetic energy, close to the regime of GRB-SN properties. The host-galaxy characteristics are consistent with GRB-SN hosts, and not with normal SN Ic hosts as well, yet this SN does not show the broad lines over extended periods of time that are typical of broad-line Type Ic SNe.

We present Spitzer observations at 3.6 and 4.5um of the supernova SN2011dh (PTF11eon) in M51 from 18 days to 625 days after explosion. The mid-infrared emission peaks at 24 days after explosion at a few times 10^7 Lsun, and decays more slowly than the visible-light bolometric luminosity. The infrared color temperature cools for the first 90 days and then is constant. Simple numerical models of a thermal echo can qualitatively reproduce the early behavior. At late times, the mid-IR light curve cannot be explained by a simple thermal echo model, suggesting additional dust heating or line emission mechanisms. We also propose that thermal echoes can serve as effective probes to uncover supernovae in heavily obscured environments, and speculate that under the right conditions, integrating the early epoch of the mid-infrared light curve may constrain the total energy in the shock breakout flash.

The intermediate Palomar Transient Factory reports our discovery of a young supernova, iPTF13bvn, in the nearby galaxy, NGC5806 (22.5Mpc). Our spectral sequence in the optical and infrared suggests a likely Type Ib classification. We identify a single, blue progenitor candidate in deep pre-explosion imaging within a 2{\sigma} error circle of 80 mas (8.7 pc). The candidate has a MB luminosity of -5.2 +/- 0.4 mag and a B-I color of 0.1+/-0.3 mag. If confirmed by future observations, this would be the first direct detection for a progenitor of a Type Ib. Fitting a power law to the early light curve, we find an extrapolated explosion date around 1.1 days before our first detection. We see no evidence of shock cooling. The pre-explosion detection limits constrain the radius of the progenitor to be smaller than a few solar radii. iPTF13bvn is also detected in cm and mm-wavelengths. Fitting a synchrotron self-absorption model to our radio data, we find a mass loading parameter of 1.3*10^12 g/cm. Assuming a wind velocity of 10^3km/s, we derive a progenitor mass loss rate of 3*10^-5Msun/yr. Our observations, taken as a whole, are consistent with a Wolf Rayet progenitor of the supernova iPTF13bvn.

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material (CSM) surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the ESO Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (~20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to non-blueshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption, and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B-V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests there are at least two distinct populations of normal, cosmologically useful SNe Ia.

We report positions, velocities and metallicities of 50 ab-type RR Lyrae (RRab) stars observed in the vicinity of the Orphan stellar stream. Using about 30 RRab stars classified as being likely members of the Orphan stream, we study the metallicity and the spatial extent of the stream. We find that RRab stars in the Orphan stream have a wide range of metallicities, from -1.5 dex to -2.7 dex. The average metallicity of the stream is -2.1 dex, identical to the value obtained by Newberg et al. (2010) using blue horizontal branch stars. We find that the most distant parts of the stream (40-50 kpc from the Sun) are about 0.3 dex more metal-poor than the closer parts (within ~30 kpc), suggesting a possible metallicity gradient along the stream's length. We have extended the previous studies and have mapped the stream up to 55 kpc from the Sun. Even after a careful search, we did not identify any more distant RRab stars that could plausibly be members of the Orphan stream. If confirmed with other tracers, this result would indicate a detection of the end of the leading arm of the stream. We have compared the distances of Orphan stream RRab stars with the best-fit orbits obtained by Newberg et al. (2010). We find that model 6 of Newberg et al. (2010) cannot explain the distances of the most remote Orphan stream RRab stars, and conclude that the best fit to distances of Orphan stream RRab stars and to the local circular velocity is provided by potentials where the total mass of the Galaxy within 60 kpc is M_{60}~2.7x10^{11} Msun, or about 60% of the mass found by previous studies. More extensive modelling that would consider non-spherical potentials and the possibility of misalignment between the stream and the orbit, is highly encouraged.

(Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.

Only a few cases of type Ic supernovae (SNe) with high-velocity ejecta have been discovered and studied. Here we present our analysis of radio and X-ray observations of a Type Ic SN, PTF12gzk. The radio emission rapidly declined less than 10 days after explosion, suggesting SN ejecta expanding at high velocity (~0.3c). The radio data also indicate that the density of the circumstellar material (CSM) around the supernova is lower by a factor of ~10 than the CSM around normal Type Ic SNe. Our observations of this rapidly declining radio SN at a distance of 58 Mpc demonstrates the potential to detect many additional radio SNe, given the new capabilities of the VLA (improved sensitivity and dynamic scheduling), that are currently missed, leading to a biased view of radio SNe Ic. Early optical discovery followed by rapid radio observations would provide a full description of the ejecta velocity distribution and CSM densities around stripped massive star explosions, as well as strong clues about the nature of their progenitor stars.

We report observations of a possible young transiting planet orbiting a previously known weak-lined T-Tauri star in the 7-10 Myr old Orion-OB1a/25-Ori region. The candidate was found as part of the Palomar Transient Factory (PTF) Orion project. It has a photometric transit period of 0.448413 +- 0.000040 days, and appears in both 2009 and 2010 PTF data. Follow-up low-precision radial velocity (RV) observations and adaptive optics imaging suggest that the star is not an eclipsing binary, and that it is unlikely that a background source is blended with the target and mimicking the observed transit. RV observations with the Hobby-Eberly and Keck telescopes yield an RV that has the same period as the photometric event, but is offset in phase from the transit center by approximately -0.22 periods. The amplitude (half range) of the RV variations is 2.4 km/s and is comparable with the expected RV amplitude that stellar spots could induce. The RV curve is likely dominated by stellar spot modulation and provides an upper limit to the projected companion mass of M_p sin i_orb < 4.8 +- 1.2 M_Jup; when combined with the orbital inclination, i orb, of the candidate planet from modeling of the transit light curve, we find an upper limit on the mass of the planetary candidate of M_p < 5.5 +- 1.4 M_Jup. This limit implies that the planet is orbiting close to, if not inside, its Roche limiting orbital radius, so that it may be undergoing active mass loss and evaporation.

Cometary activity in main-belt asteroids probes the ice content of these objects and provides clues to the history of volatiles in the inner solar system. We search the Palomar Transient Factory (PTF) survey to derive upper limits on the population size of active main-belt comets (MBCs). From data collected March 2009 through July 2012, we extracted 2 million observations of 220 thousand known main-belt objects (40% of the known population, down to 1-km diameter) and discovered 626 new objects in multi-night linked detections. We formally quantify the extendedness of a small-body observation, account for systematic variation in this metric (e.g., due to on-sky motion) and evaluate this method's robustness in identifying cometary activity using observations of 115 comets, including two known candidate MBCs and six newly-discovered non-main-belt comets (two of which were originally designated as asteroids by other surveys). We demonstrate a 66% detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100% detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidate MBCs P/2010 R2 (La Sagra) and P/2006 VW139. Using a log-constant prior, we infer 95% confidence upper limits of 33 and 22 active MBCs (per million main-belt asteroids down to 1-km diameter), for detection efficiencies of 66% and 100%, respectively. In a follow-up to this morphological search, we will perform a photometric (disk-integrated brightening) search for MBCs.

We report the discovery by the Palomar Transient Factory (PTF) of the transient source PTF11agg, which is distinguished by three primary characteristics: (1) bright, rapidly fading optical transient emission; (2) a faint, blue quiescent optical counterpart; and (3) an associated year-long, scintillating radio transient. We argue that these observed properties are inconsistent with any known class of Galactic transients, and instead suggest a cosmological origin. The detection of incoherent radio emission at such distances implies a large emitting region, from which we infer the presence of relativistic ejecta. The observed properties are all consistent with the population of long-duration gamma-ray bursts (GRBs), marking the first time such an outburst has been discovered in the distant universe independent of a high-energy trigger. We searched for possible high-energy counterparts to PTF11agg, but found no evidence for associated prompt emission. We therefore consider three possible scenarios to account for a GRB-like afterglow without a high-energy counterpart: an "untriggered" GRB (lack of satellite coverage), an "orphan" afterglow (viewing-angle effects), and a "dirty fireball" (suppressed high-energy emission). The observed optical and radio light curves appear inconsistent with even the most basic predictions for off-axis afterglow models. The simplest explanation, then, is that PTF11agg is a normal, on-axis long-duration GRB for which the associated high-energy emission was simply missed. However, we have calculated the likelihood of such a serendipitous discovery by PTF and find that it is quite small (~ 2.6%). While not definitive, we nonetheless speculate that PTF11agg may represent a new, more common (> 4 times the on-axis GRB rate at 90% confidence) class of relativistic outbursts lacking associated high-energy emission.

Owing to their utility for measurements of cosmic acceleration, Type Ia supernovae (SNe) are perhaps the best-studied class of SNe, yet the progenitor systems of these explosions largely remain a mystery. A rare subclass of SNe Ia show evidence of strong interaction with their circumstellar medium (CSM), and in particular, a hydrogen-rich CSM; we refer to them as SNe Ia-CSM. In the first systematic search for such systems, we have identified 16 SNe Ia-CSM, and here we present new spectra of 13 of them. Six SNe Ia-CSM have been well-studied previously, three were previously known but are analyzed in-depth for the first time here, and seven are new discoveries from the Palomar Transient Factory. The spectra of all SNe Ia-CSM are dominated by H{\alpha} emission (with widths of ~2000 km/s) and exhibit large H{\alpha}/H{\beta} intensity ratios (perhaps due to collisional excitation of hydrogen via the SN ejecta overtaking slower-moving CSM shells); moreover, they have an almost complete lack of He I emission. They also show possible evidence of dust formation through a decrease in the red wing of H{\alpha} 75-100 d past maximum brightness, and nearly all SNe Ia-CSM exhibit strong Na I D absorption from the host galaxy. The absolute magnitudes (uncorrected for host-galaxy extinction) of SNe Ia-CSM are found to be -21.3 <= M_R <= -19 mag, and they also seem to show ultraviolet emission at early times and strong infrared emission at late times (but no detected radio or X-ray emission). Finally, the host galaxies of SNe Ia-CSM are all late-type spirals similar to the Milky Way, or dwarf irregulars like the Large Magellanic Cloud, which implies that these objects come from a relatively young stellar population. This work represents the most detailed analysis of the SN Ia-CSM class to date.

We report on an analysis of the Extended Chandra Deep Field South (E-CDFS) region using archival data from the Very Large Array, with the goal of studying radio variability and transients at the sub-mJy level. The 49 epochs of E-CDFS observations at 1.4 GHz sample timescales from one day to 3 months. We find that only a fraction (1%) of unresolved radio sources above 40 uJy are variable at the 4-sigma level. There is no evidence that the fractional variability changes along with the known transition of radio source populations below one milliJansky. Optical identifications of the sources show that the variable radio emission is associated with the central regions of an active galactic nucleus or a star-forming galaxy. After a detailed comparison of the efficacy of various source-finding algorithms, we use the best to carry out a transient search. No transients were found. This implies that the areal density of transients with peak flux density greater than 0.21 mJy is less than 0.37 deg^-2 (at a confidence level of 95%). This result is approximately an order of magnitude below the transient rate measured at 5 GHz by Bower et al. (2007) but it is consistent with more recent upper limits from Frail et al. (2012). Our findings suggest that the radio sky at 1.4 GHz is relatively quiet. For multi-wavelength transient searches, such as the electromagnetic counterparts to gravitational waves, this frequency may be optimal for reducing the high background of false positives.

The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely.

We report the discovery of R Coronae Borealis (RCB) stars in the Andromeda galaxy (M31) using the Palomar Transient Factory (PTF). RCB stars are rare hydrogen-deficient, carbon-rich supergiant variables, most likely the merger products of two white dwarfs. These new RCBs, including two confirmed ones and two candidates, are the first to be found beyond the Milky Way and the Magellanic Clouds. All of M31 RCBs showed >1.5 mag irregular declines over timescales of weeks to months. Due to the limiting magnitude of our data (R~21-22 mag), these RCB stars have R~19.5 to 20.5 mag at maximum light, corresponding to M_R= -4 to -5, making them some of the most luminous RCBs known. Spectra of two objects show that they are warm RCBs, similar to the Milky Way RCBs RY Sgr and V854 Cen. We consider these results, derived from a pilot study of M31 variables, as an important proof-of-concept for the study of rare bright variables in nearby galaxies with the PTF or other synoptic surveys.

Some supernovae (SNe) show evidence for mass-loss events taking place prior to their explosions. Measuring their pre-outburst mass-loss rates provide essential information regarding the mechanisms that are responsible for these events. Here we present XMM-Newton and Swift X-ray observations taken after the latest, and presumably the final, outburst of SN 2009ip. We use these observations as well as new near infra-red and visible light spectra, and published radio and visible light observations to put six independent order-of-magnitude constrains on the mass-loss rate of the SN progenitor prior to the explosion. Our methods utilize: the X-ray luminosity, the bound-free absorption, the H alpha luminosity, the SN rise-time, free-free absorption, and the bolometric luminosity of the outburst detected prior to the explosion. Assuming spherical mass-loss with a wind density profile, we estimate that the effective mass-loss rate from the progenitor was between 10^-3 to 10^-2 solar masses per year, over a few years prior to the explosion, with a velocity of ~1000 km/s. This mass-loss rate corresponds to a total circum stellar matter mass of ~0.04 solar masses, within 6x10^15 cm of the SN. We note that the mass-loss rate estimate based on the H alpha luminosity is higher by an order of magnitude. This can be explained if the narrow line H alpha component is generated at radii larger than the shock radius, or if the CSM has an aspherical geometry. We discuss simple geometries which are consistent with our results.

We present first results from a new, multiyear, time domain survey of young stars in the North America Nebula complex using the Palomar Transient Factory. Our survey is providing an unprecedented view of aperiodic variability in young stars on timescales of days to years. The analyzed sample covers R_PTF = 13.5-18 and spans a range of mid-infrared color, with larger-amplitude optical variables (exceeding 0.4 mag root-mean-squared) more likely to have mid-infrared evidence for circumstellar material. This paper characterizes infrared excess stars with distinct bursts above or fades below a baseline of lower-level variability, identifying 41 examples. The light curves exhibit a remarkable diversity of amplitudes, timescales, and morphologies, with a continuum of behaviors that can not be classified into distinct groups. Among the bursters, we identify three particularly promising sources that may represent theoretically predicted short-timescale accretion instabilities. Finally, we find that fading behavior is approximately twice as common as bursting behavior on timescales of days to years, although the bursting and fading duty cycle for individual objects often varies from year to year.

Light from astronomical objects must travel through the earth's turbulent atmosphere before it can be imaged by ground-based telescopes. To enable direct imaging at maximum theoretical angular resolution, advanced techniques such as those employed by the Robo-AO adaptive-optics system must be used.

Various lines of evidence suggest that very massive stars experience extreme mass-loss episodes shortly before they explode as a supernova. Interestingly, several models predict such pre-explosion outbursts. Establishing a causal connection between these mass-loss episodes and the final supernova explosion will provide a novel way to study pre-supernova massive-star evolution. Here we report on observations of a remarkable mass-loss event detected 40 days prior to the explosion of the Type IIn supernova SN 2010mc (PTF 10tel). Our photometric and spectroscopic data suggest that this event is a result of an energetic outburst, radiating at least 6x10^47 erg of energy, and releasing about 0.01 Solar mass at typical velocities of 2000 km/s. We show that the temporal proximity of the mass-loss outburst and the supernova explosion implies a causal connection between them. Moreover, we find that the outburst luminosity and velocity are consistent with the predictions of the wave-driven pulsation model and disfavor alternative suggestions.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

Radiative transfer studies of Type Ia supernovae (SNe Ia) hold the promise of constraining both the time-dependent density profile of the SN ejecta and its stratification by element abundance which, in turn, may discriminate between different explosion mechanisms and progenitor classes. Here we present a detailed analysis of Hubble Space Telescope ultraviolet (UV) and ground-based optical spectra and light curves of the SN Ia SN 2010jn (PTF10ygu). SN 2010jn was discovered by the Palomar Transient Factory (PTF) 15 days before maximum light, allowing us to secure a time-series of four UV spectra at epochs from -11 to +5 days relative to B-band maximum. The photospheric UV spectra are excellent diagnostics of the iron-group abundances in the outer layers of the ejecta, particularly those at very early times. Using the method of 'Abundance Tomography' we have derived iron-group abundances in SN 2010jn with a precision better than in any previously studied SN Ia. Optimum fits to the data can be obtained if burned material is present even at high velocities, including significant mass fractions of iron-group elements. This is consistent with the slow decline rate (or high 'stretch') of the light curve of SN 2010jn, and consistent with the results of delayed-detonation models. Early-phase UV spectra and detailed time-dependent series of further SNe Ia offer a promising probe of the nature of the SN Ia mechanism.

We present five new outbursting AM CVn systems and one candidate discovered as part of an ongoing search for such systems using the Palomar Transient Factory (PTF). This is the first large-area, systematic search for AM CVn systems using only large-amplitude photometric variability to select candidates. Three of the confirmed systems and the candidate system were discovered as part of the PTF transient search. Two systems were found as part of a search for outbursts through the PTF photometric database. We discuss the observed characteristics of each of these systems, including the orbital periods of two systems. We also consider the position of these systems, selected in a colour-independent survey, in colour-colour space and compare to systems selected solely by their colours. We find that the colours of our newly discovered systems do not differ significantly from those of previously known systems, but significant errors preclude a definitive answer.

Robo-AO, a fully autonomous, laser guide star adaptive optics and science system, is being commissioned at Palomar Observatory's 60-inch telescope. Here we discuss the instrument, scientific goals and results of initial on-sky operation.

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be \theta_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

We present here the analysis of about 19,500 new star hours of low ecliptic latitude observations (|b| < 20 deg) obtained by the Hubble Space Telescope's FGS over a time span of more than nine years; which is an addition to the 12,000 star hours previously analyzed by Schlichting et al. (2009). Our search for stellar occultations by small Kuiper belt objects (KBOs) yielded one new candidate event corresponding to a body with a 530 +/-70m radius at a distance of about 40AU. Using bootstrap simulations, we estimate a probability of approx 5%, that this event is due to random statistical fluctuations within the new data set. Combining this new event with the single KBO occultation reported by Schlichting et al. (2009) we arrive at the following results: 1) The ecliptic latitudes of 6.6 deg and 14.4 deg of the two events are consistent with the observed inclination distribution of larger, 100km-sized KBOs. 2) Assuming that small, sub-km sized KBOs have the same ecliptic latitude distribution as their larger counterparts, we find an ecliptic surface density of KBOs with radii larger than 250m of N(r>250m) = 1.1^{+1.5}_{-0.7} x 10^7 deg^{-2}; if sub-km sized KBOs have instead a uniform ecliptic latitude distribution for -20 deg < b< 20 deg then N(r>250m) = 4.4^{+5.8}_{-2.8} x 10^6 deg^{-2}. This is the best measurement of the surface density of sub-km sized KBOs to date. 3) Assuming the KBO size distribution can be well described by a single power law given by N(>r) \propto r^{1-q}, where N(>r) is the number of KBOs with radii greater than r, and q is the power law index, we find q=3.8+/-0.2 for a KBO ecliptic latitude distribution that follows the observed distribution for larger, 100-km sized KBOs. 4) Regardless of the exact power law, our results suggest that small KBOs are numerous enough to satisfy the required supply rate for the Jupiter family comets. (Abridged)

The direct detection of binary systems in wide-field surveys is limited by the size of the stars' point-spread-functions (PSFs). A search for elongated objects can find closer companions, but is limited by the precision to which the PSF shape can be calibrated for individual stars. We have developed the BinaryFinder algorithm to search for close binaries by using precision measurements of PSF ellipticity across wide-field survey images. We show that the algorithm is capable of reliably detecting binary systems down to approximately 1/5 of the seeing limit, and can directly measure the systems' position angles, separations and contrast ratios. To verify the algorithm's performance we evaluated 100,000 objects in Palomar Transient Factory (PTF) wide-field-survey data for signs of binarity, and then used the Robo-AO robotic laser adaptive optics system to verify the parameters of 44 high-confidence targets. We show that BinaryFinder correctly predicts the presence of close companions with a <5% false-positive rate, measures the detected binaries' position angles within 2 degrees and separations within 25%, and weakly constrains their contrast ratios. When applied to the full PTF dataset, we estimate that BinaryFinder will discover and characterize ~450,000 physically-associated binary systems with separations <2 arcseconds and magnitudes brighter than R=18. New wide-field synoptic surveys with high sensitivity and sub-arcsecond angular resolution, such as LSST, will allow BinaryFinder to reliably detect millions of very faint binary systems with separations as small as 0.1 arcseconds.

Only a handful of supernovae (SNe) have been studied in multi-wavelength from radio to X-rays, starting a few days after explosion. The early detection and classification of the nearby type IIb SN2011dh/PTF11eon in M51 provides a unique opportunity to conduct such observations. We present detailed data obtained at the youngest phase ever of a core-collapse supernova (days 3 to 12 after explosion) in the radio, millimeter and X-rays; when combined with optical data, this allows us to explore the early evolution of the SN blast wave and its surroundings. Our analysis shows that the expanding supernova shockwave does not exhibit equipartition (e_e/e_B ~ 1000), and is expanding into circumstellar material that is consistent with a density profile falling like R^-2. Within modeling uncertainties we find an average velocity of the fast parts of the ejecta of 15,000 +/- 1800 km/s, contrary to previous analysis. This velocity places SN 2011dh in an intermediate blast-wave regime between the previously defined compact and extended SN IIb subtypes. Our results highlight the importance of early (~ 1 day) high-frequency observations of future events. Moreover, we show the importance of combined radio/X-ray observations for determining the microphysics ratio e_e/e_B.

We present an analysis of the Hubble diagram for 12 Type Ia supernovae (SNe Ia) observed in the near-infrared J and H bands. We select SNe exclusively from the redshift range 0.03 < z < 0.09 to reduce uncertainties coming from peculiar velocities while remaining in a cosmologically well-understood region. All of the SNe in our sample exhibit no spectral or B-band light-curve peculiarities and lie in the B-band stretch range of 0.8-1.15. Our results suggest that SNe Ia observed in the near-infrared (NIR) are the best known standard candles. We fit previously determined NIR light-curve templates to new high-precision data to derive peak magnitudes and to determine the scatter about the Hubble line. Photometry of the 12 SNe is presented in the natural system. Using a standard cosmology of (H_0, Omega_m, Lambda) = (70,0.27,0.73) we find a median J-band absolute magnitude of M_J = -18.39 with a scatter of 0.116 and a median H-band absolute magnitude of M_H = -18.36 with a scatter of 0.085. The scatter in the H band is the smallest yet measured. We search for correlations between residuals in the J- and H-band Hubble diagrams and SN properties, such as SN colour, B-band stretch and the projected distance from host-galaxy centre. The only significant correlation is between the J-band Hubble residual and the J-H pseudo-colour. We also examine how the scatter changes when fewer points in the near-infrared are used to constrain the light curve. With a single point in the H band taken anywhere from 10 days before to 15 days after B-band maximum light and a prior on the date of H-band maximum set from the date of B-band maximum, we find that we can measure distances to an accuracy of 6%. The precision of SNe Ia in the NIR provides new opportunities for precision measurements of both the expansion history of the universe and peculiar velocities of nearby galaxies.

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

We present an analysis of the maximum light, near ultraviolet (NUV; 2900-5500 A) spectra of 32 low redshift (0.001<z<0.08) Type Ia supernovae (SNe Ia), obtained with the Hubble Space Telescope (HST). We combine this spectroscopic sample with high-quality gri light curves obtained with robotic telescopes to measure photometric parameters, such as stretch, optical colour, and brightness. By comparing our data to a comparable sample of SNe Ia at intermediate-z (0.4<z<0.9), we detect modest spectral evolution (3-sigma), in the sense that our mean low-z NUV spectrum has a depressed flux compared to its intermediate-z counterpart. We also see a strongly increased dispersion about the mean with decreasing wavelength, confirming the results of earlier surveys. These trends are consistent with changes in metallicity as predicted by contemporary SN Ia spectral models. We also examine the properties of various NUV spectral diagnostics in the individual spectra. We find a general correlation between stretch and the velocity (or position) of many NUV spectral features. In particular, we observe that higher stretch SNe have larger Ca II H&K velocities, that also correlate with host galaxy stellar mass. This latter trend is probably driven by the well-established correlation between stretch and stellar mass. We find no trends between UV spectral features and optical colour. Mean spectra constructed according to whether the SN has a positive or negative Hubble residual show very little difference at NUV wavelengths, indicating that the NUV evolution and variation we identify do not directly correlate with Hubble residuals. Our work confirms and strengthens earlier conclusions regarding the complex behaviour of SNe Ia in the NUV spectral region, but suggests the correlations we find are more useful in constraining progenitor models than improving the use of SNe Ia as cosmological probes.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

We report the discovery of two new halo velocity groups (Cancer groups A and B) traced by 8 distant RR Lyrae stars and observed by the Palomar Transient Factory (PTF) survey at R.A.~129 deg, Dec~20 deg (l~205 deg, b~32 deg). Located at 92 kpc from the Galactic center (86 kpc from the Sun), these are some of the most distant substructures in the Galactic halo known to date. Follow-up spectroscopic observations with the Palomar Observatory 5.1-m Hale telescope and W. M. Keck Observatory 10-m Keck I telescope indicate that the two groups are moving away from the Galaxy at v_{gsr} = 78.0+-5.6 km/s (Cancer group A) and v_{gsr} = 16.3+-7.1 km/s (Cancer group B). The groups have velocity dispersions of \sigma_{v_{gsr}}=12.4+-5.0 km/s and \sigma_{v_{gsr}}=14.9+-6.2 km/s, and are spatially extended (about several kpc) making it very unlikely that they are bound systems, and are more likely to be debris of tidally disrupted dwarf galaxies or globular clusters. Both groups are metal-poor (median metallicities of [Fe/H] = -1.6 dex and [Fe/H] =-2.1 dex), and have a somewhat uncertain (due to small sample size) metallicity dispersion of ~0.4 dex, suggesting dwarf galaxies as progenitors. Two additional RR Lyrae stars with velocities consistent with those of the Cancer groups have been observed ~25 deg east, suggesting possible extension of the groups in that direction.

We construct a photometrically calibrated catalog of non-variable sources from the Palomar Transient Factory (PTF) observations. The first version of this catalog presented here, the PTF photometric catalog 1.0, contains calibrated R_PTF-filter magnitudes for about 21 million sources brighter than magnitude 19, over an area of about 11233 deg^2. The magnitudes are provided in the PTF photometric system, and the color of a source is required in order to convert these magnitudes into other magnitude systems. We estimate that the magnitudes in this catalog have typical accuracy of about 0.02 mag with respect to magnitudes from the Sloan Digital Sky Survey. The median repeatability of our catalog's magnitudes for stars between 15 and 16 mag, is about 0.01 mag, and it is better than 0.03 mag for 95% of the sources in this magnitude range. The main goal of this catalog is to provide reference magnitudes for photometric calibration of visible light observations. Subsequent versions of this catalog, which will be published incrementally online, will be extended to a larger sky area and will also include g_PTF-filter magnitudes, as well as variability and proper motion information.

(Abridged). The optical light curve of some SNe may be powered by the outward diffusion of the energy deposited by the explosion shock in optically thick circumstellar matter (CSM). Recently, it was shown that the radiation-mediated and -dominated shock in an optically thick wind must transform into a collisionless shock and can produce hard X-rays. The X-rays are expected to peak at late times, relative to maximum visible light. Here we report on a search, using Swift and Chandra, for X-ray emission from 28 SNe that belong to classes whose progenitors are suspected to be embedded in dense CSM (IIn/Ibn/SLSN-I). Two SNe in our sample have X-ray properties that are roughly consistent with the expectation for X-rays from a collisionless shock in optically thick CSM. Therefore, we suggest that their optical light curves are powered by shock breakout in CSM. We show that two other events were too X-ray bright during the SN maximum optical light to be explained by the shock breakout model. We conclude that the light curves of some, but not all, type-IIn/Ibn SNe are powered by shock breakout in CSM. For the rest of the SNe in our sample, including all the SLSN-I events, our X-ray limits are not deep enough and were typically obtained at too early times to conclude about their nature. We argue that the optical light curves of SNe, for which the X-ray emission peaks at late times, are likely powered by the diffusion of shock energy from a dense CSM. We comment about the possibility to detect some of these events in radio.

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

We report the discovery by the Swift hard X-ray monitor of the transient source Swift J2058.4+0516 (Sw J2058+05). Our multi-wavelength follow-up campaign uncovered a long-lived (duration >~ months), luminous X-ray (L_X,iso ~ 3 x 10^47 erg s^-1) and radio (nu L_nu,iso ~ 10^42 erg s^-1) counterpart. The associated optical emission, however, from which we measure a redshift of 1.1853, is relatively faint, and this is not due to a large amount of dust extinction in the host galaxy. Based on numerous similarities with the recently discovered GRB 110328A / Swift J164449.3+573451 (Sw J1644+57), we suggest that Sw J2058+05 may be the second member of a new class of relativistic outbursts resulting from the tidal disruption of a star by a supermassive black hole. If so, the relative rarity of these sources (compared with the expected rate of tidal disruptions) implies that either these outflows are extremely narrowly collimated (theta < 1 degree), or only a small fraction of tidal disruptions generate relativistic ejecta. Analogous to the case of long-duration gamma-ray bursts and core-collapse supernovae, we speculate that rapid spin of the black hole may be a necessary condition to generate the relativistic component. Alternatively, if powered by gas accretion (i.e., an active galactic nucleus [AGN]), Sw J2058+05 would seem to represent a new mode of variability in these sources, as the observed properties appear largely inconsistent with known classes of AGNs capable of generating relativistic jets (blazars, narrow-line Seyfert 1 galaxies).

We present the discovery of PTF 10vgv, a Type Ic supernova detected by the Palomar Transient Factory, using the Palomar 48-inch telescope (P48). R-band observations of the PTF 10vgv field with P48 probe the supernova emission from its very early phases (about two weeks before R-band maximum), and set limits on its flux in the week prior to the discovery. Our sensitive upper limits and early detections constrain the post-shock-breakout luminosity of this event. Via comparison to numerical (analytical) models, we derive an upper-limit of R \lesssim 4.5 Rsun (R \lesssim 1 Rsun) on the radius of the progenitor star, a direct indication in favor of a compact Wolf-Rayet star. Applying a similar analysis to the historical observations of SN 1994I, yields R \lesssim 1/4 Rsun for the progenitor radius of this supernova.

Colors of Trans Neptunian Objects (TNOs) are used to study the evolutionary processes of bodies in the outskirts of the Solar System, and to test theories regarding their origin. Here I describe a search for serendipitous Sloan Digital Sky Survey (SDSS) observations of known TNOs and Centaurs. I present a catalog of SDSS photometry, colors and astrometry of 388 measurements of 42 outer Solar-System objects. I find a weak evidence, at the ~2-sigma level (per trial), for a correlation between the g-r color and inclination of scattered disk objects and hot classical KBOs. I find a correlation between the g-r color and the angular momentum in the z direction of all the objects in this sample. Light curves as a function of phase angle are constructed for 13 objects. The steepness of the slopes of these light curves suggests that the coherent backscatter mechanism plays a major role in the reflectivity of outer Solar-System small objects at small phase angles. I find a weak evidence for an anti-correlation, significant at the 2-sigma confidence level (per trial), between the g-band phase angle slope parameter and the semi-major axis, as well as the aphelion distance, of these objects. I discuss the origin of this possible correlation and argue that if this correlation is real it probably indicates that "Sedna"-like objects have a different origin than other classes of TNOs. Finally, I identify several objects with large variability amplitudes (abridged).