The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background.

We report Arecibo S-band (2380 MHz; 12.6 cm) radar observations of near-Earth asteroid (3200) Phaethon during the December 2017 apparition when Phaethon passed within 0.07 au of Earth. Radar images with a resolution of 75 m per pixel reveal a roughly spheroidal shape more than 6 km in diameter at the equator with several discernible surface features hundreds of meters in extent. These include a possible crater more than 1 km across located below 30$^{\circ}$ latitude and a roughly 600-m radar-dark region near one of the poles. Overall, the radar images of Phaethon are reminiscent of those of (101955) Bennu, target of the OSIRIS-REx mission. As such, the shape of Phaethon is suspected to have an equatorial ridge similar to the top-shaped models of several other radar-observed near-Earth asteroids as well as the optical images of (162173) Ryugu returned by the Hayabusa2 spacecraft. Preliminary analysis of the radar data finds no satellites and gives no indication of a dusty coma at the time of these observations.

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

We present thermal properties and an improved shape model for potentially hazardous asteroid (162421) 2000 ET70. In addition to the radar data from 2000 ET70's apparition in 2012, our model incorporates optical lightcurves and infrared spectra that were not included in the analysis of Naidu et al. (2013, Icarus 226, 323-335). We confirm the general "clenched fist" appearance of the Naidu et al. model, but compared to their model, our best-fit model is about 10% longer along its long principal axis, nearly identical along the intermediate axis, and about 25% shorter along the short axis. We find the asteroid's dimensions to be 2.9 km $\times$ 2.2 km $\times$ 1.5 km (with relative uncertainties of about 10%, 15%, and 25%, respectively). With the available data, 2000 ET70's period and pole position are degenerate with each other. The radar and lightcurve data together constrain the pole direction to fall along an arc that is about twenty-three degrees long and eight degrees wide. Infrared spectra from the NASA InfraRed Telescope Facility (IRTF) provide an additional constraint on the pole. Thermophysical modeling, using our SHERMAN software, shows that only a subset of the pole directions, about twelve degrees of that arc, are compatible with the infrared data. Using all of the available data, we find that 2000 ET70 has a sidereal rotation period of 8.944 hours ($\pm$ 0.009 h) and a north pole direction of ecliptic coordinates $(52^{\circ}, -60^{\circ}) \pm 6^{\circ}$. The infrared data, acquired over several dates, require that the thermal properties (albedo, thermal inertia, surface roughness) must change across the asteroid's surface. By incorporating the detailed shape model and spin state into our thermal modeling, the multiple ground-based observations at different viewing geometries have allowed us to constrain the levels of the variations in the surface properties.

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.

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.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

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.

The long term study of the Sun is necessary if we are to determine the evolution of sunspot properties and thereby inform modeling of the solar dynamo, particularly on scales of a solar cycle. We aim to determine a number of sunspot properties over cycle 23 using the uniform database provided by the SOHO Michelson Doppler Imager data. We focus in particular on their distribution on the solar disk, maximum magnetic field and umbral/penumbral areas. We investigate whether the secular decrease in sunspot maximum magnetic field reported in Kitt Peak data is present also in MDI data. We have used the Sunspot Tracking And Recognition Algorithm (STARA) to detect all sunspots present in the SOHO Michelson Doppler Imager continuum data giving us 30 084 separate detections. We record information on the sunspot locations, area and magnetic field properties and corresponding information for the umbral areas detected within the sunspots, and track them through their evolution. We find the total visible umbral area is 20-40% of the total visible sunspot area at all stages of the solar cycle. We find that the number of sunspots observed follows the Solar Influences Data Centre International Sunspot Number with some interesting deviations. Finally, we use the magnetic information in our catalogue to study the long term variation of magnetic field strength within sunspot umbrae and find that it increases and decreases along with the sunspot number. However, if we were to assume a secular decreaseas was reported in the Kitt Peak data and take into account sunspots throughout the whole solar cycle we would find the maximum umbral magnetic fields to be decreasing by 23.6 \pm 3.9 Gauss per year, which is far less than has previously been observed by other studies. If we only look at the declining phase of cycle 23 we find the decrease in sunspot magnetic fields to be 70 Gauss per year.

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)

A new algorithm for automatic detection of prominences on the solar limb in 304 A EUV images is presented, and results of its application to SOHO/EIT data discussed. The detection is based on the method of moments combined with a classifier analysis aimed at discriminating between limb prominences, active regions, and the quiet corona. This classifier analysis is based on a Support Vector Machine (SVM). Using a set of 12 moments of the radial intensity profiles, the algorithm performs well in discriminating between the above three categories of limb structures, with a misclassification rate of 7%. Pixels detected as belonging to a prominence are then used as starting point to reconstruct the whole prominence by morphological image processing techniques. It is planned that a catalogue of limb prominences identified in SOHO and STEREO data using this method will be made publicly available to the scientific community.

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.

The distributions of sunspot longitude at first appearance and at disappearance display an east-west asymmetry that results from a reduction in visibility as one moves from disk centre to the limb. To first order, this is explicable in terms of simple geometrical foreshortening. However, the centre-to-limb visibility variation is much larger than that predicted by foreshortening. Sunspot visibility is also known to be affected by the Wilson effect: the apparent dish shape of the sunspot photosphere caused by the temperature-dependent variation of the geometrical position of the tau=1 layer. In this article we investigate the role of the Wilson effect on the sunspot appearance distributions, deducing a mean depth for the umbral tau=1 layer of 500 to 1500 km. This is based on the comparison of observations of sunspot longitude distribution and Monte Carlo simulations of sunspot appearance using different models for spot growth rate, growth time and depth of Wilson depression.

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

The Northern Sky Variability Survey (NSVS) is a temporal record of the sky over the optical magnitude range from 8 to 15.5. It was conducted in the course of the first generation Robotic Optical Transient Search Experiment (ROTSE-I) using a robotic system of four co-mounted unfiltered telephoto lenses equipped with CCD cameras. The survey was conducted from Los Alamos, NM, and primarily covers the entire northern sky. Some data in southern fields between declinations 0 and -38 deg is also available, although with fewer epochs and noticeably lesser quality. The NSVS contains light curves for approximately 14 million objects. With a one year baseline and typically 100-500 measurements per object, the NSVS is the most extensive record of stellar variability across the bright sky available today. In a median field, bright unsaturated stars attain a point to point photometric scatter of ~0.02 mag and position errors within 2 arcsec. At Galactic latitudes |b| < 20 deg the data quality is limited by severe blending due to ~14 arcsec pixel size. We present basic characteristics of the data set and describe data collection, analysis, and distribution. All NSVS photometric measurements are available for on-line public access from the Sky Database for Objects in Time-Domain (SkyDOT; http://skydot.lanl.gov) at LANL. Copies of the full survey photometry may also be requested on tape.

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.

The observation of a prompt optical flash from GRB990123 convincingly demonstrated the value of autonomous robotic telescope systems. Pursuing a program of rapid follow-up observations of gamma-ray bursts, the Robotic Optical Transient Search Experiment (ROTSE) has developed a next-generation instrument, ROTSE-III, that will continue the search for fast optical transients. The entire system was designed as an economical robotic facility to be installed at remote sites throughout the world. There are seven major system components: optics, optical tube assembly, CCD camera, telescope mount, enclosure, environmental sensing & protection and data acquisition. Each is described in turn in the hope that the techniques developed here will be useful in similar contexts elsewhere.

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 untriggered search for optical bursts with the ROTSE-I telephoto array. Observations were taken which monitor an effective 256 square degree field continuously over 125 hours to m_{ROTSE}=15.7. The uniquely large field, moderate limiting magnitude and fast cadence of $\sim$10 minutes permits transient searches in a new region of sensitivity. Our search reveals no candidate events. To quantify this result, we simulate potential optical bursts with peak magnitude, m_{p}, at t=10 s, which fade as f=(\frac{t}{t_{0}}) ^{\alpha_{t}}, where \alpha_t < 0. Simple estimates based on observational evidence indicate that a search of this sensitivity begins to probe the possible region occupied by GRB orphan afterglows. Our observing protocol and image sensitivity result in a broad region of high detection efficiency for light curves to the bright and slowly varying side of a boundary running from [\alpha_{t},m_{p}]=[-2.0,6.0] to [-0.3,13.2]. Within this region, the integrated rate of brief optical bursts is less than 1.1\times 10^{-8} {\rm s}^{-1} {\rm deg}^{-2}. At $\sim$22 times the observed GRB rate from BATSE, this suggests a limit on \frac{\theta_{opt}}{\theta_{\gamma}}\lesssim 5 where \theta_{opt} and \theta_{\gamma} are the optical and gamma-ray collimation angles, respectively. Several effects might explain the absence of optical bursts, and a search of the kind described here but more sensitive by about 4 magnitudes should offer a more definitive probe.

We report on the current operating status of the ROTSE-IIIa telescope, currently undergoing testing at Los Alamos National Laboratories in New Mexico. It will be shipped to Siding Spring Observatory, Australia, in first quarter 2002. ROTSE-IIIa has been in automated observing mode since early October, 2001, after completing several weeks of calibration and check-out observations. Calibrated lists of objects in ROTSE-IIIa sky patrol data are produced routinely in an automated pipeline, and we are currently automating analysis procedures to compile these lists, eliminate false detections, and automatically identify transient and variable objects. The manual application of these procedures has already led to the detection of a nova that rose over six magnitudes in two days to a maximum detected brightness of m_R~13.9 and then faded two magnitudes in two weeks. We also readily identify variable stars, includings those suspected to be variables from the Sloan Digital Sky Survey. We report on our system to allow public monitoring of the telescope operational status in real time over the WWW.

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

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

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

We report on the discovery of an outburst from a previously unknown high- latitude CV we designate ROTSE3 J015118.59-022300.1. The object was first detected on 2001 Oct 13.291 at m-ROTSE3=14.71+-0.06, after which it faded by more than 2 magnitudes over 13 d. Spectral analysis confirms a Galactic location, despite unusual properties that render classification difficult.

The X-ray nova XTE J1118+480 exhibited two outbursts in the early part of 2000. As detected by the Rossi X-ray Timing Explorer (RXTE), the first outburst began in early January and the second began in early March. Routine imaging of the northern sky by the Robotic Optical Transient Search Experiment (ROTSE) shows the optical counterpart to XTE J1118+480 during both outbursts. These data include over 60 epochs from January to June 2000. A search of the ROTSE data archives reveal no previous optical outbursts of this source in selected data between April 1998 and January 2000. While the X-ray to optical flux ratio of XTE J1118+480 was low during both outbursts, we suggest that they were full X-ray novae and not mini-outbursts based on comparison with similar sources. The ROTSE measurements taken during the March 2000 outburst also indicate a rapid rise in the optical flux that preceded the X-ray emission measured by the RXTE by approximately 10 days. Using these results, we estimate a pre-outburst accretion disk inner truncation radius of 1.2 x 10^4 Schwarzschild radii.

Gamma-ray bursts of short duration may harbor vital clues to the range of phenomena producing bursts. However, recent progress from the observation of optical counterparts has not benefitted the study of short bursts. We have searched for early optical emission from six gamma-ray bursts using the ROTSE-I telephoto array. Three of these events were of short duration, including GRB 980527 which is among the brightest short bursts yet observed. The data consist of unfiltered CCD optical images taken in response to BATSE triggers delivered via the GCN. For the first time, we have analyzed the entire 16 degree by 16 degree field covered for five of these bursts. In addition, we discuss a search for the optical counterpart to GRB 000201, a well-localized long burst. Single image sensitivities range from 13th to 14th magnitude around 10 s after the initial burst detection, and 14 - 15.8 one hour later. No new optical counterparts were discovered in this analysis suggesting short burst optical and gamma-ray fluxes are uncorrelated.

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.

In order to observe nearly simultaneous emission from Gamma-ray Bursts (GRBs), the Robotic Optical Transient Search Experiment (ROTSE) receives triggers via the GRB Coordinates Network (GCN). Since beginning operations in March, 1998, ROTSE has also taken useful data for 10 SGR events: 8 from SGR 1900+14 and 2 from SGR 1806-20. We have searched for new or variable sources in the error regions of these SGRs and no optical counterparts were observed. Limits are in the range m_ROTSE ~ 12.5 - 15.5 during the period 20 seconds to 1 hour after the observed SGR events.

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]

The Robotic Optical Transient Search Experiment (ROTSE) seeks to measure simultaneous and early afterglow optical emission from gamma-ray bursts (GRBs). A search for optical counterparts to six GRBs with localization errors of 1 square degree or better produced no detections. The earliest limiting sensitivity is m(ROTSE) > 13.1 at 10.85 seconds (5 second exposure) after the gamma-ray rise, and the best limit is m(ROTSE) > 16.0 at 62 minutes (897 second exposure). These are the most stringent limits obtained for GRB optical counterpart brightness in the first hour after the burst. Consideration of the gamma-ray fluence and peak flux for these bursts and for GRB990123 indicates that there is not a strong positive correlation between optical flux and gamma-ray emission.

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)

The ROTSE-I experiment has generated CCD photometry for the entire Northern sky in two epochs nightly since March 1998. These sky patrol data are a powerful resource for studies of astrophysical transients. As a demonstration project, we present first results of a search for periodic variable stars derived from ROTSE-I observations. Variable identification, period determination, and type classification are conducted via automatic algorithms. In a set of nine ROTSE-I sky patrol fields covering about 2000 square degrees we identify 1781 periodic variable stars with mean magnitudes between m_v=10.0 and m_v=15.5. About 90% of these objects are newly identified as variable. Examples of many familiar types are presented. All classifications for this study have been manually confirmed. The selection criteria for this analysis have been conservatively defined, and are known to be biased against some variable classes. This preliminary study includes only 5.6% of the total ROTSE-I sky coverage, suggesting that the full ROTSE-I variable catalog will include more than 32,000 periodic variable stars.

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

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

The Robotic Optical Transient Search Experiment (ROTSE) seeks to measure contemporaneous and early afterglow optical emission from gamma-ray bursts (GRBs). The ROTSE-I telescope array has been fully automated and responding to burst alerts from the GRB Coordinates Network since March 1998, taking prompt optical data for 30 bursts in its first year. We will briefly review observations of GRB990123 which revealed the first detection of an optical burst occurring during the gamma-ray emission, reaching 9th magnitude at its peak. In addition, we present here preliminary optical results for seven other gamma-ray bursts. No other optical counterparts were seen in this analysis, and the best limiting sensitivities are m(V) > 13.0 at 14.7 seconds after the gamma-ray rise, and m(V) > 16.4 at 62 minutes. These are the most stringent limits obtained for GRB optical counterpart brightness in the first hour after the burst. This analysis suggests that there is not a strong correlation between optical flux and gamma-ray emission.

The origin of gamma-ray bursts (GRBs) has been enigmatic since their discovery. The situation improved dramatically in 1997, when the rapid availability of precise coordinates for the bursts allowed the detection of faint optical and radio afterglows - optical spectra thus obtained have demonstrated conclusively that the bursts occur at cosmological distances. But, despite efforts by several groups, optical detection has not hitherto been achieved during the brief duration of a burst. Here we report the detection of bright optical emission from GRB990123 while the burst was still in progress. Our observations begin 22 seconds after the onset of the burst and show an increase in brightness by a factor of 14 during the first 25 seconds; the brightness then declines by a factor of 100, at which point (700 seconds after the burst onset) it falls below our detection threshold. The redshift of this burst, approximately 1.6, implies a peak optical luminosity of 5 times 10^{49} erg per second. Optical emission from gamma-ray bursts has been generally thought to take place at the shock fronts generated by interaction of the primary energy source with the surrounding medium, where the gamma-rays might also be produced. The lack of a significant change in the gamma-ray light curve when the optical emission develops suggests that the gamma-rays are not produced at the shock front, but closer to the site of the original explosion.

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

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

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

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

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

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.

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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