Search

We present a sample of ~5000 RR Lyrae stars selected from the recalibrated LINEAR data set and detected at heliocentric distances between 5 kpc and 30 kpc over ~8000 deg[superscript 2] of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We analyze in detail the light curve properties and Galactic distribution of the subset of ~4000 type ab RR Lyrae (RRab) stars, including a search for new halo substructures and the number density distribution as a function of Oosterhoff type. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4, but with a weaker separation than for globular cluster stars. The wide sky coverage and depth of this sample allow unique constraints for the number density distribution of halo RRab stars as a function of galactocentric distance: it can be described as an oblate ellipsoid with an axis ratio q = 0.63 and with either a single or a double power law with a power-law index in the range –2 to –3. Consistent with previous studies, we find that the Oosterhoff type II subsample has a steeper number density profile than the Oosterhoff type I subsample. Using the group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters are suggestive of tidal streams possibly originating from globular clusters. Spectroscopic follow-up of detected halo groups is encouraged., United States. Air Force (Contract FA8721-05-C-0002)

We describe photometric recalibration of data obtained by the asteroid survey LINEAR. Although LINEAR was designed for astrometric discovery of moving objects, the data set described here contains over 5 billion photometric measurements for about 25 million objects, mostly stars. We use Sloan Digital Sky Survey (SDSS) data from the overlapping ~10,000 deg[superscript 2] of sky to recalibrate LINEAR photometry and achieve errors of 0.03 mag for sources not limited by photon statistics with errors of 0.2 mag at r ~ 18. With its 200 observations per object on average, LINEAR data provide time domain information for the brightest four magnitudes of the SDSS survey. At the same time, LINEAR extends the deepest similar wide-area variability survey, the Northern Sky Variability Survey, by 3 mag. We briefly discuss the properties of about 7000 visually confirmed periodic variables, dominated by roughly equal fractions of RR Lyrae stars and eclipsing binary stars, and analyze their distribution in optical and infrared color-color diagrams. The LINEAR data set is publicly available from the SkyDOT Web site.

We present a large-scale effort of creating composite spectra of high-redshift SNe la and comparing them to low-redshift counterparts in an attempt to understand possible cosmic evolution of SNe la, which has major implications for studies of dark energy. Through the ESSENCE project, we have obtained 107 spectra of 88 high-redshift SNe la with excellent light-curve information. In addition, we have obtained 397 spectra of low-redshift SNe la through a multiple-decade effort at the Lick and Keck Observatories, and we have used 45 UV spectra obtained by HST and IUE. The low-redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the Lick composite spectra allows an excellent match of the overall SED with the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with host galaxy light than their low-redshift counterparts. This is caused by observing objects at all redshifts with similar angular slit widths, which corresponds to different projected physical distances. After correcting for the galaxy light contamination, a few marginally significant differences in the spectra remain. We have estimated the systematic errors when using current spectral templates for K-corrections to be ∼0.02 mag. The variance in the composite spectra gives an estimate of the intrinsic variance in low-redshift maximum light SN spectra of ∼3% relative flux in the optical and growing toward the UV. The difference between the maximum light low-and high-redshift spectra constrains the evolution of SN spectral features between our samples to be <10% relative flux in the rest-frame optical. Currently, galaxy contamination and the small samples of rest-frame UV spectra at low and high redshifts are the limiting factors for future studies.

We present multiepoch spectra of 13 high-redshift Type Ia supernovae (SNe Ia) drawn from the literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin and Tonry to determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time yields an apparent aging rate consistent with the 1 /(1 + z) factor (where z is the redshifi) expected in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's "tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The best-fit exponent for these models is consistent with the expected 1/(1 + z) factor.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the dark energy equation-of-state parameter, w = P/(pc 2). We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to w for a given fixed amount of telescope time. For our survey on the CTIO 4 m telescope, measuring the luminosity distances and redshifts for Supernovae at modest redshifts (z ≈ 0.5± 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find Supernovae automatically and in nearly real time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their apparent brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4 m telescope's natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 Type Ia supernovae, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for all of the Type Ia Supernovae found by ESSENCE and used in our measurement of w, presented in a companion paper by Wood-Vasey and coworkers.

The first cosmological results from the ESSENCE supernova survey (Wood-Vasey and coworkers) are extended to a wider range of cosmological models including dynamical dark energy and nonstandard cosmological models. We fold in a greater number of external data sets such as the recent Higher-z release of high-redshift supernovae (Riess and coworkers), as well as several complementary cosmological probes. Model comparison statistics such as the Bayesian and Akaike information criteria are applied to gauge the worth of models. These statistics favor models that give a good fit with fewer parameters. Based on this analysis, the preferred cosmological model is the flat cosmological constant model, where the expansion history of the universe can be adequately described with only one free parameter describing the energy content of the universe. Among the more exotic models that provide good fits to the data, we note a preference for models whose best-fit parameters reduce them to the cosmological constant model.

We present a catalog of 450 relatively high signal-to-noise ratio microlensing events observed by the MACHO collaboration between 1993 and 1999. The events are distributed throughout our fields, and as expected, they show a clear concentration toward the

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. The optical depth for "screen lensing" is essentially constant across the face of the LMC, whereas the optical depth for self-lensing shows a strong spatial dependence. We have carried out extensive simulations, based on data obtained during the first year of the project, to assess the SuperMACHO survey's ability to discriminate between these two scenarios. In our simulations we predict the expected number of observed microlensing events for various LMC models for each of our fields by adding artificial stars to the images and estimating the spatial and temporal efficiency of detecting microlensing events using Monte Carlo methods. We find that the event rate itself shows significant sensitivity to the choice of the LMC luminosity function, limiting the conclusions that 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.

We have frequency-analyzed 6391 variables classified earlier as fundamental-mode RR Lyrae (RRO) stars in the MACHO database on the Large Magellanic Cloud (LMC). The overwhelming majority (i.e., 96%) of these variables have been proved to be indeed RR0 stars, whereas the remaining ones have fallen into one of the following categories: single- and double-mode Cepheids, binaries, first-overtone and double-mode RR Lyrae stars, and nonclassified variables. Special attention has been paid to the properties of the amplitude- and phase-modulated RR0 stars (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. This sample overwhelmingly exceeds the number of Blazhko stars known in all other systems combined. The incidence rate of the Blazhko variables among the RR0 stars in the LMC is 11.9%, which is 3 times higher than their rate among the first-overtone RR Lyrae stars. No difference is found in the average brightness between the single-mode and Blazhko variables. However, the latter ones show a somewhat lower degree of skewness in their average light curves and a concomitant lower total amplitude in their modulation-free light curves. From the frequency spectra we found that variables with larger modulation amplitudes at the higher frequency side of the main pulsation component are 3 times more common than the ones showing the opposite amplitude pattern. A search 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 has failed to clearly detect such components at the ≈0.004 mag level. This means that the amplitudes of the quadruple components (if they exist) should be, on average, at least 10 times smaller than t

We present a method for solving the light curve of an eclipsing binary system that 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 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 Virginis 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 an analysis of the longest timescale microlensing events discovered by the MACHO Collaboration during a 7 year survey of the Galactic bulge. We find six 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 distances and masses of the lens objects based on a standard model of the Galactic velocity distribution. This analysis indicates that the most likely masses of five of the six lenses are greater than 1 M⊙, 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 / M⊙ = 6-3+10 and M / M⊙ = 6-3+7, respectively. The observed upper limits on the absolute brightness of main-sequence stars for these lenses are less than 1 L⊙, 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 six events may lie in the foreground Galactic disk or in the Sagittarius (Sgr) dwarf galaxy behind the bulge, but we find that bulge sources are likely to dominate our microlensing parallax event sample. Future Hubble Space Telescope 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 the Space Interferometry Mission or the Keck or VLT interferometers, as explained by Delplancke, Górski, & Richichi, will allow direct measurements of the lens masses for stellar remnant lenses as well.

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 (MV = - 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-MV 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 ∼ 3200. The SMC MACHO fields were also searched for RCB stars, but none were found.

We present the preliminary results of our astrometric study of stellar motions along the lines of sight of the Magellanic Clouds and the Galactic bulge. Using 5 years of MACHO project point-spread function photometry, we find that we can easily select stars with proper motions of 0″.03 yr-1 from these very dense fields, using the characteristic shapes of their light curves. By performing astrometry on photometrically selected, candidate, high proper motion (HPM) stars in 50 deg2, we have discovered 154 new HPM stars from ∼55 million stars monitored by the MACHO project. These new objects have proper motions as high as 0″.5 yr-1, brightnesses ranging from V ∼ 19 to V - 19, and V - R colors between 0.3 and 1.45.

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 nine 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 and first-overtone) variables, this study increased the number of known double-mode stars in the LMC to 181. During the total 6.5 yr time span of the data, 10% of the stars showed strong period changes. The size, and in general also the patterns of the period changes, exclude a simple evolutionary explanation. 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 are symmetric in their frequency spacing relative to the central component. This latter type of variables are 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 the Blazhko phenomenon by any theory depending mainly on the role of aspect angle or magnetic field. None of the current theoretical models are able to explain the observed close frequency components without invoking nonradial pulsation components in these stars.

We present the microlensing optical depth toward the Galactic bulge based on the detection of 99 events found in our Difference Image Analysis (DIA) survey. This analysis encompasses 3 yr of data, covering ̃17 million stars in ̃4 deg2, to a source-star

We present the light curves of 21 gravitational microlensing events from the first six years of the MACHO Project gravitational microlensing survey that are likely examples of lensing by binary systems. These events were manually selected from a total sample of ∼350 candidate microlensing events that were either detected by the MACHO Alert System or discovered through retrospective analyses of the MACHO database. At least 14 of these 21 events exhibit strong (caustic) features, and four of the events are well fit with lensing by large mass ratio (brown dwarf or planetary) systems, although these fits are not necessarily unique. The total binary event rate is roughly consistent with predictions based upon our knowledge of the properties of binary stars, but a precise comparison cannot be made without a determination of our binary lens event detection efficiency. Toward the Galactic bulge, we find a ratio of caustic crossing to noncaustic crossing binary lensing events of 12:4, excluding one event for which we present two fits. This suggests significant incompleteness in our ability to detect and characterize non-caustic crossing binary lensing. The distribution of mass ratios, N(q), for these binary lenses appears relatively flat. We are also able to reliably measure source-face crossing times in four of the bulge caustic crossing events, and recover from them a distribution of lens proper motions, masses, and distances consistent with a population of Galactic bulge lenses at a distance of 7 ± 1 kpc. This analysis yields two systems with companions of ∼0.05 M⊙.

MACHO Project photometry of 45 LMC first and second overtone (FO and SO, respectively) beat Cepheids that pulsate in the FO and SO has been analyzed to determine the light-curve characteristics for the SO mode of Cepheid pulsation. We predict that singly

We present photometric observations and analysis of the second microlensing event detected toward the Small Magellanic Cloud (SMC), MACHO Alert 98-SMC-1. This event was detected early enough to allow intensive observation of the light curve. These observations revealed 98-SMC-1 to be the first caustic crossing binary microlensing event toward the Magellanic Clouds to be discovered in progress. Frequent coverage of the evolving light curve allowed an accurate prediction for the date of the source crossing out of the lens caustic structure. The caustic crossing temporal width, along with the angular size of the source star, measures the proper motion of the lens with respect to the source and thus allows an estimate of the location of the lens. Lenses located in the Galactic halo would have a velocity projected to the SMC of v̂ ∼ 1500 kms-1, while an SMC lens would typically have v̂ ∼ 60 kms-1. The event light curve allows us to obtain a unique fit to the parameters of the binary lens and to estimate the proper motion of the lensing system. We have performed a joint fit to the MACHO/GMAN data presented here, including recent EROS data of this event from Afonso and collaborators. These joint data are sufficient to constrain the time t* for the lens to move an angle equal to the source angular radius: t* = 0.116 ± 0.010 days. We estimate a radius for the lensed source of R* = 1.1 ± 0.1 R⊙ from its unblended color and magnitude. This yields a projected velocity of v̂ = 76 ± 10 kms-1. Only 0.12% of halo lenses would be expected to have a v̂ value at least as small as this, while 38% of SMC lenses would be expected to have v̂ as large as this. This implies that the lensing system is more likely to reside in the SMC than in the Galactic halo. Similar observations of future Magellanic Cloud microlensing events will help to determine the contribution of MACHOS to the Galaxy's dark halo.

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 that 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 colors, positions, and light curves of microlensing events. We discuss the need for a technique to obtain the accurate microlensing timescales from blended sources and present a possible solution to this problem using the existing Hubble Space Telescope color-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 air mass is presented. Improvements to commonly used image differencing techniques are discussed.

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