We present the analysis of four hours of spectroscopic observations of NGC 6946 with the SITELLE Imaging Fourier Transform Spectrometer on the Canada-France-Hawaii Telescope, acquired to search for supernova light echoes from its ten modern supernovae. We develop a novel spectroscopic search method: identifying negatively sloped continua in the narrow-band SN3 filter as candidate highly-broadened P-Cygni profiles in the H$\alpha$ line, which would be characteristic of the spectra of supernovae ejecta. We test our methodology by looking for light echoes from any of the ten supernovae observed in NGC 6946 in the past 100 years. We find no evidence of light echoes above the survey surface brightness limit of 1$\times$10$^{-15}$erg/s/cm$^2$/arcsec$^2$.

It is increasingly suspected that the rare R Coronae Borealis (RCB) stars - hydrogen-deficient and carbon-rich supergiant stars - are the products of mergers of CO/He white-dwarf binary systems in the intermediate mass regime ($0.6<M_{Tot}<1.2 M_{\odot}$). Only 77 RCB stars are currently known in our Galaxy while up to 1000 were expected. It is necessary to find more of these peculiar and diverse stars to understand their origin and evolutionary path. We are undertaking such a dedicated search. We plan to follow up spectroscopically 2356 targets of interest that were carefully selected using the all sky 2MASS and WISE surveys. We have observed nearly 500 of these targets using optical low-resolution spectrographs. These spectra were compared to synthetic spectra from a new grid of MARCs hydrogen-deficient atmospheric models. Classical RCB stars' photospheric temperatures range mostly from 4000 K to about 8500 K, and therefore their spectra look very different showing the presence of carbon molecules C$_2$ and CN up to $\sim$6800 K and solely atomic absorption lines above that. We have put in place a series of criteria to distinguish RCB stars from other AGB carbon-rich stars. We found 45 new RCB stars, including 30 Cold ($4000<T_{eff}<6800$ K), 14 Warm ($6800<T_{eff}<8500$ K) and one hot RCB ($T_{eff}>15000$ K). Forty of these belong to the Milky Way and five are located in the Magellanic Clouds. We also confirmed that the long lasting candidate KDM 5651 is indeed a new Magellanic RCB star, increasing the total number of Magellanic Cloud RCB stars to 30. We have increased by $\sim$50\% the total number of RCB stars known, now reaching 147. We also include a list of 14 strong RCB candidates, most certainly observed during a dust obscuration phase. From the detection efficiency and success rate so far, we estimate that there should be no more than 500 RCB stars/HdC stars in the Milky Way.

R Coronae Borealis stars (RCBs) are the best candidates to be the products of white dwarf mergers in the intermediate mass range ($0.4<M_{tot}<1.1$ M$_{\odot}$). After a merger, a He envelope surrounds a CO core, and a short-lived supergiant phase starts. RCBs correspond to that phase. They are known to be hydrogen-deficient and carbon-rich supergiant stars. They are rare; They are also highly diverse in chemical abundances, brightness and photosphere temperature. To understand their evolutionary path, constrain their spatial distribution and test their formation rate, we aim to increase significantly the number of known RCBs. To do so, we need first to create a list of targets of interest rich in RCB stars, that will be subsequently followed-up spectroscopically to reveal their true nature. Fortunately, RCBs are producing dust and are therefore surrounded by warm circumstellar dust shells. Using the all-sky 2MASS and WISE infrared surveys, our goal is to discriminate RCBs among other dust-producing stars. A series of selection criteria was applied to colour-colour diagrams and to the brightness of all objects catalogued. The criteria were designed using all 102 known RCB stars. RCB SED models were also used to understand the limits of the criteria and the possible bias in terms of shell temperature selected. The analysis was complicated by the unpredictable photometric declines that RCB stars are known to undergo. With a detection efficiency of about 85\%, our analysis selected only 2356 targets among the 563 million objects catalogued in the WISE All-Sky survey. This short list of targets of interest supersedes a previously published one made from the preliminary WISE data release. This list, published in the CDS Vizier database, will allow us to increase significantly the number of known RCB stars.

We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.

Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.

The Hubble Space Telescope has imaged the central part of M87 over a 10 week span, leading to the discovery of 32 classical novae and nine fainter, likely very slow and/or symbiotic novae. In this first in a series of papers we present the M87 nova finder charts, and the light and color curves of the novae. We demonstrate that the rise and decline times, and the colors of M87 novae are uncorrelated with each other and with position in the galaxy. The spatial distribution of the M87 novae follows the light of the galaxy, suggesting that novae accreted by M87 during cannibalistic episodes are well-mixed. Conservatively using only the 32 brightest classical novae we derive a nova rate for M87: $363_{-45}^{+33}$ novae/yr. We also derive the luminosity-specific classical nova rate for this galaxy, which is $7.88_{-2.6}^{+2.3} /yr/ 10^{10}L_\odot,_{K}$. Both rates are 3-4 times higher higher than those reported for M87 in the past, and similarly higher than those reported for all other galaxies. We suggest that most previous ground-based surveys for novae in external galaxies, including M87, miss most faint, fast novae, and almost all slow novae near the centers of galaxies.

We present results of an optical search for Cepheid variable stars using the Hubble Space Telescope (HST) in 19 hosts of Type Ia supernovae (SNe Ia) and the maser-host galaxy NGC 4258, conducted as part of the SH0ES project (Supernovae and H0 for the Equation of State of dark energy). The targets include 9 newly imaged SN Ia hosts using a novel strategy based on a long-pass filter that minimizes the number of HST orbits required to detect and accurately determine Cepheid properties. We carried out a homogeneous reduction and analysis of all observations, including new universal variability searches in all SN Ia hosts, that yielded a total of 2200 variables with well-defined selection criteria -- the largest such sample identified outside the Local Group. These objects are used in a companion paper to determine the local value of H0 with a total uncertainty of 2.4%.

We present the analysis of multiband time-series data for a sample of 24 Cepheids in the field of the Large Magellanic Cloud cluster NGC1866. Very accurate BVI VLT photometry is combined with archival UBVI data, covering a large temporal window, to obtain precise mean magnitudes and periods with typical errors of 1-2% and of 1 ppm, respectively. These results represent the first accurate and homogeneous dataset for a substantial sample of Cepheid variables belonging to a cluster and hence sharing common distance, age and original chemical composition. Comparisons of the resulting multiband Period-Luminosity and Wesenheit relations to both empirical and theoretical results for the Large Magellanic Cloud are presented and discussed to derive the distance of the cluster and to constrain the mass-luminosity relation of the Cepheids. The adopted theoretical scenario is also tested by comparison with independent calibrations of the Cepheid Wesenheit zero point based on trigonometric parallaxes and Baade-Wesselink techniques. Our analysis suggests that a mild overshooting and/or a moderate mass loss can affect intermediate-mass stellar evolution in this cluster and gives a distance modulus of 18.50 +- 0.01 mag. The obtained V,I color-magnitude diagram is also analysed and compared with both synthetic models and theoretical isochrones for a range of ages and metallicities and for different efficiencies of core overshooting. As a result, we find that the age of NGC1866 is about 140 Myr, assuming Z = 0.008 and the mild efficiency of overshooting suggested by the comparison with the pulsation models.

For over a century, light echoes have been observed around variable stars and transients. The discovery of centuries-old light echoes from supernovae in the Large Magellanic Cloud has allowed the spectroscopic characterization of these events using modern instrumentation, even in the complete absence of any visual record of those events. Here we review the pivotal role the Blanco 4m telescope played in these discoveries.

We present VI-band photometry for stars in the metal-rich globular cluster NGC 6496. Our time-series data were cadenced to search for long period variables (LPVs) over a span of nearly two years, and our variability search yielded the discovery of thirteen new variable stars, of which six are LPVs, two are suspected LPVs, and five are short period eclipsing binaries. An additional star was found in the ASAS database, and we clarify its type and period. We argue that all of the eclipsing binaries are field stars, while 5-6 of the LPVs are members of NGC 6496. We compare the period-luminosity distribution of these LPVs with those of LPVs in the Large Magellanic Cloud and 47 Tucanae, and with theoretical pulsation models. We also present a VI color magnitude diagram, display the evolutionary states of the variables, and match isochrones to determine a reddening of E(B-V) ~ 0.21 mag and apparent distance modulus of ~ 15.60 mag.

We provide APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

Light echoes, light from a variable source scattered off dust, have been observed for over a century. The recent discovery of light echoes around centuries-old supernovae in the Milky Way and the Large Magellanic Cloud have allowed the spectroscopic characterization of these events, even without contemporaneous photometry and spectroscopy using modern instrumentation. Here we review the recent scientific advances using light echoes of ancient and historic transients, and focus on our latest work on SN 1987A's and Eta Carinae's light echoes.

We present follow-up optical imaging and spectroscopy of one of the light echoes of $\eta$ Carinae's 19th-century Great Eruption discovered by Rest et al. (2012). By obtaining images and spectra at the same light echo position between 2011 and 2014, we follow the evolution of the Great Eruption on a three-year timescale. We find remarkable changes in the photometric and spectroscopic evolution of the echo light. The $i$-band light curve shows a decline of $\sim 0.9$ mag in $\sim 1$ year after the peak observed in early 2011 and a flattening at later times. The spectra show a pure-absorption early G-type stellar spectrum at peak, but a few months after peak the lines of the [Ca II] triplet develop strong P-Cygni profiles and we see the appearance of [Ca II] 7291,7324 doublet in emission. These emission features and their evolution in time resemble those observed in the spectra of some Type IIn supernovae and supernova impostors. Most surprisingly, starting $\sim 300$ days after peak brightness, the spectra show strong molecular transitions of CN at $\gtrsim 6800$ \AA. The appearance of these CN features can be explained if the ejecta are strongly Nitrogen enhanced, as it is observed in modern spectroscopic studies of the bipolar Homunculus nebula. Given the spectroscopic evolution of the light echo, velocities of the main features, and detection of strong CN, we are likely seeing ejecta that contributes directly to the Homunculus nebula.

We report the results of our study of lightcurve modulation in a sample of 493 RR1 variables from the OGLE-III survey of galactic bulge fields. Each object in this list has 1500 or more I-band observations. We compare our findings with earlier studies regarding lightcurve modulation in LMC and galactic field RR1 stars. We also report the discovery of the modulated-Blazhko RR1 star OGLE-BLG-RRLYR-03825 which has a Blazhko period of 16.469 d which itself is modulated with a period of 339.2 d.

We describe numerical simulations, using the particle-in-cell (PIC) and hybrid-PIC code Lsp [T. P. Hughes et al., Phys. Rev. ST Accel. Beams 2, 110401 (1999)], of the head-on merging of two laboratory supersonic plasma jets. The goals of these experiments are to form and study astrophysically relevant collisionless shocks in the laboratory. Using the plasma jet initial conditions (density ~ 10^14-10^16 cm^(-3), temperature ~ few eV, and propagation speed ~ 20-100 km/s), large-scale simulations of jet propagation demonstrate that interactions between the two jets are essentially collisionless at the merge region. In highly resolved one- and two-dimensional simulations, we show that collisionless shocks are generated by the merging jets when immersed in applied magnetic fields (B ~ 0.1-1 kG). At expected plasma jet speeds of up to 100 km/s, our simulations do not give rise to unmagnetized collisionless shocks, which require much higher velocities. The orientation of the magnetic field and the axial and transverse density gradients of the jets have a strong effect on the nature of the interaction. We compare some of our simulation results with those of previously published PIC simulation studies of collisionless shock formation.

We report direct evidence for asymmetry in the early phases of SN 1987A via optical spectroscopy of five fields of its light echo system. The light echoes allow the first few hundred days of the explosion to be reobserved, with different position angles providing different viewing angles to the supernova. Light echo spectroscopy therefore allows a direct spectroscopic comparison of light originating from different regions of the photosphere during the early phases of SN 1987A. Gemini multi-object spectroscopy of the light echo fields shows fine-structure in the H-alpha line as a smooth function of position angle on the near-circular light echo rings. H-alpha profiles originating from the northern hemisphere of SN 1987A show an excess in redshifted emission and a blue knee, while southern hemisphere profiles show an excess of blueshifted H-alpha emission and a red knee. This fine-structure is reminiscent of the "Bochum event" originally observed for SN 1987A, but in an exaggerated form. Maximum deviation from symmetry in the H-alpha line is observed at position angles 16 and 186 degrees, consistent with the major-axis of the expanding elongated ejecta. The asymmetry signature observed in the H-alpha line smoothly diminishes as a function of viewing angle away from the poles of the elongated ejecta. We propose an asymmetric two-sided distribution of 56Ni most dominant in the southern far quadrant of SN 1987A as the most probable explanation of the observed light echo spectra. This is evidence that the asymmetry of high-velocity 56Ni in the first few hundred days after explosion is correlated to the geometry of the ejecta some 25 years later.

R Coronae Borealis stars (RCBs) are rare, hydrogen-deficient, carbon-rich supergiant variable stars that are likely the evolved merger products of pairs of CO and He white dwarfs. Only 55 RCB stars are known in our galaxy and their distribution on the sky is weighted heavily by microlensing survey field positions. A less-biased wide-area survey would provide the ability to test competing evolutionary scenarios, understand the population or populations that produce RCBs and constraint their formation rate. The ASAS-3 survey monitored the sky south of declination +28 deg since 2000 to a limiting magnitude of V=14. We searched ASAS-3 for RCB variables using a number of different methods to ensure that the probability of RCB detection was as high as possible and to reduce selection biases based on luminosity, temperature, dust production activity and shell brightness. Candidates whose light curves were visually inspected were pre-selected based on their infrared excesses due to warm dust in their circumstellar shells using the WISE and/or 2MASS catalogues, and criteria on light curve variability. We then acquired spectra of 104 stars to determine their real nature using the SSO/WiFeS spectrograph. We report 21 newly-discovered RCB stars and 2 new DY Per stars. Two previously suspected RCB candidates were also spectroscopically confirmed. Our methods allowed us to extend our detection efficiency to fainter magnitudes that would not have been easily accessible to discovery techniques based only on light curve variability. The overall detection efficiency is about 90% for RCBs with maximum light brighter than V~13. This growing sample is of great value to constrain the peculiar and disparate atmosphere composition of RCBs. Most importantly, we show that the spatial distribution and apparent magnitudes of Galactic RCB stars is consistent with RCBs being part of the Galactic bulge population.

We present the first measurements of the near-infrared (NIR), specifically the J-band, sky background in the Canadian High Arctic. There has been considerable recent interest in the development of an astronomical observatory in Ellesmere Island; initial site testing has shown promise for a world-class site. Encouragement for our study came from sky background measurements on the high Antarctic glacial plateau in winter that showed markedly lower NIR emission when compared to good mid-latitude astronomical sites due to reduced emission from OH airglow lines. This is possibly a Polar effect and may also be present in the High Arctic. To test this hypothesis, we carried out an experiment which measured the the J-band sky brightness in the High Arctic during winter. We constructed a zenith-pointing, J-band photometer, and installed it at the Polar Environment Atmospheric Research Laboratory (PEARL) near Eureka, Nunavut (latitude: 80 degrees N). We present the design of our photometer and our results from our short PEARL observing campaign in February 2012. Taken over a period of four days, our measurements indicate that the J-band sky brightness varies between 15.5-15.9 mag arcsec^2, with a measurement uncertainty of 0.15 mag dominated entirely by systematic errors present in our radiometric calibration. On our best night, we measured a fairly consistent sky brightness of 15.8 +/- 0.15 mag arcsec^2. This is not corrected for atmospheric extinction, which is typically <0.1 mag in the J-band on a good night. The measured sky brightness is comparable to an excellent mid-latitude site, but is not as dark as claimed by the Antarctic measurements. We discuss possible explanations of why we do not see as dark skies as in the Antarctic. Future winter-long sky brightness measurements are anticipated to obtain the necessary statistics to make a proper comparison with the Antarctic measurements.

Astronomical light echoes, the time-dependent light scattered by dust in the vicinity of varying objects, have been recognized for over a century. Initially, their utility was thought to be confined to mapping out the three-dimensional distribution of interstellar dust. Recently, the discovery of spectroscopically-useful light echoes around centuries-old supernovae in the Milky Way and the Large Magellanic Cloud has opened up new scientific opportunities to exploit light echoes. In this review, we describe the history of light echoes in the local Universe and cover the many new developments in both the observation of light echoes and the interpretation of the light scattered from them. Among other benefits, we highlight our new ability to spectroscopically classify outbursting objects, to view them from multiple perspectives, to obtain a spectroscopic time series of the outburst, and to establish accurate distances to the source event. We also describe the broader range of variable objects whose properties may be better understood from light echo observations. Finally, we discuss the prospects of new light echo techniques not yet realized in practice.

The source counts of galaxies discovered at sub-millimetre and millimetre wavelengths provide important information on the evolution of infrared-bright galaxies. We combine the data from six blank-field surveys carried out at 1.1 mm with AzTEC, totalling 1.6 square degrees in area with root-mean-square depths ranging from 0.4 to 1.7 mJy, and derive the strongest constraints to date on the 1.1 mm source counts at flux densities S(1100) = 1-12 mJy. Using additional data from the AzTEC Cluster Environment Survey to extend the counts to S(1100) ~ 20 mJy, we see tentative evidence for an enhancement relative to the exponential drop in the counts at S(1100) ~ 13 mJy and a smooth connection to the bright source counts at >20 mJy measured by the South Pole Telescope; this excess may be due to strong lensing effects. We compare these counts to predictions from several semi-analytical and phenomenological models and find that for most the agreement is quite good at flux densities > 4 mJy; however, we find significant discrepancies (>3sigma) between the models and the observed 1.1 mm counts at lower flux densities, and none of them are consistent with the observed turnover in the Euclidean-normalised counts at S(1100) < 2 mJy. Our new results therefore may require modifications to existing evolutionary models for low luminosity galaxies. Alternatively, the discrepancy between the measured counts at the faint end and predictions from phenomenological models could arise from limited knowledge of the spectral energy distributions of faint galaxies in the local Universe.

Instruments using arrays of many bolometers have become increasingly common in the past decade. The maps produced by such instruments typically include the filtering effects of the instrument as well as those from subsequent steps performed in the reduction of the data. Therefore interpretation of the maps is dependent upon accurately calculating the transfer function of the chosen reduction technique on the signal of interest. Many of these instruments use non-linear and iterative techniques to reduce their data because such methods can offer improved signal-to-noise over those that are purely linear, particularly for signals at scales comparable to that subtended by the array. We discuss a general approach for measuring the transfer function of principal component analysis (PCA) on point sources that are small compared to the spatial extent seen by any single bolometer within the array. The results are applied to previously released AzTEC catalogues of the COSMOS, Lockman Hole, Subaru XMM-Newton Deep Field, GOODS-North and GOODS-South fields. Source flux density and noise estimates increase by roughly +10 per cent for fields observed while AzTEC was installed at the Atacama Submillimeter Telescope Experiment and +15-25 per cent while AzTEC was installed at the James Clerk Maxwell Telescope. Detection significance is, on average, unaffected by the revised technique. The revised photometry technique will be used in subsequent AzTEC releases.

Eta Carinae (Eta Car) is one of the most massive binary stars in the Milky Way. It became the second-brightest star in the sky during its mid-19th century "Great Eruption," but then faded from view (with only naked-eye estimates of brightness). Its eruption is unique among known astronomical transients in that it exceeded the Eddington luminosity limit for 10 years. Because it is only 2.3 kpc away, spatially resolved studies of the nebula have constrained the ejected mass and velocity, indicating that in its 19th century eruption, Eta Car ejected more than 10 M_solar in an event that had 10% of the energy of a typical core-collapse supernova without destroying the star. Here we report the discovery of light echoes of Eta Carinae which appear to be from the 1838-1858 Great Eruption. Spectra of these light echoes show only absorption lines, which are blueshifted by -210 km/s, in good agreement with predicted expansion speeds. The light-echo spectra correlate best with those of G2-G5 supergiant spectra, which have effective temperatures of ~5000 K. In contrast to the class of extragalactic outbursts assumed to be analogs of Eta Car's Great Eruption, the effective temperature of its outburst is significantly cooler than allowed by standard opaque wind models. This indicates that other physical mechanisms like an energetic blast wave may have triggered and influenced the eruption.

We found the first direct evidence that the Cepheid class namesake, delta Cephei, is currently losing mass. These observations are based on data obtained with the Spitzer Space Telescope in the infrared, and with the Very Large Array in the radio. We found that delta Cephei is associated with a vast circumstellar structure, reminiscent of a bow shock. This structure is created as the wind from the star interacts with the local interstellar medium. We measure an outflow velocity of ~35.5 km/s and a mass loss rate of ~1E-7 - 1E-6 Mo/yr. The very low dust content of the outflow suggests that the wind is possibly pulsation-driven, rather than dust-driven as common for other classes of evolved stars.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies.

NSV 11154 has been confirmed as a new member of the rare hydrogen deficient R Coronae Borealis (RCB) stars based on new photometric and spectroscopic data. Using new photometry, as well as archival plates from the Harvard archive, we have constructed the historical lightcurve of NSV 11154 from 1896 to the present. The lightcurve shows the sudden, deep, irregularly spaced declines characteristic of RCB stars. The visible spectrum is typical of a cool (Teff < 5000 K) RCB star showing no hydrogen lines, strong C2 Swan bands, and no evidence of 13C. In addition, the star shows small pulsations typical of an RCB star, and an infrared excess due to circumstellar dust with a temperature of ~800 K. The distance to NSV 11154 is estimated to be ~14.5 kpc. RCB stars are very rare in the Galaxy so each additional star is important to population studies leading to a better understanding the origins of these mysterious stars. Among the known sample of RCB stars, NSV 11154 is unusual in that it lies well above the Galactic plane (5 kpc) and away from the Galactic Center which suggests that its parent population is neither thick disk nor bulge.

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

The light echo systems of historical supernovae in the Milky Way and local group galaxies provide an unprecedented opportunity to reveal the effects of asymmetry on observables, particularly optical spectra. Scattering dust at different locations on the light echo ellipsoid witnesses the supernova from different perspectives and the light consequently scattered towards Earth preserves the shape of line profile variations introduced by asymmetries in the supernova photosphere. However, the interpretation of supernova light echo spectra to date has not involved a detailed consideration of the effects of outburst duration and geometrical scattering modifications due to finite scattering dust filament dimension, inclination, and image point-spread function and spectrograph slit width. In this paper, we explore the implications of these factors and present a framework for future resolved supernova light echo spectra interpretation, and test it against Cas A and SN 1987A light echo spectra. We conclude that the full modeling of the dimensions and orientation of the scattering dust using the observed light echoes at two or more epochs is critical for the correct interpretation of light echo spectra. Indeed, without doing so one might falsely conclude that differences exist when none are actually present.

We report the first detection of asymmetry in a supernova (SN) photosphere based on SN light echo (LE) spectra of Cas A from the different perspectives of dust concentrations on its LE ellipsoid. New LEs are reported based on difference images, and optical spectra of these LEs are analyzed and compared. After properly accounting for the effects of finite dust-filament extent and inclination, we find one field where the He I and H alpha features are blueshifted by an additional ~4000 km/s relative to other spectra and to the spectra of the Type IIb SN 1993J. That same direction does not show any shift relative to other Cas A LE spectra in the Ca II near-infrared triplet feature. We compare the perspectives of the Cas A LE dust concentrations with recent three-dimensional modeling of the SN remnant (SNR) and note that the location having the blueshifted He I and H alpha features is roughly in the direction of an Fe-rich outflow and in the opposite direction of the motion of the compact object at the center of the SNR. We conclude that Cas A was an intrinsically asymmetric SN. Future LE spectroscopy of this object, and of other historical SNe, will provide additional insight into the connection of explosion mechanism to SN to SNR, as well as give crucial observational evidence regarding how stars explode.

We present the discovery of an infrared nebula around the Cepheid prototype delta Cephei and its hot companion HD 213307. Large scale (~2.1x10^4 AU) nebulosity is detected at 5.8, 8.0, 24 and 70 um. Surrounding the two stars, the 5.8 and 8.0 um emission is largely attributable to Polycyclic Aromatic Hydrocarbon (PAH) emission swept from the ISM by a wind originating from delta Cephei and/or its companion. Stochastically heated small dust grains are the most likely source of the 24 and 70 um extended emission. The 70 um emission, in particular, resembles a bow shock aligned in the direction of the proper motion of delta Cephei. This discovery supports the hypothesis that delta Cephei may be currently losing mass, at a rate in the range ~ 5x10^-9 to 6x10^-8 Mo/yr.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

A deep and detailed examination of 29 classical Cepheids with the Spitzer Space Telescope has revealed three stars with strong nearby extended emission detected in multiple bands which appears to be physically associated with the stars. RS Pup was already known to possess extended infrared emission, while the extended emission around the other two stars S Mus and delta Cep is newly discovered in our observations. Four other stars GH Lup, l Car, T Mon and X Cyg show tentative evidence for extended infrared emission. An unusual elongated extended object next to SZ Tau appears to be a background or foreground object in a chance alignment with the Cepheid. The inferred mass loss rates upper limits for S Mus and delta Cep are in the range from 1e-9 to 1e-8 Msun/yr, with the upper limit for RS Pup as high as 1e-6 Msun/yr. Mass loss during post-main-sequence evolution has been proposed as a resolution to the discrepancy between pulsational and dynamical masses of Cepheid variable stars: dust in the lost material would make itself known by the presence of an infrared bright nebula, or unresolved infrared excess. The observed frequency of infrared circumstellar emission (<24%) and the mass loss rate we estimate for our sources shows that dusty mass loss can only account for part of the Cepheid mass loss discrepancy. Nevertheless, our direct evidence that mass loss is active during the Cepheid phase is an important confirmation that these processes need to be included in evolutionary and pulsation models of these stars, and should be taken into account in the calibration of the Cepheid distance scale.

SN 2007od exhibits characteristics that have rarely been seen in a Type IIP supernova (SN). Optical V band photometry reveals a very steep brightness decline between the plateau and nebular phases of ~4.5 mag, likely due to SN 2007od containing a low mass of 56Ni. The optical spectra show an evolution from normal Type IIP with broad Halpha emission, to a complex, four component Halpha emission profile exhibiting asymmetries caused by dust extinction after day 232. This is similar to the spectral evolution of the Type IIn SN 1998S, although no early-time narrow (~200 km s-1) Halpha component was present in SN 2007od. In both SNe, the intermediate-width Halpha emission components are thought to arise in the interaction between the ejecta and its circumstellar medium (CSM). SN 2007od also shows a mid-IR excess due to new dust. The evolution of the Halpha profile and the presence of the mid-IR excess provide strong evidence that SN 2007od formed new dust before day 232. Late-time observations reveal a flattening of the visible lightcurve. This flattening is a strong indication of the presence of a light echo, which likely accounts for much of the broad, underlying Halpha component seen at late-times. We believe the multi-peaked Halpha emission is consistent with the interaction of the ejecta with a circumstellar ring or torus (for the inner components at \pm1500 km s-1), and a single blob or cloud of circumstellar material out of the plane of the CSM ring (for the outer component at -5000 km s-1). The most probable location for the formation of new dust is in the cool dense shell created by the interaction between the expanding ejecta and its CSM. Monte Carlo radiative transfer modeling of the dust emission from SN 2007od implies that up to 4x 10-4Msun of new dust has formed. This is similar to the amounts of dust formed in other CCSNe such as SNe 1999em, 2004et, and 2006jc.

We present 2323 High-Amplitude \delta-Scuti (HADS) candidates discovered in the Large Magellanic Cloud (LMC) by the SuperMACHO survey (Rest et al. 2005). Frequency analyses of these candidates reveal that several are multimode pulsators, including 119 whose largest amplitude of pulsation is in the fundamental (F) mode and 19 whose largest amplitude of pulsation is in the first overtone (FO) mode. Using Fourier decomposition of the HADS light curves, we find that the period-luminosity (PL) relation defined by the FO pulsators does not show a clear separation from the PL-relation defined by the F pulsators. This differs from other instability strip pulsators such as type c RR Lyrae. We also present evidence for a larger amplitude, subluminous population of HADS similar to that observed in Fornax (Poretti et al. 2008).

SN 2008S erupted in early 2008 in the grand design spiral galaxy NGC 6946. The progenitor was detected by Prieto et al. in Spitzer Space Telescope images taken over the four years prior to the explosion, but was not detected in deep optical images, from which they inferred a self-obscured object with a mass of about 10 Msun. We obtained Spitzer observations of SN 2008S five days after its discovery, as well as coordinated Gemini and Spitzer optical and infrared observations six months after its outburst. We have constructed radiative transfer dust models for the object before and after the outburst, using the same r^-2 density distribution of pre-existing amorphous carbon grains for all epochs and taking light-travel time effects into account for the early post-outburst epoch. We rule out silicate grains as a significant component of the dust around SN 2008S. The inner radius of the dust shell moved outwards from its pre-outburst value of 85 AU to a post-outburst value of 1250 AU, attributable to grain vaporisation by the light flash from SN 2008S. Although this caused the circumstellar extinction to decrease from Av = 15 before the outburst to 0.8 after the outburst, we estimate that less than 2% of the overall circumstellar dust mass was destroyed. The total mass-loss rate from the progenitor star is estimated to have been (0.5-1.0)x10^-4 Msun yr^-1. The derived dust mass-loss rate of 5x10^-7 Msun yr^-1 implies a total dust injection into the ISM of up to 0.01 Msun over the suggested duration of the self-obscured phase. We consider the potential contribution of objects like SN 2008S to the dust enrichment of galaxies.

Classical Cepheid variable stars have been important indicators of extragalactic distance and Galactic evolution for over a century. The Spitzer Space Telescope has opened the possibility of extending the study of Cepheids into the mid- and far-infrared, where interstellar extinction is reduced. We have obtained photometry from images of a sample of Galactic Cepheids with the IRAC and MIPS instruments on Spitzer. Here we present the first mid-infrared period-luminosity relations for Classical Cepheids in the Galaxy, and the first ever Cepheid period-luminosity relations at 24 and 70 um. We compare these relations with theoretical predictions, and with period-luminosity relations obtained in recent studies of the Large Magellanic Cloud. We find a significant period-color relation for the [3.6]-[8.0] IRAC color. Other mid-infrared colors for both Cepheids and non-variable supergiants are strongly affected by variable molecular spectral features, in particular deep CO absorption bands. We do not find strong evidence for mid-infrared excess caused by warm (~500 K) circumstellar dust. We discuss the possibility that recent detections with near-infrared interferometers of circumstellar shells around delta Cep, l Car, Polaris, Y Oph and RS Pup may be a signature of shocked gas emission in a dust-poor wind associated to pulsation-driven mass loss.

V371 Per was found to be a double-mode Cepheid with a fundamental mode period of 1.738 days, the shortest among Galactic beat Cepheids, and an unusually high period ratio of 0.731, while the other Galactic beat Cepheids have period ratios between 0.697 and 0.713. The latter suggests that the star has a metallicity [Fe/H] between -1 and -0.7. The derived distance from the Galactic Plane places it in the Thick Disk or the Halo, while all other Galactic beat Cepheids belong to the Thin Disk. There are indications from historical data that both the fundamental and first overtone periods have lengthened.

V532 Oph has been found to be a member of the rare, hydrogen-deficient R Coronae Borealis (RCB) stars from new photometric and spectroscopic data reported in this paper. The lightcurve of V532 Oph shows the sudden, deep, irregularly spaced declines characteristic of RCB stars. Its optical spectrum is typical of a warm (T(eff)~7000 K) RCB star, showing weak or absent hydrogen lines, the C2 Swan bands, and no evidence for 13C. In addition, the star shows small pulsations typical of an RCB star and an infrared excess due to circum- stellar dust. It also appears to be significantly reddened by foreground dust. The distance to V532 Oph is estimated to be 5.5-8.7 kpc. These new data show that this star was misclassified as an eclipsing binary in the General Catalog of Variable Stars. The new data presented here for V532 Oph reveal the power of high-quality, high-cadence all-sky photometric surveys, such as ASAS-3, to identify new RCB candidates on the basis of lightcurve data alone, now that they have been collecting data for durations sufficiently long to reveal multiple declines. Despite their small numbers, RCB stars may be of great importance in understanding the late stages of stellar evolution. In particular, their measured isotopic abundances imply that many, if not most, RCB stars are produced by WD mergers, which may be the low-mass counterparts of the more massive mergers thought to produce type Ia supernovae. Therefore, establishing the population of RCB stars in the Galaxy will help constrain the frequency of these WD mergers.

Due to the strong effect of systematics/trends in variable star observations, we employ the Trend Filtering Algorithm (TFA) on a subset of the MACHO database and search for variable stars. TFA has been applied successfully in planetary transit searches, where weak, short-lasting periodic dimmings are sought in the presence of noise and various systematics (due to, e.g., imperfect flat fielding, crowding, etc). These latter effects introduce colored noise in the photometric time series that can lead to a complete miss of the signal. By using a large number of available photometric time series of a given field, TFA utilizes the fact that the same types of systematics appear in several/many time series of the same field. As a result, we fit each target time series by a (least-square-sense) optimum linear combination of templates and frequency-analyze the residuals. Once a signal is found, we reconstruct the signal by employing the full model, including the signal, systematics and noise. We apply TFA on the brightest ~5300 objects from subsets of each of the MACHO Large Magellanic Cloud fields #1 and #79. We find that the Fourier frequency analysis performed on the original data detect some 60% of the objects as trend-dominated. This figure decreases essentially to zero after using TFA. Altogether, We detect 387 variables in the two fields, 183 of which would have remained undetected without using TFA. Where possible, we give preliminary classification of the variables found.

A novel method of generating artificial scintillation noise is developed and used to evaluate occultation rates and false positive rates for surveys probing the Kuiper Belt with the method of serendipitous stellar occultations. A thorough examination of survey design shows that: (1) diffraction-dominated occultations are critically (Nyquist) sampled at a rate of 2 Fsu^{-1}, corresponding to 40 s^{-1} for objects at 40 AU, (2) occultation detection rates are maximized when targets are observed at solar opposition, (3) Main Belt Asteroids will produce occultations lightcurves identical to those of Kuiper Belt Objects if target stars are observed at solar elongations of: 116 deg < epsilon < 125 deg, or 131 deg < epsilon < 141 deg, and (4) genuine KBO occultations are likely to be so rare that a detection threshold of >7-8 sigma should be adopted to ensure that viable candidate events can be disentangled from false positives.

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 report the discovery of an extensive system of scattered light echo arclets associated with the recent supernovae in the local neighbourhood of the Milky Way: Tycho (SN 1572) and Cassiopeia A. Existing work suggests that the Tycho SN was a thermonuclear explosion while the Cas A supernova was a core collapse explosion. Precise classifications according to modern nomenclature require spectra of the outburst light. In the case of ancient SNe, this can only be done with spectroscopy of their light echo, where the discovery of the light echoes from the outburst light is the first step. Adjacent light echo positions suggest that Cas A and Tycho may share common scattering dust structures. If so, it is possible to measure precise distances between historical Galactic supernovae. On-going surveys that alert on the development of bright scattered-light echo features have the potential to reveal detailed spectroscopic information for many recent Galactic supernovae, both directly visible and obscured by dust in the Galactic plane.

We report the successful identification of the type of the supernova responsible for the supernova remnant SNR 0509-675 in the Large Magellanic Cloud (LMC) using Gemini spectra of surrounding light echoes. The ability to classify outbursts associated with centuries-old remnants provides a new window into several aspects of supernova research and is likely to be successful in providing new constraints on additional LMC supernovae as well as their historical counterparts in the Milky Way Galaxy (MWG). The combined spectrum of echo light from SNR 0509-675 shows broad emission and absorption lines consistent with a supernova (SN) spectrum. We create a spectral library consisting of 26 SNe Ia and 6 SN Ib/c that are time-integrated, dust-scattered by LMC dust, and reddened by the LMC and MWG. We fit these SN templates to the observed light echo spectrum using $\chi^2$ minimization as well as correlation techniques, and we find that overluminous 91T-like SNe Ia with $\dm15<0.9$ match the observed spectrum best.

The results of a search for sub-km Kuiper Belt Objects (KBOs) with the method of serendipitous stellar occultations are reported. Photometric time series were obtained on the 1.8m telescope at the Dominion Astrophysical Observatory (DAO) in Victoria, BC, and were analyzed for the presence of occultation events. Observations were performed at 40 Hz and included a total of 5.0 star-hours for target stars in the ecliptic open cluster M35 (beta=0.9deg), and 2.1 star-hours for control stars in the off-ecliptic open cluster M34 (beta=25.7deg). To evaluate the recovery fraction of the analysis method, and thereby determine the limiting detectable size, artificial occultation events were added to simulated time series (1/f scintillation-like power-spectra), and to the real data. No viable candidate occultation events were detected. This limits the cumulative surface density of KBOs to 3.5e10 deg^{-2} (95% confidence) for KBOs brighter than m_R=35.3 (larger than ~860m in diameter, assuming a geometric albedo of 0.04 and a distance of 40 AU). An evaluation of TNO occultations reported in the literature suggests that they are unlikely to be genuine, and an overall 95%-confidence upper limit on the surface density of 2.8e9 deg^{-2} is obtained for KBOs brighter than m_R=35 (larger than ~1 km in diameter, assuming a geometric albedo of 0.04 and a distance of 40 AU) when all existing surveys are combined.

This paper has been withdrawn by the author.

This paper has been withdrawn

We present new, very late-time optical photometry and spectroscopy of the interesting Type II-P supernova, SN 2002hh, in NGC 6946. Gemini/GMOS-N has been used to acquire visible spectra at six epochs between 2004 August and 2006 July, following the evolution of the SN from age 661 to 1358 days. Few optical spectra of Type II supernovae with ages greater than one year exist. In addition, g'r'i' images were acquired at all six epochs. The spectral and photometric evolution of SN 2002hh has been very unusual. Measures of the brightness of this SN, both in the R and I bands as well as in the H-alpha emission flux, show no significant fading over an interval of nearly two years. The most straightforward explanation for this behavior is that the light being measured comes not only from the SN itself but also from an echo off of nearby dust. Echoes have been detected previously around several SNe but these echoes, at their brightest, were ~8 mag below the maximum brightness of the SN. At V~21 mag, the putative echo dominates the light of SN 2002hh and is only ~4 mag below the outburst's peak brightness. There is an estimated 6 magnitudes of total extinction in V towards SN 2002hh. The proposed explanation of a differential echo/SN absorption is inconsistent with the observed BVRI colors.

We present a set of 11 type Ia supernova (SN Ia) lightcurves with dense, pre-maximum sampling. These supernovae (SNe), in galaxies behind the Large Magellanic Cloud (LMC), were discovered by the SuperMACHO survey. The SNe span a redshift range of z = 0.11 - 0.35. Our lightcurves contain some of the earliest pre-maximum observations of SNe Ia to date. We also give a functional model that describes the SN Ia lightcurve shape (in our VR-band). Our function uses the "expanding fireball" model of Goldhaber et al. (1998) to describe the rising lightcurve immediately after explosion but constrains it to smoothly join the remainder of the lightcurve. We fit this model to a composite observed VR-band lightcurve of three SNe between redshifts of 0.135 to 0.165. These SNe have not been K-corrected or adjusted to account for reddening. In this redshift range, the observed VR-band most closely matches the rest frame V-band. Using the best fit to our functional description of the lightcurve, we find the time between explosion and observed VR-band maximum to be 17.6+-1.3(stat)+-0.07(sys) rest-frame days for a SN Ia with a VR-band Delta m_{-10} of 0.52mag. For the redshifts sampled, the observed VR-band time-of-maximum brightness should be the same as the rest-frame V-band maximum to within 1.1 rest-frame days.

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

In principle, the light from historical supernovae could still be visible as scattered-light echoes even centuries later. However, while echoes have been discovered around some nearby extragalactic supernovae well after the explosion, targeted searches have not recovered any echoes in the regions of historical Galactic supernovae. The discovery of echoes can allow us to pinpoint the supernova event both in position and age and, most importantly, allow us to acquire spectra of the echo light to type the supernova centuries after the direct light from the explosion first reached the Earth. Here we report on the discovery of three faint new variable surface brightness complexes with high apparent proper motion pointing back to well-defined positions in the Large Magellanic Cloud (LMC). These positions correspond to three of the six smallest (and likely youngest) previously catalogued supernova remnants, and are believed to be due to thermonuclear (Type Ia) supernovae. Using the distance and proper motions of these echo arcs, we estimate ages of 610 and 410 yr for the echoes #2 and #3.

We present the initial results from an in-depth re-examination of the MACHO project Large Magellanic Cloud database to identify and characterize stellar variability near the intersection of the instability strip and the main sequence. This dataset's long time-series and uniform photometry is an unprecedented resource for describing the frequency and regions of incidence of various radial and non-radial modes of excitation. The raw MACHO photometry has been investigated to identify factors responsible for most of the residual photometric variance and increase the sensitivity to small amplitudes. We present details of the search method used, the Delta Scuti variables detected thus far and discuss the implications for the completed survey.

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

Characterizing the nature and spatial distribution of the lensing objects that produce the previously measured microlensing optical depth toward the Large Magellanic Cloud (LMC) remains an open problem. We present an appraisal of the ability of the SuperMACHO Project, a next-generation microlensing survey directed toward the LMC, to discriminate between various proposed lensing populations. We consider two scenarios: lensing by a uniform foreground screen of objects and self-lensing by LMC stars. We have carried out extensive simulations, based upon data obtained during the first year of the project, to assess the SuperMACHO survey's ability to discriminate between these two scenarios. We find that the event rate itself shows significant sensitivity to the choice of the LMC luminosity function, limiting the conclusions which can be drawn from the absolute rate. If instead we determine the differential event rate across the LMC, we will decrease the impact of these systematic biases and render our conclusions more robust. With this approach the SuperMACHO Project should be able to distinguish between the two categories of lens populations. This will provide important constraints on the nature of the lensing objects and their contributions to the Galactic dark matter halo.

(abridged) Recent studies, using OGLE data for LMC Cepheids in the optical, strongly suggest that the period-luminosity (PL) relation for the Large Magellanic Cloud (LMC) Cepheids shows a break or non-linearity at a period of 10 days. In this paper we apply statistical tests, the chi-square test and the F-test, to the Cepheid data from the MACHO project to test for a non-linearity of the V- and R-band PL relations at 10 days, and extend these tests to the near infrared (JHK-band) PL relations with 2MASS data. We correct the extinction for these data by applying an extinction map towards the LMC. The statistical test we use, the F-test, is able to take account of small numbers of data points and the nature of that data on either side of the period cut at 10 days. With our data, the results we obtained imply that the VRJH-band PL relations are non-linear around a period of 10 days, while the K-band PL relation is (marginally) consistent with a single-line regression. The choice of a period of 10 days, around which this non-linearity occurs, is consistent with the results obtained when this "break" period is estimated from the data. Long period Cepheids are supplemented from the literature to increase our sample size. The photometry of these long period Cepheids is compared with our data and no trend with period is found. Our main results remain unchanged when we supplement our dataset with these long period Cepheids. By examining our data at maximum light, we also suggest arguments why errors in reddening are unlikely to be responsible for our results. The non-linearity of the mean V-band PL relation as seen in both of the OGLE and MACHO data, using different extinction maps, suggests that this non-linearity is real.

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

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

We present the first results from our next-generation microlensing survey, the SuperMACHO project. We are using the CTIO 4m Blanco telescope and the MOSAIC imager to carry out a search for microlensing toward the Large Magellanic Cloud (LMC). We plan to ascertain the nature of the population responsible for the excess microlensing rate seen by the MACHO project. Our observing strategy is optimized to measure the differential microlensing rate across the face of the LMC. We find this derivative to be relatively insensitive to the details of the LMC's internal structure but a strong discriminant between Galactic halo and LMC self lensing. In December 2003 we completed our third year of survey operations. 2003 also marked the first year of real-time microlensing alerts and photometric and spectroscopic followup. We have extracted several dozen microlensing candidates, and we present some preliminary light curves and related information. Similar to the MACHO project, we find SNe behind the LMC to be a significant contaminant - this background has not been completely removed from our current single-color candidate sample. Our follow-up strategy is optimized to discriminate between SNe and true microlensing.

The lightcurve of the Large Magellanic Cloud (LMC) variable star MACHO 81.8997.87 shows evidence for photometric variations due to both stellar pulsation, with a 2.035 day period, and eclipsing behavior, with an 800.4 day period. The primary star of the system has been identified as a first-overtone Cepheid but the nature of the secondary star has not been determined. Here we present multicolor BVI photometry of a primary eclipse of the system and fit a model to the complete lightcurve to produce an updated set of elements. These results are combined with 2MASS JHK photometry to give further insight into the identity of the companion star. We find that the companion is most consistent with a late-K or an early-M giant but also that there are a number of problems with this interpretation. The prospects for future observations of this system are also discussed.

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

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

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

We present time-series VI photometry of the metal-rich globular cluster NGC 6316. Our color-magnitude diagrams show a predominantly red-clump horizontal branch morphology with no evidence of the blue extensions seen in NGC 6388 and NGC 6441. Our data are in good agreement with published estimates of the reddening and metallicity provided the cluster has an absolute distance modulus of 15.50 mag. We have discovered more than one dozen long-period variable stars in our field of view, and argue that at least seven of them are members of NGC 6316. We have also discovered four RR Lyrae variables, although only one of them probably belongs to the cluster. The specific frequency of RR Lyrae stars in NGC 6316 is found to be less than 0.4, making it clear that this cluster is a normal, metal-rich globular and is not an analog of NGC 6388 and NGC 6441. We recommend long-term photometric monitoring of NGC 6316 to clarify the nature of the RR Lyrae member candidate as well as the pulsation classes of the long-period variables.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The MACHO collaboration has been carrying out Difference Image Analysis (DIA) since 1996 with the aim of increasing the sensitivity to the detection of gravitational microlensing. This is a preliminary report on the application of DIA to galactic bulge images in one field. We show how the DIA technique significantly increases the number of detected lensing events, by removing the positional dependence of traditional photometry schemes and lowering the microlensing event detection threshold. This technique, unlike PSF photometry, gives the unblended colours and positions of the microlensing source stars. We present a set of criteria for selecting microlensing events from objects discovered with this technique. The 16 pixel and classical microlensing events discovered with the DIA technique are presented.

This is a preliminary report on the application of Difference Image Analysis (DIA) to galactic bulge images. The aim of this analysis is to increase the sensitivity to the detection of gravitational microlensing. We discuss how the DIA technique simplifies the process of discovering microlensing events by detecting only objects which have variable flux. We illustrate how the DIA technique is not limited to detection of so called ``pixel lensing'' events, but can also be used to improve photometry for classical microlensing events by removing the effects of blending. We will present a method whereby DIA can be used to reveal the true unblended colours, positions and light curves of microlensing events. We discuss the need for a technique to obtain the accurate microlensing time scales from blended sources, and present a possible solution to this problem using the existing HST colour magnitude diagrams of the galactic bulge and LMC. The use of such a solution with both classical and pixel microlensing searches is discussed. We show that one of the major causes of systematic noise in DIA is differential refraction. A technique for removing this systematic by effectively registering images to a common airmass is presented. Improvements to commonly used image differencing techniques are discussed.

In this paper, we will review major new results regarding classical Cepheids, in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). Specifically, we discuss recent work regarding multimode Cepheids and describe new observations of a W Vir star (HV 5756) and a Cepheid which are each in eclipsing binary systems. An additional interesting pulsating supergiant in an eclipsing system is also identified. Ephemerides for eclipses for the three systems are provided.

We present time-series VI photometry of the metal-rich ([Fe/H] = -0.53) globular cluster NGC 6441. Our color-magnitude diagram shows that the extended blue horizontal branch seen in Hubble Space Telescope data exists in the outermost reaches of the cluster. The red clump slopes nearly parallel to the reddening vector. A component of this slope is due to differential reddening, but part is intrinsic. The blue horizontal branch stars are more centrally concentrated than the red clump stars. We have discovered about 50 new variable stars near NGC 6441, among them eight or more RR Lyrae stars which are very probably cluster members. Comprehensive period searches over the range 0.2-1.0 days yielded unusually long periods (0.5-0.9 days) for the fundamental pulsators compared with field RR Lyrae of the same metallicity. Three similar long-period RR Lyrae are known in other metal-rich globulars. With over ten examples in hand, it seems that a distinct sub-class of RR Lyrae is emerging. The observed properties of the horizontal branch stars are in reasonable agreement with recent models which invoke deep mixing to enhance the atmospheric helium abundance, while they conflict with models which assume high initial helium abundance. The light curves of the c-type RR Lyrae seem to have unusually long rise times and sharp minima. Reproducing these light curves in stellar pulsation models may provide another means of constraining the physical variables responsible for the anomalous blue horizontal branch extension and sloped red clump observed in NGC 6441.

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

The influence of new techniques on the discovery and characterization of pulsating variables has been enormous. In this paper, I will review the methods and results of a number of research programmes which have dramatically altered our ability to study variable stars and will likely bear fruit for many years into the future. Specifically, I will touch on results from the MACHO and EROS Projects, the HIPPARCOS mission, "flux difference" photometry and a handful of new algorithmic advances which I consider to be important.

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

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

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

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

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

The discovery of a large number of beat Cepheids in the Large Magellanic Cloud in the MACHO survey, provides an opportunity to compare the characteristics of such Cepheids over a range of metallicities. We produced a large grid of linear nonadiabatic pulsation models using the OPAL opacity tables and with compositions corresponding to those of the Milky Way, and the Large and Small Magellanic Clouds. Using the relationship between the period ratio and the main pulsation period, we are able to define a range of models which correspond to the observed beat Cepheids, and thereby constrain the physical characteristics of the LMC beat Cepheids. We are also able to make some predictions about the nature of the yet-to-be-discovered SMC beat Cepheids.

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.

We have measured the distance to the Local Group dwarf irregular galaxy IC 10 from near-infrared $JHK$ observations of Cepheids. The reddening-corrected distance modulus is $24.57\pm0.21$ mag or $0.82\pm 0.08$ Mpc. This distance lies at the short end of the range of previous estimates, and is significantly more accurate than all prior determinations. At this distance, IC 10 is comparable in total mass to the Small Magellanic Cloud, but has much more vigorous star formation.

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.

Optical/infrared photometry can double the number of proper motion measurements of Massive Compact Objects (MACHOs) relative to single band photometry. The proper motion of a MACHO can be measured by finding the ratio $q$ of the (known) radius of the source star to the Einstein radius of the MACHO, $q=\theta_s/\theta_e$. A classic method for doing this is to look for the effect on the light curve of the finite size of the source. A modification of this method proposed by Witt (1995) is to look for color changes in the light curve due to the fact that the limb darkening of the source is different in different bands. We demonstrate that the ``classical'' method is not feasible unless the MACHO actually transits the source: if the MACHO passes at say 1.5 source radii, there is still a sizable $\sim 5\%$ effect, but the light curve cannot be distinguished from point-source light curves with different parameters. However, color measurements in $V$ $(0.55\,\mu$m) and $H$ $(1.65\,\mu$m) reduce the errors by a factor $\sim 120$ and permit proper motion measurements at impact parameters of up to 2 source radii. Color maps in $V-H$ are also useful in the detection of planetary systems. Giant stars have a ``red ring'' in such maps. A planet which transits this ring gives rise to a distinctive signature which can help in the measurement of the planetary system's proper motion.

B,V CCD frames have been used to derive surface brightness profiles for NGC 3201 out to ~18 arcmin. A total of 857 radial velocities with median precision ~1 km/s for 399 member giants have been used to trace the velocity dispersion profile out to 32' (the approximate tidal radius from fits of single-mass, isotropic King-Michie models to the cluster surface brightness profiles). The median difference in radial velocity for stars on either side of an imaginary axis moved through the cluster in 1 degree steps shows a significant maximum amplitude of 1.22+/-0.25 km/s. We discuss possible explanations of this result, including: (1) cluster rotation; (2) preferential stripping of stars on prograde orbits near the limiting radius; (3) the projection of the cluster space velocity onto the plane of the sky and (4) a slight drift in the velocity zero point. It is difficult to identify the primary cause of the observed velocity field structure unambiguously, and we suspect that all of the above processes may play a role. The B,V surface brightness profiles and radial velocities have been modeled with single- & multi-mass King-Michie models and nonparametric techniques. The density and M/L profiles show good agreement over 1.5<R<10 pc, and both approaches suggest a steady rise in M/L with distance from the cluster center. Due to the low cluster luminosity, we are unable to place useful constraints on the anisotropy of the velocity dispersion profile, though the global mass-to-light ratio is well-constrained by the models as ~2.0 +/-0.2 for the multi-mass and nonparametric models, compared to ~ 1.65 +/-0.15 for models having equal-mass stars. Our best-fit, multi-mass models have mass function slopes of x~0.75 +/-0.25, consistent with findings that mass function depends on the position relative to the potential of the Galaxy.

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 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 advantages of a near-infrared variant of the Barnes-Evans method for estimating distances to Cepheid variables are described and quantified. A surface brightness-color relation for $K$ photometry and the $(V-K)_0$ color index is established using modern, high-precision angular diameters from optical interferometers. Applied to data for the galactic (cluster) Cepheid U Sgr, this method yields a distance of 0.660 $\pm$ 0.024 kpc and a true distance modulus of 9.10 $\pm$ 0.07 mag. This estimate compares with the true distance modulus of 9.37 $\pm$ 0.22 mag estimated by Gieren, Barnes, and Moffett (1993) using the classical Barnes-Evans technique. The possibility of estimating distances of LMC and SMC Cepheids directly -- without intermediate steps -- is discussed. The feasibility of determining the distance of M31 or M33 using this technique is examined and is probably within the reach of 8m-class telescopes.

In this paper we have examined the internal dynamics of the globular cluster NGC 362. A V band surface brightness profile (SBP) was constructed from CCD images, and, after it was determined that the cluster is not post core-collapse, fit with single- and multi-mass King-Michie (KM) models. The total cluster luminosity is 1.70 $\pm 0.1 \times 10^5$ L$_V$\sol. A total of 285 stellar spectra were obtained of 215 stars for radial velocity determinations. Four stars showed strong evidence for radial velocity variations and are probably members of binary systems. The true cluster binary fraction was determined from simulations to be 0.15 for circular orbits or 0.27 for orbits with a distribution function $f(e) = e$ ($e$ is eccentricity). This relatively high binary detection frequency may indicate that NGC 362 is overabundant in binaries compared to other clusters. The 208 remaining stars showed no sign of rotation. The best agreement with both the kinematic data and the SBP were for shallow mass functions $x = 0.0 - 0.5$ and intermediate amounts of anisotropy in the velocity dispersion tensor. The cluster mass is M = $2.5 - 3.5 \times 10^5$ M\sol\ for a global mass-to-light ratio of M/L$_V$ = 1.5 - 2.0 M\sol/L$_V$\sol. This low value for $x$ is in disagreement with the correlation between $x$ and the height above the Galactic disk seen for a sample of other clusters. The results are also different from the sharp turn-up in the low mass end of mass functions derived from some deep luminosity functions of three other globular clusters.

In this paper we have examined the age and internal dynamics of the young binary LMC cluster NGC 1850 using BV CCD images and echelle spectra of 52 supergiants. Isochrone fits to a BV color-magnitude diagram revealed that the primary cluster has an age of $\tau = 90 \pm 30$ Myr while the secondary member has $\tau = 6 \pm 5$ Myr. BV surface brightness profiles were constructed out to R $>$ 40 pc, and single-component King-Michie (KM) models were applied. The total cluster luminosity varied from L$_B$ = 2.60 - 2.65 $\times 10^6$ L$_B$\sol\ and L$_V$ = 1.25 - 1.35 $\times 10^6$ as the anisotropy radius varied from infinity to three times the scale radius with the isotropic models providing the best agreement with the data. Of the 52 stars with echelle spectra, a subset of 36 were used to study the cluster dynamics. The KM radial velocity distributions were fitted to these velocities yielding total cluster masses of 5.4 - 5.9 $\pm 2.4 \times 10^4$ M\sol\ corresponding to M/L$_B$ = 0.02 $\pm 0.01$ M\sol/L$_B$\sol\ or M/L$_V$ = 0.05 $\pm 0.02$ M\sol/L$_V$\sol. A rotational signal in the radial velocities has been detected at the 93\% confidence level implying a rotation axis at a position angle of 100\deg. A variety of rotating models were fit to the velocity data assuming cluster ellipticities of $\epsilon = 0.1 - 0.3$. These models provided slightly better agreement with the radial velocity data than the KM models and had masses that were systematically lower by a few percent. The preferred value for the slope of a power-law IMF is a relatively shallow, $x = 0.29 \pmm{+0.3}{-0.8}$ assuming the B-band M/L or $x = 0.71 \pmm{+0.2}{-0.4}$ for the V-band.

In this paper we investigate the internal dynamics of the LMC cluster NGC 1978 through the use of Photometric (CCD images) and kinematic (stellar radial velocities) data. We apply a variety of dynamical models to this data, including multi-mass King-Michie models and rotating and non-rotating oblate spheroid models. We discuss the cluster mass-to-light ratio and place constraints on the cluster mass function.