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

Using the most recent prototypes, design, and as-built system information, we test and quantify the capability of the Large Synoptic Survey Telescope (LSST) to discover Potentially Hazardous Asteroids (PHAs) and Near-Earth Objects (NEOs). We empirically estimate an expected upper limit to the false detection rate in LSST image differencing, using measurements on DECam data and prototype LSST software and find it to be about $450$~deg$^{-2}$. We show that this rate is already tractable with current prototype of the LSST Moving Object Processing System (MOPS) by processing a 30-day simulation consistent with measured false detection rates. We proceed to evaluate the performance of the LSST baseline survey strategy for PHAs and NEOs using a high-fidelity simulated survey pointing history. We find that LSST alone, using its baseline survey strategy, will detect $66\%$ of the PHA and $61\%$ of the NEO population objects brighter than $H=22$, with the uncertainty in the estimate of $\pm5$ percentage points. By generating and examining variations on the baseline survey strategy, we show it is possible to further improve the discovery yields. In particular, we find that extending the LSST survey by two additional years and doubling the MOPS search window increases the completeness for PHAs to $86\%$ (including those discovered by contemporaneous surveys) without jeopardizing other LSST science goals ($77\%$ for NEOs). This equates to reducing the undiscovered population of PHAs by additional $26\%$ ($15\%$ for NEOs), relative to the baseline survey.

We present new photometric data of 11 hot Jupiter transiting exoplanets (CoRoT-12b, HAT-P-5b, HAT-P-12b, HAT-P-33b, HAT-P-37b, WASP-2b, WASP-24b, WASP-60b, WASP-80b, WASP-103b, XO-3b) in order to update their planetary parameters and to constrain information about their atmospheres. These observations of CoRoT-12b, HAT-P-37b and WASP-60b are the first follow-up data since their discovery. Additionally, the first near-UV transits of WASP-80b and WASP-103b are presented. We compare the results of our analysis with previous work to search for transit timing variations (TTVs) and a wavelength dependence in the transit depth. TTVs may be evidence of a third body in the system and variations in planetary radius with wavelength can help constrain the properties of the exoplanet's atmosphere. For WASP-103b and XO-3b, we find a possible variation in the transit depths that may be evidence of scattering in their atmospheres. The B-band transit depth of HAT-P-37b is found to be smaller than its near-IR transit depth and such a variation may indicate TiO/VO absorption. These variations are detected from 2-4.6$\sigma$, so follow-up observations are needed to confirm these results. Additionally, a flat spectrum across optical wavelengths is found for 5 of the planets (HAT-P-5b, HAT-P-12b, WASP-2b, WASP-24b, WASP-80b), suggestive that clouds may be present in their atmospheres. We calculate a refined orbital period and ephemeris for all the targets, which will help with future observations. No TTVs are seen in our analysis with the exception of WASP-80b and follow-up observations are needed to confirm this possible detection.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems.

Transits of exoplanets observed in the near-UV have been used to study the scattering properties of their atmospheres and possible star-planet interactions. We observed the primary transits of 15 exoplanets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-16b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-12b, WASP-33b, WASP-36b, WASP-44b, WASP-48b, and WASP-77Ab) in the near-UV and several optical photometric bands to update their planetary parameters, ephemerides, search for a wavelength dependence in their transit depths to constrain their atmospheres, and determine if asymmetries are visible in their light curves. Here we present the first ground-based near-UV light curves for 12 of the targets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-33b, WASP-36b, WASP-48b, and WASP-77Ab). We find that none of the near-UV transits exhibit any non-spherical asymmetries, this result is consistent with recent theoretical predictions by Ben-Jaffel et al. and Turner et al. The multi-wavelength photometry indicates a constant transit depth from near-UV to optical wavelengths in 10 targets (suggestive of clouds), and a varying transit depth with wavelength in 5 targets (hinting at Rayleigh or aerosol scattering in their atmospheres). We also present the first detection of a smaller near-UV transit depth than that measured in the optical in WASP-1b and a possible opacity source that can cause such radius variations is currently unknown. WASP-36b also exhibits a smaller near-UV transit depth at 2.6$\sigma$. Further observations are encouraged to confirm the transit depth variations seen in this study.

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

The Taiwanese-American Occultation Survey (TAOS) aims to detect serendipitous occultations of stars by small (about 1 km diameter) objects in the Kuiper Belt and beyond. Such events are very rare (<0.001 events per star per year) and short in duration (about 200 ms), so many stars must be monitored at a high readout cadence. TAOS monitors typically around 500 stars simultaneously at a 5 Hz readout cadence with four telescopes located at Lulin Observatory in central Taiwan. In this paper, we report the results of the search for small Kuiper Belt Objects (KBOs) in seven years of data. No occultation events were found, resulting in a 95% c.l. upper limit on the slope of the faint end of the KBO size distribution of q = 3.34 to 3.82, depending on the surface density at the break in the size distribution at a diameter of about 90 km.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

We report the extremely high magnification (A > 1000) binary microlensing event OGLE-2007-BLG-514. We obtained good coverage around the double peak structure in the light curve via follow-up observations from different observatories. The binary lens model that includes the effects of parallax (known orbital motion of the Earth) and orbital motion of the lens yields a binary lens mass ratio of q = 0.321 +/- 0.007 and a projected separation of s = 0.072 +/- 0.001$ in units of the Einstein radius. The parallax parameters allow us to determine the lens distance D_L = 3.11 +/- 0.39 kpc and total mass M_L=1.40 +/- 0.18 M_sun; this leads to the primary and secondary components having masses of M_1 = 1.06 +/- 0.13 M_sun and M_2 = 0.34 +/- 0.04 M_sun, respectively. The parallax model indicates that the binary lens system is likely constructed by the main sequence stars. On the other hand, we used a Bayesian analysis to estimate probability distributions by the model that includes the effects of xallarap (possible orbital motion of the source around a companion) and parallax (q = 0.270 +/- 0.005, s = 0.083 +/- 0.001). The primary component of the binary lens is relatively massive with M_1 = 0.9_{-0.3}^{+4.6} M_sun and it is at a distance of D_L = 2.6_{-0.9}^{+3.8} kpc. Given the secure mass ratio measurement, the companion mass is therefore M_2 = 0.2_{-0.1}^{+1.2} M_sun. The xallarap model implies that the primary lens is likely a stellar remnant, such as a white dwarf, a neutron star or a black hole.

The mass function and statistics of binaries provide important diagnostics of the star formation process. Despite this importance, the mass function at low masses remains poorly known due to observational difficulties caused by the faintness of the objects. Here we report the microlensing discovery and characterization of a binary lens composed of very low-mass stars just above the hydrogen-burning limit. From the combined measurements of the Einstein radius and microlens parallax, we measure the masses of the binary components of $0.10\pm 0.01\ M_\odot$ and $0.09\pm 0.01\ M_\odot$. This discovery demonstrates that microlensing will provide a method to measure the mass function of all Galactic populations of very low mass binaries that is independent of the biases caused by the luminosity of the population.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

Most known extrasolar planets (exoplanets) have been discovered using the radial velocity$^{\bf 1,2}$ or transit$^{\bf 3}$ methods. Both are biased towards planets that are relatively close to their parent stars, and studies find that around 17--30% (refs 4, 5) of solar-like stars host a planet. Gravitational microlensing$^{\bf 6\rm{\bf -}\bf 9}$, on the other hand, probes planets that are further away from their stars. Recently, a population of planets that are unbound or very far from their stars was discovered by microlensing$^{\bf 10}$. These planets are at least as numerous as the stars in the Milky Way$^{\bf 10}$. Here we report a statistical analysis of microlensing data (gathered in 2002--07) that reveals the fraction of bound planets 0.5--10 AU (Sun--Earth distance) from their stars. We find that 17$_{\bf -9}^{\bf +6}$% of stars host Jupiter-mass planets (0.3--10 $\MJ$, where $\MJ {\bf = 318}$ $\Mearth$ and $\Mearth$ is Earth's mass). Cool Neptunes (10--30 $\Mearth$) and super-Earths (5--10 $\Mearth$) are even more common: their respective abundances per star are 52$_{\bf -29}^{\bf +22}$% and 62$_{\bf -37}^{\bf +35}$%. We conclude that stars are orbited by planets as a rule, rather than the exception.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

We present the analysis result of a gravitational binary-lensing event OGLE-2005-BLG-018. The light curve of the event is characterized by 2 adjacent strong features and a single weak feature separated from the strong features. The light curve exhibits noticeable deviations from the best-fit model based on standard binary parameters. To explain the deviation, we test models including various higher-order effects of the motions of the observer, source, and lens. From this, we find that it is necessary to account for the orbital motion of the lens in describing the light curve. From modeling of the light curve considering the parallax effect and Keplerian orbital motion, we are able to measure not only the physical parameters but also a complete orbital solution of the lens system. It is found that the event was produced by a binary lens located in the Galactic bulge with a distance $6.7\pm 0.3$ kpc from the Earth. The individual lens components with masses $0.9\pm 0.3\ M_\odot$ and $0.5\pm 0.1\ M_\odot$ are separated with a semi-major axis of $a=2.5 \pm 1.0$ AU and orbiting each other with a period $P=3.1 \pm 1.3$ yr. The event demonstrates that it is possible to extract detailed information about binary lens systems from well-resolved lensing light curves.

We present the first example of binary microlensing for which the parameter measurements can be verified (or contradicted) by future Doppler observations. This test is made possible by a confluence of two relatively unusual circumstances. First, the binary lens is bright enough (I=15.6) to permit Doppler measurements. Second, we measure not only the usual 7 binary-lens parameters, but also the 'microlens parallax' (which yields the binary mass) and two components of the instantaneous orbital velocity. Thus we measure, effectively, 6 'Kepler+1' parameters (two instantaneous positions, two instantaneous velocities, the binary total mass, and the mass ratio). Since Doppler observations of the brighter binary component determine 5 Kepler parameters (period, velocity amplitude, eccentricity, phase, and position of periapsis), while the same spectroscopy yields the mass of the primary, the combined Doppler + microlensing observations would be overconstrained by 6 + (5 + 1) - (7 + 1) = 4 degrees of freedom. This makes possible an extremely strong test of the microlensing solution. We also introduce a uniform microlensing notation for single and binary lenses, we define conventions, summarize all known microlensing degeneracies and extend a set of parameters to describe full Keplerian motion of the binary lenses.

The MPF mission will provide a statistical census of exoplanets with masses greater than 0.1 Earth-masses and orbital separations ranging from 0.5AU to infinity. This includes analogs to all the Solar System's planets except for Mercury, as well as most types of planets predicted by planet formation theories. Such a survey will provide results on the frequency of planets around all types of stars except those with short lifetimes. Close-in planets with separations < 0.5 AU are invisible to a space-based microlensing survey, but these can be found by Kepler. Other methods, including ground-based microlensing, cannot approach the comprehensive statistics on the mass and semi-major axis distribution of extrasolar planets that a space-based microlensing survey will provide. The terrestrial planet sensitivity of a ground-based microlensing survey is limited to the vicinity of the Einstein radius at 2-3 AU, and space-based imaging is needed to identify and determine the mass of the planetary host stars for the vast majority of planets discovered by microlensing. Thus, a space-based microlensing survey is likely to be the only way to gain a comprehensive understanding of the architecture of planetary systems, which is needed to understand planet formation and habitability. MPF can accomplish these objectives with proven technology and a cost of $333 million (excluding launch vehicle).

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.

Aims: We present a detailed analysis of OGLE 2004-BLG-482, a relatively high-magnification single-lens microlensing event which exhibits clear extended-source effects. These events are relatively rare, but they potentially contain unique information on the stellar atmosphere properties of their source star, as shown in this study. Methods: Our dense photometric coverage of the overall light curve and a proper microlensing modelling allow us to derive measurements of the OGLE 2004-BLG-482 source star's linear limb-darkening coefficients in three bands, including standard Johnson-Cousins I and R, as well as in a broad clear filter. In particular, we discuss in detail the problems of multi-band and multi-site modelling on the expected precision of our results. We also obtained high-resolution UVES spectra as part of a ToO programme at ESO VLT from which we derive the source star's precise fundamental parameters. Results: From the high-resolution UVES spectra, we find that OGLE 2004-BLG-482's source star is a red giant of MK type a bit later than M3, with Teff = 3667 +/- 150 K, log g = 2.1 +/- 1.0 and an assumed solar metallicity. This is confirmed by an OGLE calibrated colour-magnitude diagram. We then obtain from a detailed microlensing modelling of the light curve linear limb-darkening coefficients that we compare to model-atmosphere predictions available in the literature, and find a very good agreement for the I and R bands. In addition, we perform a similar analysis using an alternative description of limb darkening based on a principal component analysis of ATLAS limb-darkening profiles, and also find a very good agreement between measurements and model predictions.

[abridged] We describe the first results from the NOAO Outer Limits Survey. The survey consists of deep images of 55 0.6x0.6 degree fields at distances up to 20 degrees from the LMC/SMC, and 10 controls. The fields probe the outer structure of the Clouds, the Magellanic Stream, the Leading Arm, and the wake of the new LMC orbit. Images were taken in 5 filters on the CTIO Blanco 4-m and Mosaic2 camera, with calibration at the CTIO 0.9-m. The CRI images reach depths below the oldest LMC/SMC main sequence (MS) turnoffs, yielding probes of structure combined with ability to measure stellar ages and metallicities. M and DDO51 images allow for discrimination of LMC and SMC giant stars from foreground dwarfs, allowing us to use giants as additional probes. From photometry of 8 fields at radii of 7-19 degrees N of the LMC bar, we find MS stars associated with the LMC to 16 degrees from the LMC center, while the much rarer giants can only be convincingly detected out to 11 degrees. In one control field, we see the signature of the Milky Way globular cluster NGC 1851. The CMDs show that while at 7-degree radius LMC populations as young as 500 Myr are present, at radii more than 11 degrees only the underlying old metal-poor population remains, demonstrating the existence of a population gradient. Even at extreme distances, the dominant age is much less than "globular Cluster age." MS star counts follow an exponential decline with a scale length of 1.15 kpc, essentially the same scale length as gleaned for the inner LMC disk from prior studies. While we cannot rule out tidal structure elsewhere, detection of ordered structure to 12 disk scale lengths is unprecedented, and adds to the puzzle of the LMC's interaction history. Our results do not rule out the possible existence of an LMC stellar halo, which may only dominate the disk at yet larger radii.

We combine variability information from the MAssive Compact Halo Objects (MACHO) survey of the Large Magellanic Cloud (LMC) with infrared photometry from the Spitzer Space Telescope Surveying the Agents of a Galaxy's Evolution (SAGE) survey to create a dataset of ~30 000 variable red sources. We photometrically classify these sources as being on the first ascent of the Red Giant Branch (RGB), or as being in one of three stages along the Asymptotic Giant Branch (AGB): oxygen-rich, carbon-rich, or highly reddened with indeterminate chemistry ("extreme" AGB candidates). We present linear period-luminosity relationships for these sources using 8 separate infrared bands (J, H, K, 3.6, 4.5, 5.8, 8.0, and 24 micron) as proxies for the luminosity. We find that the wavelength dependence of the slope of the period-luminosity relationship is different for different photometrically determined classes of AGB stars. Stars photometrically classified as O-rich show the least variation of slope with wavelength, while dust enshrouded extreme AGB stars show a pronounced trend toward steeper slopes with increasing wavelength. We find that O-rich AGB stars pulsating in the fundamental mode obey a period-magnitude relation with a slope of -3.41 +/- 0.04 when magnitude is measured in the 3.6 micron band, in contrast to C-rich AGB stars, which obey a relation of slope -3.77 +/- 0.05.

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.

Gravitational microlensing is not only a successful tool for discovering distant exoplanets, but it also enables characterization of the lens and source stars involved in the lensing event. In high magnification events, the lens caustic may cross over the source disk, which allows a determination of the angular size of the source and additionally a measurement of its limb darkening. When such extended-source effects appear close to maximum magnification, the resulting light curve differs from the characteristic Paczynski point-source curve. The exact shape of the light curve close to the peak depends on the limb darkening of the source. Dense photometric coverage permits measurement of the respective limb-darkening coefficients. In the case of microlensing event OGLE 2008-BLG-290, the K giant source star reached a peak magnification of about 100. Thirteen different telescopes have covered this event in eight different photometric bands. Subsequent light-curve analysis yielded measurements of linear limb-darkening coefficients of the source in six photometric bands. The best-measured coefficients lead to an estimate of the source effective temperature of about 4700 +100-200 K. However, the photometric estimate from colour-magnitude diagrams favours a cooler temperature of 4200 +-100 K. As the limb-darkening measurements, at least in the CTIO/SMARTS2 V and I bands, are among the most accurate obtained, the above disagreement needs to be understood. A solution is proposed, which may apply to previous events where such a discrepancy also appeared.

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

The Taiwanese-American Occultation Survey (TAOS) project has collected more than a billion photometric measurements since 2005 January. These sky survey data-covering timescales from a fraction of a second to a few hundred days-are a useful source to study stellar variability. A total of 167 star fields, mostly along the ecliptic plane, have been selected for photometric monitoring with the TAOS telescopes. This paper presents our initial analysis of a search for periodic variable stars from the time-series TAOS data on one particular TAOS field, No. 151 (RA = 17$^{\rm h}30^{\rm m}6\fs$67, Dec = 27\degr17\arcmin 30\arcsec, J2000), which had been observed over 47 epochs in 2005. A total of 81 candidate variables are identified in the 3 square degree field, with magnitudes in the range 8 < R < 16. On the basis of the periodicity and shape of the lightcurves, 29 variables, 15 of which were previously unknown, are classified as RR Lyrae, Cepheid, delta Scuti, SX Phonencis, semi-regular and eclipsing binaries.

The Taiwanese-American Occultation Survey (TAOS) monitors fields of up to ~1000 stars at 5 Hz simultaneously with four small telescopes to detect occultation events from small (~1 km) Kuiper Belt Objects (KBOs). The survey presents a number of challenges, in particular the fact that the occultation events we are searching for are extremely rare and are typically manifested as slight flux drops for only one or two consecutive time series measurements. We have developed a statistical analysis technique to search the multi-telescope data set for simultaneous flux drops which provides a robust false positive rejection and calculation of event significance. In this paper, we describe in detail this statistical technique and its application to the TAOS data set.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

We analyzed data accumulated during 2005 and 2006 by the Taiwan-American Occultation Survey (TAOS) in order to detect short-period variable stars (periods of <~ 1 hour) such as delta Scuti. TAOS is designed for the detection of stellar occultation by small-size Kuiper Belt Objects (KBOs) and is operating four 50cm telescopes at an effective cadence of 5Hz. The four telescopes simultaneously monitor the same patch of the sky in order to reduce false positives. To detect short-period variables, we used the Fast Fourier Transform algorithm (FFT) inasmuch as the data points in TAOS light-curves are evenly spaced. Using FFT, we found 41 short-period variables with amplitudes smaller than a few hundredths of a magnitude and periods of about an hour, which suggest that they are low-amplitude delta Scuti stars (LADS). The light-curves of TAOS delta Scuti stars are accessible online at the Time Series Center website (http://timemachine.iic.harvard.edu)

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 the results of a search for occultation events by objects at distances between 100 and 1000 AU in lightcurves from the Taiwanese-American Occultation Survey (TAOS). We searched for consecutive, shallow flux reductions in the stellar lightcurves obtained by our survey between 7 February 2005 and 31 December 2006 with a total of $\sim4.5\times10^{9}$ three-telescope simultaneous photometric measurements. No events were detected, allowing us to set upper limits on the number density as a function of size and distance of objects in Sedna-like orbits, using simple models.

We are proposing to conduct a multicolor, synoptic infrared (IR) imaging survey of the Northern sky with a new, dedicated 6.5-meter telescope at San Pedro M\'artir (SPM) Observatory. This initiative is being developed in partnership with astronomy institutions in Mexico and the University of California. The 4-year, dedicated survey, planned to begin in 2017, will reach more than 100 times deeper than 2MASS. The Synoptic All-Sky Infrared (SASIR) Survey will reveal the missing sample of faint red dwarf stars in the local solar neighborhood, and the unprecedented sensitivity over such a wide field will result in the discovery of thousands of z ~ 7 quasars (and reaching to z > 10), allowing detailed study (in concert with JWST and Giant Segmented Mirror Telescopes) of the timing and the origin(s) of reionization. As a time-domain survey, SASIR will reveal the dynamic infrared universe, opening new phase space for discovery. Synoptic observations of over 10^6 supernovae and variable stars will provide better distance measures than optical studies alone. SASIR also provides significant synergy with other major Astro2010 facilities, improving the overall scientific return of community investments. Compared to optical-only measurements, IR colors vastly improve photometric redshifts to z ~ 4, enhancing dark energy and dark matter surveys based on weak lensing and baryon oscillations. The wide field and ToO capabilities will enable a connection of the gravitational wave and neutrino universe - with events otherwise poorly localized on the sky - to transient electromagnetic phenomena.

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 analyze OGLE-2007-BLG-050, a high magnification microlensing event (A ~ 432) whose peak occurred on 2 May, 2007, with pronounced finite-source and parallax effects. We compute planet detection efficiencies for this event in order to determine its sensitivity to the presence of planets around the lens star. Both finite-source and parallax effects permit a measurement of the angular Einstein radius \theta_E = 0.48 +/- 0.01 mas and the parallax \pi_E = 0.12 +/- 0.03, leading to an estimate of the lens mass M = 0.50 +/- 0.14 M_Sun and its distance to the observer D_L = 5.5 +/- 0.4 kpc. This is only the second determination of a reasonably precise (<30%) mass estimate for an isolated unseen object, using any method. This allows us to calculate the planetary detection efficiency in physical units (r_\perp, m_p), where r_\perp is the projected planet-star separation and m_p is the planet mass. When computing planet detection efficiency, we did not find any planetary signature and our detection efficiency results reveal significant sensitivity to Neptune-mass planets, and to a lesser extent Earth-mass planets in some configurations. Indeed, Jupiter and Neptune-mass planets are excluded with a high confidence for a large projected separation range between the planet and the lens star, respectively [0.6 - 10] and [1.4 - 4] AU, and Earth-mass planets are excluded with a 10% confidence in the lensing zone, i.e. [1.8 - 3.1] AU.

We combine all available information to constrain the nature of OGLE-2005-BLG-071Lb, the second planet discovered by microlensing and the first in a high-magnification event. These include photometric and astrometric measurements from Hubble Space Telescope, as well as constraints from higher order effects extracted from the ground-based light curve, such as microlens parallax, planetary orbital motion and finite-source effects. Our primary analysis leads to the conclusion that the host of Jovian planet OGLE-2005-BLG-071Lb is an M dwarf in the foreground disk with mass M= 0.46 +/- 0.04 Msun, distance D_l = 3.3 +/- 0.4 kpc, and thick-disk kinematics v_LSR ~ 103 km/s. From the best-fit model, the planet has mass M_p = 3.8 +/- 0.4 M_Jup, lies at a projected separation r_perp = 3.6 +/- 0.2 AU from its host and so has an equilibrium temperature of T ~ 55 K, i.e., similar to Neptune. A degenerate model less favored by \Delta\chi^2 = 2.1 (or 2.2, depending on the sign of the impact parameter) gives similar planetary mass M_p = 3.4 +/- 0.4 M_Jup with a smaller projected separation, r_\perp = 2.1 +/- 0.1 AU, and higher equilibrium temperature T ~ 71 K. These results from the primary analysis suggest that OGLE-2005-BLG-071Lb is likely to be the most massive planet yet discovered that is hosted by an M dwarf. However, the formation of such high-mass planetary companions in the outer regions of M-dwarf planetary systems is predicted to be unlikely within the core-accretion scenario. There are a number of caveats to this primary analysis, which assumes (based on real but limited evidence) that the unlensed light coincident with the source is actually due to the lens, that is, the planetary host. However, these caveats could mostly be resolved by a single astrometric measurement a few years after the event.

The Taiwanese-American Occultation Survey (TAOS) operates four fully automatic telescopes to search for occultations of stars by Kuiper Belt Objects. It is a versatile facility that is also useful for the study of initial optical GRB afterglows. This paper provides a detailed description of the TAOS multi-telescope system, control software, and high-speed imaging.

A space-based gravitational microlensing exoplanet survey will provide a statistical census of exoplanets with masses greater than 0.1 Earth-masses and orbital separations ranging from 0.5AU to infinity. This includes analogs to all the Solar System's planets except for Mercury, as well as most types of planets predicted by planet formation theories. Such a survey will provide results on the frequency of planets around all types of stars except those with short lifetimes. Close-in planets with separations < 0.5 AU are invisible to a space-based microlensing survey, but these can be found by Kepler. Other methods, including ground-based microlensing, cannot approach the comprehensive statistics on the mass and semi-major axis distribution of extrasolar planets that a space-based microlensing survey will provide. The terrestrial planet sensitivity of a ground-based microlensing survey is limited to the vicinity of the Einstein radius at 2-3 AU, and space-based imaging is needed to identify and determine the mass of the planetary host stars for the vast majority of planets discovered by microlensing. Thus, a space-based microlensing survey is likely to be the only way to gain a comprehensive understanding of the architecture of planetary systems, which is needed to understand planet formation and habitability. The proposed Microlensing Planet Finder (MPF) mission is an example of a space-based microlensing survey that can accomplish these objectives with proven technology and a cost of under $300 million (excluding launch vehicle).

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?

The star BD+29 1748 was resolved to be a close binary from its occultation by the asteroid 87 Sylvia on 2006 December 18 UT. Four telescopes were used to observe this event at two sites separated by some 80 km apart. Two flux drops were observed at one site, whereas only one flux drop was detected at the other. From the long-term variation of Sylvia, we inferred the probable shape of the shadow during the occultation, and this in turn constrains the binary parameters: the two components of BD+29 1748 have a projected separation of 0.097" to 0.110" on the sky with a position angle 104 deg to 107 deg. The asteroid was clearly resolved with a size scale ranging from 130 to 290 km, as projected onto the occultation direction. No occultation was detected for either of the two known moonlets of 87 Sylvia.

We outline a method for fitting binary-lens caustic-crossing microlensing events based on the alternative model parameterisation proposed and detailed in Cassan (2008). As an illustration of our methodology, we present an analysis of OGLE-2007-BLG-472, a double-peaked Galactic microlensing event with a source crossing the whole caustic structure in less than three days. In order to identify all possible models we conduct an extensive search of the parameter space, followed by a refinement of the parameters with a Markov Chain-Monte Carlo algorithm. We find a number of low-chi2 regions in the parameter space, which lead to several distinct competitive best models. We examine the parameters for each of them, and estimate their physical properties. We find that our fitting strategy locates several minima that are difficult to find with other modelling strategies and is therefore a more appropriate method to fit this type of events.

Results from the first two years of data from the Taiwanese-American Occultation Survey (TAOS) are presented. Stars have been monitored photometrically at 4 Hz or 5 Hz to search for occultations by small (~3 km) Kuiper Belt Objects (KBOs). No statistically significant events were found, allowing us to present an upper bound to the size distribution of KBOs with diameters 0.5 km < D < 28 km.

We present a new analysis of the long-period variables in the Large Magellanic Cloud (LMC) from the MACHO Variable Star Catalog. Three-quarters of our sample of evolved, variable stars have periodic light curves. We characterize the stars in our sample using the multiple periods found in their frequency spectra. Additionally, we use single-epoch Two Micron All Sky Survey measurements to construct the average infrared light curves for different groups of these stars. Comparison with evolutionary models shows that stars on the red giant branch (RGB) or the early asymptotic giant branch (AGB) often show non-periodic variability, but begin to pulsate with periods on the two shortest period-luminosity sequences (3 & 4) when they brighten to K_s ~ 13. The stars on the thermally pulsing AGB are more likely to pulsate with longer periods that lie on the next two P-L sequences (1 & 2), including the sequence associated with the Miras in the LMC. The Petersen diagram and its variants show that multi-periodic stars on each pair of these sequences (3 & 4, and 1 & 2) typically pulsate with periods associated only with that pair. The periods in these multi-periodic stars become longer and stronger as the star evolves. We further constrain the mechanism behind the long secondary periods (LSPs) seen in half of our sample, and find that there is a close match between the luminosity functions of the LSP stars and all of the stars in our sample, and that these star's pulsation amplitudes are relatively wavelength independent. Although this is characteristic of stellar multiplicity, the large number of these variables is problematic for that explanation.

Thirteen exo-planets have been discovered using the gravitational microlensing technique (out of which 7 have been published). These planets already demonstrate that super-Earths (with mass up to ~10 Earth masses) beyond the snow line are common and multiple planet systems are not rare. In this White Paper we introduce the basic concepts of the gravitational microlensing technique, summarise the current mode of discovery and outline future steps towards a complete census of planets including Earth-mass planets. In the near-term (over the next 5 years) we advocate a strategy of automated follow-up with existing and upgraded telescopes which will significantly increase the current planet detection efficiency. In the medium 5-10 year term, we envision an international network of wide-field 2m class telescopes to discover Earth-mass and free-floating exo-planets. In the long (10-15 year) term, we strongly advocate a space microlensing telescope which, when combined with Kepler, will provide a complete census of planets down to Earth mass at almost all separations. Such a survey could be undertaken as a science programme on Euclid, a dark energy probe with a wide-field imager which has been proposed to ESA's Cosmic Vision Programme.

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.

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 present 838 ab-type RR Lyrae stars from the Lowell Observatory Near Earth Objects Survey Phase I (LONEOS-I). These objects cover 1430 deg^2 and span distances ranging from 3-30 kpc from the Galactic Center. Object selection is based on phased, photometric data with 28-50 epochs. We use this large sample to explore the bulk properties of the stellar halo, including the spatial distribution. The period-amplitude distribution of this sample shows that the majority of these RR Lyrae stars resemble Oosterhoff type I, but there is a significant fraction (26 %) which have longer periods and appear to be Oosterhoff type II. We find that the radial distributions of these two populations have significantly different profiles (rho_{OoI} ~ R^(-2.26 +- 0.07) and rho_{OoII} ~ R^(-2.88 +- 0.11). This suggests that the stellar halo was formed by at least two distinct accretion processes and supports dual-halo models.

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 mean reddenings toward 3525 RR0 Lyrae stars from the Galactic bulge fields of the MACHO Survey. These reddenings are determined using the color at minimum $V$-band light of the RR0 Lyrae stars themselves and are found to be in general agreement with extinction estimates at the same location obtained from other methods. Using 3256 stars located in the Galactic Bulge, we derive the selective extinction coefficient $R_{V,VR}=A_V/E(V-R) = 4.3 \pm 0.2$. This value is what is expected for a standard extinction law with $R_{V,BV} = 3.1 \pm 0.3$.

We present results from a study of the orbits of eclipsing binary stars (EBs) in the Magellanic Clouds. The samples comprise 4510 EBs found in the Large Magellanic Cloud (LMC) by the MACHO project, 2474 LMC EBs found by the OGLE-II project (of which 1182 are also in the MACHO sample), 1380 in the Small Magellanic Cloud (SMC) found by the MACHO project, and 1317 SMC EBs found by the OGLE-II project (of which 677 are also in the MACHO sample); we also consider the EROS sample of 79 EBs in the bar of the LMC. Statistics of the phase differences between primary and secondary minima allow us to infer the statistics of orbital eccentricities within these samples. We confirm the well-known absence of eccentric orbit in close binary stars. We also find evidence for rapid circularization in longer period systems when one member evolves beyond the main sequence, as also found by previous studies.

We present a new sample of 4634 eclipsing binary stars in the Large Magellanic Cloud (LMC), expanding on a previous sample of 611 objects and a new sample of 1509 eclipsing binary stars in the Small Magellanic Cloud (SMC), that were identified in the light curve database of the MACHO project. We perform a cross correlation with the OGLE-II LMC sample, finding 1236 matches. A cross correlation with the OGLE-II SMC sample finds 698 matches. We then compare the LMC subsamples corresponding to center and the periphery of the LMC and find only minor differences between the two populations. These samples are sufficiently large and complete that statistical studies of the binary star populations are possible.

This paper has been withdrawn by the author.

This paper has been withdrawn

A space-based gravitational microlensing exoplanet survey will provide a statistical census of exoplanets with masses down to 0.1 Earth-masses and orbital separations ranging from 0.5AU to infinity. This includes analogs to all the Solar System's planets except for Mercury, as well as most types of planets predicted by planet formation theories. Such a survey will provide results on the frequency of planets around all types of stars except those with short lifetimes. Close-in planets with separations < 0.5 AU are invisible to a space-based microlensing survey, but these can be found by Kepler. Other methods, including ground-based microlensing, cannot approach the comprehensive statistics on the mass and semi-major axis distribution of extrasolar planets that a space-based microlensing survey will provide. The terrestrial planet sensitivity of a ground-based microlensing survey is limited to the vicinity of the Einstein radius at 2-3 AU, and space-based imaging is needed to identify and determine the mass of the planetary host stars for the vast majority of planets discovered by microlensing. Thus, a space-based microlensing survey is likely to be the only way to gain a comprehensive understanding of the nature of planetary systems, which is needed to understand planet formation and habitability. The proposed Microlensing Planet Finder (MPF) mission is an example of a space-based microlensing survey that can accomplish these objectives with proven technology and a cost that fits comfortably under the NASA Discovery Program cost cap.

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.

We present an analysis of OGLE 2004-BLG-254, a high-magnification and relatively short duration microlensing event in which the source star, a Bulge K3-giant, has been spatially resolved by a point-like lens. We have obtained dense photometric coverage of the event light curve with OGLE and PLANET telescopes, as well as a high signal-to-noise ratio spectrum taken while the source was still magnified by 20, using the UVES/VLT spectrograph. Our dense coverage of this event allows us to measure limb darkening of the source star in the I and R bands. We also compare previous measurements of linear limb-darkening coefficients involving GK-giant stars with predictions from ATLAS atmosphere models. We discuss the case of K-giants and find a disagreement between limb-darkening measurements and model predictions, which may be caused by the inadequacy of the linear limb-darkening law.

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 the favoured core-accretion model of formation of planetary systems, solid planetesimals accumulate to build up planetary cores, which then accrete nebular gas if they are sufficiently massive. Around M-dwarf stars (the most common stars in our Galaxy), this model favours the formation of Earth-mass to Neptune-mass planets with orbital radii of 1 to 10 astronomical units (AU), which is consistent with the small number of gas giant planets known to orbit M-dwarf host stars. More than 170 extrasolar planets have been discovered with a wide range of masses and orbital periods, but planets of Neptune's mass or less have not hitherto been detected at separations of more than 0.15 AU from normal stars. Here we report the discovery of a 5.5 (+5.5/-2.7) M_earth planetary companion at a separation of 2.6 (+1.5/-0.6) AU from a 0.22 (+0.21/-0.11) M_solar M-dwarf star. (We propose to name it OGLE-2005-BLG-390Lb, indicating a planetary mass companion to the lens star of the microlensing event.) The mass is lower than that of GJ876d, although the error bars overlap. Our detection suggests that such cool, sub-Neptune-mass planets may be more common than gas giant planets, as predicted by the core accretion theory.

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.

We are re-analyzing the MACHO variable star database to explore the relationships between pulsation, evolution, and mass loss in evolved stars. We will analyze the multi-periodic properties of long period variable (LPV) stars, 50% of which could not be assigned any period in the original analysis. Recent results show that the missing stars may be an important element in understanding the origin of the period-luminosity sequences observed in the LMC. Our goal is to characterize the morphology and periodic properties of these stars, and then use the stars' 2MASS colors along with theoretical isochrones to understand their mass loss and evolution. We will develop a luminosity-independent criteria to classify the LMC LPVs, which can then be applied to Galactic LPVs. This will enable a synthesis of knowledge between LPVs in the LMC and the well-studied Galactic examples.

We have identified five new R Coronae Borealis (RCB) stars in the Galactic bulge using the MACHO Project photometry database, raising the total number of known Galactic RCB stars to about 40. We have obtained spectra to confirm the identifications. The fact that four out of the five newly identified RCB stars are ``cool'' (T(eff) < 6000 K) rather than ``warm'' (T(eff) > 6000 K) suggests that the preponderance of warm RCB stars among the existing sample is a selection bias. These cool RCB stars are redder and fainter than their warm counterparts and may have been missed in surveys done with blue plates. Based on the number of new RCB stars discovered in the MACHO bulge fields, there may be ~250 RCB stars in the reddened "exclusion" zone toward the bulge.

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

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

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

We report a measurement of the shape of the source star in microlensing event MOA 2002-BLG-33. The lens for this event was a close binary whose centre-of-mass passed almost directly in front of the source star. At this time, the source star was closely bounded on all sides by a caustic of the lens. This allowed the oblateness of the source star to be constrained. We found that a/b = 1.02^{+0.04}_{-0.02} where a and b are its semi-major and semi-minor axes respectively. The angular resolution of this measurement is approximately 0.04 microarcsec. We also report HST images of the event that confirm a previous identification of the source star as an F8-G2 turn-off main-sequence star.

We use the eight year light-curve database from the MACHO (MAssive Compact Halo Objects) project together with infrared colors and magnitudes from 2MASS (the Two Micron All Sky Survey) to identify a sample of 22,000 long period variables in the Large Magellanic Cloud (referred to hereafter as LMC LPVs). A period-luminosity diagram of these stars reveals six well-defined sequences, in substantial agreement with previous analyses of samples from OGLE (Optical Gravitational Lensing Experiment). In our analysis we identify analogues to galactic LPVs in the LMC LPV sample. We find that carbon-dominated AGB stars populate only two of the sequences, one of which includes the Mira variables. The high-luminosity end of the same two sequences are also the location of the only stars with $J-K_s > 2$, indicating that they are enshrouded in dust. The unknown mechanism that drives the variability of stars in the longest period produces different morphology in the period-luminosity diagram as compared to the shortest period sequences, which are thought to be caused by pulsation. In particular, the longest period sequence extends to lower luminosity RGB stars and the luminosity function does not peak among the AGB stars. We point out several features which will constrain new models of the period-luminosity sequences.

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 present the radial velocities, metallicities and the K-band magnitudes of 74 RR Lyrae stars in the inner regions of the LMC. The intermediate resolution spectra and infrared images were obtained with FORS1 at the ESO VLT and with the SOFI infrared imager at the ESO NTT. The best 43 RR Lyrae with measured velocities yield an observed velocity dispersion of 61+-7 km/s. We obtain a true LMC RR Lyrae velocity dispersion of 53 km/s, which is higher than the velocity dispersion of any other LMC population previously measured. This is the first empirical evidence for a kinematically hot, metal-poor halo in the LMC as discussed in Minniti et al. (2003). Using Layden's (1994) modification of the Delta S method we measured the metallicity for 23 of our stars. The mean value is [Fe/H] =-1.46+-0.09 dex. The absolute magnitudes M_V and M_K of RR Lyrae stars are linear functions of metallicity. In the V band, our data agree with the Olech et al. (2003) relation, in the K band the slope is flatter. The average apparent V luminosity of 70 RR Lyrae stars is <V> = 19.45+-0.04 and the average K luminosity of 37 RR Lyrae stars is <K>=18.20+-0.06. There is no obvious relation between apparent V magnitude and LogP, while the RR Lyrae K band magnitudes show a well defined linear trend with LogP. Using the Bono et al. (2001) and Bono et al. (2003) theoretical Near-Infrared Period-Luminosity-Metallicity relations we calculate the LMC distance modulus 18.48+-0.08.

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

We report a measurement of limb darkening of a solar-like star in the very high magnification microlensing event MOA 2002-BLG-33. A 15 hour deviation from the light curve profile expected for a single lens was monitored intensively in V and I passbands by five telescopes spanning the globe. Our modelling of the light curve showed the lens to be a close binary system whose centre-of-mass passed almost directly in front of the source star. The source star was identified as an F8-G2 main sequence turn-off star. The measured stellar profiles agree with current stellar atmosphere theory to within ~4% in two passbands. The effective angular resolution of the measurements is <1 micro-arcsec. These are the first limb darkening measurements obtained by microlensing for a Solar-like star.

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.

The oldest and most metal-poor Milky Way stars form a kinematically hot halo, which motivates the two major formation scenarios for our galaxy: extended hierarchical accretion and rapid collapse. RR Lyrae stars are excellent tracers of old and metal-poor populations. We measure the kinematics of 43 RR Lyrae stars in the inner regions of the nearby Large Magellanic Cloud (LMC) galaxy. The velocity dispersion, \sigma_{true}=53\pm10 km/s, indicates that a kinematically hot metal-poor old halo also exists in the LMC. This suggests that our galaxy and smaller late-type galaxies like the LMC have similar early formation histories.

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.

Current practice in stellar evolution is to employ one-dimensional calculations that quantitatively apply only to a minority of the observed stars (single non-rotating stars, or well detached binaries). Even in these systems, astrophysicists are dependent on approximations to handle complex three-dimensional processes like convection. Understanding the structure of binary stars, like those that lead to the Type Ia supernovae used to measure the expansion of the universe, are grossly non-spherical and await a 3D treatment. To approach very large problems like multi-dimensional modeling of stars, the Lawrence Livermore National Laboratory has invested in massively parallel computers and invested even more in developing the algorithms to utilize them on complex physics problems. We have leveraged skills from across the lab to develop a 3D stellar evolution code, Djehuty (after the Egyptian god for writing and calculation) that operates efficiently on platforms with thousands of nodes, with the best available physical data (opacities, EOS, etc.). Djehuty has incorporated all basic physics for modeling stars including an accurate equation of state, radiation transport by diffusion, thermonuclear reaction rates, and hydrodynamics, and we have begun testing it in a number of applications.

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

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

The Galactic Exoplanet Survey Telescope (GEST) will observe a 2 square degree field in the Galactic bulge to search for extra-solar planets using a gravitational lensing technique. This gravitational lensing technique is the only method employing currently available technology that can detect Earth-mass planets at high signal-to-noise, and can measure the frequency of terrestrial planets as a function of Galactic position. GEST's sensitivity extends down to the mass of Mars, and it can detect hundreds of terrestrial planets with semi-major axes ranging from 0.7 AU to infinity. GEST will be the first truly comprehensive survey of the Galaxy for planets like those in our own Solar System.

We present an abundance analysis of six main sequence turnoff, subgiant, and giant branch stars toward the Galactic bulge that were observed with Keck/HIRES during microlensing events. This is an early look at the first detailed chemical analysis of main sequence stars in the Galactic bulge. Lensing events allow the effective aperture of Keck to be increased beyond its current dimensions; although, some events still stretched its spectroscopic capabilities. Future large telescopes with high resolution and high throughput spectrometers will allow the study of abundances in distant stellar populations and in less evolved stars with greater ease.

We present the photometric properties of 1279 blue variable stars within the LMC. Photometry is derived from the MACHO database. The lightcurves of the sample exhibit a variety of quasi-periodic and aperiodic outburst behavior. A characteristic feature of the photometric variation is that the objects are reddest when at maximum outburst. A subset of 102 objects were examined spectroscopically. Within this subset, 91% exhibited Balmer emission in at least one epoch, in some cases with spectacular spectral variability. The variability observed in the sample is consistent with the establishment and maintenance of the Be phenomenon.

The Hipparcos I-band calibration of horizontal-branch red clump giants as standard candles has lead to controversial results for the distance to the Large Magellanic Cloud (LMC). In an attempt to properly ascertain the corrections for interstellar extinction and clump age and metallicity, we analyze new multi-wavelength luminosity functions of the LMC red clump. Our photometry dataset in the K-band was obtained with the SOFI infrared imager at the European Southern Observatory's New Technology Telescope. In the V and I passbands, we employ data from WFPC2 onboard the Hubble Space Telescope. The LMC red clump is first identified in a K,(V-K) color-magnitude diagram. Our luminosity functions yield apparent magnitudes of K = 16.974, I = 18.206, and V = 19.233 (+- 0.009_r +- 0.02_s; random and systematic error, respectively). Compared directly to the Hipparcos red clump calibration (without a correction for age and metallicity), the LMC clump measurements imply a negative interstellar reddening correction. This unphysical result indicates a population difference between clumps. A modified calibration based on theoretical modeling yields an average reddening correction of E(B-V) = 0.089 +- 0.015_r, and a true LMC distance modulus of 18.493 +- 0.033_r +- 0.03_s. We reconcile our result with the short distance previously derived from OGLE II red clump data.

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 the first astronomical observations of the Livermore Imaging Fourier Transform Spectrograph for visible-band astronomy using the 3.5-meter Apache Point Observatory.

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 have performed a search for halo white dwarfs as high proper motion objects in a second epoch WFPC2 image of the Groth-Westphal strip. We identify 24 high proper motion objects with mu > 0.014 ''/yr. Five of these high proper motion objects are identified as strong white dwarf candidates on the basis of their position in a reduced proper motion diagram. We create a model of the Milky Way thin disk, thick disk and stellar halo and find that this sample of white dwarfs is clearly an excess above the < 2 detections expected from these known stellar populations. The origin of the excess signal is less clear. Possibly, the excess cannot be explained without invoking a fourth galactic component: a white dwarf dark halo. We present a statistical separation of our sample into the four components and estimate the corresponding local white dwarf densities using only the directly observable variables, V, V-I, and mu. For all Galactic models explored, our sample separates into about 3 disk white dwarfs and 2 halo white dwarfs. However, the further subdivision into the thin and thick disk and the stellar and dark halo, and the subsequent calculation of the local densities are sensitive to the input parameters of our model for each Galactic component. Using the lowest mean mass model for the dark halo we find a 7% white dwarf halo and six times the canonical value for the thin disk white dwarf density (at marginal statistical significance), but possible systematic errors due to uncertainty in the model parameters likely dominate these statistical error bars. The white dwarf halo can be reduced to around 1.5% of the halo dark matter by changing the initial mass function slightly. The local thin disk white dwarf density in our solution can be made consistent with the canonical value by assuming a larger thin disk scaleheight of 500 pc.

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.

The DJEHUTY project is an intensive effort at the Lawrence Livermore National Laboratory (LLNL) to produce a general purpose 3-D stellar structure and evolution code to study dynamic processes in whole stars.

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

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

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

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

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

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

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

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

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

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

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

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

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

We present UBV photometry of the eclipsing binary Harvard Variable 2274 in the Large Magellanic Cloud (LMC). The stellar parameters of this binary system were calculated by Guinan et al. (1998a) who gave both a reddening towards HV 2274 of E(B-V) = 0.083 +/- 0.006 and a distance modulus to the LMC, mu_{LMC} = 18.42 +/- 0.007. The reddening of this system was also determined by Udalski et al. (1998) who found E(B-V) = 0.149 +\- 0.015. With Udalski et al. (1998) B and V photometry Guinan et al.(1998b) obtained E(B-V) = 0.120 +\- 0.009 and mu_{LMC} = 18.30 +\- 0.07. Using our UBV photometry, we derive a reddening of E(B-V)= 0.088 +\- 0.025, consistent with the original value of Guinan et al. (1998a) and supporting a longer distance modulus to the LMC of about mu_{LMC} = 18.40 +\- 0.07. We stress the uncertainties inherent in ground-based UBV photometry and the concomitant uncertainties in determining distances based upon such photometry.

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 High Resolution Echelle Spectrometer on the Keck I telescope has been used to measure nine radial velocities having a median precision of 2 km/s for seven red giants belonging to the isolated dwarf spheroidal galaxy Andromeda II. We find a weighted mean radial velocity of -188+/-3 km/s and a central velocity dispersion of 9.3+/-2.7 km/s. There may be evidence for a radial velocity gradient across the face of the galaxy, although the significance of this result is low due to the small number of stars having measured velocities. Our best estimate for the global mass-to-light ratio of And II is M/L_V = 20.9^{+13.9}_{-10.1} in solar units. This value is similar to those of several Galactic dwarf spheroidal galaxies and is consistent with the presence of a massive dark halo in Andromeda II.

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 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 describe a concept for an imaging spectrograph for a large orbiting observatory such as NASA's proposed Next Generation Space Telescope (NGST) based on an imaging Fourier transform spectrograph (IFTS). An IFTS has several important advantages which make it an ideal instrument to pursue the scientific objectives of NGST. We review the operation of an IFTS and make a quantitative evaluation of the signal-to-noise performance of such an instrument in the context of NGST. We consider the relationship between pixel size, spectral resolution, and diameter of the beamsplitter for imaging and non-imaging Fourier transform spectrographs and give the condition required to maintain spectral modulation efficiency over the entire field of view. We give examples of scientific programs that could be performed with this facility.

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.

The bulge contains the most chemically evolved old stellar population in the Milky Way. Thanks to microlensing, it is now possible to obtain high resolution echelle spectra of bulge stars near the main sequence turn-off, and study the abundance of elements that are affected by stellar evolution, such as lithium. We have observed with the HIRES spectrograph on the Keck I 10m telescope the source star of the MACHO microlensing event 97BLG45 while it was magnified by about 1 magnitude. Here we report the detection of the Li I line at \lambda 6707.8 \AA in the echelle spectrum of this star, and measure the bulge lithium abundance for the first time: A(Li) = 2.25 \pm 0.25.

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.

We report the detection of 45 candidate microlensing events in fields toward the Galactic bulge. These come from the analysis of 24 fields containing 12.6 million stars observed for 190 days in 1993. Many of these events are of extremely high signal to noise and are remarkable examples of gravitational microlensing. The distribution of peak magnifications is shown to be consistent with the microlensing interpretation of these events. Using a sub-sample of 1.3 million ``Clump Giant" stars whose distance and detection efficiency are well known, we find 13 events and estimate the microlensing optical depth toward the Galactic Bulge as $\tau_{\rm bulge} = 3.9 {+ 1.8 \atop - 1.2} \times 10^{-6}$ averaged over an area of $\sim 12$ square degrees centered at Galactic coordinates $\ell = 2.55^\circ$ and $b = -3.64^\circ$. This is similar to the value reported by the OGLE collaboration, and is marginally higher than current theoretical models for $\tau_{\rm bulge}$. The optical depth is also seen to increase significantly for decreasing $\vert b\vert$. These results demonstrate that obtaining large numbers of microlensing events toward the Galactic bulge is feasible, and that the study of such events will have important consequences for the structure of the Galaxy and its dark halo.

We provide a status report on our search for dark matter in our Galaxy in the form of massive compact halo objects (or Machos), using gravitational microlensing of background stars. This search uses a very large format CCD camera on the dedicated 1.27m telescope at Mt.~Stromlo, Australia, and has been taking data for almost 3 years. At present, we are in the midst of analyzing our second year data for 8 million stars in the Large Magellanic Cloud. We find more microlensing events than expected from known stellar populations suggesting that Machos are indeed present in the Galactic halo, but the observed microlensing rate toward the LMC is too small to allow for a halo dominated by sub-stellar Machos. Our observations of the Galactic bulge have also yielded substantially more microlensing events than anticipated including a number of exotic ``deviant" microlensing events. The implications of these results are discussed.

The MACHO experiment is searching for dark matter in the halo of the Galaxy by monitoring more than 20 million stars in the LMC and Galactic bulge for gravitational microlensing events. The hardware consists of a 50 inch telescope, a two-color 32 megapixel ccd camera, and a network of computers. On clear nights the system generates up to 8 GB of raw data and 1 GB of reduced data. The computer system is responsible for all realtime control tasks, for data reduction, and for storing all data associated with each observation in a data base. The subject of this paper is the software system that handles these functions. It is an integrated system controlled by Petri nets that consists of multiple processes communicating via mailboxes and a bulletin board. The system is highly automated, readily extensible, and incorporates flexible error recovery capabilities. It is implemented with C++ in a Unix environment.

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.

We present the first detection of parallax effects in a gravitational microlensing event. Parallax in a gravitational microlensing event observed only from the Earth appears as a distortion of the lightcurve due to the motion of the Earth around the Sun. This distortion can be detected if the event duration is not much less than a year and if the projected velocity of the lens is not much larger than the orbital velocity of the Earth about the Sun. The event presented here has a duration of 220 days and clearly shows the distortion due to the Earth's motion. We find that the projected velocity of the lens is 75+/-5 km/s at an angle of 28+/-4 deg from the direction of increasing galactic longitude, as expected for a lens in the galactic disk. A likelihood analysis yields estimates of the distance to and mass of the lens: D_{lens} = 1.7 (+1.1/-0.7) kpc and M = 1.3 (+1.3/-0.6) Msun, suggesting that the lens is a remnant such as a white dwarf or neutron star. A less likely possibility is that the lens is a main sequence star. If so, we can add our upper limit on the observed flux from the lens to the analysis. This modifies the estimates to: D_{lens} = 2.8 (+1.1/-0.6) kpc, and M = 0.6 (+0.4/-0.2) Msun.

The MACHO collaboration reports on the analysis of our first year LMC data, 9.5 million light curves with an average of 235 observations each. Automated selection procedures give 3 events consistent with microlensing. We evaluate our experimental detection efficiency using a range of Monte- Carlo simulations. Using a `standard' halo density profile we find that a halo comprised entirely of Machos in the mass range 3 \ten{-4} to 0.06 \msun would predict > 15 detected events in this dataset; thus a standard spherical halo cannot be dominated by objects in this mass range. Assuming all three events are microlensing of halo objects and fitting a naive spherical halo model to our data yields a Macho halo fraction f =0.19+0.16-0.10, a total mass in Machos (inside 50 kpc) of 7.6+6-4 \ten{10} \msun, and a microlensing optical depth 8.8+7-5 \ten{-8} (68\% CL). Exploring a wide range of halo models we find that our constraints on the Macho fraction are quite model-dependent, but constraints on the total mass in Machos within 50 kpc are quite secure.

Results are presented from the MACHO collaboration gravitational microlensing search. The experiment and the nearly 50 microlensing events that have been detected are described. Limits on the baryonic content of the halo are given, as are estimates of the Macho contribution to the dark halo. Optical depths toward the bulge, and several unusual events such as a binary lens and a parallax event are discussed. Possible interpretations of these results are also discussed.

The MACHO Collaboration's search for baryonic dark matter via its gravitational microlensing signature has generated a massive database of time ordered photometry of millions of stars in the LMC and the bulge of the Milky Way. The search's experimental design and capabilities are reviewed and the dark matter results are briefly noted. Preliminary analysis of the approximately 39,000 variable stars discovered in the LMC database is presented and examples of periodic variables are shown. A class of aperiodically variable Be stars is described which is the closest background to microlensing which has been found. Plans for future work on variable stars using the MACHO data are described.

We have identified 4 Cepheids and 5 Miras using KPNO 4m BVRI images of an outer field in M101. The Cepheid and Mira periods range from 30 to 60 days and 350 to 800 days, respectively. We derive independent Cepheid and Mira distance moduli that agree within experimental uncertainties. We find a true distance modulus of 29.08 +- 0.13 mag.

In the past few years, the Magellanic Clouds have been the targets for several major variable star surveys. The results of these surveys are now becoming available and it is clear that a Renaissance in LMC and SMC variable star research will result. In this review, I will describe the results of such surveys and review the questions that are likely to be answered by further work. With respect to results, I will concentrate on LMC MACHO Project data, including beat Cepheids, discovery statistics, mode identification, Fourier decomposition of lightcurves, and the differences between the LMC and galactic sample.

We have monitored 8.6 million stars in the Large Magellanic Cloud for 1.1 years and have found 3 events consistent with gravitational microlensing. We place strong constraints on the Galactic halo content in the form of compact lensing objects in the mass range $10^{-4} \msun$ to $10^{-1} \msun$. Three events is fewer than expected for a standard spherical halo of objects in this mass range, but appears to exceed the number expected from known Galactic populations. Fitting a naive spherical halo model to our data yields a MACHO fraction $f = 0.20^{+0.33}_{-0.14}$, which implies a total MACHO mass (inside 50 kpc) of $8.0^{+14}_{-6}\ten{10} \msun$, and a microlensing optical depth $9^{+15}_{-7} \ten{-8}$ ($\sim 68$\% CL).

We review recent gravitational microlensing results from the EROS, MACHO, and OGLE collaborations, and present some details of the very latest MACHO results toward the Galactic Bulge. The MACHO collaboration has now discovered in excess of 40 microlensing events toward the Galactic Bulge during the 1993 observing season. A preliminary analysis of this data suggests a much higher microlensing optical depth than predicted by standard galactic models suggesting that these models will have to be revised. This may have important implications for the structure of the Galaxy and its dark halo. Also shown are MACHO data of the first microlensing event ever detected substantially before peak amplification, the first detection of parallax effects in a microlensing event, and the first caustic crossing to be resolved in a microlensing event.

We report the discovery of 45 beat Cepheids in the Large Magellanic Cloud (LMC) using the MACHO Project photometry database. The variables which are pulsating simultaneously in two radial modes are shown to break cleanly into two period-ratio groups, providing the first unambiguous evidence that the second overtone is indeed excited in real Cepheids. Thirty stars are beating in the fundamental and first overtone mode (F/1H, with a period ratio in the neighborhood of 0.72), and fifteen stars are beating in the first and second overtone (1H/2H, with a period ratio near 0.80). The F/1H period ratios are systematically higher than known Galactic beat Cepheids, indicating a metallicity dependence whose sense is in agreement with theory. Beat Cepheids in the LMC are found to select the 1H/2H mode for fundamental periods shorter than 1.25 days. We find the fraction of Cepheids excited in two modes to be about 20\% for stars with fundamental periods shorter than 2.5 days. We fail to confirm any of the proposed beat Cepheid candidates common to our sample from the surveys of Andreasen (1987) and Andreasen \& Petersen (1987). We also present finder charts and find several of the beat Cepheids to be in or near LMC clusters.

The detection of microlensing has opened the way for the development of new methods in galactic astronomy. This series of papers investigates what microlensing can teach us about the structure and shape of the dark halo. In this paper we present formulas for the microlensing rate, optical depth and event duration distributions for a simple set of axisymmetric disk-halo models. The halos are based on the "power--law models" which have simple velocity distributions. Using these models, we show that there is a large uncertainty in the predicted microlensing rate because of uncertainty in the halo parameters. For example, models which reproduce the measured galactic observables to within their errors still differ in microlensing rate towards the Magellanic Clouds by more than a factor of ten. We find that while the more easily computed optical depth correlates well with microlensing rate, the ratio of optical depth to rate can vary by a factor of two (or greater if the disk is maximal). Comparison of microlensing rates towards the Large and Small Magellanic Clouds (LMC and SMC) and M31 can be used to aid determinations of the halo flattening and rotation curve slope. For example, the ratio of microlensing rates towards the LMC and SMC is $\sim 0.7-0.8$ for E0 halos and $\sim 1.0 - 1.2$ for E7 halos (c.f. Sackett \& Gould 1993). Once the flattening has been established, the ratio of microlensing rates towards M31 and the LMC may help to distinguish between models with rising, flat or falling rotation curves. Comparison of rates along LMC and galactic bulge lines-of-sight gives useful information on the halo core radius, although this may not be so easy to extract in practice. Maximal disk models provide substantially smaller halo optical depths, shorter event durations and even larger model uncertainties.

We provide a status report on our search for dark matter in our Galaxy in the form of massive compact halo objects (MACHOs), using gravitational microlensing of background stars. This search uses a very large CCD camera on the dedicated 1.27m telescope at Mt.~Stromlo, Australia, and has been taking data for 2 years. At present, we have analysed data for 8 million stars in the Large Magellanic Cloud over 1 year, resulting in one strong candidate event and two lower-amplitude candidates. We have also analysed 5 million stars in the Galactic Bulge for 0.5 years, yielding approximately 27 microlensing events.

The MACHO project carries out regular photometric monitoring of millions of stars in the Magellanic Clouds and Galactic Bulge, to search for very rare gravitational microlensing events due to compact objects in the galactic halo and disk. A preliminary analysis of one field in the Galactic Bulge, containing {$\sim430,000$} stars observed for 190 days, reveals four stars which show clear evidence for brightenings which are time-symmetric, achromatic in our two passbands, and have shapes consistent with gravitational microlensing. This is significantly higher than the $\sim 1$ event expected from microlensing by known stars in the disk. If all four events are due to microlensing, a 95\% confidence lower limit on the optical depth towards our bulge field is $1.3 \times 10^{-6}$, and a ``best fit" value is $\tau \approx 1.6 \times 10^{-6}/\epsilon$,where $\epsilon$ is the detection efficiency of the experiment, and $\epsilon < 0.4$. If the true optical depth is close to the ``best fit" value, possible explanations include a ``maximal" disk which accounts for most of the galactic circular velocity at the solar radius, a halo which is centrally concentrated, or bulge-bulge microlensing.

There is now abundant evidence for the presence of large quantities of unseen matter surrounding normal galaxies, including our own$^{1,2}$. The nature of this `dark matter' is unknown, except that it cannot be made of normal stars, dust, or gas, as they would be easily detected. Exotic particles such as axions, massive neutrinos or other weakly interacting massive particles (collectively known as WIMPs) have been proposed)$^{3,4}$, but have yet to be detected. A less exotic alternative is normal matter in the form of bodies with masses ranging from that of a large planet to a few $ \msun$. Such objects, known collectively as massive compact halo objects$^5$ (MACHOs) might be brown dwarfs or `Jupiters' (bodies too small to produce their own energy by fusion), neutron stars, old white dwarfs, or black holes. Paczynski$^6$ suggested that MACHOs might act as gravitational microlenses, occasionally causing the apparent brightness of distant background stars temporarily to increase. We are conducting a microlensing experiment to determine whether the dark matter halo of our galaxy is made up of MACHOs. Here we report a candidate for a microlensing event, detected by monitoring the light curves of 1.8 million stars in the Large Magellanic Cloud for one year. The light curve shows no variation for most of the year of data taking, and an upward excursion lasting over 1 month, with a maximum increase of $\approx \bf 2$ mag. The most probable lens mass, inferred from the duration of the event, is $\bf \sim 0.1 \,\msun$.