This paper has been withdrawn by the author.

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

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.

The BL Lacertae object Markarian 501 was identified as a source of gamma-ray emission at the Whipple Observatory in March 1995. Here we present a flux variability analysis on several times-scales of the 233 hour data set accumulated over 213 nights (from March 1995 to July 1998) with the Whipple Observatory 10 m atmospheric Cherenkov imaging telescope. In 1995, with the exception of a single night, the flux from Markarian 501 was constant on daily and monthly time-scales and had an average flux of only 10% that of the Crab Nebula, making it the weakest VHE source detected to date. In 1996, the average flux was approximately twice the 1995 flux and showed significant month-to-month variability. No significant day-scale variations were detected. The average gamma-ray flux above ~350 GeV in the 1997 observing season rose to 1.4 times that of the Crab Nebula -- 14 times the 1995 discovery level -- allowing a search for variability on time-scales shorter than one day. Significant hour-scale variability was present in the 1997 data, with the shortest, observed on MJD 50607, having a doubling time of ~2 hours. In 1998 the average emission level decreased considerably from that of 1997 (to ~20% of the Crab Nebula flux) but two significant flaring events were observed. Thus, the emission from Markarian 501 shows large amplitude and rapid flux variability at very high energies as does Markarian 421. It also shows large mean flux level variations on year-to-year time-scales, behaviour which has not been seen from Markarian 421 so far.

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.

The close relation between ground-based TeV observations and satellite borne $\gamma$-ray measurements has been important for the understanding of blazars. The observations which involve the TeV component in blazar studies are reviewed.

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.

From considering the effect of gamma-gamma interactions on recently observed TeV gamma-ray spectra, improved limits are set to the density of extragalactic infrared (IR) photons which are robust and essentially model-independent. The resulting limits are up to two orders of magnitude more restrictive than direct observations in the 0.025-0.3eV regime. These limits are used to improve constraints on radiative neutrino decay in the mass range above 0.05eV and on Very Massive Objects (VMOs) as providing the dark matter needed to explain galaxy rotation curves.

We present the discovery of >350 GeV gamma-ray emission from the BL Lacertae object 1ES 2344+514 with the Whipple Observatory 10m gamma-ray telescope. This is the third BL Lac object detected at gamma-ray energies above 300 Gev, the other two being Markarian 421 (Mrk 421) and Mrk501. These three active galactic nuclei are all X-ray selected and have the lowest known redshifts of any BL Lac objects currently identified. The evidence for emission derives primarily from an apparent flare on December 20, 1995 when a 6 sigma excess was detected with a flux approximately 63% of the very high energy gamma-ray emission from the Crab Nebula, the standard candle for TeV gamma-ray sources. Excluding the flare, observations between October 1995 and January 1996 yield a 4 sigma detection corresponding to 11% of the VHE Crab Nebula flux. Observations spanning September 1996 to January 1997 failed to yield a significant detection of a steady flux or any flaring. For this period, the 99.9% confidence level upper limit is <8% of the Crab Nebula. The low baseline emission level and variations in nightly and yearly flux of 1ES 22344+514 are the same as the VHE emission characteristics of Mrk 421 and Mrk 501

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

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

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

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

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

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

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

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

We present multiwavelength observations of the BL Lacertae object Markarian 501 (Mrk 501) in 1997 between April 8 and April 19. Evidence of correlated variability is seen in very high energy (VHE, E > 350 GeV) gamma-ray observations taken with the Whipple Observatory gamma-ray telescope, data from the Oriented Scintillation Spectrometer Experiment of the Compton Gamma-Ray Observatory, and quicklook results from the All-Sky Monitor of the Rossi X-ray Timing Explorer while the Energetic Gamma-Ray Experiment Telescope did not detect Mrk 501. Short term optical correlations are not conclusive but the U-band flux observed with the 1.2m telescope of the Whipple Observatory was 10% higher than in March. The average energy output of Mrk 501 appears to peak in the 2 keV to 100 keV range suggesting an extension of the synchrotron emission to at least 100 keV, the highest observed in a blazar and ~100 times higher than that seen in the other TeV-emitting BL Lac object, Mrk 421. The VHE gamma-ray flux observed during this period is the highest ever detected from this object. The VHE gamma-ray energy output is somewhat lower than the 2-100 keV range but the variability amplitude is larger. The correlations seen here do not require relativistic beaming of the emission unless the VHE spectrum extends to >5 TeV.

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.

The ~70 unidentified sources of the EGRET sky survey may be one of its most important legacies. The identification of these sources at other wavelengths is critical to understanding their nature. Many have flat spectra out to 10 GeV which, if extrapolated to TeV energies, would be easily detectable relative to the steeply falling diffuse background. The Whipple Observatory gamma-ray telescope has been used to observe a number of these which were selected based on their position, intensity and spectrum and in some cases based on a possible association with a supernova remnant or pulsar. No significant emission has been detected from these sources, and upper limits are given.

The Crab Nebula has become established as the standard candle for TeV gamma-ray astronomy using the atmospheric Cherenkov technique. No evidence for variability has been seen. The spectrum of gamma rays from the Crab Nebula has been measured in the energy range 500 GeV to 8 TeV at the Whipple Observatory by the atmospheric Cherenkov imaging technique. Two methods of analysis involving independent Monte Carlo simulations and two databases of observations (1988-89 and 1995-96) were used and gave close agreement. Using the complete spectrum of the Crab Nebula, the spectrum of relativistic electrons is deduced and the spectrum of the resulting inverse Compton gamma-ray emission is in good agreement with the measured spectrum if the ambient magnetic field is about 25-30 nT.

PSR B1951+32 is a gamma-ray pulsar detected by the Energetic Gamma Ray Experiment Telescope (EGRET) and identified with the 39.5 ms radio pulsar in the supernova remnant CTB 80. The EGRET data shows no evidence for a spectral turnover. Here we report on the first observations of PSR B1951+32 beyond 30 GeV. The observations were carried out with the 10m gamma-ray telescope at the Whipple Observatory on Mt. Hopkins, Arizona. In 8.1 hours of observation we find no evidence for steady or periodic emission from PSR B1951+32 above ~260 GeV. FLux upper limits are derived and compared with model extrapolations from lower energies and the predictions of emission models.

Markarian 501 is only the second extragalactic source to be detected with high statistical certainty at TeV energies; it is similar in many ways to Markarian 421. The Whipple Observatory gamma-ray telescope has been used to observe the AGN Markarian 501 in 1996 and 1997, the years subsequent to its initial detection. The apparent variability on the one-day time-scale observed in TeV gamma rays in 1995 is confirmed and compared with the variability in Markarian 421. Observations at X-ray and optical wavelengths from 1997 are also presented.

Two active galactic nuclei have been detected at TeV energies using the atmospheric Cherenkov imaging technique. The Whipple Observatory gamma-ray telescope has been used to observe all the BL Lacertae objects in the northern hemisphere out to a redshift of 0.1. We report the tentative detection of VHE emission from a third BL Lac object, 1ES 2344+514. Progress in extending this survey out to z=0.2 will also be reported.

Based on BACODINE network notification the Whipple Observatory gamma-ray telescope has been used to search for the delayed TeV counterpart to BATSE-detected gamma-ray bursts. In the fast slew mode, any point in the sky can be reached within two minutes of the burst notification. The search strategy, necessary because of the uncertainty in burst position and limited FOV of the camera, is described.

If supernova remnants (SNRs) are the sites of cosmic-ray acceleration, the associated nuclear interactions should result in observable fluxes of TeV gamma-rays from the nearest SNRs. Measurements of the gamma-ray flux from six nearby, radio-bright, SNRs have been made with the Whipple Observatory gamma-ray telescope. No significant emission has been detected and upper limits on the $>$300 GeV flux are reported. Three of these SNRs (IC443, gamma-Cygni and W44) are spatially coincident with low latitude unidentified sources detected with EGRET. These upper limits weaken the case for the simplest models of shock acceleration and energy dependent propagation.

A search was conducted for TeV gamma-rays emitted from the direction of the ultra-high energy cosmic ray detected by the Fly's Eye Experiment with E ~ 3 x 10**20 eV. No enhancement was found at a level of 10**-10 gamma/cm**2-sec for E>350 GeV. This upper limit is consistent with theoretical estimates based on topological defects as sources of UHE cosmic rays. An upper limit was also set for the flux of TeV gamma rays from 3C147, the most prominent AGN in the error box.

Markarian 421 was the first extragalactic source to be detected with high statistical certainty at TeV energies. The Whipple Observatory gamma-ray telescope has been used to observe the Active Galactic Nucleus, Markarian 421 in 1996 and 1997. The rapid variability observed in TeV gamma rays in previous years is confirmed. Doubling times as short as 15 minutes are reported with flux levels reaching 15 photons per minute. The TeV energy spectrum is derived using two independent methods. The implications for the intergalactic infra-red medium of an observed unbroken power law spectrum up to energies of 5 TeV is discussed.

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