We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called "extra-tidal") if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area.

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

(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 Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

Quasars have long been known as intrinsically variable sources, but the physical mechanism underlying the temporal optical/UV variability is still not well understood. We propose a novel nonparametric method for modeling and forecasting the optical variability of quasars utilizing an autoencoder neural network to gain insight into the underlying processes. The autoencoder is trained with ~15,000 decade-long quasar light curves obtained by the Catalina Real-time Transient Survey selected with negligible flux contamination from the host galaxy. The autoencoder's performance in forecasting the temporal flux variation of quasars is superior to that of the damped random walk process. We find a temporal asymmetry in the optical variability and a novel relation - the amplitude of the variability asymmetry decreases as luminosity and/or black hole mass increases - is suggested with the help of autoencoded features. The characteristics of the variability asymmetry are in agreement with those from the self-organized disk instability model, which predicts that the magnitude of the variability asymmetry decreases as the ratio of the diffusion mass to inflow mass in the accretion disk increases.

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

We present extensive datasets for a class of intermediate-luminosity optical transients known as "luminous red novae" (LRNe). They show double-peaked light curves, with an initial rapid luminosity rise to a blue peak (at -13 to -15 mag), which is followed by a longer-duration red peak that sometimes is attenuated, resembling a plateau. The progenitors of three of them (NGC4490-2011OT1, M101-2015OT1, and SNhunt248), likely relatively massive blue to yellow stars, were also observed in a pre-eruptive stage when their luminosity was slowly increasing. Early spectra obtained during the first peak show a blue continuum with superposed prominent narrow Balmer lines, with P Cygni profiles. Lines of Fe II are also clearly observed, mostly in emission. During the second peak, the spectral continuum becomes much redder, Halpha is barely detected, and a forest of narrow metal lines is observed in absorption. Very late-time spectra (~6 months after blue peak) show an extremely red spectral continuum, peaking in the infrared (IR) domain. Halpha is detected in pure emission at such late phases, along with broad absorption bands due to molecular overtones (such as TiO, VO). We discuss a few alternative scenarios for LRNe. Although major instabilities of single massive stars cannot be definitely ruled out, we favour a common envelope ejection in a close binary system, with possibly a final coalescence of the two stars. The similarity between LRNe and the outburst observed a few months before the explosion of the Type IIn SN 2011ht is also discussed.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

We present deep MUSE observations of five quasars within the first Gyr of the Universe ($z\gtrsim6$), four of which display extended Ly$\alpha$ halos. After PSF-subtraction, we reveal halos surrounding two quasars for the first time, as well as confirming the presence of two more halos for which tentative detections exist in long-slit spectroscopic observations and narrow-band imaging. The four Ly$\alpha$ halos presented here are diverse in morphology and size, they each display spatial asymmetry, and none are centred on the position of the quasar. Spectra of the diffuse halos demonstrate that none are dramatically offset in velocity from the systemic redshift of the quasars ($\Delta$ v $< 200$ kms$^{-1}$), however each halo shows a broad Ly$\alpha$ line, with a velocity width of order $\sim1000$ kms$^{-1}$. Total Ly$\alpha$ luminosities range between $\sim$ $2 \times 10^{43}$ erg s$^{-1}$ and $\sim$ $2 \times 10^{44}$ erg s$^{-1}$, reaching maximum radial extents of $13 - 30$ pkpc from the quasar positions. We find larger sizes and higher Ly$\alpha$ luminosities than previous literature results at this redshift, but find no correlation between the quasar properties and the Ly$\alpha$ halo, suggesting that the detected emission is most closely related to the physical properties of the circum-galactic medium

We present a detailed study of a transient in the center of SDSS1115+0544 based on the extensive UV, optical, mid-IR light curves (LC) and spectra over 1200 days. The host galaxy is a quiescent early type galaxy at $z$ = 0.0899 with a blackhole mass of $2\times10^7M_\odot$. The transient underwent a 2.5 magnitude brightening over $\sim120$ days, reaching a peak $V$-band luminosity (extinction corrected) of $-20.9$ magnitude, then fading 0.5 magnitude over 200 days, settling into a plateau of $>600$ days. Following the optical brightening are the significant mid-IR flares at $3.4$ and $4.5\mu$m, with a peak time delay of $\sim180$ days. The mid-IR LCs are explained as the echo of UV photons by a dust medium with a radius of $5\times10^{17}$ cm, consistent with $\rm E(B-V)$ of 0.58 inferred from the spectra. This event is very energetic with an extinction corrected $L_{bol} \sim 4\times10^{44}$ erg s$^{-1}$. Optical spectra over 400 days in the plateau phase revealed newly formed broad H$\alpha, \beta$ emission with a FWHM of $\sim3750$ km s$^{-1}$ and narrow coronal lines such as [Fe VII], [Ne V]. This flare also has a steeply rising UV continuum, detected by multi-epoch $Swift$ data at $+700$ to $+900$ days post optical peak. The broad Balmer lines and the UV continuum do not show significant temporal variations. The slow evolving LCs over 1200 days, the constant Balmer lines and UV continuum at late-times rule out TDE and SN IIn as the physical model for this event. We propose that this event is a `turn-on' AGN, transitioning from a quiescent state to a type 1 AGN with a sub-Eddington accretion rate of $0.017M_\odot$/yr. This change occurred on a very short time scale of $\sim 120- 200$ days. The discovery of such a rapid `turn-on' AGN poses challenges to accretion disk theories and may indicate such event is not extremely rare.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

We present an analysis of nine years of Catalina Surveys optical photometry for 1.4 million spectroscopically confirmed SDSS galaxies. We find 717 outburst events that were not reported by ongoing transient surveys. These events have timescales ranging from weeks to years. More than two thirds of these new events are found in starforming galaxies, while such galaxies only constitute ~20% of our sample. Based on the properties of the hosts and events, we find that almost all of the new events are likely to be associated with regular supernovae. However, a small number of long-timescale events are found among the galaxies containing AGN. These events have similar properties to those recently found in the analyses of light curves of large samples of AGN. Given the lack of such events among the more than a million passive galaxies in the sample, we suggest that the long outbursts are associated with super-massive black holes or their environments.

Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400\AA\ has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine.

We present an updated catalog of 4680 northern eclipsing binaries (EBs) with Algol-type light curve morphology (i.e., with well-defined beginning and end of primary and secondary eclipses), using data from the Catalina Sky Surveys. Our work includes revised period determinations, phenomenological parameters of the light curves, and system morphology classification based on machine learning techniques. While most of the new periods are in excellent agreement with those provided in the original Catalina catalogs, improved values are now available for ~10% of the stars. A total of 3456 EBs were classified as detached and 449 as semi-detached, while 145 cannot be classified unambiguously into either subtype. The majority of the SD systems seems to be comprised of short-period Algols. By applying color criteria, we searched for K- and M-type dwarfs in these data, and present a subsample of 609 EB candidates for further investigation. We report 119 EBs (2.5% of the total sample) that show maximum quadrature light variations over long timescales, with periods bracketing the range 4.5-18 yrs and fractional luminosity variance of 0.04-0.13. We discuss possible causes for this, making use of models of variable starspot activity in our interpretation of the results.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

SNhunt151 was initially classified as a supernova (SN) impostor (nonterminal outburst of a massive star). It exhibited a slow increase in luminosity, lasting about 450 d, followed by a major brightening that reaches M_V ~ -18 mag. No source is detected to M_V > -13 mag in archival images at the position of SNhunt151 before the slow rise. Low-to-mid-resolution optical spectra obtained during the pronounced brightening show very little evolution, being dominated at all times by multicomponent Balmer emission lines, a signature of interaction between the material ejected in the new outburst and the pre-existing circumstellar medium. We also analyzed mid-infrared images from the Spitzer Space Telescope, detecting a source at the transient position in 2014 and 2015. Overall, SNhunt151 is spectroscopically a Type IIn SN, somewhat similar to SN2009ip. However, there are also some differences, such as a slow pre-discovery rise, a relatively broad light-curve peak showing a longer rise time (~ 50 d) and a slower decline, along with a negligible change in the temperature around the peak (T < 10^4 K). We suggest that SNhunt151 is the result of an outburst, or a SN explosion, within a dense circumstellar nebula, similar to those embedding some luminous blue variables like Eta Carinae and originating from past mass-loss events.

We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We provide, for the first time, robust observational constraints on the galaxy major merger fraction up to $z\approx 6$ using spectroscopic close pair counts. Deep Multi Unit Spectroscopic Explorer (MUSE) observations in the Hubble Ultra Deep Field (HUDF) and Hubble Deep Field South (HDF-S) are used to identify 113 secure close pairs of galaxies among a parent sample of 1801 galaxies spread over a large redshift range ($0.2<z<6$) and stellar masses ($10^7-10^{11} M_\odot$), thus probing about 12 Gyr of galaxy evolution. Stellar masses are estimated from spectral energy distribution (SED) fitting over the extensive UV-to-NIR HST photometry available in these deep Hubble fields, adding Spitzer IRAC bands to better constrain masses for high-redshift ($z\geqslant 3$) galaxies. These stellar masses are used to isolate a sample of 54 major close pairs with a galaxy mass ratio limit of 1:6. Among this sample, 23 pairs are identified at high redshift ($z\geqslant 3$) through their Ly$\alpha$ emission. The sample of major close pairs is divided into five redshift intervals in order to probe the evolution of the merger fraction with cosmic time. Our estimates are in very good agreement with previous close pair counts with a constant increase of the merger fraction up to $z\approx 3$ where it reaches a maximum of 20%. At higher redshift, we show that the fraction slowly decreases down to about 10% at $z\approx6$. The sample is further divided into two ranges of stellar masses using either a constant separation limit of $10^{9.5} M_\odot$ or the median value of stellar mass computed in each redshift bin. Overall, the major close pair fraction for low-mass and massive galaxies follows the same trend. These new, homogeneous, and robust estimates of the major merger fraction since $z\approx6$ are in good agreement with recent predictions of cosmological numerical simulations.

We present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3<z<6$ with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 hour, we detect 100 Ly$\alpha$ emitters and find $X_{\rm{Ly} \alpha}\gtrsim0.5$ for most of the redshift range covered, with 29 per cent of the Ly$\alpha$ sample exhibiting rest equivalent widths (rest-EWs) $\leq$ 15\AA. Adopting a range of rest-EW cuts (0 - 75\AA), we find no evidence of a dependence of $X_{\rm{Ly} \alpha}$ on either redshift or UV luminosity.

We have conducted a two-layered spectroscopic survey (1'x1' ultra deep and 3'x3' deep regions) in the Hubble Ultra Deep Field (HUDF) with the Multi Unit Spectroscopic Explorer (MUSE). The combination of a large field of view, high sensitivity, and wide wavelength coverage provides an order of magnitude improvement in spectroscopically confirmed redshifts in the HUDF; i.e., 1206 secure spectroscopic redshifts for HST continuum selected objects, which corresponds to 15% of the total (7904). The redshift distribution extends well beyond z>3 and to HST/F775W magnitudes as faint as ~30 mag (AB, 1-sigma). In addition, 132 secure redshifts were obtained for sources with no HST counterparts that were discovered in the MUSE data cubes by a blind search for emission-line features. In total, we present 1338 high quality redshifts, which is a factor of eight increase compared with the previously known spectroscopic redshifts in the same field. We assessed redshifts mainly with the spectral features [OII] at z<1.5 (473 objects) and Lya at 2.9<z<6.7 (692 objects). With respect to F775W magnitude, a 50% completeness is reached at 26.5 mag for ultra deep and 25.5 mag for deep fields, and the completeness remains >~20% up to 28-29 mag and ~27 mag, respectively. We used the determined redshifts to test continuum color selection (dropout) diagrams of high-z galaxies. The selection condition for F336W dropouts successfully captures ~80% of the targeted z~2.7 galaxies. However, for higher redshift selections (F435W, F606W, and F775W dropouts), the success rates decrease to ~20-40%. We empirically redefine the selection boundaries to make an attempt to improve them to ~60%. The revised boundaries allow bluer colors that capture Lya emitters with high Lya equivalent widths falling in the broadbands used for the color-color selection. Along with this paper, we release the redshift and line flux catalog.

Astronomy light curves are sparse, gappy, and heteroscedastic. As a result standard time series methods regularly used for financial and similar datasets are of little help and astronomers are usually left to their own instruments and techniques to classify light curves. A common approach is to derive statistical features from the time series and to use machine learning methods, generally supervised, to separate objects into a few of the standard classes. In this work, we transform the time series to two-dimensional light curve representations in order to classify them using modern deep learning techniques. In particular, we show that convolutional neural networks based classifiers work well for broad characterization and classification. We use labeled datasets of periodic variables from CRTS survey and show how this opens doors for a quick classification of diverse classes with several possible exciting extensions.

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.

A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic variability, which is not frequently observed in the CV systems. These long-term variations are possibly caused by various mechanisms, such as the precession of the accretion disk, hierarchical triple star system, magnetic field change of the companion star, and other possible mechanisms. We discuss the possible mechanisms in this study. If the long-term period is less than several tens of days, the disk precession period scenario is favored. However, the hierarchical triple star system or the variations in magnetic field strengths are most likely the predominant mechanisms for longer periods.

There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $\Delta m = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes.

We conduct a systematic search for galaxies at $z=0.1-1.5$ with [OII]$\lambda3727$, [OIII]$\lambda5007$, or H$\alpha\lambda6563$ emission lines extended over at least 30 kpc by using deep narrowband and broadband imaging in Subaru-XMM Deep Survey (SXDS) field. These extended emission-line galaxies are dubbed [OII], [OIII], or H$\alpha$ blobs. Based on a new selection method that securely selects extended emission-line galaxies, we find 77 blobs at $z=0.40-1.46$ with the isophotal area of emission lines down to $1.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ kpc$^{-2}$. Four of them are spectroscopically confirmed to be [OIII] blobs at $z=0.83$. We identify AGN activities in 8 blobs with X-ray and radio data and find that the fraction of AGN contribution increases with increasing isophotal area of the extended emission. With the Kolmogorov-Smirnov (KS) and Anderson-Darling tests, we confirm that the stellar-mass distributions of H$\alpha$ and [OII] blobs are not drawn from those of the emitters at the $>90$% confidence level in that H$\alpha$ and [OII] blobs are located at the massive end of the distributions, but cannot reject null hypothesis of being the same distributions in terms of the specific star formation rates. It is suggested that galactic-scale outflows tend to be more prominent in more massive star-forming galaxies. Exploiting our sample homogeneously selected over the large area, we derive the number densities of blobs at each epoch. The number densities of blobs decrease drastically with redshifts at the rate that is larger than that of the decrease of cosmic star formation densities.

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifiable manner whose physical origins are poorly understood. In this paper, we will review how quasars are identified by their variability and how these techniques can be improved, what physical insights into their variability can be gained from studying extreme examples of variability, and what approaches can be taken to increase the number of quasars known. These will demonstrate how astroinformatics is essential to discovering and understanding this important population.

We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40{\deg} away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distribution. The separation between the two line-of-sight components is at least 5 kpc at the edge of the VST ATLAS footprint, but appears to change along the stream, which allows us to conclude that these detections correspond to two physically independent stellar structures, rather than a mix of co-distant stellar populations within a single stream. Our discovery of a fork in the Sgr trailing arm is verified using Blue Horizontal Branch stars and our distance measurements are calibrated using RR Lyrae stars from the Catalina Real-Time Transient Survey. Comparing with numerical simulations of the Sgr dwarf disruption, the more distant of the two components in the fork matches perfectly with the track of the trailing debris. However, no obvious counterpart exists in the simulation for the closer line-of-sight component.

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10.

Context: Research on supernovae (SNe) over the past decade has confirmed that there is a distinct class of events which are much more luminous (by $\sim2$ mag) than canonical core-collapse SNe (CCSNe). These events with visual peak magnitudes $\lesssim-21$ are called superluminous SNe (SLSNe). Aims: There are a few intermediate events which have luminosities between these two classes. Here we study one such object, SN 2012aa. Methods: The optical photometric and spectroscopic follow-up observations of the event were conducted over a time span of about 120 days. Results: With V_abs at peak ~-20 mag, the SN is an intermediate-luminosity transient between regular SNe Ibc and SLSNe. It also exhibits an unusual secondary bump after the maximum in its light curve. We interpret this as a manifestation of SN-shock interaction with the CSM. If we would assume a $^{56}$Ni-powered ejecta, the bolometric light curve requires roughly 1.3 M_sun of $^{56}$Ni and an ejected mass of ~14 M_sun. This would also imply a high kinetic energy of the explosion, ~5.4$\times10^{51}$ ergs. On the other hand, the unusually broad light curve along with the secondary peak indicate the possibility of interaction with CSM. The third alternative is the presence of a central engine releasing spin energy that eventually powers the light curve over a long time. The host of the SN is a star-forming Sa/Sb/Sbc galaxy. Conclusions: Although the spectral properties and velocity evolution of SN 2012aa are comparable to those of normal SNe Ibc, its broad light curve along with a large peak luminosity distinguish it from canonical CCSNe, suggesting the event to be an intermediate-luminosity transient between CCSNe and SLSNe at least in terms of peak luminosity. We argue that SN 2012aa belongs to a subclass where CSM interaction plays a significant role in powering the SN, at least during the initial stages of evolution.

We report photometric and spectroscopic observations of the optical transient LSQ13zm. Historical data reveal the presence of an eruptive episode (that we label as `2013a') followed by a much brighter outburst (`2013b') three weeks later, that we argue to be the genuine supernova explosion. This sequence of events closely resemble those observed for SN2010mc and (in 2012) SN2009ip. The absolute magnitude reached by LSQ13zm during 2013a ($M_R=-14.87\pm0.25\,\rm{mag}$) is comparable with those of supernova impostors, while that of the 2013b event ($M_R=-18.46\pm0.21\,\rm{mag}$) is consistent with those of interacting supernovae. Our spectra reveal the presence of a dense and structured circumstellar medium, probably produced through numerous pre-supernova mass-loss events. In addition, we find evidence for high-velocity ejecta, with a fraction of gas expelled at more than 20000\kms. The spectra of LSQ13zm show remarkable similarity with those of well-studied core-collapse supernovae. From the analysis of the available photometric and spectroscopic data, we conclude that we first observed the last event of an eruptive sequence from a massive star, likely a Luminous Blue Variable, which a short time later exploded as a core-collapse supernova. The detailed analysis of archival images suggest that the host galaxy is a star-forming Blue Dwarf Compact Galaxy.

We present an update of a spectroscopic follow-up survey at low-resolution of a large number of RR Lyrae halo overdensity candidates found in the southern sky. The substructure candidates were identified in the RR Lyrae catalog of Torrealba et al. (2015) using Catalina Real-time Transient Survey (CRTS) data. Radial velocities and mean metallicities have been estimated for target stars in almost half of the original overdensities to assess their potential membership to coherent halo features.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.

The Outburst Catalogue contains a wide variety of observational properties for 722 dwarf nova-type (DN) cataclysmic variables (CVs) and 309 CVs of other types from the Catalina Real-time Transient Survey. In particular, it includes the apparent outburst and quiescent V-band magnitudes, duty cycles, limits on the recurrence time, upper- and lower-limits on the distance and absolute quiescent magnitudes, colour information, orbital parameters, and X-ray counterparts. These properties were determined by means of a classification script presented in this paper. The DN in the catalogue show a correlation between the outburst duty cycle and the orbital period (and outburst recurrence time), as well as between the quiescent absolute magnitude and the orbital period (and duty cycle). This is the largest sample of dwarf nova properties collected to date. Besides serving as a useful reference for individual systems and a means of selecting objects for targeted studies, it will prove valuable for statistical studies that aim to shed light on the formation and evolution of cataclysmic variables.

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

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.

We present new results on the evolution of the cosmic star formation rate as a function of stellar mass in the SXDS-UDS field. We make use of narrow-band selected emission line galaxies in four redshift slices between z = 1.46 and z = 0.63, and compute stellar masses by fitting a series of templates to recreate each galaxy's star formation history. We determine mass-binned luminosity functions in each redshift slice, and derive the star formation rate density (rhoSFR) as a function of mass using the [OIII] or [OII] emission lines. We calculate dust extinction and metallicity as a function of stellar mass, and investigate the effect of these corrections on the shape of the overall rhoSFR(M). We find that both these corrections are crucial for determining the shape of the rhoSFR(M), and its evolution with redshift. The fully corrected rhoSFR(M) is a relatively flat distribution, with the normalisation moving towards lower values of rhoSFR with increasing cosmic time/decreasing redshift, and requiring star formation to be truncated across all masses studied here. The peak of rhoSFR(M) is found in the 10^10.5<Msun<10^11.0 mass bin at z = 1.46. In the lower redshift slices the location of the peak is less certain, however low mass galaxies in the range 10^7.0<Msun<10^8.0 play an important part in the overall rhoSFR(M) out to at least z ~ 1.2.

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

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

We report on the search for new eclipsing white dwarf plus main-sequence (WDMS) binaries in the light curves of the Catalina surveys. We use a colour selected list of almost 2000 candidate WDMS systems from the Sloan Digital Sky Survey, specifically designed to identify WDMS systems with cool white dwarfs and/or early M type main-sequence stars. We identify a total of 17 eclipsing systems, 14 of which are new discoveries. We also find 3 candidate eclipsing systems, 2 main-sequence eclipsing binaries and 22 non-eclipsing close binaries. Our newly discovered systems generally have optical fluxes dominated by the main-sequence components, which have earlier spectral types than the majority of previously discovered eclipsing systems. We find a large number of ellipsoidally variable binaries with similar periods, near 4 hours, and spectral types M2--3, which are very close to Roche-lobe filling. We also find that the fraction of eclipsing systems is lower than found in previous studies and likely reflects a lower close binary fraction among WDMS binaries with early M-type main-sequence stars due to their enhanced angular momentum loss compared to fully convective late M type stars, hence causing them to become cataclysmic variables quicker and disappear from the WDMS sample. Our systems bring the total number of known detached, eclipsing WDMS binaries to 71.

Despite their importance, very few RR Lyrae (RRL) stars have been known to reside in binary systems. We report on a search for binary RRL in the OGLE-III Galactic bulge data. Our approach consists in the search for evidence of the light-travel time effect in so-called observed minus calculated ($O-C$) diagrams. Analysis of 1952 well-observed fundamental-mode RRL in the OGLE-III data revealed an initial sample of 29 candidates. We used the recently released OGLE-IV data to extend the baselines up to 17 years, leading to a final sample of 12 firm binary candidates. We provide $O-C$ diagrams and binary parameters for this final sample, and also discuss the properties of 8 additional candidate binaries whose parameters cannot be firmly determined at present. We also estimate that $\gtrsim 4$ per cent of the RRL reside in binary systems.

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

We use seven year's worth of observations from the Catalina Sky Survey and the Siding Spring Survey covering most of the northern and southern hemisphere at galactic latitudes higher than 20 degrees to search for serendipitously imaged moving objects in the outer solar system. These slowly moving objects would appear as stationary transients in these fast cadence asteroids surveys, so we develop methods to discover objects in the outer solar system using individual observations spaced by months, rather than spaced by hours, as is typically done. While we independently discover 8 known bright objects in the outer solar system, the faintest having $V=19.8\pm0.1$, no new objects are discovered. We find that the survey is nearly 100% efficient at detecting objects beyond 25 AU for $V\lesssim 19.1$ ($V\lesssim18.6$ in the southern hemisphere) and that the probability that there is one or more remaining outer solar system object of this brightness left to be discovered in the unsurveyed regions of the galactic plane is approximately 32%.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

Effective data visualization is a key part of the discovery process in the era of big data. It is the bridge between the quantitative content of the data and human intuition, and thus an essential component of the scientific path from data into knowledge and understanding. Visualization is also essential in the data mining process, directing the choice of the applicable algorithms, and in helping to identify and remove bad data from the analysis. However, a high complexity or a high dimensionality of modern data sets represents a critical obstacle. How do we visualize interesting structures and patterns that may exist in hyper-dimensional data spaces? A better understanding of how we can perceive and interact with multi dimensional information poses some deep questions in the field of cognition technology and human computer interaction. To this effect, we are exploring the use of immersive virtual reality platforms for scientific data visualization, both as software and inexpensive commodity hardware. These potentially powerful and innovative tools for multi dimensional data visualization can also provide an easy and natural path to a collaborative data visualization and exploration, where scientists can interact with their data and their colleagues in the same visual space. Immersion provides benefits beyond the traditional desktop visualization tools: it leads to a demonstrably better perception of a datascape geometry, more intuitive data understanding, and a better retention of the perceived relationships in the data.

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leq\delta\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 \sigma$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40\deg$.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

Over six years of operation, the Catalina Real-time Transient Survey (CRTS) has identified 1043 cataclysmic variable (CV) candidates --- the largest sample of CVs from a single survey to date. Here we provide spectroscopic identification of 85 systems fainter than g<19, including three AMCVn binaries, one helium-enriched CV, one polar and one new eclipsing CV. We analyse the outburst properties of the full sample and show that it contains a large fraction of low accretion rate CVs with long outburst recurrence times. We argue that most of the high accretion rate dwarf novae in the survey footprint have already been found and that future CRTS discoveries will be mostly low accretion rate systems. We find that CVs with white dwarf dominated spectra have significantly fewer outbursts in their CRTS light curves compared to disc-dominated CVs, reflecting the difference in their accretion rates. Comparing the CRTS sample to other samples of CVs, we estimate the overall external completeness to be 23.6 per cent, but show that as much as 56 per cent of CVs have variability amplitudes that are too small to be selected using the transient selection criteria employed by current ground-based surveys. The full table of CRTS CVs, including their outburst and spectroscopic properties examined in this paper, is provided in the online materials.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We investigate the properties of 367 ultra-short period binary candidates selected from 31,000 sources recently identified from Catalina Surveys data. Based on light curve morphology, along with WISE, SDSS and GALEX multi-colour photometry, we identify two distinct groups of binaries with periods below the 0.22 day contact binary minimum. In contrast to most recent work, we spectroscopically confirm the existence of M-dwarf+M-dwarf contact binary systems. By measuring the radial velocity variations for five of the shortest-period systems, we find examples of rare cool-white dwarf+M-dwarf binaries. Only a few such systems are currently known. Unlike warmer white dwarf systems, their UV flux and their optical colours and spectra are dominated by the M-dwarf companion. We contrast our discoveries with previous photometrically-selected ultra-short period contact binary candidates, and highlight the ongoing need for confirmation using spectra and associated radial velocity measurements. Overall, our analysis increases the number of ultra-short period contact binary candidates by more than an order of magnitude.

We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.

We present 855 cataclysmic variable candidates detected by the Catalina Real-time Transient Survey (CRTS) of which at least 137 have been spectroscopically confirmed and 705 are new discoveries. The sources were identified from the analysis of five years of data, and come from an area covering three quarters of the sky. We study the amplitude distribution of the dwarf novae CVs discovered by CRTS during outburst, and find that in quiescence they are typically two magnitudes fainter compared to the spectroscopic CV sample identified by SDSS. However, almost all CRTS CVs in the SDSS footprint have ugriz photometry. We analyse the spatial distribution of the CVs and find evidence that many of the systems lie at scale heights beyond those expected for a Galactic thin disc population. We compare the outburst rates of newly discovered CRTS CVs with the previously known CV population, and find no evidence for a difference between them. However, we find that significant evidence for a systematic difference in orbital period distribution. We discuss the CVs found below the orbital period minimum and argue that many more are yet to be identified among the full CRTS CV sample. We cross-match the CVs with archival X-ray catalogs and find that most of the systems are dwarf novae rather than magnetic CVs.

We present optical photometry and spectra of the super luminous type II/IIn supernova CSS121015:004244+132827 (z=0.2868) spanning epochs from -30 days (rest frame) to more than 200 days after maximum. CSS121015 is one of the more luminous supernova ever found and one of the best observed. The photometric evolution is characterized by a relatively fast rise to maximum (~40 days in the SN rest frame), and by a linear post-maximum decline. The light curve shows no sign of a break to an exponential tail. A broad Halpha is first detected at ~ +40 days (rest-frame). Narrow, barely-resolved Balmer and [O III] 5007 A lines, with decreasing strength, are visible along the entire spectral evolution. The spectra are very similar to other super luminous supernovae (SLSNe) with hydrogen in their spectrum, and also to SN 2005gj, sometimes considered a type Ia interacting with H-rich CSM. The spectra are also similar to a subsample of H-deficient SLSNe. We propose that the properties of CSS121015 are consistent with the interaction of the ejecta with a massive, extended, opaque shell, lost by the progenitor decades before the final explosion, although a magnetar powered model cannot be excluded. Based on the similarity of CSS121015 with other SLSNe (with and without H), we suggest that the shocked-shell scenario should be seriously considered as a plausible model for both types of SLSN.

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.

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

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

Using imaging from the Pan-STARRS1 survey, we identify a precursor outburst at epochs 287 and 170 days prior to the reported explosion of the purported Type IIn supernova (SN) 2011ht. In the Pan-STARRS data, a source coincident with SN 2011ht is detected exclusively in the \zps\ and \yps-bands. An absolute magnitude of M$_z\simeq$-11.8 suggests that this was an outburst of the progenitor star. Unfiltered, archival Catalina Real Time Transient survey images also reveal a coincident source from at least 258 to 138 days before the main event. We suggest that the outburst is likely to be an intrinsically red eruption, although we cannot conclusively exclude a series of erratic outbursts which were observed only in the redder bands by chance. This is only the fourth detection of an outburst prior to a claimed SN, and lends credence to the possibility that many more interacting transients have pre-explosion outbursts, which have been missed by current surveys.

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

This paper presents a new period finding method based on conditional entropy that is both efficient and accurate. We demonstrate its applicability on simulated and real data. We find that it has comparable performance to other information-based techniques with simulated data but is superior with real data, both for finding periods and just identifying periodic behaviour. In particular, it is robust against common aliasing issues found with other period-finding algorithms.

We report on multifrequency observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036 performed from radio to gamma-rays during 2008 August-2012 November by Fermi-LAT, Swift (XRT and UVOT), OVRO, VLBA, and VLA. No significant variability has been observed in gamma-rays, with 0.1-100 GeV flux that ranged between (3-7)x10^-8 ph/cm^2/s using 3-month time bins. The photon index of the LAT spectrum (Gamma=2.60+/-0.06) and the apparent isotropic gamma-ray luminosity, L(0.1-100 GeV)= 7.8x10^45 erg/s, over 51 months are typical of a flat spectrum radio quasar. The radio spectral variability and the one-sided structure, in addition to the observed gamma-ray luminosity, suggest a relativistic jet with a high Doppler factor. In contrast to SBS 0846+513, the VLBA at 15 GHz did not observe superluminal motion for PKS 1502+036. Despite having the optical characteristics typical of a narrow-line Seyfert 1 galaxy, radio and gamma-ray properties of PKS 1502+036 are found to be similar to those of a blazar at the low end of the black hole mass distribution for blazars. This is in agreement with what has been found in the case of the other gamma-ray emitting narrow-line Seyfert 1 SBS 0846+513.

We present new results on the cosmic star formation history in the SXDS-UDS field out to z=1.6. We compile narrow-band data from the Subaru Telescope and the Visible and Infrared Survey Telescope for Astronomy (VISTA) in conjunction with broad-band data from the SXDS and UDS, to make a selection of 5725 emission-line galaxies in 12 redshift slices, spanning 10 Gyr of cosmic time. We determine photometric redshifts for the sample using 11-band photometry, and use a spectroscopically confirmed subset to fine tune the resultant redshift distribution. We use the maximum-likelihood technique to determine luminosity functions in each redshift slice and model the selection effects inherent in any narrow-band selection statistically, to obviate the retrospective corrections ordinarily required. The deep narrow-band data are sensitive to very low star formation rates (SFRs), and allow an accurate evaluation of the faint end slope of the Schechter function, alpha. We find that alpha is particularly sensitive to the assumed faintest broadband magnitude of a galaxy capable of hosting an emission line, and propose that this limit should be empirically motivated. For this analysis we base our threshold on the limiting observed equivalent widths of emission lines in the local Universe. We compute the characteristic SFR of galaxies in each redshift slice, and the integrated SFR density, rho_{SFR}. We find our results to be in good agreement with the literature and parametrize the evolution of the SFR density as rho_{SFR} proportional to (1 + z)^{4.58} confirming a steep decline in star formation activity since z = 1.6.

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

We present the analysis of 1,207 RR Lyrae found in photometry taken by the Catalina Survey's Mount Lemmon telescope. By combining accurate distances for these stars with measurements for ~14,000 type-AB RR Lyrae from the Catalina Schmid telescope, we reveal an extended association that reaches Galactocentric distances beyond 100 kpc and overlaps the Sagittarius streams system. This result confirms earlier evidence for the existence of an outer halo tidal stream resulting from a disrupted stellar system. By comparing the RR Lyrae source density with that expected based on halo models, we find the detection has ~8 sigma significance. We investigate the distances, radial velocities, metallicities, and period-amplitude distribution of the RR Lyrae. We find that both radial velocities and distances are inconsistent with current models of the Sagittarius stream. We also find tentative evidence for a division in source metallicities for the most distant sources. Following prior analyses, we compare the locations and distances of the RR Lyrae with photometrically selected candidate horizontal branch stars and find supporting evidence that this structure spans at least 60 deg of the sky. We investigate the prospects of an association between the stream and unusual globular cluster NGC 2419.

There are many occasions when one does not have complete information in order to classify objects into different classes, and yet it is important to do the best one can since other decisions depend on that. In astronomy, especially time-domain astronomy, this situation is common when a transient is detected and one wishes to determine what it is in order to decide if one must follow it. We propose to use the Difference Boosting Neural Network (DBNN) which can boost differences between feature vectors of different objects in order to differentiate between them. We apply it to the publicly available data of the Catalina Real-Time Transient Survey (CRTS) and present preliminary results. We also describe another use with a stellar spectral library to identify spectra based on a few features. The technique itself is more general and can be applied to a varied class of problems.

We present the analysis of 12227 type-ab RR Lyrae found among the 200 million public lightcurves in the Catalina Surveys Data Release 1 (CSDR1). These stars span the largest volume of the Milky Way ever surveyed with RR Lyrae, covering ~20,000 square degrees of the sky (0 < RA < 360, -22 < Dec < 65 deg) to heliocentric distances of up to 60kpc. Each of the RR Lyrae are observed between 60 and 419 times over a six-year period. Using period finding and Fourier fitting techniques we determine periods and apparent magnitudes for each source. We find that the periods at generally accurate to sigma = 0.002% by comparison with 2842 previously known RR Lyrae and 100 RR Lyrae observed in overlapping survey fields. We photometrically calibrate the light curves using 445 Landolt standard stars and show that the resulting magnitudes are accurate to ~0.05 mags using SDSS data for ~1000 blue horizontal branch stars and 7788 of the RR Lyrae. By combining Catalina photometry with SDSS spectroscopy, we analyze the radial velocity and metallicity distributions for > 1500 of the RR Lyrae. Using the accurate distances derived for the RR Lyrae, we show the paths of the Sagittarius tidal streams crossing the sky at heliocentric distances from 20 to 60 kpc. By selecting samples of Galactic halo RR Lyrae, we compare their velocity, metallicity, and distance with predictions from a recent detailed N-body model of the Sagittarius system. We find that there are some significant differences between the distances and structures predicted and our observations.

Recent observations by Mauerhan et al. have shown the unprecedented transition of the previously identified luminous blue variable (and supernova impostor) SN 2009ip to a real Type IIn supernova (SN) explosion. We present high-cadence optical imaging of SN 2009ip obtained between 2012 UT Sep. 23.6 and Oct. 9.6, using 0.3-0.4 meter aperture telescopes from the Coral Towers Observatory in Cairns, Australia. The light curves show well-defined phases, including very rapid brightening early on (0.5 mag in 6 hr observed during the night of Sep. 24), a transition to a much slower rise between Sep. 25 and Sep. 28, and a plateau/peak around Oct. 7. These changes are coincident with the reported spectroscopic changes that, most likely, mark the start of a strong interaction between the fast SN ejecta and a dense circumstellar medium formed during the earlier LBV eruptions. In the 16-day observing period SN 2009ip brightened by 3.7 mag from I = 17.4 mag on Sep. 23.6 (MI = -14.2) to I = 13.7 mag (MI = -17.9) on Oct. 9.6, radiating ~2x10^49 erg in the optical wavelength range. Currently, SN 2009ip is more luminous than most Type IIP SNe and comparable to other Type IIn SNe.

We analyse the Catalina Real-time Transient Survey light curves of 835 spectroscopically confirmed white dwarf plus main-sequence binaries from the Sloan Digital Sky Survey (SDSS) with g<19, in search of new eclipsing systems. We identify 29 eclipsing systems, 12 of which were previously unknown. This brings the total number of eclipsing white dwarf plus main-sequence binaries to 49. Our set of new eclipsing systems contains two with periods of 1.9 and 2.3 days making them the longest period eclipsing white dwarf binaries known. We also identify one system which shows very large ellipsoidal modulation (almost 0.3 magnitudes), implying that the system is both very close to Roche-lobe overflow and at high inclination. However, our follow up photometry failed to firmly detect an eclipse meaning that this system either contains a cool white dwarf and hence the eclipse is very shallow and undetectable in our red-sensitive photometry or that it is non-eclipsing. Radial velocity measurements for the main-sequence stars in three of our newly identified eclipsing systems imply that their white dwarf masses are lower than those inferred from modelling their SDSS spectra. 13 non-eclipsing post common envelope binaries were also identified, from either reflection or ellipsoidal modulation effects. The white dwarfs in our newly discovered eclipsing systems span a wide range of parameters, including; low mass (~0.3Msol), very hot (80,000K) and a DC white dwarf. The spectral types of the main-sequence stars range from M2 to M6. This makes our sample ideal for testing white dwarf and low-mass star mass-radius relationships as well as close binary evolution.

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This presents some unusual challenges: the data are sparse, heterogeneous and incomplete; evolving in time; and most of the relevant information comes not from the data stream itself, but from a variety of archival data and contextual information (spatial, temporal, and multi-wavelength). We are exploring a variety of novel techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. The current surveys are already overwhelming our ability to effectively follow all of the potentially interesting events, and these challenges will grow by orders of magnitude over the next decade as the more ambitious sky surveys get under way. While we focus on an application in a specific domain (astrophysics), these challenges are more broadly relevant for event or anomaly detection and knowledge discovery in massive data streams.

Astronomy has been at the forefront of the development of the techniques and methodologies of data intensive science for over a decade with large sky surveys and distributed efforts such as the Virtual Observatory. However, it faces a new data deluge with the next generation of synoptic sky surveys which are opening up the time domain for discovery and exploration. This brings both new scientific opportunities and fresh challenges, in terms of data rates from robotic telescopes and exponential complexity in linked data, but also for data mining algorithms used in classification and decision making. In this paper, we describe how an informatics-based approach-part of the so-called "fourth paradigm" of scientific discovery-is emerging to deal with these. We review our experiences with the Palomar-Quest and Catalina Real-Time Transient Sky Surveys; in particular, addressing the issue of the heterogeneity of data associated with transient astronomical events (and other sensor networks) and how to manage and analyze it.

We present time-resolved optical spectroscopy of the dwarf nova CSS100603:112253-111037. Its optical spectrum is rich in helium, with broad, double-peaked emission lines produced in an accretion disc. We measure a line flux ratio HeI5876/H_alpha = 1.49 +/- 0.04, a much higher ratio than is typically observed in dwarf novae. The orbital period, as derived from the radial velocity of the line wings, is 65.233 +/- 0.015 minutes. In combination with the previously measured superhump period, this implies an extreme mass ratio of M_2/M_1 = 0.017 +/- 0.004. The H_alpha and HeI6678 emission lines additionally have a narrow central spike, as is often seen in the spectra of AM CVn type stars. Comparing its properties with CVs, AM CVn systems and hydrogen binaries below the CV period minimum, we argue that CSS100603:112253-111037 is the first compelling example of an AM CVn system forming via the evolved CV channel. With the addition of this system, evolved cataclysmic variables (CVs) now account for seven per cent of all known semi-detached white dwarf binaries with Porb < 76 min. Two recently discovered binaries may further increase this figure. Although the selection bias of this sample is not yet well defined, these systems support the evolved CV model as a possible formation channel for ultracompact accreting binaries. The orbital periods of the three ultracompact hydrogen accreting binaries overlap with those of the long period AM CVn stars, but there are currently no known systems in the period range 67 - 76 minutes.

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

Sky surveys represent a fundamental data basis for astronomy. We use them to map in a systematic way the universe and its constituents, and to discover new types of objects or phenomena. We review the subject, with an emphasis on the wide-field imaging surveys, placing them in a broader scientific and historical context. Surveys are the largest data generators in astronomy, propelled by the advances in information and computation technology, and have transformed the ways in which astronomy is done. We describe the variety and the general properties of surveys, the ways in which they may be quantified and compared, and offer some figures of merit that can be used to compare their scientific discovery potential. Surveys enable a very wide range of science; that is perhaps their key unifying characteristic. As new domains of the observable parameter space open up thanks to the advances in technology, surveys are often the initial step in their exploration. Science can be done with the survey data alone or a combination of different surveys, or with a targeted follow-up of potentially interesting selected sources. Surveys can be used to generate large, statistical samples of objects that can be studied as populations, or as tracers of larger structures. They can be also used to discover or generate samples of rare or unusual objects, and may lead to discoveries of some previously unknown types. We discuss a general framework of parameter spaces that can be used for an assessment and comparison of different surveys, and the strategies for their scientific exploration. As we move into the Petascale regime, an effective processing and scientific exploitation of such large data sets and data streams poses many challenges, some of which may be addressed in the framework of Virtual Observatory and Astroinformatics, with a broader application of data mining and knowledge discovery technologies.

We present phase-resolved spectroscopy of two new short period low accretion rate magnetic binaries, SDSSJ125044.42+154957.3 (Porb = 86 min) and SDSSJ151415.65+074446.5 (Porb = 89 min). Both systems were previously identified as magnetic white dwarfs from the Zeeman splitting of the Balmer absorption lines in their optical spectra. Their spectral energy distributions exhibit a large near-infrared excess, which we interpret as a combination of cyclotron emission and possibly a late type companion star. No absorption features from the companion are seen in our optical spectra. We derive the orbital periods from a narrow, variable H_alpha emission line which we show to originate on the companion star. The high radial velocity amplitude measured in both systems suggests a high orbital inclination, but we find no evidence for eclipses in our data. The two new systems resemble the polar EF Eri in its prolonged low state and also SDSSJ121209.31+013627.7, a known magnetic white dwarf plus possible brown dwarf binary, which was also recovered by our method.

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

The Catalina Real-time Transient Survey (CRTS) currently covers 33,000 deg^2 of the sky in search of transient astrophysical events, with time baselines ranging from 10 minutes to ~7 years. Data provided by the Catalina Sky Survey provides an unequaled baseline against which >4,000 unique optical transient events have been discovered and openly published in real-time. Here we highlight some of the discoveries of CRTS.

The time domain community wants robust and reliable tools to enable production of and subscription to community-endorsed event notification packets (VOEvent). The VAO Transient Facility (VTF) is being designed to be the premier brokering service for the community, both collecting and disseminating observations about time-critical astronomical transients but also supporting annotations and the application of intelligent machine-learning to those observations. This distinguishes two types of activity associated with the facility: core infrastructure and user services. In this paper, we will review the prior art in both areas and describe the planned capabilities of the VTF. In particular, we will focus on scalability and quality-of-service issues required by the next generation of sky surveys, such as LSST and SKA.

Synoptic sky surveys are becoming the largest data generators in astronomy, and they are opening a new research frontier, that touches essentially every field of astronomy. Opening of the time domain to a systematic exploration will strengthen our understanding of a number of interesting known phenomena, and may lead to the discoveries of as yet unknown ones. We describe some lessons learned over the past decade, and offer some ideas that may guide strategic considerations in planning and execution of the future synoptic sky surveys.

Exploration of the time domain - variable and transient objects and phenomena - is rapidly becoming a vibrant research frontier, touching on essentially every field of astronomy and astrophysics, from the Solar system to cosmology. Time domain astronomy is being enabled by the advent of the new generation of synoptic sky surveys that cover large areas on the sky repeatedly, and generating massive data streams. Their scientific exploration poses many challenges, driven mainly by the need for a real-time discovery, classification, and follow-up of the interesting events. Here we describe the Catalina Real-Time Transient Survey (CRTS), that discovers and publishes transient events at optical wavelengths in real time, thus benefiting the entire community. We describe some of the scientific results to date, and then focus on the challenges of the automated classification and prioritization of transient events. CRTS represents a scientific and a technological testbed and precursor for the larger surveys in the future, including the Large Synoptic Survey Telescope (LSST) and the Square Kilometer Array (SKA).

We describe the development of a system for an automated, iterative, real-time classification of transient events discovered in synoptic sky surveys. The system under development incorporates a number of Machine Learning techniques, mostly using Bayesian approaches, due to the sparse nature, heterogeneity, and variable incompleteness of the available data. The classifications are improved iteratively as the new measurements are obtained. One novel feature is the development of an automated follow-up recommendation engine, that suggest those measurements that would be the most advantageous in terms of resolving classification ambiguities and/or characterization of the astrophysically most interesting objects, given a set of available follow-up assets and their cost functions. This illustrates the symbiotic relationship of astronomy and applied computer science through the emerging discipline of AstroInformatics.

We present the discovery and follow-up observations of SN 2008jb, a core-collapse supernova in the dwarf irregular galaxy ESO 302-14 at 9.6 Mpc. This transient was missed by galaxy-targeted surveys and was only found in archival optical images obtained by CRTS and ASAS. It was detected shortly after explosion and reached a bright optical maximum, Vmax = 13.6 mag (M_Vmax ~ -16.5). The shape of the light curve shows a plateau of 100 days, followed by a drop of 1.4 mag in V-band to a decline with the approximate Co 56 decay slope, consistent with 0.04 Msun of Ni 56 synthesized in the explosion. A spectrum obtained 2 years after explosion shows a broad, boxy Halpha emission line, which is unusual for type IIP supernovae. We detect the supernova in archival Spitzer and WISE images obtained 8-14 months after explosion, which show clear signs of warm dust emission. The dwarf irregular host galaxy has a low gas-phase oxygen abundance, 12 + log(O/H) = 8.2 (~1/5 Solar), similar to those of the SMC and the hosts of long gamma-ray bursts and luminous core-collapse supernovae. We study the host environment using GALEX FUV, R-band, and Halpha images and find that the supernova occurred in a large star-formation complex. The morphology of the Halpha emission appears as a large shell (R = 350 pc) surrounding the FUV and optical emission. We estimate an age of ~9 Myr and a total mass of ~2 x 10^5 Msun for the star-formation complex. These properties are consistent with the expanding Halpha supershells observed in well-studied nearby dwarf galaxies, which are tell-tale signs of feedback from the cumulative effect of massive star winds and supernovae. The age estimated for the star-forming region suggests a relatively high-mass progenitor star with initial mass of ~20 Msun. We discuss the implications of these findings in the study of core-collapse supernova progenitors. (Abridged)

Using archival SDSS multi-epoch imaging data (Stripe 82), we have searched for the tidal disruption of stars by super-massive black holes in non-active galaxies. Two candidate tidal disruption events (TDEs) are identified. They have optical black-body temperatures 2 10^4 K and observed peak luminosities M_g=-18.3 and -20.4; their cooling rates are very low, qualitatively consistent with expectations for tidal disruption flares. Their properties are examined using i) SDSS imaging to compare them to other flares observed in the search, ii) UV emission measured by GALEX and iii) spectra of the hosts and of one of the flares. Our pipeline excludes optically identifiable AGN hosts, and our variability monitoring over 9 years provides strong evidence that these are not flares in hidden AGNs. The spectra and color evolution of the flares are unlike any SN observed to date, their strong late-time UV emission is particularly distinctive, and they are nuclear at high resolution, arguing against their being first cases of a previously-unobserved class of SNe or more extreme examples of known SN types. Taken together, the observed properties are difficult to reconcile with a SN or AGN-flare explanation, although an entirely new process specific to the inner few-hundred parsecs of non-active galaxies cannot be excluded. Our observed rate and method show the feasibility of obtaining a candidate TDE sample of hundreds of events and O(1) purity, using geometric resolution and host and flare color alone. A by-product of this work is quantification of the power-spectrum of extreme flares in AGNs.

We present high-speed ULTRACAM and SOFI photometry and X-shooter spectroscopy of the recently discovered 94 minute orbital period eclipsing white dwarf / main-sequence binary SDSS J085746.18+034255.3 (CSS 03170) and use these observations to measure the system parameters. We detect a shallow secondary eclipse and hence are able to determine an orbital inclination of 85.5 +/- 0.2 deg. The white dwarf has a mass of 0.51 +/- 0.05 Msun and a radius of 0.0247 +/- 0.0008 Rsun. With a temperature of 35,300 +/- 400K the white dwarf is highly over-inflated if it has a carbon-oxygen core, however if it has a helium core then its mass and radius are consistent with evolutionary models. Therefore, the white dwarf in SDSS J085746.18+034255.3 is most likely a helium core white dwarf with a mass close to the upper limit expected from evolution. The main-sequence star is an M8 dwarf with a mass of 0.09 +/- 0.01 Msun and a radius of 0.110 +/- 0.004 Rsun placing it close to the hydrogen burning limit. The system emerged from a common envelope ~20 million years ago and will reach a semi-detached configuration in ~400 million years, becoming a cataclysmic variable with a period of 66 minutes, below the period minimum.

In course of the search for the optical identifications associated with ROSAT X-ray sources we have found a highly variable object with the very unusual long-term behavior, color indices and high X-ray-to-optical flux ratio. We report the archival photometric light curve from the Catalina Sky Survey, optical spectroscopy from RTT150 and time-resolved photometry from Astrotel-Caucasus telescope. The object appears to be the magnetic cataclysmic variable (polar) with orbital period of P=3.20 hr.

Using Liverpool Telescope+RISE photometry we identify the 2.78 hour period binary star CSS 41177 as a detached eclipsing double white dwarf binary with a 21,100K primary star and a 10,500K secondary star. This makes CSS 41177 only the second known eclipsing double white dwarf binary after NLTT 11748. The 2 minute long primary eclipse is 40% deep and the secondary eclipse 10% deep. From Gemini+GMOS spectroscopy we measure the radial velocities of both components of the binary from the H{\alpha} absorption line cores. These measurements, combined with the light curve information, yield white dwarf masses of M1 = 0.283\pm0.064Msun and M2 = 0.274\pm0.034Msun, making them both helium core white dwarfs. As an eclipsing, double-lined spectroscopic binary CSS 41177 is ideally suited to measuring precise, model-independent masses and radii. The two white dwarfs will merge in roughly 1.1 Gyr to form a single sdB star.

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

Catalina Real-Time Transient Survey (CRTS) is a synoptic sky survey uses data streams from 3 wide-field telescopes in Arizona and Australia, covering the total area of ~30,000 deg2, down to the limiting magnitudes ~ 20 - 21 mag per exposure, with time baselines from 10 min to 6 years (and growing); there are now typically ~ 200 - 300 exposures per pointing, and coadded images reach deeper than 23 mag. The basic goal of CRTS is a systematic exploration and characterization of the faint, variable sky. The survey has detected ~ 3,000 high-amplitude transients to date, including ~ 1,000 supernovae, hundreds of CVs (the majority of them previously uncatalogued), and hundreds of blazars / OVV AGN, highly variable and flare stars, etc. CRTS has a complete open data philosophy: all transients are published immediately electronically, with no proprietary period at all, and all of the data (images, light curves) will be publicly available in the near future, thus benefiting the entire astronomical community. CRTS is a scientific and technological testbed and precursor for the grander synoptic sky surveys to come.

Although white dwarfs are believed to be the end point of most stellar evolution, unlike main sequence stars, they have not yet been the subject of dedicated time-domain surveys for exoplanets. We discuss how their size and distinctive colour make them excellent targets for wide-field searches for exoplanets. In particular, we note that planets of Earth-size can give rise to multi-magnitude eclipses of massive white dwarfs. Such a large signal is almost unmistakable and would be detectable even with very low-precision photometry. For objects of smaller size, the high accuracy photometry currently being used to detect Super-Earth and smaller planets transiting Sun-sized stars, is capable of revealing minor planets down to R~100km as they transit white dwarfs. Such observations can be used to test current evidence for asteroid-size objects being the cause for dust rings which have recently been observed for a number of white dwarfs. No other current exoplanet search method is capable of detecting such exo-asteroids. As an initial test of this search strategy, we combine synoptic data from the Catalina Sky Survey with multi-colour photometry and spectra from the Sloan Digital Sky Survey to search ~12,000 white dwarf lightcurves for eclipsing events. We find 20 new eclipsing white dwarf binary systems with low-mass companions. This doubles the number of known eclipsing white dwarfs and is expected to enable the determination of accurate white dwarf radii. Three of the discoveries have radii consistent with substellar systems and show no evidence of flux from the eclipsing object in their SDSS optical spectra, or near-IR data.

We report the discovery of a six-month-long mid-infrared transient, SDWFS-MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO Deep Wide-Field Survey Bootes field. The transient, located in a z=0.19 low luminosity (M_[4.5]~-18.6 mag, L/L_MilkyWay~0.01) metal-poor (12+log(O/H)~7.8) irregular galaxy, peaked at a mid-infrared absolute magnitude of M_[4.5]~-24.2 in the 4.5 micron Spitzer/IRAC band and emitted a total energy of at least 10^51 ergs. The optical emission was likely fainter than the mid-infrared, although our constraints on the optical emission are poor because the transient peaked when the source was "behind" the Sun. The Spitzer data are consistent with emission by a modified black body with a temperature of ~1350 K. We rule out a number of scenarios for the origin of the transient such as a Galactic star, AGN activity, GRB, tidal disruption of a star by a black hole and gravitational lensing. The most plausible scenario is a supernova exploding inside a massive, optically thick circumstellar medium, composed of multiple shells of previously ejected material. If the proposed scenario is correct, then a significant fraction (~10%) of the most luminous supernova may be self-enshrouded by dust not only before but also after the supernova occurs. The spectral energy distribution of the progenitor of such a supernova would be a slightly cooler version of eta Carina, peaking at 20-30 microns.

With orbital periods of the order of tens of minutes or less, the AM Canum Venaticorum stars are ultracompact, hydrogen deficient binaries with the shortest periods of any binary subclass, and are expected to be among the strongest gravitational wave sources in the sky. To date, the only known eclipsing source of this type is the P = 28 min binary SDSS J0926+3624. We present multiband, high time resolution light curves of this system, collected with WHT/ULTRACAM in 2006 and 2009. We supplement these data with additional observations made with LT/RISE, XMM_Newton and the Catalina Real-Time Transient Survey. From light curve models we determine the mass ratio to be q = M2 / M1 = 0.041 +/- 0.002 and the inclination to be 82.6 +/- 0.3 deg. We calculate the mass of the primary white dwarf to be 0.85 +/- 0.04 solar masses and the donor to be 0.035 +/- 0.003 solar masses, implying a partially degenerate state for this component. We observe superhump variations that are characteristic of an elliptical, precessing accretion disc. Our determination of the superhump period excess is in agreement with the established relationship between this parameter and the mass ratio, and is the most precise calibration of this relationship at low q. We also observe a quasi-periodic oscillation in the 2006 data, and we examine the outbursting behaviour of the system over a 4.5 year period.

This paper describes the aims, objectives and first results of the observational program for the study of distant core-collapse supernovae (SNe) with redshifts z < 0.3. This work is done within the framework of an international cooperation program on the SNe monitoring at the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences, and other telescopes. We study both the early phases of events (SN type determination, redshift estimation, and a search for manifestations of a wind envelope), and the nebular phase (the effects of explosion asymmetry). The SNe, associated with cosmic gamma-ray bursts are of particular interest. An interpretation of our observational data along with the data obtained on other telescopes is used to test the existing theoretical models of both the SN explosion, and the surrounding circumstellar medium. In 2009 we observed 30 objects; the spectra were obtained for 12 of them. We determined the types, phases after maximum, and redshifts for five SNe (SN 2009db, SN 2009dy, SN 2009dw, SN 2009ew, SN 2009ji). Based on the obtained photometric data a discovery of two more SNe was confirmed (SN 2009bx and SN 2009cb). A study of two type II supernovae in the nebular phase (SN 2008gz and SN 2008in) is finalized, four more objects (SN 2008iy, SN 2009ay, SN 2009bw, SN 2009de) are currently monitored.

Stars with initial masses 10 M_{solar} < M_{initial} < 100 M_{solar} fuse progressively heavier elements in their centres, up to inert iron. The core then gravitationally collapses to a neutron star or a black hole, leading to an explosion -- an iron-core-collapse supernova (SN). In contrast, extremely massive stars (M_{initial} > 140 M_{solar}), if such exist, have oxygen cores which exceed M_{core} = 50 M_{solar}. There, high temperatures are reached at relatively low densities. Conversion of energetic, pressure-supporting photons into electron-positron pairs occurs prior to oxygen ignition, and leads to a violent contraction that triggers a catastrophic nuclear explosion. Tremendous energies (>~ 10^{52} erg) are released, completely unbinding the star in a pair-instability SN (PISN), with no compact remnant. Transitional objects with 100 M_{solar} < M_{initial} < 140 M_{solar}, which end up as iron-core-collapse supernovae following violent mass ejections, perhaps due to short instances of the pair instability, may have been identified. However, genuine PISNe, perhaps common in the early Universe, have not been observed to date. Here, we present our discovery of SN 2007bi, a luminous, slowly evolving supernova located within a dwarf galaxy (~1% the size of the Milky Way). We measure the exploding core mass to be likely ~100 M_{solar}, in which case theory unambiguously predicts a PISN outcome. We show that >3 M_{solar} of radioactive 56Ni were synthesized, and that our observations are well fit by PISN models. A PISN explosion in the local Universe indicates that nearby dwarf galaxies probably host extremely massive stars, above the apparent Galactic limit, perhaps resulting from star formation processes similar to those that created the first stars in the Universe.

We report on the discovery and initial observations of the energetic type IIn supernova (SN), 2008fz. The optical energy emitted by SN 2008fz (based on the light curve over a 88 day period), is possibly the most ever observed for a supernova (1.4 x 10^51 erg). The event was more luminous than the type IIn SN 2006gy, but exhibited same smooth, slowly evolving light curve. As is characteristic of type IIn SN, the early spectra of 2008fz initially exhibited narrow Balmer lines which were replaced by a broader component at later times. The spectra also show a blue continuum with no signs of Ca or Na absorption, suggesting that there is little extinction due to intragalatic dust in the host or circumstellar material. No host galaxy is identified in prior coadded images reaching R ~ 22. From the supernova's redshift, z=0.133, we place an upper limit on the host of M_R=-17. The presence of the SN within such a faint host follows the majority of recently discovered highly luminous SN. A possible reason for this occurrence is the very high star formation rate occurring in low-mass galaxies in combination with the low metallicity environment, which makes the production of very massive stars possible. We determine the peak absolute magnitude of the event to be M_V = -22.3 from the initial photometry and the redshift distance, placing it among the most luminous supernovae discovered.

By cross matching blue objects from SDSS with GALEX and the astrometric catalogues USNO-B1.0, GSC2.3 and CMC14, 64 new dwarf nova candidates with one or more observed outbursts have been identified. 14 of these systems are confirmed as cataclysmic variables through existing and follow-up spectroscopy. A study of the amplitude distribution and an estimate of the outburst frequency of these new dwarf novae and those discovered by the Catalina Real-time Transient Survey (CRTS) indicates that besides systems that are faint because they are farther away, there also exists a population of intrinsically faint dwarf novae with rare outbursts.

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables.

The Palomar Transient Factory (PTF) is a wide-field experiment designed to investigate the optical transient and variable sky on time scales from minutes to years. PTF uses the CFH12k mosaic camera, with a field of view of 7.9 deg^2 and a plate scale of 1 asec/pixel, mounted on the the Palomar Observatory 48-inch Samuel Oschin Telescope. The PTF operation strategy is devised to probe the existing gaps in the transient phase space and to search for theoretically predicted, but not yet detected, phenomena, such as fallback supernovae, macronovae, .Ia supernovae and the orphan afterglows of gamma-ray bursts. PTF will also discover many new members of known source classes, from cataclysmic variables in their various avatars to supernovae and active galactic nuclei, and will provide important insights into understanding galactic dynamics (through RR Lyrae stars) and the Solar system (asteroids and near-Earth objects). The lessons that can be learned from PTF will be essential for the preparation of future large synoptic sky surveys like the Large Synoptic Survey Telescope. In this paper we present the scientific motivation for PTF and describe in detail the goals and expectations for this experiment.

Skyalert.org is a web application to collect and disseminate observations about time-critical astronomical transients, and to add annotations and intelligent machine-learning to those observations. The information is "pushed" to subscribers, who may be either humans (email, text message etc) or they may be machines that control telescopes. Subscribers can prepare precise "trigger rules" to decide which events should reach them and their robots, rules that may be based on sky position, or on the specific vocabulary of parameters that define a particular type of observation. Our twin thrusts are automation of process, and discrimination of interesting events.

We present Hubble Space Telescope (HST) WFPC2 photometry of 13 microlensed source stars from the 5.7 year Large Magellanic Cloud (LMC) survey conducted by the MACHO Project. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. None of these sources is coincident with a background galaxy, which rules out the possibility that the MACHO LMC microlensing sample is contaminated with misidentified supernovae or AGN in galaxies behind the LMC. This supports the conclusion that the MACHO LMC microlensing sample has only a small amount of contamination due to non-microlensing forms of variability. We compare the WFPC2 source star magnitudes with the lensed flux predictions derived from microlensing fits to the light curve data. In most cases the source star brightness is accurately predicted. Finally, we develop a statistic which constrains the location of the Large Magellanic Cloud (LMC) microlensing source stars with respect to the distributions of stars and dust in the LMC and compare this to the predictions of various models of LMC microlensing. This test excludes at > 90% confidence level models where more than 80% of the source stars lie behind the LMC. Exotic models that attempt to explain the excess LMC microlensing optical depth seen by MACHO with a population of background sources are disfavored or excluded by this test. Models in which most of the lenses reside in a halo or spheroid distribution associated with either the Milky Way or the LMC are consistent which these data, but LMC halo or spheroid models are favored by the combined MACHO and EROS microlensing results.

Digital synoptic sky surveys pose several new object classification challenges. In surveys where real-time detection and classification of transient events is a science driver, there is a need for an effective elimination of instrument-related artifacts which can masquerade as transient sources in the detection pipeline, e.g., unremoved large cosmic rays, saturation trails, reflections, crosstalk artifacts, etc. We have implemented such an Artifact Filter, using a supervised neural network, for the real-time processing pipeline in the Palomar-Quest (PQ) survey. After the training phase, for each object it takes as input a set of measured morphological parameters and returns the probability of it being a real object. Despite the relatively low number of training cases for many kinds of artifacts, the overall artifact classification rate is around 90%, with no genuine transients misclassified during our real-time scans. Another question is how to assign an optimal star-galaxy classification in a multi-pass survey, where seeing and other conditions change between different epochs, potentially producing inconsistent classifications for the same object. We have implemented a star/galaxy multipass classifier that makes use of external and a priori knowledge to find the optimal classification from the individually derived ones. Both these techniques can be applied to other, similar surveys and data sets.

Exploration of time domain is now a vibrant area of research in astronomy, driven by the advent of digital synoptic sky surveys. While panoramic surveys can detect variable or transient events, typically some follow-up observations are needed; for short-lived phenomena, a rapid response is essential. Ability to automatically classify and prioritize transient events for follow-up studies becomes critical as the data rates increase. We have been developing such methods using the data streams from the Palomar-Quest survey, the Catalina Sky Survey and others, using the VOEventNet framework. The goal is to automatically classify transient events, using the new measurements, combined with archival data (previous and multi-wavelength measurements), and contextual information (e.g., Galactic or ecliptic latitude, presence of a possible host galaxy nearby, etc.); and to iterate them dynamically as the follow-up data come in (e.g., light curves or colors). We have been investigating Bayesian methodologies for classification, as well as discriminated follow-up to optimize the use of available resources, including Naive Bayesian approach, and the non-parametric Gaussian process regression. We will also be deploying variants of the traditional machine learning techniques such as Neural Nets and Support Vector Machines on datasets of reliably classified transients as they build up.

We report on the results from the first six months of the Catalina Real-time Transient Survey (CRTS). In order to search for optical transients with timescales of minutes to years, the CRTS analyses data from the Catalina Sky Survey which repeatedly covers twenty six thousand of square degrees on the sky. The CRTS provides a public stream of transients that are bright enough to be followed up using small telescopes. Since the beginning of the survey, all CRTS transients have been made available to astronomers around the world in real-time using HTML tables, RSS feeds and VOEvents. As part of our public outreach program the detections are now also available in KML through Google Sky. The initial discoveries include over 350 unique optical transients rising more than two magnitudes from past measurements. Sixty two of these are classified as supernovae, based on light curves, prior deep imaging and spectroscopic data. Seventy seven are due to cataclysmic variables (only 13 previously known), while an additional 100 transients were too infrequently sampled to distinguish between faint CVs and SNe. The remaining optical transients include AGN, Blazars, high proper motions stars, highly variable stars (such as UV Ceti stars) and transients of an unknown nature. Our results suggest that there is a large population of SNe missed by many current supernova surveys because of selection biases. These objects appear to be associated with faint host galaxies. We also discuss the unexpected discovery of white dwarf binary systems through dramatic eclipses.

There is an increasing number of large, digital, synoptic sky surveys, in which repeated observations are obtained over large areas of the sky in multiple epochs. Likewise, there is a growth in the number of (often automated or robotic) follow-up facilities with varied capabilities in terms of instruments, depth, cadence, wavelengths, etc., most of which are geared toward some specific astrophysical phenomenon. As the number of detected transient events grows, an automated, probabilistic classification of the detected variables and transients becomes increasingly important, so that an optimal use can be made of follow-up facilities, without unnecessary duplication of effort. We describe a methodology now under development for a prototype event classification system; it involves Bayesian and Machine Learning classifiers, automated incorporation of feedback from follow-up observations, and discriminated or directed follow-up requests. This type of methodology may be essential for the massive synoptic sky surveys in the future.

We describe briefly the Palomar-Quest (PQ) digital synoptic sky survey, including its parameters, data processing, status, and plans. Exploration of the time domain is now the central scientific and technological focus of the survey. To this end, we have developed a real-time pipeline for detection of transient sources. We describe some of the early results, and lessons learned which may be useful for other, similar projects, and time-domain astronomy in general. Finally, we discuss some issues and challenges posed by the real-time analysis and scientific exploitation of massive data streams from modern synoptic sky surveys.

This paper has been withdrawn

We review some of the recent developments and challenges posed by the data analysis in modern digital sky surveys, which are representative of the information-rich astronomy in the context of Virtual Observatory. Illustrative examples include the problems of an automated star-galaxy classification in complex and heterogeneous panoramic imaging data sets, and an automated, iterative, dynamical classification of transient events detected in synoptic sky surveys. These problems offer good opportunities for productive collaborations between astronomers and applied computer scientists and statisticians, and are representative of the kind of challenges now present in all data-intensive fields. We discuss briefly some emergent types of scalable scientific data analysis systems with a broad applicability.

We present the results from the discovery and study of ~3000 chromospherically active giant and subgiant stars toward the Galactic bulge. We find that these stars are predominantly RS CVn binaries with rotation periods between 10 and 100 days. We discover that the average rotational period of these stars decreases with their distance from the Galactic plane. We find that the primary stars in the RS CVn systems are predominantly first-ascent giants. Our research also suggests that, if these stars have spot cycles like the sun, then the cycle period must be longer than 10 years on average. We confirm that the amplitude of the spot-induced modulations observed in the light curves of these objects is generally larger at minimum light than at maximum. Furthermore, we confirm that the amplitudes of the modulations due to stellar spots generally increases as the observed change in average brightness increases. We find no evidence for a relationship between a CA star's brightness and its rotational period. However, the average period does increase with colour for stars with periods >~30 days.

We present a measurement of the systemic proper motion of the Large Magellanic Cloud (LMC) from astrometry with the High Resolution Camera (HRC) of the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). We observed LMC fields centered on 21 background QSOs that were discovered from their optical variability in the MACHO database. The QSOs are distributed homogeneously behind the central few degrees of the LMC. With 2 epochs of HRC data and a ~2 year baseline we determine the proper motion of the LMC to better than 5% accuracy: mu_W = -2.03 +/- 0.08 mas/yr; mu_N = 0.44 +/- 0.05 mas/yr. This is the most accurate proper motion measurement for any Milky Way satellite thus far. When combined with HI data from the Magellanic Stream this should provide new constraints on both the mass distribution of the Galactic Halo and models of the Stream.

Using 7 years of MACHO survey data, we present a new determination of the optical depth to microlensing towards the Galactic bulge. We select the sample of 62 microlensing events (60 unique) on clump giant sources and perform a detailed efficiency analysis. We use only the clump giant sources because these are bright bulge stars and are not as strongly affected by blending as other events. Using a subsample of 42 clump events concentrated in just 4.5 square degrees, we find tau = 2.17^{+0.47}_{-0.38} x 10^{-6} at (l,b) = (1.50, -2.68), somewhat smaller than found in most previous MACHO studies, but in excellent agreement with recent theoretical predictions. We also present the optical depth in each of the 19 fields in which we detected events, and find limits on optical depth for fields with no events. The errors in optical depth in individual fields are dominated by Poisson noise. We measure optical depth gradients (1.06 +/- 0.71) x 10^{-6} deg^{-1} and (0.29 +/- 0.43) x 10^{-6} deg^{-1} in the galactic latitude b and longitude l directions, respectively. Finally, we discuss the possibility of anomalous duration distribution of events in the field 104 centered on (l,b) = (3.11, -3.01) as well as investigate spatial clustering of events in all fields.

We present a catalog of 450 high signal-to-noise microlensing events observed by the MACHO collaboration between 1993 and 1999. The events are distributed throughout our fields and, as expected, they show clear concentration toward the Galactic center. No optical depth is given for this sample since no blending efficiency calculation has been performed, and we find evidence for substantial blending. In a companion paper we give optical depths for the sub-sample of events on clump giant source stars, where blending is not a significant effect. Several events with sources that may belong to the Sagittarius dwarf galaxy are identified. For these events even relatively low dispersion spectra could suffice to classify these events as either consistent with Sagittarius membership or as non-Sagittarius sources. Several unusual events, such as microlensing of periodic variable source stars, binary lens events, and an event showing extended source effects are identified. We also identify a number of contaminating background events as cataclysmic variable stars.

We present the results from an analysis of Hubble Space Telescope High Resolution Camera data for the Large Magellanic Cloud microlensing event MACHO-LMC-5. By determining the parallax and proper motion of this object we find that the lens is an M dwarf star at a distance of 578(+65/-53)pc with a proper motion of 21.39 +/- 0.04 mas/yr. Based on the kinematics and location of this star it more likely to be part of the Galactic thick disk than thin disk population. We confirm that the microlensing event LMC-5 is a jerk-parallax event.

Fourier coefficents have been derived for the $V$ and $R$ light curves of 785 overtone RR Lyrae variables in 16 MACHO fields near the bar of the LMC. The $\phi_{31}$ and $R_{21}$ coefficients have been compared with those of the first overtone RR Lyrae variables in the Galactic globular clusters NGC 6441, M107, M5, M3, M2, $\omega$ Centauri and M68. The results indicate that many of the LMC variables have properties similar to the ones in M2, M3, M5 and the Oosterhoff type I variables in $\omega$ Cen, but they are different from the Oosterhoff type II variables in $\omega$ Cen. Equations derived from hydrodynamic pulsation models have been used to calculate the luminosity and temperature for the 330 bona fide first-overtone variables. The results indicate that they have $\log L$ in the range 1.6 to $1.8\lsun$ and $\log T_{eff}$ between 3.85 and 3.87. Based on these temperatures, a mean color excess $E(V-R) =0.08$ mag, equivalent to $E(B-V)=0.14$ mag, has been estimated for these 330 stars. The 80 M5-like variables (selected according to their location in the $\phi_{31}-\log P$ plot) are used to determine a LMC distance. After correcting for the effects of extinction and crowding, a mean apparent magnitude $<V_0>=18.99 \pm 0.02$ (statistical) $\pm 0.16$ (systematic) has been estimated for these 80 stars. Combining this with a mean absolute magnitude $M_V=0.56\pm 0.06$ for M5-like stars derived from Baade-Wesselink analyses, main sequence fitting, Fourier parameters and the trigonometric parallax of RR Lyrae, we derive an LMC distance modulus $\mu=18.43\pm 0.06$ (statistical) $\pm 0.16$ (systematic) mag. The large systematic error arises from the difficulties of correcting for interstellar extinction and for crowding.

We have frequency analyzed 6391 variables classified earlier as fundamental mode RR Lyrae (RR0) stars in the MACHO database on the Large Magellanic Cloud (LMC). The overwhelming majority of these variables have been proved to be indeed RR0 stars, whereas the remaining ones have fallen in one of the following categories: single- and double-mode Cepheids, binaries, first overtone and double-mode RR Lyrae stars and non-classified variables. Special attention has been paid to the Blazhko stars. We found altogether 731 Blazhko variables showing either a doublet or an equidistant triplet pattern at the main pulsation component in their frequency spectra. The incidence rate of the Blazhko variables among the RR0 stars in the LMC is 11.9%, which is three times higher than their rate among the first overtone RR Lyrae stars. From the frequency spectra we found that variables with larger modulation amplitudes at the higher frequency side of the main pulsation component are three times more frequent than the ones showing the opposite amplitude pattern. A search made for a modulation component with the Blazhko period in the average brightness of the individual variables showed the existence of such a modulation with an overall amplitude of ~0.006 mag. On the other hand, a similar search for quadruple modulation patterns around the main pulsation component have failed to clearly detect such components at the ~0.004 mag level. This means that the amplitudes of the quadruple components should be, on the average, at least ten times smaller than those of the triplet components. This finding and the existence of Blazhko variables with highly asymmetric modulation amplitudes not only question the validity of the magnetic oblique rotator model, but also puts stringent constraints on models based on mode coupling theories.

We describe a few recent microlensing results from the MACHO Collaboration. The aim of the MACHO Project was the identification and quantitative description of dark and luminous matter in the Milky Way using microlensing toward the Magellanic Clouds and Galactic bulge. We start with a discussion of the HST follow-up observations of the microlensing events toward the LMC detected in the first 5 years of the experiment. Using color-magnitude diagrams we attempt to distinguish between two possible locations of the microlensing sources: 1) in the LMC or 2) behind the LMC. We conclude that unless the extinction is extremely patchy, it is very unlikely that most of the LMC events have sources behind the LMC. During an examination of the HST images of the 13 LMC events we found a very red object next to the source star of event LMC-5. Astrometry, microlensing parallax fit, and a spectrum suggest that in this case we directly image the lens - a low-mass disk star. Then we focus on the majority of events observed by the MACHO Project, which are detected toward the Galactic bulge. We argue that the microlensing optical depth toward the bulge is best measured using events that have clump giant sources, which are almost unaffected by blending. From this sample we derive a low optical depth toward the Galactic bulge of (1.4 +/- 0.3) x 10^{-6}, in good agreement with other observational constraints and with theoretical models. The presence of many long-duration events among the bulge candidates allows us to investigate the microlensing parallax effect. Events with the strongest parallax signal are probably due to massive remnants. Events MACHO-96-BLG-5 and MACHO-98-BLG-6 might have been caused by the 6-solar-mass black holes.

We present a new method for differentiating between planetary transits and eclipsing binaries based on the presence of the ellipsoidal light variations. These variations can be used to detect stellar secondaries with masses ~0.2 M_sun orbiting sun-like stars at a photometric accuracy level which has already been achieved in transit surveys. By removing candidates exhibiting this effect it is possible to greatly reduce the number of objects requiring spectroscopic follow up with large telescopes. Unlike the usual candidate selection method, which are primarily based on the estimated radius of the orbiting object, this technique is not biased against bona-fide planets and brown dwarfs with large radii, because the amplitude of the effect depends on the transiting object's mass and orbital distance. In many binary systems, where a candidate planetary transit is actually due to the partial eclipse of two normal stars, the presence of flux variations due to the gravity darkening effect will show the true nature of these systems. We show that many of the recent OGLE-III photometric transit candidates exhibit the presence of significant variations in their light curves and are likely to be due to stellar secondaries. We find that the light curves of white dwarf transits will generally not mimic those of small planets because of significant gravitationally induced flux variations. We discuss the relative merits of methods used to detect transit candidates which are due to stellar blends rather than planets. We outline how photometric observations taken in two bands can be used to detect the presence of stellar blends.

The recent detections of a large population of faint submillimetre sources, an excess halo gamma-ray background, and the extreme scattering events observed for extragalactic radio sources have been explained as being due to baryonic dark matter in the form of small, dark, gas clouds. In this paper we present the results of a search for the transient stellar obscurations such clouds are expected to cause. We examine the Macho project light curves of 48 x 10^6 stars toward the Galactic bulge, LMC and SMC for the presence of dark cloud extinction events. We find no evidence for the existence of a population of dark gas clouds with Av > 0.2 and masses between ~ 10^-4 and 10^-2 M_solar in the Galactic disk or halo. However, it is possible that such dark cloud populations could exist if they are clustered in regions away from the observed lines of sight.

We present 47 spectroscopically-confirmed quasars discovered behind the Magellanic Clouds identified via photometric variability in the MACHO database. Thirty-eight quasars lie behind the Large Magellanic Cloud and nine behind the Small Magellanic Cloud, more than tripling the number of quasars previously known in this region. The quasars cover the redshift interval 0.2 < z < 2.8 and apparent mean magnitudes 16.6 < V < 20.1. We discuss the details of quasar candidate selection based on time variability in the MACHO database and present results of spectroscopic follow-up observations. Our follow-up detection efficiency was 20%; the primary contaminants were emission-line Be stars in Magellanic Clouds. For the 47 quasars discovered behind the Magellanic Clouds plus an additional 12 objects previously identified in this region, we present 7.5-year MACHO V- and R-band lightcurves with average sampling times of 2-10 days.

We present an analysis of the longest timescale microlensing events discovered by the MACHO Collaboration during a 7 year survey of the Galactic bulge. We find 6 events that exhibit very strong microlensing parallax signals due, in part, to accurate photometric data from the GMAN and MPS collaborations. The microlensing parallax fit parameters are used in a likelihood analysis, which is able to estimate the distance and masses of the lens objects based upon a standard model of the Galactic velocity distribution. This analysis indicates that the most likely masses of 5 of the 6 lenses are > 1 Msun, which suggests that a substantial fraction of the Galactic lenses may be massive stellar remnants. This could explain the observed excess of long timescale microlensing events. The lenses for events MACHO-96-BLG-5 and MACHO-98-BLG-6 are the most massive, with mass estimates of M/Msun = 6 +10/-3 and M/Msun = 6 +7/-3, respectively. The observed upper limits on the absolute brightness of main sequence stars for these lenses are < 1 Lsun, so both lenses are black hole candidates. The black hole interpretation is also favored by a likelihood analysis with a Bayesian prior using a conventional model for the lens mass function. We consider the possibility that the source stars for some of these 6 events may lie in the foreground or background of the Galactic bulge, but we find that this is unlikely. Future HST observations of these events can either confirm the black hole lens hypothesis or detect the lens stars and provide a direct measurement of their masses. Future observations of similar events by SIM or the Keck or VLTI interferometers will allow direct measurements of the lens masses for stellar remnant lenses as well.

We present a method for solving the lightcurve of an eclipsing binary system which contains a Cepheid variable as one of its components as well as the solutions for three eclipsing Cepheids in the Large Magellanic Cloud (LMC). A geometric model is constructed in which the component stars are assumed to be spherical and on circular orbits. The emergent system flux is computed as a function of time, with the intrinsic variations in temperature and radius of the Cepheid treated self-consistently. Fitting the adopted model to photometric observations, incorporating data from multiple bandpasses, yields a single parameter set best describing the system. This method is applied to three eclipsing Cepheid systems from the MACHO Project LMC database: MACHO ID's 6.6454.5, 78.6338.24 and 81.8997.87. A best-fit value is obtained for each system's orbital period and inclination and for the relative radius, color and limb-darkening coefficients of each star. Pulsation periods and parameterizations of the intrinsic color variations of the Cepheids are also obtained and the amplitude of the radial pulsation of each Cepheid is measured directly. The system 6.6454.5 is found to contain a 4.97-day Cepheid, which cannot be definitely classified as Type I or Type II, with an unexpectedly brighter companion. The system 78.6338.24 consists of a 17.7-day, W Vir Class Type II Cepheid with a smaller, dimmer companion. The system 81.8997.87 contains an intermediate-mass, 2.03-day overtone Cepheid with a dimmer, red giant secondary.

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. This class of events allows one to obtain robust conclusions because relatively bright clump stars are not strongly affected by blending. We discuss: 1) the selection of `giant' events, 2) the distribution of event durations, 3) the anomalous character of event durations and optical depth in the MACHO field 104 centered on (l,b) = (3.1, -3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9, -3.8), and present a map of the spatial distribution of the optical depth. When field 104 is removed from the sample, the optical depth drops to (1.4 +/- 0.3) x 10^{-6}, which is in excellent agreement with infrared-based models of the central Galactic region.

We report the discovery of eight new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The discovery of these new stars increases the number of known RCB stars in the LMC to thirteen. We have also discovered four stars similar to the Galactic variable DY Per. These stars decline much more slowly and are cooler than the RCB stars. The absolute luminosities of the Galactic RCB stars are unknown since there is no direct measurement of the distance to any Galactic RCB star. Hence, the importance of the LMC RCB stars. We find a much larger range of absolute magnitudes (M(V) = -2.5 to -5 mag) than inferred from the small pre-MACHO sample of LMC RCB stars. It is likely that there is a temperature - M(V)relationship with the cooler stars being intrinsically fainter. Cool (~5000 K) RCB stars are much more common than previously thought based on the Galactic RCB star sample. Using the fairly complete sample of RCB stars discovered in the MACHO fields, we have estimated the likely number of RCB stars in the Galaxy to be ~3,200. The SMC MACHO fields were also searched for RCB stars but none were found.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

We attempt to determine whether the MACHO microlensing source stars are drawn from the average population of the LMC or from a population behind the LMC by examining the HST color-magnitude diagram (CMD) of microlensing source stars. We present WFPC2 HST photometry of eight MACHO microlensing source stars and the surrounding fields in the LMC. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. We consider in detail a model for the background population of source stars based on that presented by Zhao, Graff & Guhathakurta. In this model, the source stars have an additional reddening <E(B-V)> = 0.13 mag and a slightly larger distance modulus <Delta mu> ~ 0.3 mag than the average LMC population. We also investigate a series of source star models, varying the relative fraction of source stars drawn from the average and background populations and the displacement of the background population from the LMC. Due to the small number of analyzed events the distribution of probabilities of different models is rather flat. A shallow maximum occurs at a fraction s_LMC ~ 0.8 of the source stars in the LMC. This is consistent with the interpretation that a significant fraction of observed microlensing events are due to lenses in the Milky Way halo, but does not definitively exclude other models.

We present photometry and analysis of the microlensing alert MACHO 96-LMC-2. The ~3% photometry provided by the Global Microlensing Alert Network follow--up effort reveals a periodic modulation in the lightcurve. We attribute this to binarity of the lensed source. Microlensing fits to a rotating binary source magnified by a single lens converge on two minima, separated by delta chi^2 ~ 1. The most significant fit X1 predicts a primary which contributes ~100% of the light, a dark secondary, and an orbital period (T) of 9.2 days. The second fit X2 yields a binary source with two stars of roughly equal mass and luminosity, and T = 21.2 days. The lensed object appears to lie on the upper LMC main sequence. We estimate the mass of the primary component of the binary system, M ~2 M_sun. For the preferred model X1, we explore the range of dark companions by assuming 0.1 M_sun and 1.4 M_sun objects in models X1a and X1b, respectively. We find lens velocities projected to the LMC in these models of v^hat_X1a = 18.3 +/- 3.1 km/s and v^hat_X1b = 188 +/- 32 k/ms. In both these cases, a likelihood analysis suggests an LMC lens is preferred over a Galactic halo lens, although only marginally so in model X1b. We also find v^hat_X2 = 39.6 +/- 6.1 k/ms, where the likelihood for the lens location is strongly dominated by the LMC disk. In all cases, the lens mass is consistent with that of an M-dwarf. The LMC self-lensing rate contributed by 96-LMC-2 is consistent with model self-lensing rates. (Abridged)

We report on a search for long duration microlensing events towards the Large Magellanic Cloud. We find none, and therefore put limits on the contribution of high mass objects to the Galactic dark matter. At 95% confidence level we exclude objects in the mass range 0.3 solar masses to 30.0 solar masses from contributing more than 4 times 10^11 solar masses to the Galactic halo. Combined with earlier results, this means that objects with masses under 30 solar masses cannot make up the entire dark matter halo if the halo is of typical size. For a typical dark halo, objects with masses under 10 solar masses contribute less than 40% of the dark matter.

We present a ~ 5-yr optical light curve of the recurrent Be/X-ray transient A0538-66 obtained as a by-product of the MACHO Project. These data reveal both a long-term modulation at P = 420.8 +/- 0.8 d and a short-term modulation at 16.6510 +/- 0.0022 d which, within errors, confirms the previously found orbital period. Furthermore, the orbital activity is only seen at certain phases of the 421 d cycle suggesting that the long-term modulation is related to variations in the Be star envelope.

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. In particular, we discuss: 1) the selection of `giant' events, 2) distribution of impact parameters, 3) distribution of event durations, 4) the concentration of long duration events in MACHO field 104 centered on (l,b) = (3.1,-3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9,-3.8). We discuss future work and prospects for building a coherent model of the Galaxy.

Using the exceptional long-term monitoring capabilities of the MACHO project, we present here the optical history of LMC X-2 for a continuous 6-yr period. These data were used to investigate the previously claimed periodicities for this source of 8.15 h and 12.54 d : we find upper amplitude limits of 0.10 mag and 0.09 mag, respectively.

More than 1300 variables classified provisionally as first overtone RR Lyrae pulsators in the MACHO variable star database of the Large Magellanic Cloud (LMC) have been subjected to standard frequency analysis. Based on the remnant power in the prewhitened spectra, we found 70% of the total population to be monoperiodic. The remaining 30% (411 stars) are classified as one of 9 types according to their frequency spectra. Several types of RR Lyrae pulsational behavior are clearly identified here for the first time. Together with the earlier discovered double-mode (fundamental & first overtone) variables this study increased the number of the known double-mode stars in the LMC to 181. During the total 6.5yr time span of the data, 10% of the stars show strong period changes. We also discovered two additional types of multifrequency pulsators with low occurrence rates of 2% for each. In the first type there remains one closely spaced component after prewhitening by the main pulsation frequency. In the second type the number of remnant components is two, they are also closely spaced, and, in addition, they are symmetric in their frequency spacing relative to the central component. This latter type of variables is associated with their relatives among the fundamental pulsators, known as Blazhko variables. Their high frequency (~20%) among the fundamental mode variables versus the low occurrence rate of their first overtone counterparts makes it more difficult to explain Blazhko phenomenon by any theory depending mainly on the role of aspect angle or magnetic field. Current theoretical models invoke nonradial pulsation components in these stars.

We present the microlensing optical depth towards the Galactic bulge based on the detection of 99 events found in our Difference Image Analysis (DIA) survey. This analysis encompasses three years of data, covering ~ 17 million stars in ~ 4 deg^2, to a source star baseline magnitude limit of V = 23. The DIA technique improves the quality of photometry in crowded fields, and allows us to detect more microlensing events with faint source stars. We find this method increases the number of detection events by 85% compared with the standard analysis technique. DIA light curves of the events are presented and the microlensing fit parameters are given. The total microlensing optical depth is estimated to be tau_(total)= 2.43^(+0.39/-0.38) x 10^(-6) averaged over 8 fields centered at l=2.68 and b=-3.35. For the bulge component we find tau_(bulge)=3.23^(+0.52/-0.50) x 10^(-6) assuming a 25% stellar contribution from disk sources. These optical depths are in good agreement with the past determinations of the MACHO Alcock et al. (1997) and OGLE Udalski et al. (1994) groups, and are higher than predicted by contemporary Galactic models. We show that our observed event timescale distribution is consistent with the distribution expected from normal mass stars, if we adopt the stellar mass function of Scalo (1986) as our lens mass function. However, we note that as there is still disagreement about the exact form of the stellar mass function, there is uncertainty in this conclusion. Based on our event timescale distribution we find no evidence for the existence of a large population of brown dwarfs in the direction of the Galactic bulge.

The MACHO project is a search for dark matter in the form of massive compact halo objects (MACHOs). The project has photometrically monitored tens of millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge in search of rare gravitational microlensing events caused by these otherwise invisible objects. In 5.7 years of observations toward the LMC some 13-17 microlensing events have been observed by the MACHO survey, allowing powerful statements to be made about the nature of the dark population in the halo of our Galaxy. A critical component of these statements is an accurate determination of the survey's detection efficiency. The detection efficiency is a complicated function of temporal sampling, stellar crowding (the luminosity function), image quality, photometry, time-series analysis, and criteria used to select the microlensing candidates. Such a complex interdependence is most naturally solved using a Monte Carlo approach. Here we describe the details of the Monte Carlo used to calculate the efficiency presented in the MACHO 5.7-year LMC results. Here we correct several shortcomings of past determinations, including (1) adding fainter source stars (2.5 magnitudes below our faintest detected "stars"), (2) an up-to-date luminosity function for the LMC, (3) better sampling of real images in both stellar density and observing conditions, (4) an improved scheme for adding artificial microlensing onto a random sample of real lightcurves, and many other improvements. [Abridged]

We present a 9 million star color-magnitude diagram (9M CMD) of the LMC bar. The 9M CMD reveals a complex superposition of different age and metallicity stellar populations, with important stellar evolutionary phases occurring over 3 orders of magnitude in number density. First, we count the non-variable supergiants, the associated Cepheids, and measure the effective temperatures defining the instability strip. Lifetime predictions of stellar evolution theory are tested, with implications for the origin of low-luminosity Cepheids. The highly-evolved AGB stars have a bimodal distribution in brightness, which we interpret as discrete old populations (>1 Gyr). The faint AGB may be metal-poor and very old. We identify the clusters NGC 411 and M3 as templates for the admixture of old stellar populations. However, there are indications that the old and metal-poor field population has a red HB morphology: the RR Lyraes lie on the red edge of the instability strip, the AGB-bump is very red, and the ratio of AGB-bump stars to RR Lyraes is quite large. If the HB second parameter is age, the old and metal-poor field population likely formed after the oldest clusters. Lifetime predictions of stellar evolution theory lead us to associate a significant fraction of the red HB clump giants with the same old and metal-poor population producing the RR Lyraes and the AGB-bump. In this case, compared to the age-dependent luminosity predictions of stellar evolution theory, the red HB clump is too bright relative to the RR Lyraes and AGB-bump. Last, the surface density profile of RR Lyraes is fit by an exponential, favoring a disk-like rather than spheroidal distribution. We conclude that the age of the LMC disk is probably similar to the age of the Galactic disk. (ABRIDGED)

We report on our search for microlensing towards the Large Magellanic Cloud (LMC). Analysis of 5.7 years of photometry on 11.9 million stars in the LMC reveals 13 - 17 microlensing events. This is significantly more than the $\sim$ 2 to 4 events expected from lensing by known stellar populations. The timescales ($\that$) of the events range from 34 to 230 days. We estimate the microlensing optical depth towards the LMC from events with $2 < \that < 400$ days to be 1.2 ^{+0.4}_ {-0.3} \ten{-7}$, with an additional 20% to 30% of systematic error. The spatial distribution of events is mildly inconsistent with LMC/LMC disk self-lensing, but is consistent with an extended lens distribution such as a Milky Way or LMC halo. Interpreted in the context of a Galactic dark matter halo, consisting partially of compact objects, a maximum likelihood analysis gives a MACHO halo fraction of 20% for a typical halo model with a 95% confidence interval of 8% to 50%. A 100% MACHO halo is ruled out at the 95% C.L. for all except our most extreme halo model. Interpreted as a Galactic halo population, the most likely MACHO mass is between $ 0.15 \msun$ and $ 0.9 \msun$, depending on the halo model, and the total mass in MACHOs out to 50 kpc is found to be 9+4-3 10^{10} msun, independent of the halo model. These results are marginally consistent with our previous results, but are lower by about a factor of two. Besides a larger data set, this work also includes an improved efficiency determination, improved likelihood analysis, and more thorough testing of systematic errors, especially with respect to the treatment of potential backgrounds to microlensing, such as supernovae in galaxies behind the LMC. [Abridged]

The MACHO Project is a microlensing survey that monitors the brightnesses of 60 million stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud, and Galactic bulge. Our database presently contains about 80 billion photometric measurements, a significant fraction of all astronomical photometry. We describe the calibration of MACHO two-color photometry and transformation to the standard Kron-Cousins V and R system. Calibrated MACHO photometry may be properly compared with all other observations on the Kron-Cousins standard system, enhancing the astrophysical value of these data. For 9 million stars in the LMC bar, independent photometric measurements of 20,000 stars with V < 18 mag in field-overlap regions demonstrate an internal precision sigV = 0.021, sigR = 0.019, and sig(V-R) = 0.028 mag. The accuracy of the zero-point in this calibration is estimated to be +-0.035 mag for stars with colors in the range -0.1 < (V-R) < 1.2 mag. A comparison of calibrated MACHO photometry with published photometric sequences and new Hubble Space Telescope observations shows agreement. The current calibration zero-point uncertainty for the remainder of the MACHO photometry database is estimated to be +- 0.10 mag in V or R and +-0.04 mag in (V-R). We describe the first application of calibrated MACHO photometry data: the construction of a color-magnitude diagram used to calculate our experimental sensitivity to detect microlensing in the LMC.

The MACHO collaboration has been carrying out Difference Image Analysis (DIA) since 1996 with the aim of increasing the sensitivity to the detection of gravitational microlensing. This is a preliminary report on the application of DIA to galactic bulge images in one field. We show how the DIA technique significantly increases the number of detected lensing events, by removing the positional dependence of traditional photometry schemes and lowering the microlensing event detection threshold. This technique, unlike PSF photometry, gives the unblended colours and positions of the microlensing source stars. We present a set of criteria for selecting microlensing events from objects discovered with this technique. The 16 pixel and classical microlensing events discovered with the DIA technique are presented.

This is a preliminary report on the application of Difference Image Analysis (DIA) to galactic bulge images. The aim of this analysis is to increase the sensitivity to the detection of gravitational microlensing. We discuss how the DIA technique simplifies the process of discovering microlensing events by detecting only objects which have variable flux. We illustrate how the DIA technique is not limited to detection of so called ``pixel lensing'' events, but can also be used to improve photometry for classical microlensing events by removing the effects of blending. We will present a method whereby DIA can be used to reveal the true unblended colours, positions and light curves of microlensing events. We discuss the need for a technique to obtain the accurate microlensing time scales from blended sources, and present a possible solution to this problem using the existing HST colour magnitude diagrams of the galactic bulge and LMC. The use of such a solution with both classical and pixel microlensing searches is discussed. We show that one of the major causes of systematic noise in DIA is differential refraction. A technique for removing this systematic by effectively registering images to a common airmass is presented. Improvements to commonly used image differencing techniques are discussed.

We present photometric observations and analysis of the second microlensing event detected towards the Small Magellanic Cloud (SMC), MACHO Alert 98-SMC-1. This event was detected early enough to allow intensive observation of the lightcurve. These observations revealed 98-SMC-1 to be the first caustic crossing, binary microlensing event towards the Magellanic Clouds to be discovered in progress. Frequent coverage of the evolving lightcurve allowed an accurate prediction for the date of the source crossing out of the lens caustic structure. The caustic crossing temporal width, along with the angular size of the source star, measures the proper motion of the lens with respect to the source, and thus allows an estimate of the location of the lens. Lenses located in the Galactic halo would have a velocity projected to the SMC of v^hat ~1500 km/s, while an SMC lens would typically have v^hat ~60 km/s. We have performed a joint fit to the MACHO/GMAN data presented here, including recent EROS data of this event. These joint data are sufficient to constrain the time for the lens to move an angle equal to the source angular radius; 0.116 +/- 0.010 days. We estimate a radius for the lensed source of 1.4 +/- 0.1 R_sun. This yields a projected velocity of v^hat = 84 +/- 9 km/s. Only 0.15% of halo lenses would be expected to have a v^hat value at least as small as this, while 31% of SMC lenses would be expected to have v^hat as large as this. This implies that the lensing system is more likely to reside in the SMC than in the Galactic halo.