Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $\sigma_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $\alpha=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history.

We explore the long-term evolution of a bias-free orbital representation of the cometary nuclei (with diameters above 2 km) of the Kuiper belt, using the so-called L7 synthetic model from CFEPS, which consists of three dynamical sub-populations: the Classical, the Resonant, and the Scattering. The dynamical evolution of belt particles is studied under the gravitational influence of the Sun and the four giant planets, as well as of the 34 largest known trans-Neptunian objects (TNOs with $H_V < 4$). Here we indistinctly call Dwarf Planets (DPs) to the full sample of 34 large TNOs. Over a 1 Gyr timescale, we analyze the secular influence of the DPs over Kuiper belt disk particles and their contribution to the injection rate of new visible Jupiter Family Comets (JFCs). We found that DPs globally increase the number of JFCs by 12.6% when compared with the comets produced by the giant planets alone. When considering each population separately, we find that the increment produced by DPs is 17%, 12%, and 3% for the Classical, Resonant, and Scattering populations, respectively. Given the rate of escapes from the Kuiper belt, we find upper limits to the number of objects in each population required to maintain the JFCs in steady-state; the results are $55.9\times10^6$, $78.5\times10^6$, and $274.3\times10^6$ for the Scattering, Resonant, and Classical populations, respectively. Finally, we find that the Plutinos are the most important source of comets which were originally in a resonant configuration, where the presence of Pluto alone enhances by 10% the number of JFCs.

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) is acquiring near-simultaneous $g$, $r$, and $J$ photometry of unprecedented precision with the Gemini North Telescope, targeting nearly a hundred trans-Neptunian objects (TNOs) brighter than $m_r=23.6$ mag discovered in the Outer Solar System Origins Survey. Combining the optical and near-infrared photometry with the well-characterized detection efficiency of the Col-OSSOS target sample will provide the first flux-limited compositional dynamical map of the outer Solar System. In this paper, we describe our observing strategy and detail the data reduction processes we employ, including techniques to mitigate the impact of rotational variability. We present optical and near-infrared colors for 35 TNOs. We find two taxonomic groups for the dynamically excited TNOs, the neutral and red classes, which divide at $g-r \simeq 0.75$. Based on simple albedo and orbital distribution assumptions, we find that the neutral class outnumbers the red class, with a ratio of 4:1 and potentially as high as 11:1. Including in our analysis constraints from the cold classical objects, which are known to exhibit unique albedos and $r-z$ colors, we find that within our measurement uncertainty, our observations are consistent with the primordial Solar System protoplanetesimal disk being neutral-class-dominated, with two major compositional divisions in $grJ$ color space.

We present variability measurements and partial light curves of Trans-Neptunian Objects (TNOs) from a two-night pilot study using Hyper Suprime-Cam (HSC) on the Subaru Telescope (Maunakea, Hawai'i, USA). Subaru's large aperture (8-m) and HSC's large field of view (1.77 square degrees) allow us to obtain measurements of multiple objects with a range of magnitudes in each telescope pointing. We observed 65 objects with m_r = 22.6--25.5 mag in just six pointings, allowing 20--24 visits of each pointing over the two nights. Our sample, all discovered in the recent Outer Solar System Origins Survey (OSSOS), span absolute magnitudes H_r = 6.2--10.8 mag and thus investigates smaller objects than previous light curve projects have typically studied. Our data supports the existence of a correlation between light curve amplitude and absolute magnitude seen in other works, but does not support a correlation between amplitude and orbital inclination. Our sample includes a number of objects from different dynamical populations within the trans-Neptunian region, but we do not find any relationship between variability and dynamical class. We were only able to estimate periods for 12 objects in the sample and found that a longer baseline of observations is required for reliable period analysis. We find that 31 objects (just under half of our sample) have variability greater than 0.4 magnitudes during all of the observations; in smaller 1.25 hr, 1.85 hr and 2.45 hr windows, the median variability is 0.13, 0.16 and 0.19 mags, respectively. The fact that variability on this scale is common for small TNOs has important implications for discovery surveys (such as OSSOS or the Large Synoptic Survey Telescope) and color measurements.

Both physical and dynamical properties must be considered to constrain the origins of the dynamically excited distant Solar System populations. We present high-precision (g-r) colors for 25 small (Hr>5) dynamically excited Trans-Neptunian Objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey (Col-OSSOS). We combine our dataset with previously published measurements and consider a set of 229 colors of outer Solar System objects on dynamically excited orbits. The overall color distribution is bimodal and can be decomposed into two distinct classes, termed `gray' and `red', that each has a normal color distribution. The two color classes have different inclination distributions: red objects have lower inclinations than the gray ones. This trend holds for all dynamically excited TNO populations. Even in the worst-case scenario, biases in the discovery surveys cannot account for this trend: it is intrinsic to the TNO population. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs, our finding solves the conundrum of centaurs being the only outer Solar System population identified so far to exhibit this property (Tegler et al. 2016). The different orbital distributions of the gray and red dynamically excited TNOs provide strong evidence that their colors are due to different formation locations in a disk of planetesimals with a compositional gradient.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

The recent discovery by Pan-STARRS1 of 1I/2017 U1 (`Oumuamua), on an unbound and hyperbolic orbit, offers a rare opportunity to explore the planetary formation processes of other stars, and the effect of the interstellar environment on a planetesimal surface. 1I/`Oumuamua's close encounter with the inner Solar System in 2017 October was a unique chance to make observations matching those used to characterize the small-body populations of our own Solar System. We present near-simultaneous g$^\prime$, r$^\prime$, and J photometry and colors of 1I/`Oumuamua from the 8.1-m Frederick C. Gillett Gemini North Telescope, and $gri$ photometry from the 4.2 m William Herschel Telescope. Our g$^\prime$r$^\prime$J observations are directly comparable to those from the high-precision Colours of the Outer Solar System Origins Survey (Col-OSSOS), which offer unique diagnostic information for distinguishing between outer Solar System surfaces. The J-band data also provide the highest signal-to-noise measurements made of 1I/`Oumuamua in the near-infrared. Substantial, correlated near-infrared and optical variability is present, with the same trend in both near-infrared and optical. Our observations are consistent with 1I/`Oumuamua rotating with a double-peaked period of $8.10 \pm 0.42$ hours and being a highly elongated body with an axial ratio of at least 5.3:1, implying that it has significant internal cohesion. The color of the first interstellar planetesimal is at the neutral end of the range of Solar System $g-r$ and $r-J$ solar-reflectance colors: it is like that of some dynamically excited objects in the Kuiper belt and the less-red Jupiter Trojans.

Several different classes of trans-Neptunian objects (TNOs) have been identified based on their optical and near-infrared colors. As part of the Colours of the Outer Solar System Origins Survey, we have obtained $g$, $r$, and $z$ band photometry of 26 TNOs using Subaru and Gemini Observatories. Previous color surveys have not utilized $z$ band reflectance, and the inclusion of this band reveals significant surface reflectance variations between sub-populations. The colors of TNOs in $g-r$ and $r-z$ show obvious structure, and appear consistent with the previously measured bi-modality in $g-r$. The distribution of colors of the two dynamically excited surface types can be modeled using the two-component mixing models from Fraser \& Brown (2012). With the combination of $g-r$ and $r-z$, the dynamically excited classes can be separated cleanly into red and neutral surface classes. In $g - r$ and $r - z$, the two dynamically excited surface groups are also clearly distinct from the cold classical TNO surfaces, which are red, with $g-r\gtrsim$0.85 and $r-z\lesssim$0.6, while all dynamically excited objects with similar $g-r$ colors exhibit redder $r-z$ colors. The $z$ band photometry makes it possible for the first time to differentiate the red excited TNO surfaces from the red cold classical TNO surfaces. The discovery of different $r-z$ colors for these cold classical TNOs makes it possible to search for cold classical surfaces in other regions of the Kuiper belt and to completely separate cold classical TNOs from the dynamically excited population, which overlaps in orbital parameter space.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is currently the deepest wide- field survey in progress. The 8.2 m aperture of Subaru telescope is very powerful in detect- ing faint/small moving objects, including near-Earth objects, asteroids, centaurs and Tran- Neptunian objects (TNOs). However, the cadence and dithering pattern of the HSC-SSP are not designed for detecting moving objects, making it difficult to do so systematically. In this paper, we introduce a new pipeline for detecting moving objects (specifically TNOs) in a non-dedicated survey. The HSC-SSP catalogs are re-arranged into the HEALPix architecture. Then, the stationary detections and false positive are removed with a machine learning al- gorithm to produce a list of moving object candidates. An orbit linking algorithm and visual inspections are executed to generate the final list of detected TNOs. The preliminary results of a search for TNOs using this new pipeline on data from the first HSC-SSP data release (Mar 2014 to Nov 2015) are also presented.

The cold classical Kuiper belt objects have low inclinations and eccentricities and are the only Kuiper belt population suspected to have formed in situ. Compared with the dynamically excited populations, which exhibit a broad range of colours and a low binary fraction of ~10% cold classical Kuiper belt objects typically have red optical colours with ~30% of the population found in binary pairs; the origin of these differences remains unclear. We report the detection of a population of blue-coloured, tenuously bound binaries residing among the cold classical Kuiper belt objects. Here we show that widely separated binaries could have survived push-out into the cold classical region during the early phases of Neptune's migration. The blue binaries may be contaminants, originating at ~38 au, and could provide a unique probe of the formative conditions in a region now nearly devoid of objects. The idea that the blue objects, which are predominantly binary, are the products of push-out requires that the planetesimals formed entirely as multiples. Plausible formation routes include planetesimal formation via pebble accretion and subsequent binary production through dynamic friction and binary formation during the collapse of a cloud of solids.

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

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.

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

We suggest a new web-based approach for browsing and visualizing data produced by a network of telescopes, such as those of the ongoing TAOS and the forthcoming TAOS II projects. We propose a modern client-side technology and we present two examples based on two software packages developed for different kinds of server- side database approaches. In spite our examples are specific for the browsing of TAOS light curves, the software is coded in a way to be suitable for the use in several types of astronomical projects.

A preliminary data analysis of the stellar light curves obtained by the robotic telescopes of the TAOS project is presented. We selected a data run relative to one of the stellar fields observed by three of the four TAOS telescopes, and we investigate the common trend and the correlation between the light curves. We propose two ways to remove these trends and show the preliminary results. A project aimed at flagging interesting behaviors, such as stellar variability, and to set up an automated follow-up with the San Pedro M\'artir Facilities is on the way.

We present the analysis on our new limits of the dark matter (DM) halo consisting of primordial black holes (PBHs) or massive compact halo objects (MACHOs). We present a search of the first two years of publicly available Kepler mission data for potential signatures of gravitational microlensing caused by these objects, as well as an extensive analysis of the astrophysical sources of background error. These include variable stars, flare events, and comets or asteroids which are moving through the Kepler field. We discuss the potential of detecting comets using the Kepler lightcurves, presenting measurements of two known comets and one unidentified object, most likely an asteroid or comet. After removing the background events with statistical cuts, we find no microlensing candidates. We therefore present our Monte Carlo efficiency calculation in order to constrain the PBH DM with masses in the range of 2 x 10^-9 solar masses to 10^-7 solar masses. We find that PBHs in this mass range cannot make up the entirety of the DM, thus closing a full order of magnitude in the allowed mass range for PBH DM.

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.

We present the study on GRB 071112C X-ray and optical light curves. In these two wavelength ranges, we have found different temporal properties. The R-band light curve showed an initial rise followed by a single power-law decay, while the X-ray light curve was described by a single power-law decay plus a flare-like feature. Our analysis shows that the observed temporal evolution cannot be described by the external shock model in which the X-ray and optical emission are produced by the same emission mechanism. No significant color changes in multi-band light curves and a reasonable value of the initial Lorentz factor ({\Gamma}0 = 275 \pm 20) in a uniform ISM support the afterglow onset scenario as the correct interpretation for the early R-band rise. The result suggests the optical flux is dominated by afterglow. Our further investigations show that the X-ray flux could be created by an additional feature related to energy injection and X-ray afterglow. Different theoretical interpretations indicate the additional feature in X-ray can be explained by either late internal dissipation or local inverse-Compton scattering in the external shock.

If the Dark Matter consists of primordial black holes (PBHs), we show that gravitational lensing of stars being monitored by NASA's Kepler search for extra-solar planets can cause significant numbers of detectable microlensing events. A search through the roughly 150,000 lightcurves would result in large numbers of detectable events for PBHs in the mass range $5 \ten{-10}\msun$ to $\aten{-4}\msun$. Non-detection of these events would close almost two orders of magnitude of the mass window for PBH dark matter. The microlensing rate is higher than previously noticed due to a combination of the exceptional photometric precision of the Kepler mission and the increase in cross section due to the large angular sizes of the relatively nearby Kepler field stars. We also present a new formalism for calculating optical depth and microlensing rates in the presence of large finite-source effects.

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

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 conducted a search for occultations of bright stars by Kuiper Belt Objects (KBOs) to estimate the density of sub-km KBOs in the sky. We report here the first results of this occultation survey of the outer solar system conducted in June 2007 and June/July 2008 at the MMT Observatory using Megacam, the large MMT optical imager. We used Megacam in a novel shutterless continuous--readout mode to achieve high precision photometry at 200 Hz. We present an analysis of 220 star hours at signal-to-noise ratio of 25 or greater. The survey efficiency is greater than 10% for occultations by KBOs of diameter d>=0.7 km, and we report no detections in our dataset. We set a new 95% confidence level upper limit for the surface density \Sigma_N(d) of KBOs larger than 1 km: \Sigma_N(d>=1 km) <= 2.0e8 deg^-2, and for KBOs larger than 0.7 km \Sigma_N(d>= 0.7 km) <= 4.8e8 deg^-2.

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.

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.

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.

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.

This paper has been withdrawn by the author.

This paper has been withdrawn

The serendipitous detection of stellar occultations by Outer Solar System objects is a powerful method for ascertaining the small end ($r \lesssim 15$ km) of the size distribution of Kuiper Belt Objects and may potentially allow the exploration of objects as far out as the Oort Cloud. The design and implementation of an occultation survey is aided by a detailed understanding of how diffraction and observational parameters affect the detection of occultation events. In this study, stellar occultations are simulated, accounting for diffraction effects, finite source sizes, finite bandwidths, stellar spectra, sampling, and signal-to-noise. Finally, the possibility of detecting small Outer Solar System objects from the Kuiper Belt all the way out to the Oort Cloud is explored for three photometric systems: a proposed space telescope, Whipple (Kaplan et al. 2003), the Taiwanese-American Occultation Survey (Lehner et al. 2006), and the Multi Mirror Telescope (Bianco 2007).

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.

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

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

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

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

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

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

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

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

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

The status of dark matter searches with inorganic scintillator detectors at Boulby mine is reviewed and the results of tests with a CsI(Tl) crystal are presented. The objectives of the latter experiment were to study anomalous fast events previously observed and to identify ways to remove this background. Clear indications were found that these events were due to surface contamination of crystals by alphas, probably from radon decay. A new array of unencapsulated NaI(Tl) crystals immersed either in liquid paraffin or pure nitrogen atmosphere is under construction at Boulby. Such an approach allows complete control of the surface of the crystals and the ability to remove any surface contamination. First data from the unencapsulated NaI(Tl) do not show the presence of anomalous fast events.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

A low pressure time projection chamber for the detection of WIMPs is discussed. Discrimination against Compton electron background in such a device should be very good, and directional information about the recoil atoms would be obtainable. If a full 3-D reconstruction of the recoil tracks can be achieved, Monte Carlo studies indicate that a WIMP signal could be identified with high confidence from as few as 30 detected WIMP-nucleus scattering events.

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

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

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.

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

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.