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 a new analysis of the Jupiter+Saturn analog system, OGLE-2006-BLG-109Lb,c, which was the first double planet system discovered with the gravitational microlensing method. This is the only multi-planet system discovered by any method with measured masses for the star and both planets. In addition to the signatures of two planets, this event also exhibits a microlensing parallax signature and finite source effects that provide a direct measure of the masses of the star and planets, and the expected brightness of the host star is confirmed by Keck AO imaging, yielding masses of M_* = 0.51(+0.05-0.04) M_sun, M_b = 231+-19 M_earth, M_c = 86+-7 M_earth. The Saturn-analog planet in this system had a planetary light curve deviation that lasted for 11 days, and as a result, the effects of the orbital motion are visible in the microlensing light curve. We find that four of the six orbital parameters are tightly constrained and that a fifth parameter, the orbital acceleration, is weakly constrained. No orbital information is available for the Jupiter-analog planet, but its presence helps to constrain the orbital motion of the Saturn-analog planet. Assuming co-planar orbits, we find an orbital eccentricity of eccentricity = 0.15 (+0.17-0.10) and an orbital inclination of i = 64 (+4-7) deg. The 95% confidence level lower limit on the inclination of i > 49 deg. implies that this planetary system can be detected and studied via radial velocity measurements using a telescope of >30m aperture.

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.

Cosmology with Type Ia supernovae heretofore has required extensive spectroscopic follow-up to establish a redshift. Though tolerable at the present discovery rate, the next generation of ground-based all-sky survey instruments will render this approach unsustainable. Photometry-based redshift determination is a viable alternative, but introduces non-negligible errors that ultimately degrade the ability to discriminate between competing cosmologies. We present a strictly template-based photometric redshift estimator and compute redshift reconstruction errors in the presence of photometry and statistical errors. With reasonable assumptions for a cadence and supernovae distribution, these redshift errors are combined with systematic errors and propagated using the Fisher matrix formalism to derive lower bounds on the joint errors in $\Omega_w$ and $\Omega_w'$ relevant to the next generation of ground-based all-sky survey.

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.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of ~0.71 and ~0.27 times the mass of Jupiter and orbital separations of ~2.3 and ~4.6 astronomical units orbiting a primary star of mass ~0.50 solar masses at a distance of ~1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.

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 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.

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.

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.

(abridged) Recent studies, using OGLE data for LMC Cepheids in the optical, strongly suggest that the period-luminosity (PL) relation for the Large Magellanic Cloud (LMC) Cepheids shows a break or non-linearity at a period of 10 days. In this paper we apply statistical tests, the chi-square test and the F-test, to the Cepheid data from the MACHO project to test for a non-linearity of the V- and R-band PL relations at 10 days, and extend these tests to the near infrared (JHK-band) PL relations with 2MASS data. We correct the extinction for these data by applying an extinction map towards the LMC. The statistical test we use, the F-test, is able to take account of small numbers of data points and the nature of that data on either side of the period cut at 10 days. With our data, the results we obtained imply that the VRJH-band PL relations are non-linear around a period of 10 days, while the K-band PL relation is (marginally) consistent with a single-line regression. The choice of a period of 10 days, around which this non-linearity occurs, is consistent with the results obtained when this "break" period is estimated from the data. Long period Cepheids are supplemented from the literature to increase our sample size. The photometry of these long period Cepheids is compared with our data and no trend with period is found. Our main results remain unchanged when we supplement our dataset with these long period Cepheids. By examining our data at maximum light, we also suggest arguments why errors in reddening are unlikely to be responsible for our results. The non-linearity of the mean V-band PL relation as seen in both of the OGLE and MACHO data, using different extinction maps, suggests that this non-linearity is real.

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.

In this paper, we recalibrate the Cepheid distance to some nearby galaxies observed by the HST Key Project (KP) and the Sandage-Tammann-Saha (STS) group. We use much of the KP methodology in our analysis but apply new techniques, based on Fourier methods to estimate the mean magnitudes to published Cepheid data. We also apply different calibrating PL relations to estimate Cepheid distances, and investigate the sensitivity of the distance moduli to the adopted calibrating PL relation. We re-determine the OGLE LMC PL relations using a more conservative approach and also study the effect of using Galactic PL relations on the distance scale. For the KP galaxies, we find good agreement with an average discrepancy of -0.002 and 0.075mag. when using the LMC and Galaxy, respectively, as a calibrating PL relation. However, the Cepheid distance to STS galaxies are shorter by -0.178mag. and -0.108mag. on average again when using the LMC and Galaxy PL relations, respectively. We also calculate the distance to these galaxies using PL relations at maximum light and find very good agreement with mean light PL distances. However, after correcting for metallicity effects, the difference between the distance moduli obtained using the two sets of calibrating PL relations becomes negligible. This suggests that Cepheids in the LMC and Galaxy do follow different PL relations, at least at the median periods of the target galaxies.

Fourier decomposition is a well established technique used in stellar pulsation. However the quality of reconstructed light curves using this method is reduced when the observed data have uneven phase coverage. We use simulated annealing techniques together with Fourier decomposition to improve the quality of the Fourier decomposition for many OGLE LMC fundamental mode Cepheids. This method restricts the range that Fourier amplitudes can take. The ranges are specified by well sampled Cepheids in the Galaxy and Magellanic Clouds. We also apply this method to reconstructing Cepheid light curves observed by the HST. These typically consist of 12 V band and 4 I band points. We employ a direct Fourier fit to the 12 V band points using the simulated annealing method mentioned above and explicitly derive and use Fourier interrelations to reconstruct the I band light curve. We discuss advantages and drawbacks of this method when applied to HST Cepheid data over existing template methods. Application of this method to reconstruct the light curves of Cepheids observed in NGC 4258 shows that the derived Cepheid distance ($\mu_0=29.38\pm0.06$, random error) is consistent with its geometrical distance ($\mu_0=29.28\pm0.09$) derived from observations of water maser.

We present results from the application of a global photometric calibration (GPC) procedure to calibration data from the first 2 years of The Two Micron All Sky Survey (2MASS). The GPC algorithm uses photometry of both primary standards and moderately bright `tracer' stars in 35 2MASS calibration fields. During the first two years of the Survey, each standard was observed on approximately 50 nights, with about 900 individual measurements. Based on the photometry of primary standard stars and secondary tracer stars and under the assumption that the nightly zeropoint drift is linear, GPC ties together all calibration fields and all survey nights simultaneously, producing a globally optimized solution. Calibration solutions for the Northern and Southern hemisphere observatories are found separately, and are tested for global consistency based on common fields near the celestial equator. Several results from the GPC are presented, including establishing candidate secondary standards, monitoring of near-infrared atmospheric extinction coefficients, and verification of global validity of the standards. The solution gives long-term averages of the atmospheric extinction coefficients, A_J=0.096, A_H=0.026, A_{K_s}=0.066 (North) and A_J=0.092, A_H=0.031, A_{K_s}=0.065 (South), with formal error of 0.001. The residuals show small seasonal variations, most likely due to changing atmospheric content of water vapor. Extension of the GPC to approximately 100 field stars in each of the 35 calibration fields yields a catalog of more than two thousand photometric standards ranging from 10th to 14th magnitude, with photometry that is globally consistent to $\sim 1%$.

We derive structural parameters and evidence for extended tidal debris from star count and preliminary standard candle analyses of the Large Magellanic Cloud based on Two Micron All Sky Survey (2MASS) data. The full-sky coverage and low extinction in K_s presents an ideal sample for structural analysis of the LMC. The star count surface densities and deprojected inclination for both young and older populations are consistent with previous work. We use the full areal coverage and large LMC diameter to Galactrocentric distance ratio to infer the same value for the disk inclination based on perspective. A standard candle analysis based on a sample of carbon long-period variables (LPV) in a narrow color range, 1.6<J-K_s<1.7 allows us to probe the three-dimensional structure of the LMC along the line of sight. The intrinsic brightness distribution of carbon LPVs in selected fields implies that $\sigma_M\simlt 0.2^m$ for this color cut. The sample provides a {\it direct} determination of the LMC disk inclination: $42.3^\circ\pm 7.2^\circ$. Distinct features in the photometric distribution suggest several distinct populations. We interpret this as the presence of an extended stellar component of the LMC, which may be as thick as 14 kpc, and intervening tidal debris at roughly 15 kpc from the LMC.

We present a morphological analysis of the feature-rich 2MASS LMC color-magnitude diagram, identifying Galactic and LMC populations and estimating the density of LMC populations alone. We also present the projected spatial distributions of various stellar populations. Major populations are identified based on matching morphological features of the CMD with expected positions of known populations, isochrone fits, and analysis of the projected spatial distributions. The LMC populations along the first-ascent RGB and AGB are quantified. We find the RGB tip at $K_s=12.3\pm0.1$. Preliminary isochrone analysis is done for giant populations in the bar and the outer regions of the Cloud. We find no significant differences in metallicities and ages between the fields. The observed LMC giant branch is well-fit by published tracks in the CIT/CTIO system with a distance modulus of $\mu=18.5\pm0.1$, reddening $E_{B-V}=0.15-0.20$, metallicity $Z=0.004^{+0.002}_{-0.001}$ and age 3-13 Gyr. Analysis of deep 2MASS engineering data with six times the standard exposure produces similar estimates.

Previous work reported a bar signature in color-selected IRAS variable stars. Here, we estimate the source density of these variables while consistently accounting for spatial incompleteness in data using a likelihood approach. The existence of the bar is confirmed with shoulder of 4 kpc and position angle of 19 +/- 1 degrees. The ratio of non-axisymmetric to axisymmetric components gives similar estimate for the bar size 3.3 +/- 0.1 kpc and position angle 24 +/- 2 degrees. We estimate a scale length 4.00 +/- 0.55 kpc for the IRAS variable population, suggesting that these stars represent the old disk population. We use this density reconstruction to estimate the optical depth to microlensing for the large-scale bar in the Galactic disk. We find an enhancement over an equivalent axisymmetric disk by 30% but still too small to account for the MACHO result. In addition, we predict a significant asymmetry at positive and negative longitudes along lines of sight through the end of the bar with optical depths comparable to that in Baade's window. An infrared microlensing survey may be a sensitive tool for detecting or constraining structural asymmetries. More generally, this is a pilot study for Bayesian star count analyses. Bayesian approach allows the assessment of prior probabilities to the unknown parameters of the model; the resulting likelihood function is straightforwardly modified to incorporate all available data.