Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis. Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm. Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $\rho_b=0.98_{-0.14}^{+0.10}\;\text{g cm}^{-3}$ and $\rho_c=0.494_{-0.077}^{+0.066}\;\text{g cm}^{-3}$.

We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package.

We use the LOFAR Two-metre Sky Survey (LoTSS) Data Release I to identify the groups of galaxies (and individual galaxies) from the Hickson Compact Groups and Magnitude Limited Compact Groups samples that emit at the frequency of 150\,MHz, characterise their radio emission (extended or limited to the galaxies), and compare new results to earlier observations and theoretical predictions. The detection of 73 systems (and 7 more -- probably) out of 120, of which as many as 17 show the presence of extended radio structures, confirms the previous hypothesis of the common character of the magnetic field inside galaxy groups and its detectability. In order to investigate the future potential of low-frequency radio studies of galaxy groups, we also present a more detailed insight into four radio-emitting systems, for which the strength of the magnetic field inside their intergalactic medium (IGM) is calculated. The estimated values are comparable to that found inside star-forming galaxies, suggesting a dynamical and evolutionary importance of the magnetic field in galaxy groups.

We present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive \textit{ab initio} simulations of hydrogen-helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal compressibility, and the sound velocity. Important transport properties such as the electrical and thermal conductivity, opacity, and shear viscosity are also calculated. Further results for associated quantities including magnetic and thermal diffusivity, kinematic shear viscosity, as well as the static Love number $k_2$ and the equidistance are presented. In comparison to Jupiter-mass planets, the behavior inside massive giant planets and brown dwarfs is stronger dominated by degenerate matter. We discuss the implications on possible dynamics and magnetic fields of those massive objects. The consistent data set compiled here may serve as starting point to obtain material and transport properties for other substellar H-He objects with masses above one Jovian mass and finally may be used as input for dynamo simulations.

The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

During its four years of photometric observations, the Kepler space telescope detected thousands of exoplanets and exoplanet candidates. One of Kepler's greatest heritages has been the confirmation and characterization of hundreds of multi-planet systems via Transit Timing Variations (TTVs). However, there are many interesting candidate systems displaying TTVs on such long time scales that the existing Kepler observations are of insufficient length to confirm and characterize them by means of this technique. To continue with Kepler's unique work we have organized the "Kepler Object of Interest Network" (KOINet). The goals of KOINet are, among others, to complete the TTV curves of systems where Kepler did not cover the interaction timescales well. KOINet has been operational since March, 2014. Here we show some promising first results obtained from analyzing seven primary transits of KOI-0410.01, KOI-0525.01, KOI-0760.01, and KOI-0902.01 in addition to Kepler data, acquired during the first and second observing seasons of KOINet. While carefully choosing the targets we set demanding constraints about timing precision (at least 1 minute) and photometric precision (as good as 1 part per thousand) that were achieved by means of our observing strategies and data analysis techniques. For KOI-0410.01, new transit data revealed a turn-over of its TTVs. We carried out an in-depth study of the system, that is identified in the NASA's Data Validation Report as false positive. Among others, we investigated a gravitationally-bound hierarchical triple star system, and a planet-star system. While the simultaneous transit fitting of ground and space-based data allowed for a planet solution, we could not fully reject the three-star scenario. New data, already scheduled in the upcoming 2018 observing season, will set tighter constraints on the nature of the system.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

We conducted a 12-month monitoring campaign of 33 T Tauri stars (TTS) in Taurus. Our goal was to monitor objects that possess a disk but have a weak Halpha line, a common accretion tracer for young stars, to determine whether they host a passive circumstellar disk. We used medium-resolution optical spectroscopy to assess the objects' accretion status and to measure the Halpha line. We found no convincing example of passive disks; only transition disk and debris disk systems in our sample are non-accreting. Among accretors, we find no example of flickering accretion, leading to an upper limit of 2.2% on the duty cycle of accretion gaps assuming that all accreting TTS experience such events. Combining literature results with our observations, we find that the reliability of traditional Halpha-based criteria to test for accretion is high but imperfect, particularly for low-mass TTS. We find a significant correlation between stellar mass and the full width at 10 per cent of the peak (W10%) of the Halpha line that does not seem to be related to variations in free-fall velocity. Finally, our data reveal a positive correlation between the Halpha equivalent width and its W10%, indicative of a systematic modulation in the line profile whereby the high-velocity wings of the line are proportionally more enhanced than its core when the line luminosity increases. We argue that this supports the hypothesis that the mass accretion rate on the central star is correlated with the Halpha W10% through a common physical mechanism.

We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than $K_p \approx 13$, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline, and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online.

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%.

We present a study of the statistical flare rates of M dwarfs (dMs) with close white dwarf (WD) companions (WD+dM; typical separations < 1 au). Our previous analysis demonstrated that dMs with close WD companions are more magnetically active than their field counterparts. One likely implication of having a close binary companion is increased stellar rotation through disk-disruption, tidal effects, and/or angular momentum exchange; increased stellar rotation has long been associated with an increase in stellar activity. Previous studies show a strong correlation between dMs that are magnetically active (showing H{\alpha} in emission) and the frequency of stellar flare rates. We examine the difference between the flare rates observed in close WD+dM binary systems and field dMs. Our sample consists of a subset of 181 close WD+dM pairs from Morgan et al. (2012) observed in the Sloan Digital Sky Survey Stripe 82, where we obtain multi-epoch observations in the Sloan ugriz-bands. We find an increase in the overall flaring fraction in the close WD+dM pairs (0.09$\pm$0.03%) compared to the field dMs (0.0108$\pm$0.0007%; Kowalski et al. 2009) and a lower flaring fraction for active WD+dMs (0.05$\pm$0.03%) compared to active dMs (0.28$\pm$0.05%; Kowalski et al. 2009). We discuss how our results constrain both the single and binary dM flare rates. Our results also constrain dM multiplicity, our knowledge of the Galactic transient background, and may be important for the habitability of attending planets around dMs with close companions.

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.

We have measured electron impact ionization (EII) for Fe 7+ from the ionization threshold up to 1200 eV. The measurements were performed using the TSR heavy ion storage ring. The ions were stored long enough prior to measurement to remove most metastables, resulting in a beam of 94% ground state ions. Comparing with the previously recommended atomic data, we find that the Arnaud & Raymond (1992) cross section is up to about 40\% larger than our measurement, with the largest discrepancies below about 400~eV. The cross section of Dere (2007) agrees to within 10%, which is about the magnitude of the experimental uncertainties. The remaining discrepancies between measurement and the most recent theory are likely due to shortcomings in the theoretical treatment of the excitation-autoionization contribution.

We provide software with a graphical user interface to calculate the phenomenology of a wide class of dark energy models featuring multiple scalar fields. The user chooses a subclass of models and, if desired, initial conditions, or else a range of initial parameters for Monte Carlo. The code calculates the energy density of components in the universe, the equation of state of dark energy, and the linear growth of density perturbations, all as a function of redshift and scale factor. The output also includes an approximate conversion into the average equation of state, as well as the common $(w_0, w_a)$ parametrization. The code is available here: http://github.com/kahinton/Dark-Energy-UI-and-MC

So far the highly unstable phase of luminous blue variables (LBVs) has not been understood well. It is still uncertain why and which massive stars enter this phase. Investigating the variabilities by looking for a possible regular or even (semi-)periodic behaviour could give a hint at the underlying mechanism for these variations and might answer the question of where these variabilities originate. Finding out more about the LBV phase also means understanding massive stars better in general, which have (e.g. by enriching the ISM with heavy elements, providing ionising radiation and kinetic energy) a strong and significant influence on the ISM, hence also on their host galaxy. Photometric and spectroscopic data were taken for the LBV Var C in M33 to investigate its recent status. In addition, scanned historic plates, archival data, and data from the literature were gathered to trace Var C's behaviour in the past. Its long-term variability and periodicity was investigated. Our investigation of the variability indicates possible (semi-)periodic behaviour with a period of 42.3 years for Var C. That Var C's light curve covers a time span of more than 100 years means that more than two full periods of the cycle are visible. The critical historic maximum around 1905 is less strong but discernible even with the currently rare historic data. The semi-periodic and secular structure of the light curve is similar to the one of LMC R71. Both light curves hint at a new aspect in the evolution of LBVs.

Two of the most commonly used tools to constrain the primordial non-Gaussianity are the bispectrum and the Minkowski functionals of CMB temperature anisotropies. These two measures of non-Gaussianity in principle provide distinct (though correlated) information, but in the past constraints from them have only been loosely compared, and not statistically combined. In this work we evaluate, for the first time, the covariance matrix between the local non-Gaussianity coefficient fnl estimated through the bispectrum and Minkowski functionals. We find that the estimators are positively correlated, with corerlation coefficient r ~ 0.3. Using the WMAP7 data to combine the two measures and accounting for the point-source systematics, we find the combined constraint fnl=37+/-28, which has a ~20% smaller error than either of the individual constraints.

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN.

We present new equations of state (EOS) for hydrogen and helium covering a wide range of temperatures from 60 K to 10$^7$ K and densities from $10^{-10}$ g/cm$^3$ to $10^3$ g/cm$^3$. They include an extended set of ab initio EOS data for the strongly correlated quantum regime with an accurate connection to data derived from other approaches for the neighboring regions. We compare linear-mixing isotherms based on our EOS tables with available real-mixture data. A first important astrophysical application of this new EOS data is the calculation of interior models for Jupiter and the comparison with recent results. Secondly, mass-radius relations are calculated for Brown Dwarfs which we compare with predictions derived from the widely used EOS of Saumon, Chabrier and van Horn. Furthermore, we calculate interior models for typical Brown Dwarfs with different masses, namely Corot-3b, Gliese-229b and Corot-15b, and the Giant Planet KOI-889b. The predictions for the central pressures and densities differ by up to 10$\%$ dependent on the EOS used. Our EOS tables are made available in the supplemental material of this paper.

Spacecraft data reveal a very Earth-like Jovian magnetic field. This is surprising since numerical simulations have shown that the vastly different interiors of terrestrial and gas planets can strongly affect the internal dynamo process. Here we present the first numerical dynamo that manages to match the structure and strength of the observed magnetic field by embracing the newest models for Jupiter's interior. Simulated dynamo action primarily occurs in the deep high electrical conductivity region while zonal flows are dynamically constrained to a strong equatorial jet in the outer envelope of low conductivity. Our model reproduces the structure and strength of the observed global magnetic field and predicts that secondary dynamo action associated to the equatorial jet produces banded magnetic features likely observable by the Juno mission. Secular variation in our model scales to about 2000 nT per year and should also be observable during the one year nominal mission duration.

We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation.

We present the fiducial main sequence stellar locus traced by 10 photometric colors observed by SDSS, 2MASS, and WISE. Median colors are determined using 1,052,793 stars with r-band extinction less than 0.125. We use this locus to measure the dust extinction curve relative to the r-band, which is consistent with previous measurements in the SDSS and 2MASS bands. The WISE band extinction coefficients are larger than predicted by standard extinction models. Using 13 lines of sight, we find variations in the extinction curve in H, Ks, and WISE bandpasses. Relative extinction decreases towards Galactic anti-center, in agreement with prior studies. Relative extinction increases with Galactic latitude, in contrast to previous observations. This indicates a universal mid-IR extinction law does not exist due to variations in dust grain size and chemistry with Galactocentric position. A preliminary search for outliers due to warm circumstellar dust is also presented, using stars with high signal-to-noise in the W3-band. We find 199 such outliers, identified by excess emission in Ks-W3. Inspection of SDSS images for these outliers reveals a large number of contaminants due to nearby galaxies. Six sources appear to be genuine dust candidates, yielding a fraction of systems with infrared excess of 0.12$\pm$0.05%.

This paper describes the data release of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey conducted between 2005 and 2007. Light curves, spectra, classifications, and ancillary data are presented for 10,258 variable and transient sources discovered through repeat ugriz imaging of SDSS Stripe 82, a 300 deg2 area along the celestial equator. This data release is comprised of all transient sources brighter than r~22.5 mag with no history of variability prior to 2004. Dedicated spectroscopic observations were performed on a subset of 889 transients, as well as spectra for thousands of transient host galaxies using the SDSS-III BOSS spectrographs. Photometric classifications are provided for the candidates with good multi-color light curves that were not observed spectroscopically. From these observations, 4607 transients are either spectroscopically confirmed, or likely to be, supernovae, making this the largest sample of supernova candidates ever compiled. We present a new method for SN host-galaxy identification and derive host-galaxy properties including stellar masses, star-formation rates, and the average stellar population ages from our SDSS multi-band photometry. We derive SALT2 distance moduli for a total of 1443 SN Ia with spectroscopic redshifts as well as photometric redshifts for a further 677 purely-photometric SN Ia candidates. Using the spectroscopically confirmed subset of the three-year SDSS-II SN Ia sample and assuming a flat Lambda-CDM cosmology, we determine Omega_M = 0.315 +/- 0.093 (statistical error only) and detect a non-zero cosmological constant at 5.7 sigmas.

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($\tau$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources.

We present and discuss an extensive data set for the non-Blazhko ab-type RR Lyrae star SDSSJ015450+001501, including optical SDSS ugriz light curves and spectroscopic data, LINEAR and CSS unfiltered optical light curves, and infrared 2MASS JHKs and WISE W1 and W2 light curves. Most notably, light curves obtained by 2MASS include close to 9000 photometric measures collected over 3.3 years and provide exceedingly precise view of near-IR variability. These data demonstrate that static atmosphere models are insufficient to explain multi-band photometric light curve behavior and present strong constraints for non-linear pulsation models for RR Lyrae stars. It is a challenge to modelers to produce theoretical light curves that can explain data presented here, which we make publicly available.

We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available.

We have measured dissociative recombination of HCl+ with electrons using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. We present the measured absolute merged beams recombination rate coefficient for collision energies from 0 to 4.5 eV. We have also developed a new method for deriving the cross section from the measurements. Our approach does not suffer from approximations made by previously used methods. The cross section was transformed to a plasma rate coefficient for the electron temperature range from T=10 to 5000 K. We show that the previously used HCl+ DR data underestimate the plasma rate coefficient by a factor of 1.5 at T=10 K and overestimate it by a factor of 3.0 at T=300 K. We also find that the new data may partly explain existing discrepancies between observed abundances of chlorine-bearing molecules and their astrochemical models.

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (\rho_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and \rho_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable.

The Green Bank Telescope (GBT) is the largest fully steerable radio telescope in the world and is one of our greatest tools for discovering and studying radio pulsars. Over the last decade, the GBT has successfully found over 100 new pulsars through large-area surveys. Here I discuss the two most recent---the GBT 350 MHz Drift-scan survey and the Green Bank North Celestial Cap survey. The primary science goal of both surveys is to find interesting individual pulsars, including young pulsars, rotating radio transients, exotic binary systems, and especially bright millisecond pulsars (MSPs) suitable for inclusion in Pulsar Timing Arrays, which are trying to directly detect gravitational waves. These two surveys have combined to discover 85 pulsars to date, among which are 14 MSPs and many unique and fascinating systems. I present highlights from these surveys and discuss future plans. I also discuss recent results from targeted GBT pulsar searches of globular clusters and Fermi sources.

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise.

We present follow-up observations and analysis of the recently discovered short period low-mass eclipsing binary, SDSS J001641-000925. With an orbital period of 0.19856 days, this system has one of the shortest known periods for an M dwarf binary system. Medium-resolution spectroscopy and multi-band photometry for the system are presented. Markov chain Monte Carlo modeling of the light curves and radial velocities yields estimated masses for the stars of M1 = 0.54 +/- 0.07 Msun and M2 = 0.34 +/- 0.04 Msun, and radii of R1 = 0.68 +/- 0.03 Rsun and R2 = 0.58 +/- 0.03 Rsun respectively. This solution places both components above the critical Roche overfill limit, providing strong evidence that SDSS J001641-000925 is the first verified M-dwarf contact binary system. Within the follow-up spectroscopy we find signatures of non-solid body rotation velocities, which we interpret as evidence for mass transfer or loss within the system. In addition, our photometry samples the system over 9 years, and we find strong evidence for period decay at the rate of dP/dt ~8 s/yr. Both of these signatures raise the intriguing possibility that the system is in over-contact, and actively losing angular momentum, likely through mass loss. This places SDSS J001641-000925 as not just the first M-dwarf over-contact binary, but one of the few systems of any spectral type known to be actively undergoing coalescence. Further study SDSS J001641-000925 is on-going to verify the nature of the system, which may prove to be a unique astrophysical laboratory.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2\sigma\ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (\rho$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

We use data from the WMAP temperature maps to constrain a scale-dependent generalization of the popular 'local' model for primordial non-Gaussianity. In the model where the parameter fNL is allowed to run with scale k, fNL(k) = fNL* (k/k_piv)^n, we constrain the running to be n = 0.30(+1.9)(-1.2) at 95% confidence, marginalized over the amplitude fNL*. The constraints depend somewhat on the prior probabilities assigned to the two parameters. In the near future, constraints from a combination of Planck and large-scale structure surveys are expected to improve this limit by about an order of magnitude and usefully constrain classes of inflationary models.

We forecast combined future constraints from the cosmic microwave background and large-scale structure on the models of primordial non-Gaussianity. We study the generalized local model of non-Gaussianity, where the parameter f_NL is promoted to a function of scale, and present the principal component analysis applicable to an arbitrary form of f_NL(k). We emphasize the complementarity between the CMB and LSS by using Planck, DES and BigBOSS surveys as examples, forecast constraints on the power-law f_NL(k) model, and introduce the figure of merit for measurements of scale-dependent non-Gaussianity.

We present R-Band light curves of Type II supernovae (SNe) from the Caltech Core Collapse Project (CCCP). With the exception of interacting (Type IIn) SNe and rare events with long rise times, we find that most light curve shapes belong to one of three distinct classes: plateau, slowly declining and rapidly declining events. The last class is composed solely of Type IIb SNe which present similar light curve shapes to those of SNe Ib, suggesting, perhaps, similar progenitor channels. We do not find any intermediate light curves, implying that these subclasses are unlikely to reflect variance of continuous parameters, but rather might result from physically distinct progenitor systems, strengthening the suggestion of a binary origin for at least some stripped SNe. We find a large plateau luminosity range for SNe IIP, while the plateau lengths seem rather uniform at approximately 100 days. As analysis of additional CCCP data goes on and larger samples are collected, demographic studies of core collapse SNe will likely continue to provide new constraints on progenitor scenarios.

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

We propose a robust principal component analysis (PCA) framework for the exploitation of multi-band photometric measurements in large surveys. Period search results are improved using the time series of the first principal component due to its optimized signal-to-noise ratio.The presence of correlated excess variations in the multivariate time series enables the detection of weaker variability. Furthermore, the direction of the largest variance differs for certain types of variable stars. This can be used as an efficient attribute for classification. The application of the method to a subsample of Sloan Digital Sky Survey Stripe 82 data yielded 132 high-amplitude Delta Scuti variables. We found also 129 new RR Lyrae variables, complementary to the catalogue of Sesar et al., 2010, extending the halo area mapped by Stripe 82 RR Lyrae stars towards the Galactic bulge. The sample comprises also 25 multiperiodic or Blazhko RR Lyrae stars.

The amount and distribution of heavy elements in Jupiter gives indications on the process of its formation and evolution. Core mass and metallicity predictions however depend on the equations of state used, and on model assumptions. We present an improved ab initio hydrogen equation of state, H-REOS.2 and compute the internal structure and thermal evolution of Jupiter within the standard three-layer approach. The advance over our previous Jupiter models with H-REOS.1 by Nettelmann et al.(2008) is that the new models are also consistent with the observed 2 or more times solar heavy element abundances in Jupiter's atmosphere. Such models have a rock core mass Mcore=0-8 ME, total mass of heavy elements MZ=28-32 ME, a deep internal layer boundary at 4 or more Mbar, and a cooling time of 4.4-5.0 Gyrs when assuming homogeneous evolution. We also calculate two-layer models in the manner of Militzer et al.(2008) and find a comparable large core of 16-21 ME, out of which ~11 ME is helium, but a significantly higher envelope metallicity of 4.5 times solar. According to our preferred three-layer models, neither the characteristic frequency (nu0 ~156 microHz) nor the normalized moment of inertia (~0.276) are sensitive to the core mass but accurate measurements could well help to rule out some classes of models.

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

We review current methods for building PSF-matching kernels for the purposes of image subtraction or coaddition. Such methods use a linear decomposition of the kernel on a series of basis functions. The correct choice of these basis functions is fundamental to the efficiency and effectiveness of the matching - the chosen bases should represent the underlying signal using a reasonably small number of shapes, and/or have a minimum number of user-adjustable tuning parameters. We examine methods whose bases comprise multiple Gauss-Hermite polynomials, as well as a form free basis composed of delta-functions. Kernels derived from delta-functions are unsurprisingly shown to be more expressive; they are able to take more general shapes and perform better in situations where sum-of-Gaussian methods are known to fail. However, due to its many degrees of freedom (the maximum number allowed by the kernel size) this basis tends to overfit the problem, and yields noisy kernels having large variance. We introduce a new technique to regularize these delta-function kernel solutions, which bridges the gap between the generality of delta-function kernels, and the compactness of sum-of-Gaussian kernels. Through this regularization we are able to create general kernel solutions that represent the intrinsic shape of the PSF-matching kernel with only one degree of freedom, the strength of the regularization lambda. The role of lambda is effectively to exchange variance in the resulting difference image with variance in the kernel itself. We examine considerations in choosing the value of lambda, including statistical risk estimators and the ability of the solution to predict solutions for adjacent areas. Both of these suggest moderate strengths of lambda between 0.1 and 1.0, although this optimization is likely dataset dependent.

We present the first rates of flares from M dwarf stars in both red optical and near infrared (NIR) filters. We have studied ~50,000 M dwarfs from the SDSS Stripe 82 area, and 1,321 M dwarfs from the 2MASS Calibration Scan Point Source Working Database that overlap SDSS imaging fields. We assign photometric spectral types from M0 to M6 using (r-i) and (i-z) colors for every star in our sample. Stripe 82 stars each have 50-100 epochs of data, while 2MASS Calibration stars have ~1900 epochs. From these data we estimate the observed rates and theoretical detection thresholds for flares in eight photometric bands as a function of spectral type. Optical flare rates are found to be in agreement with previous studies, while the frequency per hour of NIR flare detections is found to be more than two orders of magnitude lower. An excess of small amplitude flux increases in all bands exhibits a power-law distribution, which we interpret as the result of flares below our detection thresholds. In order to investigate the recovery efficiency for flares in each filter, we extend a two-component flare model into the NIR. Quiescent M0-M6 spectral templates were used with the model to predict the photometric response of flares from u to Ks. We determine that red optical filters are sensitive to flares with u-band amplitudes >2 mag, and NIR filters to flares with delta u>4.5 mag. Our model predicts that M0 stars have the best color-contrast for J-band detections, but M4-M6 stars should show the highest rate of NIR flares with amplitudes of delta J ~0.01 mag. Characterizing flare rates and photometric variations at longer wavelengths is important for predicting the signatures of M dwarf variability in next-generation surveys, and we discuss their impact on surveys such as LSST.

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements.

We describe photometric recalibration of data obtained by the asteroid survey LINEAR. Although LINEAR was designed for astrometric discovery of moving objects, the dataset described here contains over 5 billion photometric measurements for about 25 million objects, mostly stars. We use SDSS data from the overlapping ~10,000 deg^2 of sky to recalibrate LINEAR photometry, and achieve errors of 0.03 mag for sources not limited by photon statistics, with errors of 0.2 mag at r~18. With its 200 observations per object on average, LINEAR data provide time domain information for the brightest 4 magnitudes of SDSS survey. At the same time, LINEAR extends the deepest similar wide-area variability survey, the Northern Sky Variability Survey, by 3 mag. We briefly discuss the properties of about 7,000 visually confirmed periodic variables, dominated by roughly equal fractions of RR Lyrae stars and eclipsing binary stars, and analyze their distribution in optical and infra-red color-color diagrams. The LINEAR dataset is publicly available from the SkyDOT website (http://skydot.lanl.gov).

The luminosity changes of most types of variable stars are correlated in the different wavelengths, and these correlations may be exploited for several purposes: for variability detection, for distinction of microvariability from noise, for period search or for classification. Principal component analysis is a simple and well-developed statistical tool to analyze correlated data. We will discuss its use on variable objects of Stripe 82 of the Sloan Digital Sky Survey, with the aim of identifying new RR Lyrae and SX Phoenicis-type candidates. The application is not straightforward because of different noise levels in the different bands, the presence of outliers that can be confused with real extreme observations, under- or overestimated errors and the dependence of errors on the magnitudes. These particularities require robust methods to be applied together with the principal component analysis. The results show that PCA is a valuable aid in variability analysis with multi-band data.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

We present a catalog of periodic stellar variability in the "Stripe 82" region of the Sloan Digital Sky Survey (SDSS). After aggregating and recalibrating catalog-level data from the survey, we ran a period-finding algorithm (Supersmoother) on all point-source lightcurves. We used color selection to identify systems that are likely to contain low-mass stars, in particular M dwarfs and white dwarfs. In total, we found 207 candidates, the vast majority of which appear to be in eclipsing binary systems. The catalog described in this paper includes 42 candidate M dwarf / white dwarf pairs, 4 white-dwarf pairs, 59 systems whose colors indicate they are composed of 2 M dwarfs and whose lightcurve shapes suggest they are in detached eclipsing binaries, and 28 M dwarf systems whose lightcurve shapes suggest they are in contact binaries. We find no detached systems with periods longer than 3 days, thus the majority of our sources are likely to have experienced orbital spin-up and enhanced magnetic activity. Indeed, twenty-six of twenty-seven M dwarf systems that we have spectra for show signs of chromospheric magnetic activity, far higher than the 24% seen in field stars of the same spectral type. We also find binaries composed of stars that bracket the expected boundary between partially and fully convective interiors, which will allow the measurement of the stellar mass-radius relationship across this transition. The majority of our contact systems have short orbital periods, with small variance (0.02 days) in the sample near the observed cutoff of 0.22 days. The accumulation of these stars at short orbital period suggests that the process of angular momentum loss, leading to period evolution, becomes less efficient at short periods. (Abridged)

We generalize the local model of primordial non-Gaussianity by promoting the parameter fNL to a general scale-dependent function fNL(k). We calculate the resulting bispectrum and the effect on the bias of dark matter halos, and thus the extent to which fNL(k) can be measured from the large-scale structure observations. By calculating the principal components of fNL(k), we identify scales where this form of non-Gaussianity is best constrained and estimate the overlap with previously studied local and equilateral non-Gaussian models.

Time-resolved databases with large spatial coverage are quickly becoming a standard tool for all types of astronomical studies. We report preliminary results from our search for stellar flares in the 2MASS calibration fields. A sample of 4343 M dwarfs, spatially matched between the SDSS and the 2MASS calibration fields, each with hundreds to thousands of epochs in near infrared bandpasses, is analyzed using a modified Welch-Stetson index to characterize the variability. A Monte Carlo model was used to assess the noise of the variability index. We find significnat residuals above the noise with power-law slopes of -3.37 and -4.05 for our JH and HKs distributions respectively. This is evidence for flares being observed from M dwarfs in infrared photometry.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

We present three transits of GJ 1214b, observed as part of the Apache Point Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). We used APOSTLE r-band lightcurves in conjunction with previously gathered data of GJ 1214b to re-derive system parameters. By using parameters such as transit duration and ingress/egress length we are able to reduce the degeneracies between parameters in the fitted transit model, which is a preferred condition for Markov Chain Monte Carlo techniques typically used to quantify uncertainties in measured parameters. The joint analysis of this multi-wavelength dataset confirms earlier estimates of system parameters including planetary orbital period, the planet-to-star radius ratio and stellar density. We fit the photometric spectralenergy distribution of GJ 1214 to derive stellar luminosity, which we then use to derive its absolute mass and radius. From these derived stellar properties and previously published radial velocity data we were able to refine estimates of the absolute parameters for the planet GJ 1214b. Transit times derived from our study show no evidence for strong transit timing variations. Some lightcurves we present show features that we believe are due to stellar activity. During the first night we observed a rise in the out-of-eclipse flux of GJ 1214 with a characteristic fast-rise exponential decay shape commonly associated with stellar flares. On the second night we observed a minor brightening during transit, which we believe might have been caused by the planet obscuring a star-spot on the stellar disk.

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

Large, multi-frequency imaging surveys, such as the Large Synaptic Survey Telescope (LSST), need to do near-real time analysis of very large datasets. This raises a host of statistical and computational problems where standard methods do not work. In this paper, we study a proposed method for combining stacks of images into a single summary image, sometimes referred to as a template. This task is commonly referred to as image coaddition. In part, we focus on a method proposed in previous work, which outlines a procedure for combining stacks of images in an online fashion in the Fourier domain. We evaluate this method by comparing it to two straightforward methods through the use of various criteria and simulations. Note that the goal is not to propose these comparison methods for use in their own right, but to ensure that additional complexity also provides substantially improved performance.

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Koz{\l}owski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate based on data from the Sloan Digital Sky Survey II (SDSS-II) Supernova Survey. The adopted sample of supernovae (SNe) includes 516 SNe Ia at redshift z \lesssim 0.3, of which 270 (52%) are spectroscopically identified as SNe Ia. The remaining 246 SNe Ia were identified through their light curves; 113 of these objects have spectroscopic redshifts from spectra of their host galaxy, and 133 have photometric redshifts estimated from the SN light curves. Based on consideration of 87 spectroscopically confirmed non-Ia SNe discovered by the SDSS-II SN Survey, we estimate that 2.04+1.61-0.95 % of the photometric SNe Ia may be misidentified. The sample of SNe Ia used in this measurement represents an order of magnitude increase in the statistics for SN Ia rate measurements in the redshift range covered by the SDSS-II Supernova Survey. If we assume a SN Ia rate that is constant at low redshift (z < 0.15), then the SN observations can be used to infer a value of the SN rate of rV = (2.69+0.34+0.21-0.30-0.01) x10^{-5} SNe yr^{-1} Mpc-3 (H0 /(70 km s^{-1} Mpc^{-1}))^{3} at a mean redshift of ~ 0.12, based on 79 SNe Ia of which 72 are spectroscopically confirmed. However, the large sample of SNe Ia included in this study allows us to place constraints on the redshift dependence of the SN Ia rate based on the SDSS-II Supernova Survey data alone. Fitting a power-law model of the SN rate evolution, r_V(z) = A_p x ((1 + z)/(1 + z0))^{\nu}, over the redshift range 0.0 < z < 0.3 with z0 = 0.21, results in A_p = (3.43+0.15-0.15) x 10^{-5} SNe yr^{-1} Mpc-3 (H0 /(70 km s^{-1} Mpc^{-1}))^{3} and \nu = 2.04+0.90-0.89.

ABRIDGED We present measurements of the Type Ia supernova (SN) rate in galaxy clusters based on data from the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. The cluster SN Ia rate is determined from 9 SN events in a set of 71 C4 clusters at z <0.17 and 27 SN events in 492 maxBCG clusters at 0.1 < z < 0.3$. We find values for the cluster SN Ia rate of $({0.37}^{+0.17+0.01}_{-0.12-0.01}) \mathrm{SNu}r h^{2}$ and $({0.55}^{+0.13+0.02}_{-0.11-0.01}) \mathrm{SNu}r h^{2}$ ($\mathrm{SNu}x = 10^{-12} L_{x\sun}^{-1} \mathrm{yr}^{-1}$) in C4 and maxBCG clusters, respectively, where the quoted errors are statistical and systematic, respectively. The SN rate for early-type galaxies is found to be $({0.31}^{+0.18+0.01}_{-0.12-0.01}) \mathrm{SNu}r h^{2}$ and $({0.49}^{+0.15+0.02}_{-0.11-0.01})$ $\mathrm{SNu}r h^{2}$ in C4 and maxBCG clusters, respectively. The SN rate for the brightest cluster galaxies (BCG) is found to be $({2.04}^{+1.99+0.07}_{-1.11-0.04}) \mathrm{SNu}r h^{2}$ and $({0.36}^{+0.84+0.01}_{-0.30-0.01}) \mathrm{SNu}r h^{2}$ in C4 and maxBCG clusters. The ratio of the SN Ia rate in cluster early-type galaxies to that of the SN Ia rate in field early-type galaxies is ${1.94}^{+1.31+0.043}_{-0.91-0.015}$ and ${3.02}^{+1.31+0.062}_{-1.03-0.048}$, for C4 and maxBCG clusters. The SN rate in galaxy clusters as a function of redshift...shows only weak dependence on redshift. Combining our current measurements with previous measurements, we fit the cluster SN Ia rate data to a linear function of redshift, and find $r_{L} = $ $[(0.49^{+0.15}_{-0.14}) +$ $(0.91^{+0.85}_{-0.81}) \times z]$ $\mathrm{SNu}B$ $h^{2}$. A comparison of the radial distribution of SNe in cluster to field early-type galaxies shows possible evidence for an enhancement of the SN rate in the cores of cluster early-type galaxies... we estimate the fraction of cluster SNe that are hostless to be $(9.4^+8._3-5.1)%$.

We present 2323 High-Amplitude \delta-Scuti (HADS) candidates discovered in the Large Magellanic Cloud (LMC) by the SuperMACHO survey (Rest et al. 2005). Frequency analyses of these candidates reveal that several are multimode pulsators, including 119 whose largest amplitude of pulsation is in the fundamental (F) mode and 19 whose largest amplitude of pulsation is in the first overtone (FO) mode. Using Fourier decomposition of the HADS light curves, we find that the period-luminosity (PL) relation defined by the FO pulsators does not show a clear separation from the PL-relation defined by the F pulsators. This differs from other instability strip pulsators such as type c RR Lyrae. We also present evidence for a larger amplitude, subluminous population of HADS similar to that observed in Fornax (Poretti et al. 2008).

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

We present an improved analysis of halo substructure traced by RR Lyrae stars in the SDSS stripe 82 region. With the addition of SDSS-II data, a revised selection method based on new ugriz light curve templates results in a sample of 483 RR Lyrae stars that is essentially free of contamination. The main result from our first study persists: the spatial distribution of halo stars at galactocentric distances 5--100 kpc is highly inhomogeneous. At least 20% of halo stars within 30 kpc from the Galactic center can be statistically associated with substructure. We present strong direct evidence, based on both RR Lyrae stars and main sequence stars, that the halo stellar number density profile significantly steepens beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars are associated with substructure. By using a novel method that simultaneously combines data for RR Lyrae and main sequence stars, and using photometric metallicity estimates for main sequence stars derived from deep co-added u-band data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing arm) towards R.A.2h-3h and Dec~0 deg to be 0.3 dex higher ([Fe/H]=-1.2) than that of surrounding halo field stars. Together with a similar result for another major halo substructure, the Monoceros stream, these results support theoretical predictions that an early forming, smooth inner halo, is metal poor compared to high surface brightness material that have been accreted onto a later-forming outer halo. The mean metallicity of stars in the outer halo that are not associated with detectable clumps may still be more metal-poor than the bulk of inner-halo stars, as has been argued from other data sets.

We present an analysis of the luminosity distances of Type Ia Supernovae from the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey in conjunction with other intermediate redshift (z<0.4) cosmological measurements including redshift-space distortions from the Two-degree Field Galaxy Redshift Survey (2dFGRS), the Integrated Sachs-Wolfe (ISW) effect seen by the SDSS, and the latest Baryon Acoustic Oscillation (BAO) distance scale from both the SDSS and 2dFGRS. We have analysed the SDSS-II SN data alone using a variety of "model-independent" methods and find evidence for an accelerating universe at >97% level from this single dataset. We find good agreement between the supernova and BAO distance measurements, both consistent with a Lambda-dominated CDM cosmology, as demonstrated through an analysis of the distance duality relationship between the luminosity (d_L) and angular diameter (d_A) distance measures. We then use these data to estimate w within this restricted redshift range (z<0.4). Our most stringent result comes from the combination of all our intermediate-redshift data (SDSS-II SNe, BAO, ISW and redshift-space distortions), giving w = -0.81 +0.16 -0.18(stat) +/- 0.15(sys) and Omega_M=0.22 +0.09 -0.08 assuming a flat universe. This value of w, and associated errors, only change slightly if curvature is allowed to vary, consistent with constraints from the Cosmic Microwave Background. We also consider more limited combinations of the geometrical (SN, BAO) and dynamical (ISW, redshift-space distortions) probes.

We use the new SNe Ia discovered by the SDSS-II Supernova Survey together with additional supernova datasets as well as observations of the cosmic microwave background and baryon acoustic oscillations to constrain cosmological models. This complements the analysis presented by Kessler et al. in that we discuss and rank a number of the most popular non-standard cosmology scenarios. When this combined data-set is analyzed using the MLCS2k2 light-curve fitter, we find that more exotic models for cosmic acceleration provide a better fit to the data than the Lambda-CDM model. For example, the flat DGP model is ranked higher by our information criteria tests than the standard model. When the dataset is instead analyzed using the SALT-II light-curve fitter, the standard cosmological constant model fares best. Our investigation also includes inhomogeneous Lemaitre-Tolman-Bondi (LTB) models. While our LTB models can be made to fit the supernova data as well as any other model, the extra parameters they require are not supported by our information criteria analysis.

We present measurements of the Hubble diagram for 103 Type Ia supernovae (SNe) with redshifts 0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift SN Ia surveys. We combine the SDSS-II measurements with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy Survey, the Hubble Space Telescope, and a compilation of nearby SN Ia measurements. Combining the SN Hubble diagram with measurements of Baryon Acoustic Oscillations from the SDSS Luminous Red Galaxy sample and with CMB temperature anisotropy measurements from WMAP, we estimate the cosmological parameters w and Omega_M, assuming a spatially flat cosmological model (FwCDM) with constant dark energy equation of state parameter, w. For the FwCDM model and the combined sample of 288 SNe Ia, we find w = -0.76 +- 0.07(stat) +- 0.11(syst), Omega_M = 0.306 +- 0.019(stat) +- 0.023(syst) using MLCS2k2 and w = -0.96 +- 0.06(stat) +- 0.12(syst), Omega_M = 0.265 +- 0.016(stat) +- 0.025(syst) using the SALT-II fitter. We trace the discrepancy between these results to a difference in the rest-frame UV model combined with a different luminosity correction from color variations; these differences mostly affect the distance estimates for the SNLS and HST supernovae. We present detailed discussions of systematic errors for both light-curve methods and find that they both show data-model discrepancies in rest-frame $U$-band. For the SALT-II approach, we also see strong evidence for redshift-dependence of the color-luminosity parameter (beta). Restricting the analysis to the 136 SNe Ia in the Nearby+SDSS-II samples, we find much better agreement between the two analysis methods but with larger uncertainties.

We present ugriz light curves for 146 spectroscopically confirmed or spectroscopically probable Type Ia supernovae from the 2005 season of the SDSS-II Supernova survey. The light curves have been constructed using a photometric technique that we call scene modelling, which is described in detail here; the major feature is that supernova brightnesses are extracted from a stack of images without spatial resampling or convolution of the image data. This procedure produces accurate photometry along with accurate estimates of the statistical uncertainty, and can be used to derive photometry taken with multiple telescopes. We discuss various tests of this technique that demonstrate its capabilities. We also describe the methodology used for the calibration of the photometry, and present calibrated magnitudes and fluxes for all of the spectroscopic SNe Ia from the 2005 season.

We present a flare rate analysis of 50,130 M dwarf light curves in SDSS Stripe 82. We identified 271 flares using a customized variability index to search ~2.5 million photometric observations for flux increases in the u- and g-bands. Every image of a flaring observation was examined by eye and with a PSF-matching and image subtraction tool to guard against false positives. Flaring is found to be strongly correlated with the appearance of H-alpha in emission in the quiet spectrum. Of the 99 flare stars that have spectra, we classify 8 as relatively inactive. The flaring fraction is found to increase strongly in stars with redder colors during quiescence, which can be attributed to the increasing flare visibility and increasing active fraction for redder stars. The flaring fraction is strongly correlated with |Z| distance such that most stars that flare are within 300 pc of the Galactic plane. We derive flare u-band luminosities and find that the most luminous flares occur on the earlier-type M dwarfs. Our best estimate of the lower limit on the flaring rate (averaged over Stripe 82) for flares with \Delta u \ge 0.7 magnitudes on stars with u < 22 is 1.3 flares hour^-1 square degree^-1 but can vary significantly with the line-of-sight.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0<z<4. Aside from cosmology, SNe Ia regulate galactic and cluster chemical evolution, inform stellar evolution, and are laboratories for extreme physics. Essential probes of SNe Ia in these contexts include spectroscopy from the UV to the IR, X-ray cluster and SN remnant observations, spectropolarimetry, and advanced theoretical studies. While there are an abundance of discovery facilities planned, there is a deficit of follow-up resources. Living in the systematics era demands deep understanding rather than larger statistics. NOAO ReSTAR initiative to build 2-4m telescopes would provide necessary follow-up capability. Finally, to fully exploit LSST, well-matched wide-field spectroscopic capabilities are desirable.

The next decade of survey astronomy has the potential to transform our knowledge of variable stars. Stellar variability underpins our knowledge of the cosmological distance ladder, and provides direct tests of stellar formation and evolution theory. Variable stars can also be used to probe the fundamental physics of gravity and degenerate material in ways that are otherwise impossible in the laboratory. The computational and engineering advances of the past decade have made large-scale, time-domain surveys an immediate reality. Some surveys proposed for the next decade promise to gather more data than in the prior cumulative history of astronomy. The actual implementation of these surveys will have broad implications for the types of science that will be enabled. We examine the design considerations for an optimal time-domain photometric survey dedicated to variable star science, including: observing cadence, wavelength coverage, photometric and astrometric accuracy, single-epoch and cumulative depth, overall sky coverage, and data access by the broader astronomical community. The best surveys must combine aspects from each of these considerations to fully realize the potential for the next decade of time-domain science.

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.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

We report the discovery of a minor planet (2006 SQ372) on an orbit with a perihelion of 24 AU and a semimajor axis of 796 AU. Dynamical simulations show that this is a transient orbit and is unstable on a timescale of 200 Myrs. Falling near the upper semimajor axis range of the scattered disk and the lower semimajor axis range of the Oort Cloud, previous membership in either class is possible. By modeling the production of similar orbits from the Oort Cloud as well as from the scattered disk, we find that the Oort Cloud produces 16 times as many objects on SQ372-like orbits as the scattered disk. Given this result, we believe this to be the most distant long-period comet ever discovered. Furthermore, our simulation results also indicate that 2000 OO67 has had a similar dynamical history. Unaffected by the "Jupiter-Saturn Barrier," these two objects are most likely long-period comets from the inner Oort Cloud.

Gliese 436 is an M dwarf with a mass of 0.45 Msun and hosts the extrasolar planet GL 436b [3, 6, 7, 2], which is currently the least massive transiting planet with a mass of ~23.17 Mearth [10], and the only planet known to transit an M dwarf. GL 436b represents the first transiting detection of the class of extrasolar planets known as "Hot Neptunes" that have masses within a few times that of Neptune's mass (~17 Mearth) and orbital semimajor axis <0.1 AU about the host star. Unlike most other known transiting extrasolar planets, GL 436b has a high eccentricity (e~0.16). This brings to light a new parameter space for habitability zones of extrasolar planets with host star masses much smaller than typical stars of roughly a solar mass. This unique system is an ideal candidate for orbital perturbation and transit-time variation (TTV) studies to detect smaller, possibly Earth-mass planets in the system. In April 2008 we began a long-term intensive campaign to obtain complete high-precision light curves using the Apache Point Observatory's 3.5-meter telescope, NMSU's 1-meter telescope (located at APO), and Sommers Bausch Observatory's 24" telescope. These light curves are being analyzed together, along with amateur and other professional astronomer observations. Results of our analysis are discussed. Continued measurements over the next few years are needed to determine if additional planets reside in the system, and to study the impact of other manifestations on the light curves, such as star spots and active regions.

Significant progress in the description of quasar variability has been recently made by employing SDSS and POSS data. Common to most studies is a fundamental assumption that photometric observations at two epochs for a large number of quasars will reveal the same statistical properties as well-sampled light curves for individual objects. We critically test this assumption using light curves for a sample of $\sim$2,600 spectroscopically confirmed quasars observed about 50 times on average over 8 years by the SDSS stripe 82 survey. We find that the dependence of the mean structure function computed for individual quasars on luminosity, rest-frame wavelength and time is qualitatively and quantitatively similar to the behavior of the structure function derived from two-epoch observations of a much larger sample. We also reproduce the result that the variability properties of radio and X-ray selected subsamples are different. However, the scatter of the variability structure function for fixed values of luminosity, rest-frame wavelength and time is similar to the scatter induced by the variance of these quantities in the analyzed sample. Hence, our results suggest that, although the statistical properties of quasar variability inferred using two-epoch data capture some underlying physics, there is significant additional information that can be extracted from well-sampled light curves for individual objects.

The Large Synoptic Survey Telescope (LSST) will be a large, wide-field ground-based system designed to obtain, starting in 2015, multiple images of the sky that is visible from Cerro Pachon in Northern Chile. About 90% of the observing time will be devoted to a deep-wide-fast survey mode which will observe a 20,000 deg$^2$ region about 1000 times during the anticipated 10 years of operations (distributed over six bands, $ugrizy$). Each 30-second long visit will deliver 5$\sigma$ depth for point sources of $r\sim24.5$ on average. The co-added map will be about 3 magnitudes deeper, and will include 10 billion galaxies and a similar number of stars. We discuss various measurements that will be automatically performed for these 20 billion sources, and how they can be used for classification and determination of source physical and other properties. We provide a few classification examples based on SDSS data, such as color classification of stars, color-spatial proximity search for wide-angle binary stars, orbital-color classification of asteroid families, and the recognition of main Galaxy components based on the distribution of stars in the position-metallicity-kinematics space. Guided by these examples, we anticipate that two grand classification challenges for LSST will be 1) rapid and robust classification of sources detected in difference images, and 2) {\it simultaneous} treatment of diverse astrometric and photometric time series measurements for an unprecedentedly large number of objects.

We present ground-based observations of the transiting Neptune-mass planet Gl 436b obtained with the 3.5-meter telescope at Apache Point Observatory and other supporting telescopes. Included in this is an observed transit in early 2005, over two years before the earliest reported transit detection. We have compiled all available transit data to date and perform a uniform modeling using the JKTEBOP code. We do not detect any transit timing variations of amplitude greater than ~1 minute over the ~3.3 year baseline. We do however find possible evidence for a self-consistent trend of increasing orbital inclination, transit width, and transit depth, which supports the supposition that Gl 436b is being perturbed by another planet of < 12 Earth Masses in a non-resonant orbit.

We present the discovery of 130.87 day periodic near-infrared flux variability for the Class II T Tauri star WL 4 (= 2MASS J16271848-2429059, ISO-Oph 128). Our data are from the 2MASS Calibration Point Source Working Database, and constitute 1580 observations in J, H and Ks of a field in Rho Ophiuchus used to calibrate the 2MASS All-Sky Survey. We identify a light curve for WL 4 with eclipse amplitudes of ~0.4 mag lasting more than one-quarter the period, and color variations in J-H and H-Ks of ~0.1 mag. The long period cannot be explained by stellar rotation. We propose that WL 4 is a triple YSO system, with an inner binary orbital period of 130.87 days. We posulate that we are observing each component of the inner binary alternately being eclipsed by a circum-binary disk with respect to our line of sight. This system will be useful in investigating terrestrial zone YSO disk properties and dynamics at ~1 Myr.

We present a measurement of the rate of type Ia supernovae (SNe Ia) from the first of three seasons of data from the SDSS-II Supernova Survey. For this measurement, we include 17 SNe Ia at redshift $z\le0.12$. Assuming a flat cosmology with $\Omega_m = 0.3=1-\Omega_\Lambda$, we find a volumetric SN Ia rate of $[2.93^{+0.17}_{-0.04}({\rm systematic})^{+0.90}_{-0.71}({\rm statistical})] \times 10^{-5} {\rm SNe} {\rm Mpc}^{-3} h_{70}^3 {\rm year}^{-1}$, at a volume-weighted mean redshift of 0.09. This result is consistent with previous measurements of the SN Ia rate in a similar redshift range. The systematic errors are well controlled, resulting in the most precise measurement of the SN Ia rate in this redshift range. We use a maximum likelihood method to fit SN rate models to the SDSS-II Supernova Survey data in combination with other rate measurements, thereby constraining models for the redshift-evolution of the SN Ia rate. Fitting the combined data to a simple power-law evolution of the volumetric SN Ia rate, $r_V \propto (1+z)^{\beta}$, we obtain a value of $\beta = 1.5 \pm 0.6$, i.e. the SN Ia rate is determined to be an increasing function of redshift at the $\sim 2.5 \sigma$ level. Fitting the results to a model in which the volumetric SN rate, $r_V=A\rho(t)+B\dot \rho(t)$, where $\rho(t)$ is the stellar mass density and $\dot \rho(t)$ is the star formation rate, we find $A = (2.8 \pm 1.2) \times 10^{-14} \mathrm{SNe} \mathrm{M}_{\sun}^{-1} \mathrm{year}^{-1}$, $B = (9.3^{+3.4}_{-3.1})\times 10^{-4} \mathrm{SNe} \mathrm{M}_{\sun}^{-1}$.

We report the discovery of an extensive system of scattered light echo arclets associated with the recent supernovae in the local neighbourhood of the Milky Way: Tycho (SN 1572) and Cassiopeia A. Existing work suggests that the Tycho SN was a thermonuclear explosion while the Cas A supernova was a core collapse explosion. Precise classifications according to modern nomenclature require spectra of the outburst light. In the case of ancient SNe, this can only be done with spectroscopy of their light echo, where the discovery of the light echoes from the outburst light is the first step. Adjacent light echo positions suggest that Cas A and Tycho may share common scattering dust structures. If so, it is possible to measure precise distances between historical Galactic supernovae. On-going surveys that alert on the development of bright scattered-light echo features have the potential to reveal detailed spectroscopic information for many recent Galactic supernovae, both directly visible and obscured by dust in the Galactic plane.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

We present multiepoch spectra of 13 high-redshift Type Ia supernovae (SNe Ia) drawn from the literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin & Tonry to determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time yields an apparent aging rate consistent with the 1/(1+z) factor (where z is the redshift) expected in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's "tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The best-fit exponent for these models is consistent with the expected 1/(1+z) factor.

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

AM CVn systems are a select group of ultracompact binaries with the shortest orbital periods of any known binary subclass; mass-transfer is likely from a low-mass (partially-)degenerate secondary onto a white dwarf primary, driven by gravitational radiation. In the past few years, the Sloan Digital Sky Survey (SDSS) has provided five new AM CVns. Here we report on two further candidates selected from more recent SDSS data. SDSS J1208+3550 is similar to the earlier SDSS discoveries, recognized as an AM CVn via its distinctive spectrum which is dominated by helium emission. From the expanded SDSS Data Release 6 (DR6) spectroscopic area, we provide an updated surface density estimate for such AM CVns of order 10^{-3.1} to 10^{-2.5} per deg^2 for 15<g<20.5. In addition, we present another new candidate AM CVn, SDSS J2047+0008, that was discovered in the course of followup of SDSS-II supernova candidates. It shows nova-like outbursts in multi-epoch imaging data; in contrast to the other SDSS AM CVn discoveries, its (outburst) spectrum is dominated by helium absorption lines, reminiscent of KL Dra and 2003aw. The variability selection of SDSS J2047+0008 from the 300 deg^2 of SDSS Stripe 82 presages further AM CVn discoveries in future deep, multicolor, and time-domain surveys such as LSST. The new additions bring the total SDSS yield to seven AM CVns thus far, a substantial contribution to this rare subclass, versus the dozen previously known.

We report the successful identification of the type of the supernova responsible for the supernova remnant SNR 0509-675 in the Large Magellanic Cloud (LMC) using Gemini spectra of surrounding light echoes. The ability to classify outbursts associated with centuries-old remnants provides a new window into several aspects of supernova research and is likely to be successful in providing new constraints on additional LMC supernovae as well as their historical counterparts in the Milky Way Galaxy (MWG). The combined spectrum of echo light from SNR 0509-675 shows broad emission and absorption lines consistent with a supernova (SN) spectrum. We create a spectral library consisting of 26 SNe Ia and 6 SN Ib/c that are time-integrated, dust-scattered by LMC dust, and reddened by the LMC and MWG. We fit these SN templates to the observed light echo spectrum using $\chi^2$ minimization as well as correlation techniques, and we find that overluminous 91T-like SNe Ia with $\dm15<0.9$ match the observed spectrum best.

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Database. Additional photometry from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we derive an orbital period of 2.6390157 +/- 0.0000016 days. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses (M_1 = 0.66 +/- 0.03 M_sun; M_2 = 0.62 +/- 0.03 M_sun) and radii (R_1 = 0.64 +/- 0.08 R_sun; R_2 = 0.61 +/- 0.09 R_sun) of the components, which are consistent with empirical mass-radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of H-alpha emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.

I provide an incomplete inventory of the astronomical variability that will be found by next-generation time-domain astronomical surveys. These phenomena span the distance range from near-Earth satellites to the farthest Gamma Ray Bursts. The surveys that detect these transients will issue alerts to the greater astronomical community; this decision process must be extremely robust to avoid a slew of ``false'' alerts, and to maintain the community's trust in the surveys. I review the functionality required of both the surveys and the telescope networks that will be following them up, and the role of VOEvents in this process. Finally, I offer some ideas about object and event classification, which will be explored more thoroughly by other articles in these proceedings.

We have developed an end-to-end photometric data processing pipeline to compare current photometric algorithms commonly used on ground-based imaging data. This testbed is exceedingly adaptable, and enables us to perform many research and development tasks, including image subtraction and co-addition, object detection and measurements, the production of photometric catalogs, and the creation and stocking of database tables with time-series information. This testing has been undertaken to evaluate existing photometry algorithms for consideration by a next-generation image processing pipeline for the Large Synoptic Survey Telescope (LSST). We outline the results of our tests for four packages: The Sloan Digital Sky Survey's (SDSS) Photo package, Daophot and Allframe, DoPhot, and two versions of Source Extractor (SExtractor). The ability of these algorithms to perform point-source photometry, astrometry, shape measurements, star-galaxy separation, and to measure objects at low signal-to-noise is quantified. We also perform a detailed crowded field comparison of Daophot and Allframe, and profile the speed and memory requirements in detail for SExtractor. We find that both Daophot and Photo are able to perform aperture photometry to high enough precision to meet LSST's science requirements, and less adequately at PSF-fitting photometry. Photo performs the best at simultaneous point and extended-source shape and brightness measurements. SExtractor is the fastest algorithm, and recent upgrades in the software yield high-quality centroid and shape measurements with little bias towards faint magnitudes. Allframe yields the best photometric results in crowded fields.

The Sloan Digital Sky Survey-II Supernova Survey has identified a large number of new transient sources in a 300 sq. deg. region along the celestial equator during its first two seasons of a three-season campaign. Multi-band (ugriz) light curves were measured for most of the sources, which include solar system objects, Galactic variable stars, active galactic nuclei, supernovae (SNe), and other astronomical transients. The imaging survey is augmented by an extensive spectroscopic follow-up program to identify SNe, measure their redshifts, and study the physical conditions of the explosions and their environment through spectroscopic diagnostics. During the survey, light curves are rapidly evaluated to provide an initial photometric type of the SNe, and a selected sample of sources are targeted for spectroscopic observations. In the first two seasons, 476 sources were selected for spectroscopic observations, of which 403 were identified as SNe. For the Type Ia SNe, the main driver for the Survey, our photometric typing and targeting efficiency is 90%. Only 6% of the photometric SN Ia candidates were spectroscopically classified as non-SN Ia instead, and the remaining 4% resulted in low signal-to-noise, unclassified spectra. This paper describes the search algorithm and the software, and the real-time processing of the SDSS imaging data. We also present the details of the supernova candidate selection procedures and strategies for follow-up spectroscopic and imaging observations of the discovered sources.

We present the first large-scale effort of creating composite spectra of high-redshift type Ia supernovae (SNe Ia) and comparing them to low-redshift counterparts. Through the ESSENCE project, we have obtained 107 spectra of 88 high-redshift SNe Ia with excellent light-curve information. In addition, we have obtained 397 spectra of low-redshift SNe through a multiple-decade effort at Lick and Keck Observatories, and we have used 45 UV spectra obtained by HST/IUE. The low-redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the Lick composite spectra allows a nearly perfect match of the overall spectral-energy distribution with the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with host-galaxy light than their low-redshift counterparts. This is caused by observing objects at all redshifts with the same slit width, which corresponds to different projected distances. After correcting for the galaxy-light contamination, subtle differences in the spectra remain. We have estimated the systematic errors when using current spectral templates for K-corrections to be ~0.02 mag. The variance in the composite spectra give an estimate of the intrinsic variance in low-redshift maximum-light SN spectra of ~3% in the optical and growing toward the UV. The difference between the maximum light low and high-redshift spectra constrain SN evolution between our samples to be < 10% in the rest-frame optical.

This paper has been withdrawn by the author.

This paper has been withdrawn

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

The first cosmological results from the ESSENCE supernova survey (Wood-Vasey et al. 2007) are extended to a wider range of cosmological models including dynamical dark energy and non-standard cosmological models. We fold in a greater number of external data sets such as the recent Higher-z release of high-redshift supernovae (Riess et al. 2007) as well as several complementary cosmological probes. Model comparison statistics such as the Bayesian and Akaike information criteria are applied to gauge the worth of models. These statistics favor models that give a good fit with fewer parameters. Based on this analysis, the preferred cosmological model is the flat cosmological constant model, where the expansion history of the universe can be adequately described with only one free parameter describing the energy content of the universe. Among the more exotic models that provide good fits to the data, we note a preference for models whose best-fit parameters reduce them to the cosmological constant model.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

We present extensive u'g'r'i'BVRIYJHKs photometry and optical spectroscopy of SN 2005hk. These data reveal that SN 2005hk was nearly identical in its observed properties to SN 2002cx, which has been called ``the most peculiar known type Ia supernova.'' Both supernovae exhibited high ionization SN 1991T-like pre-maximum spectra, yet low peak luminosities like SN 1991bg. The spectra reveal that SN 2005hk, like SN 2002cx, exhibited expansion velocities that were roughly half those of typical type Ia supernovae. The R and I light curves of both supernovae were also peculiar in not displaying the secondary maximum observed for normal type Ia supernovae. Our YJH photometry of SN 2005hk reveals the same peculiarity in the near-infrared. By combining our optical and near-infrared photometry of SN 2005hk with published ultraviolet light curves obtained with the Swift satellite, we are able to construct a bolometric light curve from ~10 days before to ~60 days after B maximum. The shape and unusually low peak luminosity of this light curve, plus the low expansion velocities and absence of a secondary maximum at red and near-infrared wavelengths, are all in reasonable agreement with model calculations of a 3D deflagration which produces ~0.25 M_sun of 56Ni.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the equation of state parameter of the dark energy. We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to the dark energy equation of state parameter w=P/rho c^2 for a given fixed amount of telescope time. For our survey on the CTIO 4m telescope, measuring the luminosity distances and redshifts for supernovae at modest redshifts (z~0.5 +- 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find supernovae automatically and in near real-time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4m natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 type Ia SNe, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for the all type Ia supernovae found by ESSENCE and used in our measurement of w, presented in Wood-Vasey et al, 2007.

We present constraints on the dark energy equation-of-state parameter, w=P/(rho c^2), using 60 Type Ia supernovae (SNe Ia) from the ESSENCE supernova survey. We derive a set of constraints on the nature of the dark energy assuming a flat Universe. By including constraints on (Omega_M, w) from baryon acoustic oscillations, we obtain a value for a static equation-of-state parameter w=-1.05^{+0.13}_{-0.12} (stat; 1 sigma) +- 0.11 (sys) and Omega_M=0.274^{+0.033}_{-0.020} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.96. These results are consistent with those reported by the SuperNova Legacy Survey in a similar program measuring supernova distances and redshifts. We evaluate sources of systematic error that afflict supernova observations and present Monte Carlo simulations that explore these effects. Currently, the largest systematic currently with the potential to affect our measurements is the treatment of extinction due to dust in the supernova host galaxies. Combining our set of ESSENCE SNe Ia with the SuperNova Legacy Survey SNe Ia, we obtain a joint constraint of w=-1.07^{+0.09}_{-0.09} (stat; 1 sigma) +- 0.12 (sys), Omega_M=0.267^{+0.028}_{-0.018} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.91. The current SN Ia data are fully consistent with a cosmological constant.

We present a set of 11 type Ia supernova (SN Ia) lightcurves with dense, pre-maximum sampling. These supernovae (SNe), in galaxies behind the Large Magellanic Cloud (LMC), were discovered by the SuperMACHO survey. The SNe span a redshift range of z = 0.11 - 0.35. Our lightcurves contain some of the earliest pre-maximum observations of SNe Ia to date. We also give a functional model that describes the SN Ia lightcurve shape (in our VR-band). Our function uses the "expanding fireball" model of Goldhaber et al. (1998) to describe the rising lightcurve immediately after explosion but constrains it to smoothly join the remainder of the lightcurve. We fit this model to a composite observed VR-band lightcurve of three SNe between redshifts of 0.135 to 0.165. These SNe have not been K-corrected or adjusted to account for reddening. In this redshift range, the observed VR-band most closely matches the rest frame V-band. Using the best fit to our functional description of the lightcurve, we find the time between explosion and observed VR-band maximum to be 17.6+-1.3(stat)+-0.07(sys) rest-frame days for a SN Ia with a VR-band Delta m_{-10} of 0.52mag. For the redshifts sampled, the observed VR-band time-of-maximum brightness should be the same as the rest-frame V-band maximum to within 1.1 rest-frame days.

We present ground-based and Hubble Space Telescope optical observations of the X-ray flash (XRF) 020903, covering 300 days. The afterglow showed a very rapid rise in the first day, followed by a relatively slow decay in the next few days. There was a clear bump in the light curve after ~25 days, accompanied by a drastic change in the spectral energy distribution. The light curve and the spectral energy distribution are naturally interpreted as the emergence -- and subsequent decay -- of a supernova (SN), similar to SN 1998bw. At peak luminosity, the SN is estimated to be 0.8 +/- 0.1 mag fainter than SN1998bw. This argues in favor of the existence of a supernova associated with this X-ray flash. A spectrum obtained 35 days after the burst shows emission lines from the host galaxy. We use this spectrum to put an upper limit on the oxygen abundance of the host at [O/H] < -0.6 dex. We also discuss a possible trend between the softness of several bursts and the early behavior of the optical afterglow, in the sense that XRFs and X-ray rich GRBs seem to have a plateau phase or even a rising light curve. This can be naturally explained in models where XRFs are similar to GRBs but seen off the jet axis.

In principle, the light from historical supernovae could still be visible as scattered-light echoes even centuries later. However, while echoes have been discovered around some nearby extragalactic supernovae well after the explosion, targeted searches have not recovered any echoes in the regions of historical Galactic supernovae. The discovery of echoes can allow us to pinpoint the supernova event both in position and age and, most importantly, allow us to acquire spectra of the echo light to type the supernova centuries after the direct light from the explosion first reached the Earth. Here we report on the discovery of three faint new variable surface brightness complexes with high apparent proper motion pointing back to well-defined positions in the Large Magellanic Cloud (LMC). These positions correspond to three of the six smallest (and likely youngest) previously catalogued supernova remnants, and are believed to be due to thermonuclear (Type Ia) supernovae. Using the distance and proper motions of these echo arcs, we estimate ages of 610 and 410 yr for the echoes #2 and #3.

Using archival data of low-redshift (z < 0.01) Type Ia supernovae (SN Ia) and recent observations of high-redshift (0.16 < z <0.64; Matheson et al. 2005) SN Ia, we study the "uniformity'' of the spectroscopic properties of nearby and distant SN Ia. We find no difference in the measures we describe here. In this paper, we base our analysis solely on line-profile morphology, focusing on measurements of the velocity location of maximum absorption (vabs) and peak emission (vpeak). We find that the evolution of vabs and vpeak for our sample lines (Ca II 3945, Si II 6355, and S II 5454, 5640) is similar for both the low- and high-redshift samples. We find that vabs for the weak S II 5454, 5640 lines, and vpeak for S II 5454, can be used to identify fast-declining [dm15 > 1.7] SN Ia, which are also subluminous. In addition, we give the first direct evidence in two high-z SN Ia spectra of a double-absorption feature in Ca II 3945, an event also observed, though infrequently, in low-redshift SN Ia spectra (6/22 SN Ia in our local sample). We report for the first time the unambiguous and systematic intrinsic blueshift of peak emission of optical P-Cygni line profiles in Type Ia spectra, by as much as 8000 km/s. All the high-z SN Ia analyzed in this paper were discovered and followed up by the ESSENCE collaboration, and are now publicly available.

The proper usage of Type Ia supernovae (SNe Ia) as distance indicators has revolutionized cosmology, and added a new dominant component to the energy density of the Universe, dark energy. Following the discovery and confirmation era, the currently ongoing SNe Ia surveys aim to determine the properties of the dark energy. ESSENCE is a five year ground-based supernova survey aimed at finding and characterizing 200 SNe Ia in the redshift domain z=[0.2-0.8]. The goal of the project is to put constraints on the equation of state parameter, w, of the dark energy with an accuracy of <10%. This paper presents these ongoing efforts in the context of the current developments in observational cosmology.

Characterizing the nature and spatial distribution of the lensing objects that produce the previously measured microlensing optical depth toward the Large Magellanic Cloud (LMC) remains an open problem. We present an appraisal of the ability of the SuperMACHO Project, a next-generation microlensing survey directed toward the LMC, to discriminate between various proposed lensing populations. We consider two scenarios: lensing by a uniform foreground screen of objects and self-lensing by LMC stars. We have carried out extensive simulations, based upon data obtained during the first year of the project, to assess the SuperMACHO survey's ability to discriminate between these two scenarios. We find that the event rate itself shows significant sensitivity to the choice of the LMC luminosity function, limiting the conclusions which can be drawn from the absolute rate. If instead we determine the differential event rate across the LMC, we will decrease the impact of these systematic biases and render our conclusions more robust. With this approach the SuperMACHO Project should be able to distinguish between the two categories of lens populations. This will provide important constraints on the nature of the lensing objects and their contributions to the Galactic dark matter halo.

We present broad-band light curves of nine supernovae ranging in redshift from 0.5 to 0.8. The supernovae were discovered as part of the ESSENCE project, and the light curves are a combination of Cerro Tololo 4-m and Hubble Space Telescope (HST) photometry. On the basis of spectra and/or light-curve fitting, eight of these objects are definitely Type Ia supernovae, while the classification of one is problematic. The ESSENCE project is a five-year endeavor to discover about 200 high-redshift Type Ia supernovae, with the goal of tightly constraining the time average of the equation-of-state parameter [w = p/(rho c^2)] of the "dark energy." To help minimize our systematic errors, all of our ground-based photometry is obtained with the same telescope and instrument. In 2003 the highest-redshift subset of ESSENCE supernovae was selected for detailed study with HST. Here we present the first photometric results of the survey. We find that all but one of the ESSENCE SNe have slowly declining light curves, and the sample is not representative of the low-redshift set of ESSENCE Type Ia supernovae. This is unlikely to be a sign of evolution in the population. We attribute the decline-rate distribution of HST events to a selection bias at the high-redshift edge of our sample and find that such a bias will infect other magnitude-limited SN Ia searches unless appropriate precautions are taken.

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.

AM CVn systems are a rare (about a dozen previously known) class of cataclysmic variables, arguably encompassing the shortest orbital periods (down to about 10 minutes) of any known binaries. Both binary components are thought to be degenerate (or partially so), likely with mass-transfer from a helium-rich donor onto a white dwarf, driven by gravitational radiation. Although rare, AM CVn systems are of high interest as possible SN Ia progenitors, and because they are predicted to be common sources of gravity waves in upcoming experiments such as LISA. We have identified four new AM CVn candidates from the Sloan Digital Sky Survey (SDSS) spectral database. All four show hallmark spectroscopic characteristics of the AM CVn class: each is devoid of hydrogen features, and instead shows a spectrum dominated by helium. All four show double-peaked emission, indicative of helium-dominated accretion disks. Limited time-series CCD photometric follow-on data have been obtained for three of the new candidates from the ARC 3.5m; most notably, a 28.3 minute binary period with sharp, deep eclipses is discovered in one case, SDSS J0926+3624. This is the first confirmed eclipsing AM CVn, and our data allow initial estimates of binary parameters for this ultracompact system. The four new SDSS objects also provide a substantial expansion of the currently critically-small sample of AM CVn systems.

We present previously unpublished photometry of three Large Magellanic Cloud (LMC) microlensing events and show that the new photometry confirms the microlensing interpretation of these events. These events were discovered by the MACHO Project alert system and were also recovered by the analysis of the 5.7 year MACHO data set. This new photometry provides a substantial increase in the signal-to-noise ratio over the previously published photometry and in all three cases, the gravitational microlensing interpretation of these events is strengthened. The new data consist of MACHO-Global Microlensing Alert Network (GMAN) follow-up images from the CTIO 0.9 telescope plus difference imaging photometry of the original MACHO data from the 1.3m "Great Melbourne" telescope at Mt. Stromlo. We also combine microlensing light curve fitting with photometry from high resolution HST images of the source stars to provide further confirmation of these events and to show that the microlensing interpretation of event MACHO-LMC-23 is questionable. Finally, we compare our results with the analysis of Belokurov, Evans & Le Du who have attempted to classify candidate microlensing events with a neural network method, and we find that their results are contradicted by the new data and more powerful light curve fitting analysis for each of the four events considered in this paper. The failure of the Belokurov, Evans & Le Du method is likely to be due to their use of a set of insensitive statistics to feed their neural networks.

In preparation for the Supernova Survey of the Sloan Digital Sky Survey (SDSS) II, a proposed 3-year extension to the SDSS, we have conducted an early engineering and science run during the fall of 2004, which consisted of approximately 20 scheduled nights of repeated imaging of half of the southern equatorial stripe. Transient supernova-like events were detected in near real-time and photometric measurements were made in the five SDSS filter bandpasses with a cadence of ~2 days. Candidate type Ia supernovae (SNe) were pre-selected based on their colors, light curve shape, and the properties of the host galaxy. Follow-up spectroscopic observations were performed with the Astrophysical Research Consortium 3.5m telescope and the 9.2m Hobby-Eberly Telescope to confirm their types and measure the redshifts. The 2004 campaign resulted in 22 spectroscopically confirmed SNe, which includes 16 type Ia, 5 type II, and 1 type Ib/c. These SN Ia will help fill in the sparsely sampled redshift interval of z = 0.05 - 0.35, the so-called 'redshift desert', in the Hubble diagram. Detailed investigation of the spectral properties of these moderate-redshift SNe Ia will also provide a bridge between local SNe and high-redshift objects, and will help us understand the systematics for future cosmological applications that require high photometric precision. Finally, the large survey volume also provides the opportunity to select unusual supernovae for spectroscopic study that are poorly sampled in other surveys. We report on some of the early results from this program and discuss potential future applications.

We present the results of spectroscopic observations of targets discovered during the first two years of the ESSENCE project. The goal of ESSENCE is to use a sample of ~200 Type Ia supernovae (SNe Ia) at moderate redshifts (0.2 < z < 0.8) to place constraints on the equation of state of the Universe. Spectroscopy not only provides the redshifts of the objects, but also confirms that some of the discoveries are indeed SNe Ia. This confirmation is critical to the project, as techniques developed to determine luminosity distances to SNe Ia depend upon the knowledge that the objects at high redshift are the same as the ones at low redshift. We describe the methods of target selection and prioritization, the telescopes and detectors, and the software used to identify objects. The redshifts deduced from spectral matching of high-redshift SNe Ia with low-redshift SNe Ia are consistent with those determined from host-galaxy spectra. We show that the high-redshift SNe Ia match well with low-redshift templates. We include all spectra obtained by the ESSENCE project, including 52 SNe Ia, 5 core-collapse SNe, 12 active galactic nuclei, 19 galaxies, 4 possibly variable stars, and 16 objects with uncertain identifications.

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 first results from our next-generation microlensing survey, the SuperMACHO project. We are using the CTIO 4m Blanco telescope and the MOSAIC imager to carry out a search for microlensing toward the Large Magellanic Cloud (LMC). We plan to ascertain the nature of the population responsible for the excess microlensing rate seen by the MACHO project. Our observing strategy is optimized to measure the differential microlensing rate across the face of the LMC. We find this derivative to be relatively insensitive to the details of the LMC's internal structure but a strong discriminant between Galactic halo and LMC self lensing. In December 2003 we completed our third year of survey operations. 2003 also marked the first year of real-time microlensing alerts and photometric and spectroscopic followup. We have extracted several dozen microlensing candidates, and we present some preliminary light curves and related information. Similar to the MACHO project, we find SNe behind the LMC to be a significant contaminant - this background has not been completely removed from our current single-color candidate sample. Our follow-up strategy is optimized to discriminate between SNe and true microlensing.

We report on the methodology and first results from the Deep Lens Survey transient search. We utilize image subtraction on survey data to yield all sources of optical variability down to 24th magnitude. Images are analyzed immediately after acquisition, at the telescope and in near-real time, to allow for followup in the case of time-critical events. All classes of transients are posted to the web upon detection. Our observing strategy allows sensitivity to variability over several decades in timescale. The DLS is the first survey to classify and report all types of photometric and astrometric variability detected, including solar system objects, variable stars, supernovae, and short timescale phenomena. Three unusual optical transient events were detected, flaring on thousand-second timescales. All three events were seen in the B passband, suggesting blue color indices for the phenomena. One event (OT 20020115) is determined to be from a flaring Galactic dwarf star of spectral type dM4. From the remaining two events, we find an overall rate of \eta = 1.4 events deg-2 day-1 on thousand-second timescales, with a 95% confidence limit of \eta < 4.3. One of these events (OT 20010326) originated from a compact precursor in the field of galaxy cluster Abell 1836, and its nature is uncertain. For the second (OT 20030305) we find strong evidence for an extended extragalactic host. A dearth of such events in the R passband yields an upper 95% confidence limit on short timescale astronomical variability between 19.5 < R < 23.4 of \eta_R < 5.2. We report also on our ensemble of astrometrically variable objects, as well as an example of photometric variability with an undetected precursor.

We report the weak lensing discovery, spectroscopic confirmation, and weak lensing tomography of a massive cluster of galaxies at $z=0.68$, demonstrating that shear selection of clusters works at redshifts high enough to be cosmologically interesting. The mass estimate from weak lensing, $11.1 +- 2.8 x 10^{14} (r/Mpc)$ solar masses within projected radius r, agrees with that derived from the spectroscopy ($\sigma_v = 980 km s^{-1}$), and with the position of an arc which is likely to be a strongly lensed background galaxy. The redshift estimate from weak lensing tomography is consistent with the spectroscopy, demonstrating the feasibility of baryon-unbiased mass surveys. This tomographic technique will be able to roughly identify the redshifts of any dark clusters which may appear in shear-selected samples, up to z ~ 1.

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 present a (V-R)-based reddening map of about 43 square degrees of the Galactic bulge/bar. The map is constructed using template image photometry from the MACHO microlensing survey, contains 9717 resolution elements, and is based on (V-R)-color averages of the entire color-magnitude diagrams (CMDs) in 4 by 4 arc-minute tiles. The conversion from the observed color to the reddening follows from an assumption that CMDs of all bulge fields would look similar in the absence of extinction. Consequently, the difference in observed color between various fields originates from varying contribution of the disk extinction summed along different lines of sight. We check that our (V-R) colors correlate very well with infrared and optical reddening maps. We show that a dusty disk obeying a cosec|b| extinction law, E(V-R) = 0.0274 cosec|b|, provides a good approximation to the extinction toward the MACHO bulge/bar fields. The large-scale (V-R)-color and visual extinction map presented here is publicly available in the electronic edition of the Journal and on the World Wide Web.

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.

The Deep Lens Survey (DLS) is a deep BVRz' imaging survey of seven 2x2 degree fields, with all data to be made public. The primary scientific driver is weak gravitational lensing, but the survey is also designed to enable a wide array of other astrophysical investigations. A unique feature of this survey is the search for transient phenomena. We subtract multiple exposures of a field, detect differences, classify, and release transients on the Web within about an hour of observation. Here we summarize the scientific goals of the DLS, field and filter selection, observing techniques and current status, data reduction, data products and release, and transient detections. Finally, we discuss some lessons which might apply to future large surveys such as LSST.

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 re-analysis of the longest timescale gravitational microlensing event discovered to date: MACHO-99-BLG-22/OGLE-1999-BUL-32, which was discovered by both the MACHO and OGLE microlensing alert systems. Our analysis of this microlensing parallax event includes a likelihood analysis of the lens position based upon a standard model of the Galactic velocity distribution, and this implies that the lens could be a black hole of ~100 M_solar at a distance of a few hundred parsecs in the Galactic disk or a massive stellar remnant (black hole or neutron star) in the Galactic bulge. Our new analysis includes data from the MACHO, GMAN, and MPS collaborations in addition to the OGLE data used in a previous analysis by Mao et al (2002). The crucial feature that distinguishes our analysis from that of Mao et al is an accurate constraint on the direction of lens motion and an analysis of the implications of this direction.

We present UBVRIz lightcurves of the Type Ia SN 1999ee and the Type Ib/c SN 1999ex, both located in the galaxy IC 5179. SN 1999ee has an extremely well sampled lightcurve spanning from 10 days before Bmax through 53 days after peak. Near maximum we find systematic differences ~0.05 mag in photometry measured with two different telescopes, even though the photometry is reduced to the same local standards around the supernova using the specific color terms for each instrumental system. We use models for our bandpasses and spectrophotometry of SN 1999ee to derive magnitude corrections (S-corrections) and remedy this problem. This exercise demonstrates the need of accurately characterizing the instrumental system before great photometric accuracies of Type Ia supernovae can be claimed. It also shows that this effect can have important astrophysical consequences since a small systematic shift of 0.02 mag in the B-V color can introduce a 0.08 mag error in the extinction corrected peak B magnitudes of a supernova and thus lead to biased cosmological parameters. The data for the Type Ib/c SN 1999ex present us with the first ever observed shock breakout of a supernova of this class. These observations show that shock breakout occurred 18 days before Bmax and support the idea that Type Ib/c supernovae are due to core collapse of massive stars rather than thermonuclear disruption of white dwarfs.

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 observations of microlensing event MACHO-98-BLG-35 which reached a peak magnification factor of almost 80. These observations by the Microlensing Planet Search (MPS) and the MOA Collaborations place strong constraints on the possible planetary system of the lens star and show intriguing evidence for a low mass planet with a mass fraction $4\times 10^{-5} \leq \epsilon \leq 2\times 10^{-4}$. A giant planet with $\epsilon = 10^{-3}$ is excluded from 95% of the region between 0.4 and 2.5 $R_E$ from the lens star, where $R_E$ is the Einstein ring radius of the lens. This exclusion region is more extensive than the generic "lensing zone" which is $0.6 - 1.6 R_E$. For smaller mass planets, we can exclude 57% of the "lensing zone" for $\epsilon = 10^{-4}$ and 14% of the lensing zone for $\epsilon = 10^{-5}$. The mass fraction $\epsilon = 10^{-5}$ corresponds to an Earth mass planet for a lensing star of mass $\sim 0.3 \msun$. A number of similar events will provide statistically significant constraints on the prevalence of Earth mass planets. In order to put our limits in more familiar terms, we have compared our results to those expected for a Solar System clone averaging over possible lens system distances and orientations. We find that such a system is ruled out at the 90% confidence level. A copy of the Solar System with Jupiter replaced by a second Saturn mass planet can be ruled out at 70% confidence. Our low mass planetary signal (few Earth masses to Neptune mass) is significant at the $4.5\sigma$ confidence level. If this planetary interpretation is correct, the MACHO-98-BLG-35 lens system constitutes the first detection of a low mass planet orbiting an ordinary star without gas giant planets.

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 study of extra-solar planetary systems has emerged as a new discipline of observational astronomy in the past few years with the discovery of a number of extra-solar planets. The properties of most of these extra-solar planets were not anticipated by theoretical work on the formation of planetary systems. Here we report observations and light curve modeling of gravitational microlensing event MACHO-97-BLG-41, which indicates that the lens system consists of a planet orbiting a binary star system. According to this model, the mass ratio of the binary star system is 3.8:1 and the stars are most likely to be a late K dwarf and an M dwarf with a separation of about 1.8 AU. A planet of about 3 Jupiter masses orbits this system at a distance of about 7 AU. If our interpretation of this light curve is correct, it represents the first discovery of a planet orbiting a binary star system and the first detection of a Jovian planet via the gravitational microlensing technique. It suggests that giant planets may be common in short period binary star systems.

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.

In this paper, we will review major new results regarding classical Cepheids, in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). Specifically, we discuss recent work regarding multimode Cepheids and describe new observations of a W Vir star (HV 5756) and a Cepheid which are each in eclipsing binary systems. An additional interesting pulsating supergiant in an eclipsing system is also identified. Ephemerides for eclipses for the three systems are provided.

We present the observations of the binary lensing event MACHO-98-SMC-1 conducted at the Mt.~Stromlo 74" telescope by the Microlensing Planet Search (MPS) collaboration. The MPS data constrain the first caustic crossing to have occurred after 1998 June 5.55 UT and thus directly rule out one of the two fits presented by the PLANET collaboration (model II). This substantially reduces the uncertainty in the the relative proper motion estimations of the lens object. We perform joint binary microlensing fits of the MPS data together with the publicly available data from the EROS, MACHO/GMAN and OGLE collaborations. We also study the binary lens fit parameters previously published by the PLANET and MACHO/GMAN collaborations by using them as initial values for $\chi^2$ minimization. Fits based on the PLANET model I appear to be in conflict with the GMAN-CTIO data. From our best fit, we find that the lens system has a proper motion of $\mu = 1.3\pm 0.2 \kmsk$ with respect to the source, which implies that the lens system is most likely to be located in the Small Magellanic Cloud strengthening the conclusion of previous reports.

We present an analysis of the period distribution of $\sim 1800$ Cepheids in the Large Magellanic Cloud, based on data obtained by the MACHO microlensing experiment and on a previous catalogue by Payne-Gaposchkin. Using stellar evolution and pulsation models, we construct theoretical period-frequency distributions that are compared to the observations. These models reveal that a significant burst of star formation has occurred recently in the LMC ($\sim 1.15\times 10^8$ years). We also show that during the last $\sim 10^8$ years, the main center of star formation has been propagating from SE to NW along the bar. We find that the evolutionary masses of Cepheids are still smaller than pulsation masses by $\sim 7$ % and that the red edge of the Cepheid instability strip could be slightly bluer than indicated by theory. There are $\sim 600$ Cepheids with periods below $\sim 2.5$ days cannot be explained by evolution theory. We suggest that they are anomalous Cepheids; a number of these stars are double-mode Cepheids.

MACHO Project photometry of 45 LMC FO/SO beat Cepheids which pulsate in the first and second overtone (FO and SOo, respectively) has been analysed to determine the lightcurve characteristics for the SO mode of Cepheid pulsation. We predict that singly-periodic SO Cepheids will have nearly sinusoidal lightcurves; that we will only be able to discern SO Cepheids from fundamental (F) and (FO) Cepheids for P <= 1.4 days; and that the SO distribution will overlap the short-period edge of the LMC FO Cepheid period-luminosity relation (when both are plotted as a function of photometric period). We also report the discovery of one SO Cepheid candidate, MACHO*05:03:39.6$-$70:04:32, with a photometric period of 0.775961 +/- 0.000019 days and an instrumental amplitude of 0.047 +/- 0.009 mag in V.

Using the novel data set of the MACHO project, we show that three ``Blazhko Effect'' RR Lyrae stars show nearly-pure amplitude modulation of a single pulsation mode. This is strong observational evidence that the ``Oblique Pulsator Model'' is the correct solution to this 90-year-old problem.

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.

We report the first discovery of a gravitational microlensing candidate towards a new population of source stars, the Small Magellanic Cloud (SMC). The candidate event's light curve shows no variation for 3 years before an upward excursion lasting $ \sim 217$ days that peaks around January 11, 1997 at a magnification of $ \sim 2.1$. Microlensing events towards the Large Magellanic Cloud and the Galactic bulge have allowed important conclusions to be reached on the stellar and dark matter content of the Milky Way. The SMC gives a new line-of-sight through the Milky Way, and is expected to prove useful in determining the flattening of the Galactic halo.

We review the pulsating stars contained in the top 24 fields of the MACHO bulge database, with special emphasis on the red semiregular stars. Based on period, amplitude and color cuts, we have selected a sample of 2000 semiregular variables with $15<P<100$ days. Their color-magnitude diagram is presented, and period-luminosity relation is studied, as well as their spatial distribution. We find that they follow the bar, unlike the RR Lyrae in these fields.

We report the discovery of 20 1H/2H and 7 F/1H beat Cepheids in the SMC by the MACHO Project. We utilize the 20 1H/2H stars to determine lightcurve shape for the SMC second-overtone (2H) mode of Cepheid pulsation. We predict, similar to the findings of Alcock et al. (1997, ApJ, submitted), that 2H Cepheids will have nearly or purely sinusoidal light variations; that the P-L relation for 2H Cepheids will not be distinguishable from the P-L relation for 1H Cepheids within photometric accuracy; and that 2H stars may be discernable from F and 1H stars using the amplitude-period diagram and Fourier parameter progressions for periods P < 0.7 days, our current sample 2H period limit.

We report the discovery of RV Tauri stars in the Large Magellanic Cloud. In light and colour curve behaviour, the RV Tauri stars appear to be a direct extension of the type II Cepheids to longer periods. A single period- luminosity-colour relationship is seen to describe both the type II Cepheids and the RV Tauri stars in the LMC. We derive the relation: $V_{\circ} = 17.89 (\pm 0.20) - 2.95 (\pm 0.12) \log_{10}P + 5.49 (\pm 0.35) \bar{(V-R)_{\circ}}$, valid for type II Cepheids and RV Tauri stars in the period range $0.9 < \log_{10}P < 1.75$. Assuming a distance modulus to the Large Magellanic Cloud of 18.5, the relation in terms of the absolute luminosities becomes: $M_{V} = -0.61 (\pm 0.20) - 2.95 (\pm 0.12) \log_{10}P + 5.49 (\pm 0.35) \bar{(V-R)_{\circ}}$.

The MACHO Project monitors millions of stars in the Large Magellanic Cloud, the Small Magellanic Cloud and the bulge of the Milky Way searching for the gravitational microlensing signature of baryonic dark matter. This Project has yielded surprising results. An analysis of two years of data monitoring the Large Magellanic Cloud points to {$\sim 50%$} of the mass of the Milky Way's halo in compact objects of {$\sim 0.5 M_{\odot}$}. An analysis of one year of monitoring the bulge has yielded more microlensing than predicted without the invocation of a massive bar or significant disk dark matter. The huge database of light curves created by this search is yielding information on extremely rare types of astrophysical variability as well as providing temporal detail for the study of well known variable astrophysical phenomena. The variable star catalog created from this database is previewed and example light curves are presented.

We describe the search for delta Scuti stars in the MACHO database of bulge fields. Concentrating on a sample of high amplitude delta Scutis, we examine the light curves and pulsation modes. We also discuss their spatial distribution and evolutionary status using mean colors and absolute magnitudes.

Mean colors and magnitudes of RR Lyrae stars in 24 fields towards the Galactic bulge from the MACHO database are presented. Accurate mean reddenings are computed for these fields on the basis of the mean colors. The distribution along the line of sight of the RR Lyrae population is examined on the basis of the mean magnitudes, and it is shown that the bulk of the RR Lyrae population is not barred. Only the RR Lyrae in the inner fields closer to the Galactic center (l<4, b>-4) show evidence for a bar. The red giant clump stars in the MACHO fields, however, clearly show a barred distribution, confirming the results of previous studies. Given the different spatial distribution, the RR Lyrae and the clump giants trace two different populations. The RR Lyrae would represent the inner extension of the Galactic halo in these fields.

The recent suggestion that the microlensing events observed towards the Large Magellanic Cloud are due to an intervening Sgr-like dwarf galaxy is examined. A search for foreground RR Lyrae in the MACHO photometry database yields 20 stars whose distance distribution follow the expected halo density profile. Cepheid and red giant branch clump stars in the MACHO database are consistent with membership in the LMC. There is also no evidence in the literature for a distinct kinematic population, for intervening gas, or for the turn-off of such a hypothetical galaxy. We conclude that if the lenses are in a foreground galaxy, it must be a particularly dark galaxy.

We measure the zero point of the Stanek (1996) extinction map by comparing the observed (V-K) colors of 20 RR Lyrae stars (type ab) found in the MACHO survey with their intrinsic (V-K)_0 colors as a function of period as determined from nearby RR Lyrae stars. We find that the zero point of the Stanek map should be changed by \Delta A_V = -0.11 +/- 0.05 mag, in excellent agreement with the recent measurement of Gould, Popowski & Terndrup (1997) using K giants.

We present analysis of MACHO Alert 95-30, a dramatic gravitational microlensing event towards the Galactic bulge whose peak magnification departs significantly from the standard point-source microlensing model. Alert 95-30 was observed in real-time by the Global Microlensing Alert Network (GMAN), which obtained densely sampled photometric and spectroscopic data throughout the event. We interpret the light-curve ``fine structure'' as indicating transit of the lens across the extended face of the source star. This signifies resolution of a star several kpc distant. We find a lens angular impact parameter theta_{min}/theta_{source} = 0.715 +/- 0.003. This information, along with the radius and distance of the source, provides an additional constraint on the lensing system. Spectroscopic and photometric data indicate the source is an M4 III star of radius 61 +/- 12 Rsun, located on the far side of the bulge at 9 kpc. We derive a lens angular velocity, relative to the source, of 21.5 +/- 4.9 km/s/kpc, where the error is dominated by uncertainty in the source radius. Likelihood analysis yields a median lens mass of 0.67{+2.53}{-0.46} Msun, located with 80% probability in the Galactic bulge at a distance of 6.93{+1.56}{-2.25} kpc. If the lens is a main-sequence star, we can include constraints on the lens luminosity. This modifies our estimates to M_lens = 0.53{+0.52}{-0.35} Msun and D_lens = 6.57{+0.99}{-2.25} kpc. Spectra taken during the event show that the absorption line equivalent widths of H alpha and the TiO bands near 6700 A vary, as predicted for microlensing of an extended source. This is most likely due to center-to-limb variation in the stellar spectral lines. These data demonstrate the feasibility of using microlensing limb crossings as a tool to probe stellar atmospheres directly.

We present the lightcurves of two microlensing events from the MACHO Project data that are likely to be due to lenses with masses similar to Jupiter's mass. Although the MACHO Project survey data are not sufficient to definitively establish the identification of planetary mass lenses in these cases, observations by microlensing follow-up networks such as GMAN and PLANET should be able to definitively determine the planetary nature of similar events which may occur in the near future.

The MACHO project is searching for dark matter in the form of massive compact halo objects (Machos), by monitoring the brightness of millions of stars in the Magellanic Clouds to search for gravitational microlensing events. Analysis of our first 2.3 years of data for 8.5 million stars in the LMC yields 8 candidate microlensing events, well in excess of the $\approx 1$ event expected from lensing by known low-mass stars. The event timescales range from 34 to 145 days, and the estimated optical depth is $\sim$ 2 x 10^{-7}, about half of that expected from a `standard' halo. Likelihood analysis indicates the typical lens mass is $0.5^{+0.3}_{-0.2} \Msun$, suggesting they may be old white dwarfs.

The approximately 1500 Cepheid variables in 22 LMC fields of the MACHO Project survey have been calibrated and analysed. In this paper, we report improved period ratios for a total of 73 beat Cepheids and provide a first look at the Fourier decomposition parameters for both singly- and doubly-periodic Cepheids. We also note an unusual amplitude-changing pulsator.

The MACHO Project has catalogued $\sim 8000$ RR Lyrae stars in the Large Magellanic Cloud, $\sim 1800$ in the Galactic bulge, and $\sim$ 50 in the Sgr dwarf galaxy. These variables are excellent distance indicators, and are used as tools to study the structure of the Large Magellanic Cloud and the bulge. The large datasets also probe uncommon pulsation modes. A number of double-mode RR Lyrae stars (RRd) are found in the Large Magellanic Cloud sample. These stars provide important clues for understanding the formation and evolution of the inner Galaxy, the Large Magellanic Cloud and the Sgr dwarf galaxy. A large number of second overtone pulsators (RRe) are found in the LMC and bulge. Finally, the RR Lyrae belonging to the Sgr dwarf yield an accurate distance to this galaxy. Their presence also alerts us of the very interesting possibility of distant sources for bulge microlensing events.

We report the discovery of 73 double-mode RR Lyrae (RRd) stars in fields near the bar of the LMC. The stars are detected among the MACHO database of short-period variables that currently contains about 7900 RR Lyrae stars. Fundamental periods (P_0) for these stars are found in the range 0.46-0.55 days and first overtone-to-fundamental period ratios are found to be in the range 0.742 < P_1/P_0 < 0.748. A significant fraction of our current sample have period ratios smaller than any previously discovered RRd variables. We present mean magnitudes, colors, and lightcurve properties for all LMC RRd stars detected to date. The range in period ratios is unexpectedly large. We present a determination of absolute magnitudes for these stars based primarily on pulsation theory and the assumption that all observed stars are at the fundamental blue edge (FBE) of the instability strip. Comparison of the calibrated MACHO V and R_KC photometry with these derived absolute magnitudes yields an absorption-corrected distance modulus to the LMC of 18.57 +/- 0.19 mag which is in good agreement with that found (18.5) through comparison of galactic and LMC Cepheids. Adopting this luminosity calibration, we derive an increase in the distance modulus, and thus a reduction in the age found via isochrone fitting for M15 of about 33% and discuss the implications for cosmology.

The MACHO Project is a search for dark matter in the form of massive compact halo objects (Machos). Photometric monitoring of millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge is used to search for gravitational microlensing events caused by these otherwise invisible objects. Analysis of the first 2.1 years of photometry of 8.5 million stars in the LMC reveals 8 candidate microlensing events. This is substantially more than the number expected ($\sim 1.1$) from lensing by known stellar populations. The timescales ($\that$) of the events range from 34 to 145 days. We estimate the total microlensing optical depth towards the LMC from events with $2 < \that < 200$ days to be $\tau_2^{200} = 2.9 ^{+1.4}_{-0.9} \ten{-7}$ based upon our 8 event sample. This exceeds the optical depth, $\tau_{\rm backgnd} = 0.5 \ten{-7}$, expected from known stars, and the difference is to be compared with the optical depth predicted for a ``standard" halo composed entirely of Machos: $\tau_{halo} = 4.7\ten{-7}$. Likelihood analysis gives a fairly model independent estimate of the halo mass in Machos within 50 kpc of $2.0^{+1.2}_{-0.7} \ten{11} \msun$, about half of the ``standard halo" value. We also find a most probable Macho mass of $0.5^{+0.3}_{-0.2}\msun$, although this value is strongly model dependent. Additionally, the absence of short duration events places stringent upper limits on the contribution of low-mass Machos: objects from $10^{-4} \msun$ to $0.03 \msun$ contribute $\simlt 20\%$ of the ``standard" dark halo.

The MACHO Project is searching for galactic dark matter in the form of massive compact halo objects (Machos). Millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge are photometrically monitored in an attempt to detect rare gravitational microlensing events caused by otherwise invisible Machos. Analysis of two years of photometry on 8.5 million stars in the LMC reveals 8 candidate microlensing events, far more than the $\sim1$ event expected from lensing by low-mass stars in known galactic populations. From these eight events we estimate the optical depth towards the LMC from events with $2 < \that < 200$ days to be $\tau_2^{200} \approx 2.9 ^{+1.4}_{-0.9} \ten{-7}$. This exceeds the optical depth of $0.5\ten{-7}$ expected from known stars and is to be compared with an optical depth of $4.7\ten{-7}$ predicted for a ``standard'' halo composed entirely of Machos. The total mass in this lensing population is $\approx 2^{+1.2}_{-0.7} \ten{11} \msun$ (within 50 kpc from the Galactic center). Event timescales yield a most probable Macho mass of $0.5^{+0.3}_{-0.2}\msun$, although this value is quite model dependent.

The MACHO collaboration has recently analyzed 2.1 years of photometric data for about 8.5 million stars in the Large Magellanic Cloud (LMC). This analysis has revealed 8 candidate microlensing events and a total microlensing optical depth of $\tau_{meas} = 2.9 +1.4/-0.9 \times 10^{-7}$. This significantly exceeds the number of events (1.1) and the microlensing optical depth predicted from known stellar populations: $\tau_{back} = 5.4\times 10^{-8}$, but it is consistent with models in which about half of the standard dark halo mass is composed of Machos of mass $\sim 0.5 \msun$. One of these 8 events appears to be a binary lensing event with a caustic crossing that is partially resolved which allows us to estimate the distance to the lenses. If the source star is not a short period binary star, then we show that the lens system is very likely to reside in the LMC. However, if we assume that the optical depth for LMC-LMC lensing is large enough to account for our entire lensing signal, then the binary event does not appear to be consistent with lensing of a single LMC source star by a binary residing in the LMC. Thus, while the binary lens may indeed reside in the LMC, there is no indication that most of the lenses reside in the LMC.

We report the discovery of 30 type a,b RR Lyrae (RRab) which are likely members of the Sagittarius (Sgr) dwarf galaxy. Accurate positions, periods, amplitudes and magnitudes are presented. Their distances are determined with respect to RRab in the Galactic bulge found also in the MACHO 1993 data. For R$_{\odot} = 8$ kpc, the mean distance to these stars is $D = 22 \pm 1$ kpc, smaller than previous determinations for this galaxy. This indicates that Sgr has an elongated main body extending for more than 10 kpc, which is inclined along the line of sight, with its northern part (in Galactic coordinates) closer to us. The size and shape of Sgr give clues about the past history of this galaxy. If the shape of Sgr follows the direction of its orbit, the observed spatial orientation suggests that Sgr is moving away from the Galactic plane. Also, Sgr stars may be the sources of some of the microlensing events seen towards the bulge.

We report the discovery of two new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The identification of both stars has been confirmed spectroscopically. One is a cool RCB star (T_eff about 5000 K) characterized by very strong Swan bands of C_2 and violet bands of CN, and weak or absent Balmer lines, G-band and 12C-13C bands. The second star is an example of a hot RCB star of which only 3 were previously known to exist in the Galaxy and none in the LMC. Its spectrum is characterized by several C II lines in emission. Both stars have shown deep declines of Delta V > 4 mag in brightness. The new stars are significantly fainter at maximum light than the three previously known LMC RCB stars. The amount of reddening toward these stars is somewhat uncertain but both seem to have absolute magnitudes, M_V, about half a magnitude fainter than the other three stars. Estimates of M_Bol find that the hot RCB star lies in the range of the other three stars while the cool RCB star is fainter. The two cool LMC RCB stars are the faintest at M_Bol. The discovery of these two new stars brings to five the number of known RCB stars in the LMC and demonstrates the utility of the MACHO photometric database for the discovery of new RCB stars.

The MACHO project has been monitoring about ten million stars in the Large Magellanic Cloud in the search for gravitational microlensing events caused by massive compact halo objects (Machos) in the halo of the Milky Way. In our standard analysis, we have searched this data set for well sampled, long duration microlensing lightcurves, detected several microlensing events consistent with Machos in the 0.1 < m < 1.0 M_sun mass range, and set limits on the abundance of objects with masses 1e-5 < m < 0.1 M_sun. In this paper, we present a different type of analysis involving the search for very short time scale brightenings of stars which is used to set strong limits on the abundance of lower mass Machos. Our analysis of the first two years of data toward the LMC indicates that Machos with masses in the range 2.5e-7 < m < 5.2e-4 M_sun cannot make up the entire mass of a standard spherical dark halo. Combining these results with those from the standard analysis, we find that the halo dark matter may not be comprised of objects with masses 2.5e-7 < m < 8.1e-2 M_sun.

Real-time detection of microlensing has moved from proof of concept in 1994 to a steady stream of events this year. Global dissemination of these events by the MACHO and OGLE collaborations has made possible intensive photometric and spectroscopic followup from widely dispersed sites confirming the microlensing hypothesis. Improved photometry and increased temporal resolution from followup observations greatly increases the possibility of detecting deviations from the standard point-source, point-lens, inertial motion microlensing model. These deviations are crucial in understanding individual lensing systems by breaking the degeneracy between lens mass, position and velocity. We report here on GMAN (Global Microlensing Alert Network), the coordinated followup of MACHO alerts.