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.

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

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 an analysis of 436 lines of sight with extinction data covering wavelength range from near-infrared (NIR) to ultraviolet (UV). We use J, H, K photometry from 2MASS database, the IR intrinsic colors from Wegner (1994), and UV extinction data from Wegner (2002). We exclude 19 lines of sight (4%) from the original sample because of suspected photometric problems. We derive total to selective extinction ratios (R_V) based on the Cardelli, Clayton, and Mathis (1989; CCM) law, which is typically used to fit the extinction data for both diffuse and dense interstellar medium. We conclude that CCM law is able to fit well most of the extinction curves in our sample (71%), and we present a catalog of R_V and visual extinction (A_V) values for those cases. We divide the remaining lines of sight with peculiar extinction into two groups according to two main behaviors: a) the NIR or/and UV wavelength regions cannot be reproduced by CCM formula (14% of the entire sample), b) the NIR and UV extinction data taken separately are best fit by CCM laws with significantly different values of R_V (10% of the entire sample). We present examples of such curves. We caution that some peculiarities of the extinction curves may not be intrinsic but simply caused by faulty data. The study of the intrinsically peculiar cases could help us to learn about the physical processes that affect dust in the interstellar medium, e.g., formation of mantles on the surface of grains, evaporation, growing or shattering.

We present the analysis of a large sample of lines of sight with extinction curves covering wavelength range from near-infrared (NIR) to ultraviolet (UV). We derive total to selective extinction ratios based on the Cardelli, Clayton and Mathis (1989, CCM) law, which is typically used to fit the extinction data both for diffuse and dense interstellar medium. We conclude that the CCM law is able to fit most of the extinction curves in our sample. We divide the remaining lines of sight with peculiar extinction into two groups according to two main behaviors: a) the optical/IR or/and UV wavelength region cannot be reproduced by the CCM formula; b) the optical/NIR and UV extinction data are best fit by the CCM law with different values of R_V. We present examples of such curves. The study of both types of peculiar cases can help us to learn about the physical processes that affect dust in the interstellar medium, e.g., formation of mantles on the surface of grains, evaporation, growing or shattering.

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

We present determinations of the total to selective extinction ratio R_V and visual extinction A_V values for Milky Way stars using ultraviolet color excesses. We extend the analysis of Gnacinski and Sikorski (1999) by using non-equal weights derived from observational errors. We present a detailed discussion of various statistical errors. In addition, we estimate the level of systematic errors by considering different normalization of the extinction curve adopted by Wegner (2002). Our catalog of 782 R_V and A_V values and their errors is available in the electronic form on the World Wide Web.

It is commonly accepted that quasar redshifts have a cosmological character and that most of the quasars are at Gigaparsec distances. However, there are some cases where several quasars with completely different redshifts and a nearby active galaxy are aligned in a certain way or occupy a very small patch on the sky, which is claimed by some authors to be unlikely to happen by chance. Is there a small subset of quasars with non-cosmological redshifts? For quasars apparently associated with galaxies, we consider two scenarios for the origin of their redshift: 1. a standard, cosmological scenario, 2. a velocity-induced Doppler shift of a nearby object's spectrum (local, ejection scenario). We argue for a simple astrometric test which can distinguish between these two sources of quasar redshifts by constraining their proper motions. We give the predictions for the maximum possible proper motions of a quasar for the cosmological and local scenarios of the origin of their redshifts. We apply these theoretical results to the Bukhmastova (2001) catalog, which contains more than 8000 close QSO-galaxy associations. In the standard interpretation of quasar redshifts, their typical proper motions are a fraction of micro arc-second, and beyond the reach of planned astrometric missions like GAIA and SIM. On the other hand, the quasars ejected from local AGNs at velocities close to the speed of light would have proper motions 5-6 orders of magnitude larger, which would easily be measurable with future astrometric missions. The distributions of proper motions for the cosmological and local scenarios are very well separated. Moreover, the division corresponds nicely to the expected accuracy from GAIA and SIM.

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 present a (V-R)-based reddening map of about 43 square degrees of the Galactic bulge/bar. The map is constructed using template image photometry from the MACHO microlensing survey, contains 9717 resolution elements, and is based on (V-R)-color averages of the entire color-magnitude diagrams (CMDs) in 4 by 4 arc-minute tiles. The conversion from the observed color to the reddening follows from an assumption that CMDs of all bulge fields would look similar in the absence of extinction. Consequently, the difference in observed color between various fields originates from varying contribution of the disk extinction summed along different lines of sight. We check that our (V-R) colors correlate very well with infrared and optical reddening maps. We show that a dusty disk obeying a cosec|b| extinction law, E(V-R) = 0.0274 cosec|b|, provides a good approximation to the extinction toward the MACHO bulge/bar fields. The large-scale (V-R)-color and visual extinction map presented here is publicly available in the electronic edition of the Journal and on the World Wide Web.

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

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

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

I make a new evaluation of the microlensing optical depth toward the Galactic bar from Difference Image Analysis (DIA) of the MACHO Collaboration. First, I present supplementary evidence that MACHO field 104 located at (l,b) = (3.11,-3.01) is anomalous in terms of the event duration distribution. I argue that both the event durations and the very high optical depth of field 104 are not representative and, therefore, exclude this field as an outlier. In addition, I eliminate field 159 at (l,b) = (6.35,-4.40) based mainly on its separate location, but also on unexplained statistical properties of the event durations. The remaining six DIA fields form a very homogeneous and spatially compact set that is very suitable for averaging. The weighting of the optical depth values for these six DIA fields results in a total optical depth tau_{tot} = 2.01^{+0.34}_{-0.32} x 10^{-6} at (l,b) = (2.22,-3.18). If a fraction of all sources, f_{disk}, assumed to be in the disk, does not contribute to microlensing, then the optical depth toward the sources in the bar is tau_{bar} = 2.23^{+0.38}_{-0.35} x 10^{-6} 0.9/(1-f_{disk}). Both tau_{tot} and tau_{bar} are substantially lower than the original estimates of Alcock et al. Most of the change in the DIA-based optical depths comes from a more appropriate statistical treatment of the results in individual fields and not from the removal of fields 104 and 159. When taken together with tau_{bar} = (1.4 +/- 0.3) x 10^{-6} at (l,b) = (3.9,-3.8) as derived from clump giants, this new result suggests that the conclusions from microlensing experiments are in reasonable agreement with the expectations from infrared-based Galactic models.

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 entertain the idea that robust theoretical expectations can become a tool in removing hidden observational or data-reduction biases. We illustrate this approach for a specific problem associated with gravitational microlensing. Using the fact that a group is more than just a collection of individuals, we derive formulae for correcting the distribution of the dimensionless impact parameters of events, u_min. We refer to the case when undetected biases in the u_min distribution can be alleviated by multiplication of impact parameters of all events by a common constant factor. We show that in this case the general maximum likelihood problem of solving an infinite number of equations reduces to two constraints, and we find an analytic solution. Under the above assumptions, this solution represents a state in which the ``entropy'' of a microlensing ensemble is at its maximum, that is, the distribution of u_min resembles a specific, theoretically expected, box-like distribution to the highest possible extent. We also show that this technique does not allow one to correct the parameters of individual events on the event by event basis independently from each other.

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.

I explore the consequences of making the RR Lyrae and clump giant distance scales consistent in the solar neighborhood, Galactic bulge and Large Magellanic Cloud (LMC). I employ two major assumptions: 1) that the absolute magnitude - metallicity, M_V(RR) - [Fe/H], relation for RR Lyrae stars is universal, and 2) that absolute I-magnitudes of clump giants, M_I(RC), in Baade's Window are known (e.g., can be inferred from the local Hipparcos-based calibration or theoretical modeling). A comparison between the solar neighborhood and Baade's Window sets M_V(RR) at [Fe/H] = -1.6 in the range (0.59 +/- 0.05, 0.70 +/- 0.05), somewhat brighter than the statistical parallax solution. More luminous RR Lyrae stars imply younger ages of globular cluster, which would be in better agreement with the conclusions from the currently favored stellar evolution and cosmological models. A comparison between Baade's Window and the LMC sets the M_I^{LMC}(RC) in the range (-0.33 +/- 0.09, -0.53 +/- 0.09). The distance modulus to the LMC, mu^{LMC}, is between 18.24 +/- 0.08 and 18.44 +/- 0.07. Unlike M_I^{LMC}(RC), this range in mu^{LMC} does NOT depend on the adopted value of the dereddened LMC clump magnitude, I_0^{LMC}(RC). I argue that the currently available information is insufficient to select the correct distance scale with high confidence.

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]

I demonstrate that the two unexpected results in the local Universe: anomalous intrinsic (V-I)_0 colors of RR Lyrae stars and clump giants in the Galactic center, and very short distances to Magellanic Clouds inferred from clump giants, can be at least partially resolved with a modified coefficient of selective extinction A_V/E(V-I). With this modification, I find a new clump-giant distance modulus to the Large Magellanic Cloud, mu_{LMC} = 18.27 +/- 0.07, which is 0.09 larger than the Udalski (1998b) result. When distance estimates from the red clump, RR Lyrae stars and the eclipsing binary HV2274 are combined, one obtains mu_{LMC} = 18.31 +/- 0.04 (internal).

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

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

I demonstrate that the two unexpected results in the local Universe: 1) anomalous intrinsic (V-I)_0 colors of the clump giants and RR Lyrae stars in the Galactic center, and 2) very short distances to the Magellanic Clouds (LMC, SMC) as inferred from clump giants, are connected with each other. The (V-I)_0 anomaly is partially resolved by using the photometry from the phase-II of the Optical Gravitational Lensing Experiment (OGLE) rather than phase-I. The need for V- or I-magnitude-based change in the bulge (V-I)_0 is one option to explain the remaining color discrepancy. Such change may originate in a coefficient of selective extinction A_V/E(V-I) smaller than typically assumed. Application of the (V-I)_0 correction (independent of its source) doubles the slope of the absolute magnitude - metallicity relation for clump giants, so that M_I(RC) = -0.23 + 0.19[Fe/H]. Consequently, the estimates of the clump distances to the LMC and SMC are affected. Udalski's (1998c) distance modulus of mu_{LMC} = 18.18 +/- 0.06 increases to mu_{LMC} = 18.27 +/- 0.07. The distance modulus to the SMC increases by 0.12 to mu_{SMC} = 18.77 +/- 0.08. I argue that a more comprehensive assessment of the metallicity effect on M_I(RC) is needed.

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

We compare (V-I)_0-(V-K)_0 color-color and (V-I)_0-log P period-color diagrams for Baade's Window and local RRab Lyrae stars. We find that for a fixed log P the Baade's Window RR Lyrae stars are ~0.17 magnitudes redder in (V-I)_0 than the local RR Lyrae stars. We also show that there is no such effect observed in (V-K)_0. We argue that an extinction misestimate towards Baade's Window is not a plausible explanation of the discrepancy. Unlike Baade's Window RR Lyrae stars, the local ones follow a black-body color-color relation and are well approximated by theoretical models. We test two parameters, metallicity and surface gravity, and find that their effects are too small to explain the (V-I)_0 discrepancy between the two groups of stars. We do not provide any explanation for the anomalous (V-I)_0 behavior of the Baade's Window RR Lyrae stars. We note that a similar effect for clump giant stars has been recently reported by Paczynski and we caution that RR Lyrae stars and clump giants, often used as standard candles, can be subject to the same type of systematics.

We review seven methods of measuring the absolute magnitude M_V of RR Lyrae stars in light of the Hipparcos mission and other recent developments. We focus on identifying possible systematic errors and rank the methods by relative immunity to such errors. For the three most robust methods, statistical parallax, trigonometric parallax, and cluster kinematics, we find M_V (at [Fe/H] = -1.6) of 0.77 +/- 0.13, 0.71 +/- 0.15, 0.67 +/- 0.10. These methods cluster consistently around 0.71 +/- 0.07. We find that Baade-Wesselink and theoretical models both yield a broad range of possible values (0.45-0.70 and 0.45-0.65) due to systematic uncertainties in the temperature scale and input physics. Main-sequence fitting gives a much brighter M_V = 0.45 +/- 0.04 but this may be due to a difference in the metallicity scales of the cluster giants and the calibrating subdwarfs. White-dwarf cooling-sequence fitting gives 0.67 +/- 0.13 and is potentially very robust, but at present is too new to be fully tested for systematics. If the three most robust methods are combined with Walker's mean measurement for 6 LMC clusters, V_{0,LMC} = 18.98 +/- 0.03 at [Fe/H] = -1.9, then mu_{LMC} = 18.33 +/- 0.08.

We sing the praises of the central limit theorem. Having previously removed all other possible causes of significant systematic error in the statistical parallax determination of RR Lyrae absolute magnitudes, we investigate systematic errors from two final sources of input data: apparent magnitudes and extinctions. We find corrections due to each of ~0.05 mag, i.e., ~1/2 the statistical error. However, these are of opposite sign and so roughly cancel. The apparent magnitude system that we previously adopted from Layden et al. was calibrated to the photometry of Clube & Dawe. Using Hipparcos photometry we show that the Clube & Dawe system is ~0.06 mag too bright. Extinctions were previously pinned to the HI-based map of Burstein & Heiles. We argue that A_V should rather be based on new COBE/IRAS dust-emission map of Schlegel, Finkbeiner & Davis. This change increases the mean A_V by ~0.05 mag. We find M_V=0.77 +/- 0.13 at [Fe/H]=-1.60 for a pure sample of 147 halo RR Lyraes, or M_V=0.80 +/- 0.11 at [Fe/H]=-1.71 if we incorporate kinematic information from 716 non-kinematically selected non-RR Lyrae stars from Beers & Sommer-Larsen. These are 2 and 3 sigma fainter than recent determinations of M_V from main sequence fitting of clusters using Hipparcos measurements of subdwarfs by Reid and Gratton et al. Since statistical parallax is being cleared of systematic errors and since the chance of a >2 sigma statistical fluctuation is <1/20, we conclude that these brighter determinations may be in error. In the course of three papers, we have corrected 6 systematic errors whose absolute values total 0.20 mag. Had these, contrary to the expectation of the central limit theorem, all lined up one way, they could have resolved the conflict in favor of the brighter determinations. In fact, the net change was only 0.06 mag.

In this paper we investigate whether a misestimate of proper motions could have been a source of substantial systematic errors in the statistical parallax determination of the absolute magnitude of RR Lyrae stars. In an earlier paper, we showed that the statistical parallax method is extremely robust and rather insensitive to various systematic effects. The main potential problem with this method would therefore arise from systematically bad observational inputs, primarily radial velocities and proper motions. In that paper, we demonstrated that the radial velocities have not been systematically misestimated. Here we turn our attention to proper motions. We compare three different catalogs of proper motions --- Lick, Hipparcos and the one compiled by Wan et al. (WMJ). We find that the WMJ catalog is too heterogeneous to be a reliable source. We analyze the sample of 165 halo RR Lyrae stars with either Lick or Hipparcos proper motions. For the stars with both Lick and Hipparcos proper motions we use the weighted means of reported values. Various possible biases are investigated through vigorous Monte Carlo simulations and we evaluate small corrections due to Malmquist bias, anisotropic positions of the stars on the sky, and non-Gaussian distribution of stellar velocities. The mean RR Lyrae absolute magnitude is M_V = 0.74 +/- 0.12 at the mean metallicity of the sample <[Fe/H]> = -1.60, only 0.01 mag brighter than the value obtained in the previous study which did not incorporate Hipparcos proper motions. The faint absolute magnitudes of RR Lyrae stars confirmed by this analysis gives strong support to the short distance scale.

We present a mathematical analysis of the statistical parallax method. The method yields physical insight into the maximum-likelihood determinations of the luminosity and velocity distribution and enables us to conduct a vigorous Monte Carlo investigation into various systematic effects. We apply our analytic formalism to the RR Lyrae sample of Layden et al. The velocity distribution of RR Lyrae stars is highly non-Gaussian, with kurtoses K_\pi = 2.04, K_\theta = 3.22 and K_z = 4.28 in the three principal directions, but this has almost no effect on either the best fit or the uncertainty of the luminosity determination. Indeed, our principal result is that the statistical parallax method is extremely robust in the face of all systematic effects that we considered. Our analysis, applied to the Layden et al. RR Lyrae sample, strictly confirms the majority of their results. The mean RR Lyrae absolute magnitude is M_V = 0.75 +/- 0.13 at the mean metallicity of the sample <[Fe/H]> = -1.61, compared to M_V = 0.71 +/- 0.12 obtained by Layden et al. Most of the difference is due to Malmquist bias which was not considered in previous studies. We also analyze a semi-independent non-kinematically selected sample of stars with metallicities at the [Fe/H] <= 1.5 taken from Layden et al. and Beers & Sommer-Larsen and obtain M_V = 0.79 +/- 0.12 at <[Fe/H]> = -1.79. Additionally, this analysis yields measurements of the radial bulk motion (4 +/- 10 km/s) and vertical bulk motion (0 +/- 6 km/s) of the halo relative to the Local Standard of Rest.

We have detected seven stars with a common proper motion which are located within 2.5 arcminutes of the globular cluster NGC 6522 in the Baade's Window field of the Galactic bulge. We argue that these stars are members of the cluster, and derive a weighted mean proper motion and heliocentric radial velocity of mu_l = 1.4 +/- 0.2 mas/yr, mu_b = -6.2 +/- 0.2 mas/yr, v = -28.5 +/- 6.5 km/s. We rederive the distance to NGC 6522 (0.91 +/- 0.04 R_0, where R_0 is the Galactocentric distance) and metallicity ([Fe/H] = -1.28 +/- 0.12), making use of recent revisions in the foreground extinction toward the cluster (A_V = 1.42 +/- 0.05). We find the spatial velocity of the cluster and conclude that the cluster stays close to the Galactic center, and may have experienced significant bulge/disk shocking during its lifetime.

The non-Euclidean geometry of spacetime induces an anisotropy in the apparent correlation function of high-redshift quasars. This effect can constrain the cosmological constant \Lambda independent of any assumptions about evolution of luminosities, sizes, or clustering. We examine the prospects for distinguishing between low-density (\Omega_0=0.1-0.4) cosmological models with flat and open space geometry using the quasar samples anticipated from the 2dF and Sloan redshift surveys. We show that even these large quasar surveys are likely to reside in the "sparse sampling" regime, so that measurement errors in the correlation function obey Poisson statistics. As a result: (a) one can devise a simple maximum-likelihood scheme for estimating clustering parameters, (b) one can generate Monte Carlo realizations of correlation function measurements without creating artificial quasar distributions, and (c) for fixed quasar number, a deeper survey over a smaller area has greater statistical power than a shallow, large-area survey. Adopting recent (quite uncertain) estimates of the quasar correlation length, we find that the 2dF and Sloan samples can provide clear discrimination between flat and open geometries for \Omega <= 0.2 but only marginal discrimination for \Omega = 0.4. Clear discrimination is possible for \Omega = 0.4 if the true quasar correlation length is a factor of two larger, and a high-density survey of 30,000 quasars in 200 square degrees would provide clear discrimination even for the lower correlation length.

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 measure the zero point of the Stanek (1996) extinction map by comparing the observed (V-K) colors of 206 K giant stars with their intrinsic (V-K)_0 colors as derived from their H\beta indices. We find that the zero point of the Stanek map should be changed by \Delta A_V = -0.10 +/- 0.06 mag, obtaining as a bonus a three-fold reduction of the previous statistical error. The most direct way to test for systematic errors in this determination would be to conduct a parallel measurement based on the (V-K) colors of RR Lyraes (type ab).

The distance modulus of a galaxy determined from Cepheids differs from its distance modulus determined from Tully-Fisher by an amount that is proportional to the galaxy's line width, W: Delta mu ~ (-1.73\pm 0.46)log(W). While a miscalibration of the slope of Tully-Fisher could in principle produce this effect, we argue that such a miscalibration is very unlikely. The other possible explanation is that the inferred Cepheid luminosity is correlated with the rotation speed v (and hence W) of its parent galaxy: L ~ v^{0.7}. This proportionality is superposed on the well-established relation between period and luminosity.) Such a dependence would be expected if Cepheid luminosity is correlated with metallicity, since galaxies with deeper potential wells tend to retain more of their metals. It would induce a discrepancy of 0.83 mag between H_0 determinations that are calibrated using M31 (log W=2.71) and the H_0 determination from SNIa, which is calibrated in galaxies with mean line width <log W>=2.23. A discrepancy of 0.83 mag corresponds to the difference between 80 km/s/Mpc and 55 km/s/Mpc, very similar to the actual discrepancy reported in the current literature.