In this catalog we compile information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) as well as all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae found in 2017, totaling 308 supernovae. We also present UV through near-IR magnitudes gathered from public databases of all host galaxies for the supernovae in the sample. We perform statistical analyses of our full bright supernova sample, which now contains 949 supernovae discovered since 2014 May 1, including supernovae from our previous catalogs. This is the fourth of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team, and this work presents updated data and measurements, including light curves, redshifts, classifications, and host galaxy identifications, that supersede information contained in any previous publications.

Hydrogen-poor superluminous supernovae (SLSN-I) are a class of rare and energetic explosions discovered in untargeted transient surveys in the past decade. The progenitor stars and the physical mechanism behind their large radiated energies ($\sim10^{51}$ erg) are both debated, with one class of models primarily requiring a large rotational energy, while the other requires very massive progenitors to either convert kinetic energy into radiation via interaction with circumstellar material (CSM), or engender a pair-instability explosion. Observing the structure of the CSM around SLSN-I offers a powerful test of some scenarios, though direct observations are scarce. Here, we present a series of spectroscopic observations of the SLSN-I iPTF16eh, which reveal both absorption and time- and frequency-variable emission in the Mg II resonance doublet. We show that these observations are naturally explained as a resonance scattering light echo from a circumstellar shell. Modeling the evolution of the emission, we find a shell radius of 0.1 pc and velocity of 3300 km s$^{-1}$, implying the shell was ejected three decades prior to the supernova explosion. These properties match theoretical predictions of pulsational pair-instability shell ejections, and imply the progenitor had a He core mass of $\sim 50-55~{\rm M}_{\odot}$, corresponding to an initial mass of $\sim 115~{\rm M}_{\odot}$.

Wide-field surveys are discovering a growing number of rare transients whose physical origin is not yet well understood. Here, we present optical and UV data and analysis of iPTF16asu, a luminous, rapidly-evolving, high velocity, stripped-envelope supernova. With a rest-frame rise-time of just 4 days and a peak absolute magnitude of $M_{\rm g}=-20.4$ mag, the light curve of iPTF16asu is faster and more luminous than previous rapid transients. The spectra of iPTF16asu show a featureless, blue continuum near peak that develops into a Type Ic-BL spectrum on the decline. We show that while the late-time light curve could plausibly be powered by $^{56}$Ni decay, the early emission requires a different energy source. Non-detections in the X-ray and radio strongly constrain any associated gamma-ray burst to be low-luminosity. We suggest that the early emission may have been powered by either a rapidly spinning-down magnetar, or by shock breakout in an extended envelope of a very energetic explosion. In either scenario a central engine is required, making iPTF16asu an intriguing transition object between superluminous supernovae, Type Ic-BL supernovae, and low-energy gamma-ray bursts.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

We present ground-based and \textit{Swift} observations of iPTF16fnl, a likely tidal disruption event (TDE) discovered by the intermediate Palomar Transient Factory (iPTF) survey at 66.6 Mpc. The lightcurve of the object peaked at absolute $M_g=-17.2$ mag. The maximum bolometric luminosity (from optical and UV) was $L_p~\simeq~(1.0\,\pm\,0.15) \times 10^{43}$ erg/s, an order of magnitude fainter than any other optical TDE discovered so far. The luminosity in the first 60 days is consistent with an exponential decay, with $L \propto e^{-(t-t_0)/\tau}$, where $t_0$=~57631.0 (MJD) and $\tau\simeq 15$ days. The X-ray shows a marginal detection at $L_X=2.4^{1.9}_{-1.1}\times 10^{39}$ erg/s (\textit{Swift} X-ray Telescope). No radio counterpart was detected down to 3$\sigma$, providing upper limits for monochromatic radio luminosity of $\nu L_{\nu} < 2.3\times10^{36}$ erg/s and $\nu L_{\nu}<1.7\times 10^{37}$ erg/s (VLA, 6.1 and 22 GHz). The blackbody temperature, obtained from combined \textit{Swift} UV and optical photometry, shows a constant value of 19,000 K. The transient spectrum at peak is characterized by broad He II and H$\alpha$ emission lines, with an FWHM of about 14,000 km/s and 10,000 km/s respectively. He I lines are also detected at $\lambda\lambda$ 5875 and 6678. The spectrum of the host is dominated by strong Balmer absorption lines, which are consistent with a post-starburst (E+A) galaxy with an age of $\sim$650 Myr and solar metallicity. The characteristics of iPTF16fnl make it an outlier on both luminosity and decay timescales, as compared to other optically selected TDEs. The discovery of such a faint optical event suggests a higher rate of tidal disruptions, as low luminosity events may have gone unnoticed in previous searches.

We report the discovery by the intermediate Palomar Transient Factory (iPTF) of a candidate tidal disruption event (TDE) iPTF16axa at $z=0.108$, and present its broadband photometric and spectroscopic evolution from 3 months of follow-up observations with ground-based telescopes and Swift. The light curve is well fitted with a $t^{-5/3}$ decay, and we constrain the rise-time to peak to be $<$49 rest-frame days after disruption, which is roughly consistent with the fallback timescale expected for the $\sim 5\times$10$^{6}$ $M_\odot$ black hole inferred from the stellar velocity dispersion of the host galaxy. The UV and optical spectral energy distribution (SED) is well described by a constant blackbody temperature of T$\sim$ 3$\times$10$^4$ K over the monitoring period, with an observed peak luminosity of 1.1$\times$10$^{44}$ erg s$^{-1}$. The optical spectra are characterized by a strong blue continuum and broad HeII and H$\alpha$ lines characteristic of TDEs. We compare the photometric and spectroscopic signatures of iPTF16axa with 11 TDE candidates in the literature with well-sampled optical light curves. Based on a single-temperature fit to the optical and near-UV photometry, most of these TDE candidates have peak luminosities confined between log(L [erg s$^{-1}$]) = 43.4-44.4, with constant temperatures of a few $\times 10^{4}$ K during their power-law declines, implying blackbody radii on the order of ten times the tidal disruption radius, that decrease monotonically with time. For TDE candidates with hydrogen and helium emission, the high helium-to-hydrogen ratios suggest that the emission arises from high-density gas, where nebular arguments break down. We find no correlation between the peak luminosity and the black hole mass, contrary to the expectations for TDEs to have $\dot{M} \propto M_{\rm BH}^{-1/2}$.

With the advent of new wide-field, high-cadence optical transient surveys, our understanding of the diversity of core-collapse supernovae has grown tremendously in the last decade. However, the pre-supernova evolution of massive stars, that sets the physical backdrop to these violent events, is theoretically not well understood and difficult to probe observationally. Here we report the discovery of the supernova iPTF13dqy = SN 2013fs, a mere ~3 hr after explosion. Our rapid follow-up observations, which include multiwavelength photometry and extremely early (beginning at ~6 hr post-explosion) spectra, map the distribution of material in the immediate environment (<~ 10^15 cm) of the exploding star and establish that it was surrounded by circumstellar material (CSM) that was ejected during the final ~1 yr prior to explosion at a high rate, around 10^-3 solar masses per year. The complete disappearance of flash-ionised emission lines within the first several days requires that the dense CSM be confined to within <~ 10^15 cm, consistent with radio non-detections at 70--100 days. The observations indicate that iPTF13dqy was a regular Type II SN; thus, the finding that the probable red supergiant (RSG) progenitor of this common explosion ejected material at a highly elevated rate just prior to its demise suggests that pre-supernova instabilities may be common among exploding massive stars.

We present a radio-quiet quasar at z=0.237 discovered "turning on" by the intermediate Palomar Transient Factory (iPTF). The transient, iPTF 16bco, was detected by iPTF in the nucleus of a galaxy with an archival SDSS spectrum with weak narrow-line emission characteristic of a low-ionization emission line region (LINER). Our follow-up spectra show the dramatic appearance of broad Balmer lines and a power-law continuum characteristic of a luminous (L_bol~10^45 erg/s) type 1 quasar 12 years later. Our photometric monitoring with PTF from 2009-2012, and serendipitous X-ray observations from the XMM-Newton Slew Survey in 2011 and 2015, constrain the change of state to have occurred less than 500 days before the iPTF detection. An enhanced broad Halpha to [OIII]5007 line ratio in the type 1 state relative to other changing-look quasars also is suggestive of the most rapid change of state yet observed in a quasar. We argue that the >10 increase in Eddington ratio inferred from the brightening in UV and X-ray continuum flux is more likely due to an intrinsic change in the accretion rate of a pre-existing accretion disk, than an external mechanism such as variable obscuration, microlensing, or the tidal disruption of a star. However, further monitoring will be helpful in better constraining the mechanism driving this change of state. The rapid "turn on" of the quasar is much shorter than the viscous infall timescale of an accretion disk, and requires a disk instability that can develop around a ~10^8 M_sun black hole on timescales less than a year.

We present basic statistics for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) during its first year-and-a-half of operations, spanning 2013 and 2014. We also present the same information for all other bright ($m_V\leq17$), spectroscopically confirmed supernovae discovered from 2014 May 1 through the end of 2014, providing a comparison to the ASAS-SN sample starting from the point where ASAS-SN became operational in both hemispheres. In addition, we present collected redshifts and near-UV through IR magnitudes, where available, for all host galaxies of the bright supernovae in both samples. This work represents a comprehensive catalog of bright supernovae and their hosts from multiple professional and amateur sources, allowing for population studies that were not previously possible because the all-sky emphasis of ASAS-SN redresses many previously existing biases. In particular, ASAS-SN systematically finds bright supernovae closer to the centers of host galaxies than either other professional surveys or amateurs, a remarkable result given ASAS-SN's poorer angular resolution. This is the first of a series of yearly papers on bright supernovae and their hosts that will be released by the ASAS-SN team.

We present the results of optical, near-infrared, and mid-infrared observations of M101 OT2015-1 (PSN J14021678+5426205), a luminous red transient in the Pinwheel galaxy (M101), spanning a total of 16 years. The lightcurve showed two distinct peaks with absolute magnitudes $M_r\leq-12.4$ and $M_r \simeq-12$, on 2014 November 11 and 2015 February 17, respectively. The spectral energy distributions during the second maximum show a cool outburst temperature of $\approx$3700 K and low expansion velocities ($\approx-$300 \kms) for the H I, Ca II, Ba II and K I lines. From archival data spanning 15 to 8 years before the outburst, we find a single source consistent with the optically discovered transient which we attribute to being the progenitor; it has properties consistent with being an F-type yellow supergiant with $L$~$\sim$~8.7~$\times\ 10^4$ \Lsun, $T_{\rm{eff}}\approx$7000~K and an estimated mass of $\rm{M1}= 18\pm 1$ \Msun. This star has likely just finished the H burning phase in the core, started expanding, and is now crossing the Hertzsprung gap. Based on the combination of observed properties, we argue that the progenitor is a binary system, with the more evolved system overfilling the Roche lobe. Comparison with binary evolution models suggests that the outburst was an extremely rare phenomenon, likely associated with the ejection of the common envelope. The initial mass of the binary progenitor system fills the gap between the merger candidates V838 Mon (5$-$10 \Msun) and NGC~4490-OT~(30~\Msun).

Type Ic supernovae (SNe Ic) arise from the core-collapse of H (and He) poor stars, which could be either single WR stars or lower-mass stars stripped of their envelope by a companion. Their light curves are radioactively powered and usually show a fast rise to peak ($\sim$10-15 d), without any early (first few days) emission bumps (with the exception of broad-lined SNe Ic) as sometimes seen for other types of stripped-envelope SNe (e.g., Type IIb SN 1993J and Type Ib SN 2008D). We have studied iPTF15dtg, a spectroscopically normal SN Ic with an early excess in the optical light curves followed by a long ($\sim$30 d) rise to the main peak. It is the first spectroscopically-normal double-peaked SN Ic observed. We aim to determine the properties of this explosion and of its progenitor star. Optical photometry and spectroscopy of iPTF15dtg was obtained with multiple telescopes. The resulting light curves and spectral sequence are analyzed and modelled with hydrodynamical and analytical models, with particular focus on the early emission. Results. iPTF15dtg is a slow rising SN Ic, similar to SN 2011bm. Hydrodynamical modelling of the bolometric properties reveals a large ejecta mass ($\sim$10 $M_{\odot}$) and strong $^{56}$Ni mixing. The luminous early emission can be reproduced if we account for the presence of an extended ($\sim$500 R$_{\odot}$), low-mass ($\sim$0.045 M$_{\odot}$) envelope around the progenitor star. Alternative scenarios for the early peak, such as the interaction with a companion, a shock-breakout (SBO) cooling tail from the progenitor surface, or a magnetar-driven SBO are not favored. The large ejecta mass and the presence of H and He free extended material around the star suggest that the progenitor of iPTF15dtg was a massive ($\gtrsim$ 35 M$_{\odot}$) WR star suffering strong mass loss.

Temporal variability of narrow absorption lines in high-resolution spectra of Type Ia supernovae (SNe Ia) is studied to search for circumstellar matter. Time series which resolve the profiles of absorption lines such as Na I D or Ca II H&K are expected to reveal variations due to photoionisation and subsequent recombination of the gases. The presence, composition, and geometry of circumstellar matter may hint at the elusive progenitor system of SNe Ia and could also affect the observed reddening law. To date, there are few known cases of time-varying Na I D absorption in SNe Ia, all of which occurred during relatively late phases of the supernova evolution. Photoionisation, however, is predicted to occur during the early phases of SNe Ia, when the supernova peaks in the ultraviolet. We therefore attempt to observe early-time absorption-line variations by obtaining high-resolution spectra of SNe before maximum light. We have obtained photometry and high-resolution spectroscopy of SNe Ia 2013gh and iPTF 13dge, to search for absorption- line variations. Furthermore, we study interstellar absorption features in relation to the observed photometric colours of the SNe. Results. Both SNe display deep Na I D and Ca II H&K absorption features. Furthermore, small but significant variations are detected in a feature of the Na I D profile of SN 2013gh. The variations are consistent with either geometric effects of rapidly moving or patchy gas clouds or photoionisation of Na I gas at R \approx 1019 cm from the explosion. Our analysis indicates that it is necessary to focus on early phases to detect photoionisation effects of gases in the circumstellar medium of SNe Ia. Different absorbers such as Na I and Ca II can be used to probe for matter at different distances from the SNe.

Studies of high-redshift gamma-ray bursts (GRBs) provide important information about the early Universe such as the rates of stellar collapsars and mergers, the metallicity content, constraints on the re-ionization period, and probes of the Hubble expansion. Rapid selection of high-z candidates from GRB samples reported in real time by dedicated space missions such as Swift is the key to identifying the most distant bursts before the optical afterglow becomes too dim to warrant a good spectrum. Here we introduce "machine-z", a redshift prediction algorithm and a "high-z" classifier for Swift GRBs based on machine learning. Our method relies exclusively on canonical data commonly available within the first few hours after the GRB trigger. Using a sample of 284 bursts with measured redshifts, we trained a randomized ensemble of decision trees (random forest) to perform both regression and classification. Cross-validated performance studies show that the correlation coefficient between machine-z predictions and the true redshift is nearly 0.6. At the same time our high-z classifier can achieve 80% recall of true high-redshift bursts, while incurring a false positive rate of 20%. With 40% false positive rate the classifier can achieve ~100% recall. The most reliable selection of high-redshift GRBs is obtained by combining predictions from both the high-z classifier and the machine-z regressor.

Supernova (SN) 1987A was a peculiar H-rich event with a long-rising (LR) light curve (LC), stemming from a compact blue supergiant star (BSG). Only a few similar events have been presented in the literature. We present new data for a sample of 6 LR Type II SNe (SNe II), 3 of which were discovered and observed by the Palomar Transient Factory (PTF) and 3 observed by the Caltech Core-Collapse Project (CCCP). Our aim is to enlarge the family of LR SNe II, characterizing their properties. Spectra, LCs, and host-galaxies (HG) of these SNe are presented. Comparisons with known SN 1987A-like events are shown, with emphasis on the absolute magnitudes, colors, expansion velocities, and HG metallicities. Bolometric properties are derived from the multiband LC. By modeling the early-time LCs with scaling relations derived from the SuperNova Explosion Code (SNEC) models of MESA progenitor stars, we estimate the progenitor radii of these SNe and other progenitor parameters. We present PTF12kso, a LR SN II with the largest amount of 56Ni mass for this class. PTF09gpn and PTF12kso are found at the lowest HG metallicities for this SN group. The variety of early LC luminosities depends on the wide range of progenitor radii, from a few tens of solar radii (SN 2005ci) up to thousands (SN 2004ek) with intermediate cases between 100 (PTF09gpn) and 300 solar radii (SN 2004em). We confirm that LR SNe II with LC shapes closely resembling that of SN 1987A generally arise from BSGs. However, some of them likely have progenitors with larger radii (~300 solar radii, typical of yellow supergiants) and can thus be regarded as intermediate cases between normal SNe IIP and SN 1987A-like SNe. Some extended red supergiant (RSG) stars such as the progenitor of SN 2004ek can also produce LR SNe II if they synthesized a large amount of 56Ni. Low HG metallicity is confirmed as a characteristic of BSG SNe.

We present ground-based and Swift photometric and spectroscopic observations of the candidate tidal disruption event (TDE) ASASSN-14li, found at the center of PGC 043234 ($d\simeq90$ Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN). The source had a peak bolometric luminosity of $L\simeq10^{44}$ ergs s$^{-1}$ and a total integrated energy of $E\simeq7\times10^{50}$ ergs radiated over the $\sim6$ months of observations presented. The UV/optical emission of the source is well-fit by a blackbody with roughly constant temperature of $T\sim35,000$ K, while the luminosity declines by roughly a factor of 16 over this time. The optical/UV luminosity decline is broadly consistent with an exponential decline, $L\propto e^{-t/t_0}$, with $t_0\simeq60$ days. ASASSN-14li also exhibits soft X-ray emission comparable in luminosity to the optical and UV emission but declining at a slower rate, and the X-ray emission now dominates. Spectra of the source show broad Balmer and helium lines in emission as well as strong blue continuum emission at all epochs. We use the discoveries of ASASSN-14li and ASASSN-14ae to estimate the TDE rate implied by ASAS-SN, finding an average rate of $r \simeq 4.1 \times 10^{-5}~{\rm yr}^{-1}$ per galaxy with a 90% confidence interval of $(2.2 - 17.0) \times 10^{-5}~{\rm yr}^{-1}$ per galaxy. ASAS-SN found roughly 1 TDE for every 70 Type Ia supernovae in 2014, a rate that is much higher than that of other surveys.

Temporal and spectral information extracted from a stream of photons received from astronomical sources is the foundation on which we build understanding of various objects and processes in the Universe. Typically astronomers fit a number of models separately to light curves and spectra to extract relevant features. These features are then used to classify, identify, and understand the nature of the sources. However, these feature extraction methods may not be optimally sensitive to unknown properties of light curves and spectra. One can use the raw light curves and spectra as features to train classifiers, but this typically increases the dimensionality of the problem, often by several orders of magnitude. We overcome this problem by integrating light curves and spectra to create an abstract image and using wavelet analysis to extract important features from the image. Such features incorporate both temporal and spectral properties of the astronomical data. Classification is then performed on those abstract features. In order to demonstrate this technique, we have used gamma-ray burst (GRB) data from the NASA's Swift mission to classify GRBs into high- and low-redshift groups. Reliable selection of high-redshift GRBs is of considerable interest in astrophysics and cosmology.

We report the discovery of ASASSN-15lh (SN 2015L), which we interpret as the most luminous supernova yet found. At redshift z = 0.2326, ASASSN-15lh reached an absolute magnitude of M_{u,AB} = -23.5+/-0.1 and bolometric luminosity L_bol = (2.2+/-0.2)x 10^45 ergs s^-1, which is more than twice as luminous as any previously known supernova. It has several major features characteristic of the hydrogen-poor super-luminous supernovae (SLSNe-I), whose energy sources and progenitors are currently poorly understood. In contrast to most previously known SLSNe-I that reside in star-forming dwarf galaxies, ASASSN-15lh appears to be hosted by a luminous galaxy (M_K ~ -25.5) with little star formation. In the 4 months since first detection, ASASSN-15lh radiated (1.1+/- 0.2)x10^52 ergs, challenging the magnetar model for its engine.

Extreme coronal-line emitter (ECLE) SDSSJ095209.56+214313.3, known by its strong, fading, high ionization lines, has been a long standing candidate for a tidal disruption event, however a supernova origin has not yet been ruled out. Here we add several new pieces of information to the puzzle of the nature of the transient that powered its variable coronal lines: 1) an optical light curve from the Lincoln Near Earth Asteroid Research (LINEAR) survey that serendipitously catches the optical flare, and 2) late-time observations of the host galaxy with the Swift Ultraviolet and Optical Telescope (UVOT) and X-ray telescope (XRT) and the ground-based Mercator telescope. The well-sampled, $\sim10$-year long, unfiltered LINEAR light curve constrains the onset of the flare to a precision of $\pm5$ days and enables us to place a lower limit on the peak optical magnitude. Difference imaging allows us to estimate the location of the flare in proximity of the host galaxy core. Comparison of the \textsl{GALEX} data (early 2006) with the recently acquired Swift UVOT (June 2015) and Mercator observations (April 2015) demonstrate a decrease in the UV flux over a $\sim 10$ year period, confirming that the flare was UV-bright. The long-lived UV-bright emission, detected 1.8 rest-frame years after the start of the flare, strongly disfavors a SN origin. These new data allow us to conclude that the flare was indeed powered by the tidal disruption of a star by a supermassive black hole and that TDEs are in fact capable of powering the enigmatic class of ECLEs.

Supernovae (SNe) embedded in dense circumstellar material (CSM) may show prominent emission lines in their early-time spectra ($\leq 10$ days after the explosion), owing to recombination of the CSM ionized by the shock-breakout flash. From such spectra ("flash spectroscopy"), we can measure various physical properties of the CSM, as well as the mass-loss rate of the progenitor during the year prior to its explosion. Searching through the Palomar Transient Factory (PTF and iPTF) SN spectroscopy databases from 2009 through 2014, we found 12 Type II SNe showing flash-ionized (FI) signatures in their first spectra. All are younger than 10 days. These events constitute 14\% of all 84 SNe in our sample having a spectrum within 10 days from explosion, and 18\% of SNe~II observed at ages $<5$ days, thereby setting lower limits on the fraction of FI events. We classified as "blue/featureless" (BF) those events having a first spectrum which is similar to that of a black body, without any emission or absorption signatures. It is possible that some BF events had FI signatures at an earlier phase than observed, or that they lack dense CSM around the progenitor. Within 2 days after explosion, 8 out of 11 SNe in our sample are either BF events or show FI signatures. Interestingly, we found that 19 out of 21 SNe brighter than an absolute magnitude $M_R=-18.2$ belong to the FI or BF groups, and that all FI events peaked above $M_R=-17.6$ mag, significantly brighter than average SNe~II.

Gamma-ray bursts (GRBs) are detectable out to very large distances and as such are potentially powerful cosmological probes. Historically, the angular distribution of GRBs provided important information about their origin and physical properties. As a general population, GRBs are distributed isotropically across the sky. However, there are published reports that once binned by duration or redshift, GRBs display significant clustering. We have studied the redshift- and duration-dependent clustering of GRBs using proximity measures and kernel density estimation. Utilizing bursts detected by BATSE, Fermi/GBM and Swift/BAT, we found marginal evidence for clustering in very short duration GRBs lasting less than 100 ms. Our analysis provides little evidence for significant redshift-dependent clustering of GRBs.

We present the results of the spectroscopic and photometric monitoring campaign of ASASSN-15ed. The transient was discovered quite young by the All Sky Automated Survey for SuperNovae (ASAS-SN). Amateur astronomers allowed us to sample the photometric SN evolution around maximum light, which we estimate to have occurred on JD = 2457087.4 +- 0.6 in the r-band. Its apparent r-band magnitude at maximum was r = 16.91 +- 0.10, providing an absolute magnitude M(r) ~ -20.04 +- 0.20, which is slightly more luminous than the typical magnitudes estimated for Type Ibn SNe. The post-peak evolution was well monitored, and the decline rate (being in most bands around 0.1 mag/d during the first 25 d after maximum) is marginally slower than the average decline rates of SNe Ibn during the same time interval. The object was initially classified as a Type Ibn SN because early-time spectra were characterized by a blue continuum with superimposed narrow P-Cygni lines of He I, suggesting the presence of a slowly moving (1200-1500 km/s), He-rich circumstellar medium. Later on, broad P-Cygni He I lines became prominent. The inferred velocities, as measured from the minimum of the broad absorption components, were between 6000 and 7000 km/s. As we attribute these broad features to the SN ejecta, this is the first time we have observed the transition of a Type Ibn SN to a Type Ib SN.

We report here the discovery by the Intermediate Palomar Transient Factory (iPTF) of iPTF14yb, a luminous ($M_{r}\approx-27.8$ mag), cosmological (redshift 1.9733), rapidly fading optical transient. We demonstrate, based on probabilistic arguments and a comparison with the broader population, that iPTF14yb is the optical afterglow of the long-duration gamma-ray burst GRB 140226A. This marks the first unambiguous discovery of a GRB afterglow prior to (and thus entirely independent of) an associated high-energy trigger. We estimate the rate of iPTF14yb-like sources (i.e., cosmologically distant relativistic explosions) based on iPTF observations, inferring an all-sky value of $\Re_{\mathrm{rel}}=610$ yr$^{-1}$ (68% confidence interval of 110-2000 yr$^{-1}$). Our derived rate is consistent (within the large uncertainty) with the all-sky rate of on-axis GRBs derived by the Swift satellite. Finally, we briefly discuss the implications of the nondetection to date of bona fide "orphan" afterglows (i.e., those lacking detectable high-energy emission) on GRB beaming and the degree of baryon loading in these relativistic jets.

We investigate the effect that the absorption of high-energy (above 100 MeV) photons produced in GRB afterglow shocks has on the light-curves and spectra of Fermi-LAT afterglows. Afterglows produced by the interaction of a relativistic outflow with a wind-like medium peak when the blast-wave deceleration sets in, and the afterglow spectrum could be hardening before that peak, as the optical thickness to pair-formation is decreasing. In contrast, in afterglows produced in the interaction with a homogeneous medium, the optical thickness to pair-formation should increase and yield a light-curve peak when it reaches unity, followed by a fast light-curve decay, accompanied by a spectral softening. If energy is injected in the blast-wave, then the accelerated increase of the optical thickness yields a convex afterglow light-curve. Other features, such as a double-peak light-curve or a broad hump, can arise from the evolution of the optical thickness to photon-photon absorption. Fast decays and convex light-curves are seen in a few LAT afterglows, but the expected spectral softening is rarely seen in (and difficult to measure with) LAT observations. Furthermore, for the effects of photon-photon attenuation to shape the high-energy afterglow light-curve without attenuating it too much, the ejecta initial Lorentz factor must be in a relatively narrow range (50-200), which reduces the chance of observing those effects.

The complex multiwavelength emission of GRB afterglow 130427A (monitored in the radio up to 10 days, in the optical and X-ray until 50 days, and at GeV energies until 1 day) can be accounted for by a hybrid reverse-forward shock synchrotron model, with inverse-Compton emerging only above a few GeV. The high ratio of the early optical to late radio flux requires that the ambient medium is a wind and that the forward-shock synchrotron spectrum peaks in the optical at about 10 ks. The latter has two consequences: the wind must be very tenuous and the optical emission before 10 ks must arise from the reverse-shock, as suggested also by the bright optical flash that Raptor has monitored during the prompt emission phase (<100 s). The VLA radio emission is from the reverse-shock, the Swift X-ray emission is mostly from the forward-shock, but the both shocks give comparable contributions to the Fermi GeV emission. The weak wind implies a large blast-wave radius (8 t_{day}^{1/2} pc), which requires a very tenuous circumstellar medium, suggesting that the massive stellar progenitor of GRB 130427A resided in a super-bubble.

The optical light that is generated simultaneously with the x-rays and gamma-rays during a gamma-ray burst (GRB) provides clues about the nature of the explosions that occur as massive stars collapse to form black holes. We report on the bright optical flash and fading afterglow from the powerful burst GRB 130427A and present a comparison with the properties of the gamma-ray emission that show correlation of the optical and >100 MeV photon flux light curves during the first 7,000 seconds. We attribute this correlation to co-generation in an external shock. The simultaneous, multi-color, optical observations are best explained at early times by reverse shock emission generated in the relativistic burst ejecta as it collides with surrounding material and at late times by a forward shock traversing the circumburst environment. The link between optical afterglow and >100 MeV emission suggests that nearby early peaked afterglows will be the best candidates for studying GRB emission at GeV/TeV energies.

We discuss the results of the analysis of multi-wavelength data for the afterglows of GRB 081007 and GRB 090424, two bursts detected by Swift. One of them, GRB 081007, also shows a spectroscopically confirmed supernova, SN 2008hw, which resembles SN 1998bw in its absorption features, while the maximum luminosity is only about half as large as that of SN 1998bw. Bright optical flashes have been detected in both events, which allows us to derive solid constraints on the circumburst-matter density profile. This is particularly interesting in the case of GRB 081007, whose afterglow is found to be propagating into a constant-density medium, yielding yet another example of a GRB clearly associated with a massive star progenitor which did not sculpt the surroundings with its stellar wind. There is no supernova component detected in the afterglow of GRB 090424, likely due to the brightness of the host galaxy, comparable to the Milky Way. We show that the afterglow data are consistent with the presence of both forward- and reverse-shock emission powered by relativistic outflows expanding into the interstellar medium. The absence of optical peaks due to the forward shock strongly suggests that the reverse shock regions should be mildly magnetized. The initial Lorentz factor of outflow of GRB 081007 is estimated to be \Gamma ~ 200, while for GRB 090424 a lower limit of \Gamma > 170 is derived. We also discuss the prompt emission of GRB 081007, which consists of just a single pulse. We argue that neither the external forward-shock model nor the shock-breakout model can account for the prompt emission data and suggest that the single-pulse-like prompt emission may be due to magnetic energy dissipation of a Poynting-flux dominated outflow or to a dissipative photosphere.

We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available.

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. 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. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the 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. Using a 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 is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

The peaks of 30 optical afterglows and 14 X-ray light-curves display a good anticorrelation of the peak flux with the peak epoch: F_p ~ t_p^{-2.0} in the optical, F_p ~ t_p^{-1.6} in the X-ray, the distributions of the peak epochs being consistent with each other. We investigate the ability of two forward-shock models for afterglow light-curve peaks -- an observer location outside the initial jet aperture and the onset of the forward-shock deceleration -- to account for those peak correlations. For both models, the slope of the F_p - t_p relation depends only on the slope of the afterglow spectrum. We find that only a conical jet seen off-aperture and interacting with a wind-like medium can account for both the X-ray peak relation, given the average X-ray spectral slope beta_x = 1.0, and for the larger slope of the optical peak relation. However, any conclusion about the origin of the peak flux - peak epoch correlation is, at best, tentative, because the current sample of X-ray peaks is too small to allow a reliable measurement of the F_p - t_p relation slope and because more than one mechanism and/or one afterglow parameter may be driving that correlation.

We consider density estimators based on the nearest neighbors method applied to discrete point distibutions in spaces of arbitrary dimensionality. If the density is constant, the volume of a hypersphere centered at a random location is proportional to the expected number of points falling within the hypersphere radius. The distance to the $N$-th nearest neighbor alone is then a sufficient statistic for the density. In the non-uniform case the proportionality is distorted. We model this distortion by normalizing hypersphere volumes to the largest one and expressing the resulting distribution in terms of the Legendre polynomials. Using Monte Carlo simulations we show that this approach can be used to effectively address the tradeoff between smoothing bias and estimator variance for sparsely sampled distributions.

We report new mid-eclipse times of the two close binaries NSVS14256825 and HS0705+6700, harboring an sdB primary and a low-mass main-sequence secondary. Both objects display clear variations in the measured orbital period, which can be explained by the action of a third object orbiting the binary. If this interpretation is correct, the third object in NSVS14256825 is a giant planet with a mass of roughly 12 M_Jup. For HS0705+6700, we provide evidence that strengthens the case for the suggested periodic nature of the eclipse time variation and reduces the uncertainties in the parameters of the brown dwarf implied by that model. The derived period is 8.4 yr and the mass is 31 M_Jup, if the orbit is coplanar with the binary. This research is part of the PlanetFinders project, an ongoing collaboration between professional astronomers and student groups at high schools.

VOEvent defines the content and meaning of a standard information packet for representing, transmitting, publishing and archiving information about a transient celestial event, with the implication that timely follow-up is of interest. The objective is to motivate the observation of targets-of-opportunity, to drive robotic telescopes, to trigger archive searches, and to alert the community. VOEvent is focused on the reporting of photon events, but events mediated by disparate phenomena such as neutrinos, gravitational waves, and solar or atmospheric particle bursts may also be reported. Structured data is used, rather than natural language, so that automated systems can effectively interpret VOEvent packets. Each packet may contain zero or more of the "who, what, where, when & how" of a detected event, but in addition, may contain a hypothesis (a "why") regarding the nature of the underlying physical cause of the event. Citations to previous VOEvents may be used to place each event in its correct context. Proper curation is encouraged throughout each event's life cycle from discovery through successive follow-ups. VOEvent packets gain persistent identifiers and are typically stored in databases reached via registries. VOEvent packets may therefore reference other packets in various ways. Packets are encouraged to be small and to be processed quickly. This standard does not define a transport layer or the design of clients, repositories, publishers or brokers; it does not cover policy issues such as who can publish, who can build a registry of events, who can subscribe to a particular registry, nor the intellectual property issues.

We report the discovery of a six-month-long mid-infrared transient, SDWFS-MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO Deep Wide-Field Survey Bootes field. The transient, located in a z=0.19 low luminosity (M_[4.5]~-18.6 mag, L/L_MilkyWay~0.01) metal-poor (12+log(O/H)~7.8) irregular galaxy, peaked at a mid-infrared absolute magnitude of M_[4.5]~-24.2 in the 4.5 micron Spitzer/IRAC band and emitted a total energy of at least 10^51 ergs. The optical emission was likely fainter than the mid-infrared, although our constraints on the optical emission are poor because the transient peaked when the source was "behind" the Sun. The Spitzer data are consistent with emission by a modified black body with a temperature of ~1350 K. We rule out a number of scenarios for the origin of the transient such as a Galactic star, AGN activity, GRB, tidal disruption of a star by a black hole and gravitational lensing. The most plausible scenario is a supernova exploding inside a massive, optically thick circumstellar medium, composed of multiple shells of previously ejected material. If the proposed scenario is correct, then a significant fraction (~10%) of the most luminous supernova may be self-enshrouded by dust not only before but also after the supernova occurs. The spectral energy distribution of the progenitor of such a supernova would be a slightly cooler version of eta Carina, peaking at 20-30 microns.

We present observations and analysis of the broadband afterglow of Swift GRB 071025. Using optical and infrared (RIYJHK) photometry, we derive a photometric redshift of 4.4 < z < 5.2; at this redshift our simultaneous multicolour observations begin at ~30 s after the GRB trigger in the host frame and during the initial rising phase of the afterglow. We associate the light curve peak at 580 s in the observer frame with the formation of the forward shock, giving an estimate of the initial Lorentz factor Gamma_0 ~ 200. The red spectral energy distribution (even in regions not affected by the Lyman-alpha break) provides secure evidence of a large dust column. However, the inferred extinction curve shows a prominent flat component between 2000-3000 Angstroms in the rest-frame, inconsistent with any locally observed template but well-fit by models of dust formed by supernovae. Time-dependent fits to the extinction profile reveal no evidence of dust destruction and limit the decrease in the extinction column to Delta A_3000 < 0.54 mag after t = 50 s in the rest frame. Our observations provide evidence of a transition in dust properties at z~5, in agreement with studies of high-z quasars, and suggest that SN-formed dust continues to dominate the opacity of typical galaxies at this redshift.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

On 19 March 2008, the northern sky was the stage of a spectacular optical transient that for a few seconds remained visible to the naked eye. The transient was associated with GRB 080319B, a gamma-ray burst at a luminosity distance of about 6 Gpc (standard cosmology), making it the most luminous optical object ever recorded by human kind. We present comprehensive sky monitoring and multi-color optical follow-up observations of GRB 080319B collected by the RAPTOR telescope network covering the development of the explosion and the afterglow before, during, and after the burst. The extremely bright prompt optical emission revealed features that are normally not detectable. The optical and gamma-ray variability during the explosion are correlated, but the optical flux is much greater than can be reconciled with single emission mechanism and a flat gamma-ray spectrum. This extreme optical behavior is best understood as synchrotron self-Compton model (SSC). After a gradual onset of the gamma-ray emission, there is an abrupt rise of the prompt optical flux suggesting that variable self-absorption dominates the early optical light curve. Our simultaneous multi-color optical light curves following the flash show spectral evolution consistent with a rapidly decaying red component due to large angle emission and the emergence of a blue forward shock component from interaction with the surrounding environment. While providing little support for the reverse shock that dominates the early afterglow, these observations strengthen the case for the universal role of the SSC mechanism in generating gamma-ray bursts.

We present a direct detection of the gravitational lens that caused the microlensing event MACHO-95-BLG-37. This is the first fully resolved microlensing system involving a source in the Galactic bulge, and the second such system in general. The lens and source are clearly resolved in images taken with the High Resolution Channel of the Advanced Camera for Surveys on board the Hubble Space Telescope (HST) ~9 years after the microlensing event. The presently available data are not sufficient for the final, unambiguous identification of the gravitational lens and the microlensed source. While the light curve models combined with the high resolution photometry for individual objects indicate that the source is red and the lens is blue, the color-magnitude diagram for the line of sight and the observed proper motions strongly support the opposite case. The first scenario points to a metal-poor lens with mass M = ~0.6 M_Sun at the distance D_l = ~4 kpc. In the second scenario the lens could be a main-sequence star with M = 0.8 - 0.9 M_Sun about half-way to the Galactic bulge or in the foreground disk, depending on the extinction.

Gravitational microlensing surveys target very dense stellar fields in the local group. As a consequence the microlensed source stars are often blended with nearby unresolved stars. The presence of `blending' is a cause of major uncertainty when determining the lensing properties of events towards the Galactic centre. After demonstrating empirical cases of blending we utilize Monte Carlo simulations to probe the effects of blending. We generate artificial microlensing events using an HST luminosity function convolved to typical ground-based seeing, adopting a range of values for the stellar density and seeing. We find that a significant fraction of bright events are blended, contrary to the oft-quoted assumption that bright events should be free from blending. We probe the effect that this erroneous assumption has on both the observed event timescale distribution and the optical depth, using realistic detection criteria relevent to the different surveys. Importantly, under this assumption the latter quantity appears to be reasonably unaffected across our adopted values for seeing and density. The timescale distribution is however biased towards smaller values, even for the least dense fields. The dominant source of blending is from lensing of faint source stars, rather than lensing of bright source stars blended with nearby fainter stars. We also explore other issues, such as the centroid motion of blended events and the phenomena of `negative' blending. Furthermore, we breifly note that blending can affect the determination of the centre of the red clump giant region from an observed luminosity function. This has implications for a variety of studies, e.g. mapping extinction towards the bulge and attempts to constrain the parameters of the Galactic bar through red clump giant number counts. (Abridged)

We report on observations of a gamma-ray burst (GRB 061126) with an extremely bright (R ~ 12 mag at peak) early-time optical afterglow. The optical afterglow is already fading as a power law 22 seconds after the trigger, with no detectable prompt contribution in our first exposure, which was coincident with a large prompt-emission gamma-ray pulse. The optical--infrared photometric spectral energy distribution is an excellent fit to a power law, but it exhibits a moderate red-to-blue evolution in the spectral index at about 500 s after the burst. This color change is contemporaneous with a switch from a relatively fast decay to slower decay. The rapidly decaying early afterglow is broadly consistent with synchrotron emission from a reverse shock, but a bright forward-shock component predicted by the intermediate- to late-time X-ray observations under the assumptions of standard afterglow models is not observed. Indeed, despite its remarkable early-time brightness, this burst would qualify as a dark burst at later times on the basis of its nearly flat optical-to-X-ray spectral index. Our photometric spectral energy distribution provides no evidence of host-galaxy extinction, requiring either large quantities of grey dust in the host system (at redshift 1.1588 +/- 0.0006, based upon our late-time Keck spectroscopy) or separate physical origins for the X-ray and optical afterglows.

We use data from the Northern Sky Variability Survey (NSVS), obtained from the first generation Robotic Optical Transient Search Experiment (ROTSE-I), to identify and study RR Lyrae variable stars in the solar neighborhood. We initially identified 1197 RRab (RR0) candidate stars brighter than the ROTSE median magnitude V = 14. Periods, amplitudes, and mean V magnitudes are determined for a subset of 1188 RRab stars with well defined light curves. Metallicities are determined for 589 stars by the Fourier parameter method and by the relationship between period, amplitude, and [Fe/H]. We comment upon the difficulties of clearly classifying RRc (RR1) variables in the NSVS dataset. Distances to the RRab stars are calculated using an adopted luminosity-metallicity relation with corrections for interstellar extinction. The 589 RRab stars in our final sample are used to study the properties of the RRab population within 5 kpc of the Sun. The Bailey diagram of period versus amplitude shows that the largest component of this sample belongs to Oosterhoff type I. Metal-rich ([Fe/H] > -1) RRab stars appear to be associated with the Galactic disk. Our metal-rich RRab sample may include a thin disk as well as a thick disk population, although the uncertainties are too large to establish this. There is some evidence among the metal-rich RRab stars for a decline in scale height with increasing [Fe/H], as was found by Layden (1995). The distribution of RRab stars with -1 < [Fe/H] < -1.25 indicates that within this metallicity range the RRab stars are a mixture of stars belonging to halo and disk populations.

The taxonomy of optical emission detected during the critical first few minutes after the onset of a gamma-ray burst (GRB) defines two broad classes: prompt optical emission correlated with prompt gamma-ray emission, and early optical afterglow emission uncorrelated with the gamma-ray emission. The standard theoretical interpretation attributes prompt emission to internal shocks in the ultra-relativistic outflow generated by the internal engine; early afterglow emission is attributed to shocks generated by interaction with the surrounding medium. Here we report on observations of a bright GRB that, for the first time, clearly show the temporal relationship and relative strength of the two optical components. The observations indicate that early afterglow emission can be understood as reverberation of the energy input measured by prompt emission. Measurements of the early afterglow reverberations therefore probe the structure of the environment around the burst, whereas the subsequent response to late-time impulsive energy releases reveals how earlier flaring episodes have altered the jet and environment parameters. Many GRBs are generated by the death of massive stars that were born and died before the Universe was ten per cent of its current age, so GRB afterglow reverberations provide clues about the environments around some of the first stars.

We present a proper motion mini-survey of 35 fields in the vicinity of Baade window, (l, b) = (1 deg, -4 deg), sampling roughly a 5 x 2.5 deg region of the Galactic bar. Our second epoch observations collected with the ACS/HRC instrument on board the Hubble Space Telescope were combined with the archival WFPC2/PC images. The resulting time baselines are in the range of 4 - 8 years. Precise proper motions of 15,863 stars were determined in the reference frame defined by the mean motion of stars with magnitudes between I_F814W = 16.5 - 21.5 along the line of sight. We clearly detect small gradients in proper motion dispersions (sigma_l, sigma_b) ~ (3.0, 2.5) mas/yr, and in the amount of anisotropy (sigma_l/sigma_b ~ 1.2). Both the longitude dispersion sigma_l and its ratio to the vertical motion sigma_b increase toward the Galactic plane. The decline of the anisotropy ratio sigma_l/sigma_b toward the minor axis of the bulge is mostly due to increasing sigma_b. We also find, for the first time, a significant negative covariance term in the transverse velocity field sigma_lb/(sigma_l*sigma_b) ~ -0.10. Our results extend by a factor of ~15 the number of the Galactic bar fields with good proper motion dispersions.

The RAPid Telescopes for Optical Response (RAPTOR) system at Los Alamos National Laboratory observed GRB 060206 starting 48.1 minutes after gamma-ray emission triggered the Burst Alert Telescope (BAT) on-board the Swift satellite. The afterglow light curve measured by RAPTOR shows a spectacular re-brightening by ~1 mag about 1 h after the trigger and peaks at R ~ 16.4 mag. Shortly after the onset of the explosive re-brightening the OT doubled its flux on a time-scale of about 4 minutes. The total R-band fluence received from GRB 060206 during this episode is 2.3e-9 erg/cm2. In the rest frame of the burst (z = 4.045) this yields an isotropic equivalent energy release of ~0.7e50 erg in just a narrow UV band 130 +/- 22 nm. We discuss the implications of RAPTOR observations for untriggered searches for fast optical transients and studies of GRB environments at high redshift.

We present a measurement of the microlensing optical depth toward the Galactic Bulge based on 4 years of the OGLE-II survey using Red Clump Giant (RCG). Using 32 events we find tau=2.55_{-0.46}^{+0.57}* 10^{-6} at (l,b)=(1.16, -2.75). Taking into account the measured gradient along the Galactic latitude b, tau = [ (4.48+/- 2.37) + (0.78+/- 0.84)* b]* 10^{-6}, this value is consistent with previous measurements using RCG sources and recent theoretical predictions. We determine the microlensing parameters and select events using a model light curve with the flux blending. We find that ~38% of the OGLE-II events which appear to have RCG sources are actually due to much fainter stars blended with a bright companion. We show explicitly that model fits without blending result in similar tau estimates through partial cancellation of contributions from higher detection efficiency, underestimated time-scales and larger number of selected events. This approach, however, leads to biased time-scale distributions and event rates. Consequently, microlensing studies should carefully consider source confusion effects even for bright stars.

The RAPid Telescopes for Optical Response (RAPTOR) system at Los Alamos National Laboratory observed GRB 050319 starting 25.4 seconds after gamma-ray emission triggered the Burst Alert Telescope (BAT) on-board the Swift satellite. Our well sampled light curve of the early optical afterglow is composed of 32 points (derived from 70 exposures) that measure the flux decay during the first hour after the GRB. The GRB 050319 light curve measured by RAPTOR can be described as a relatively gradual flux decline (power-law index alpha = -0.37) with a transition, at about 400 s after the GRB, to a faster flux decay (alpha = -0.91). The addition of other available measurements to the RAPTOR light curve suggests that another emission component emerged after 10^4 s. We hypothesize that the early afterglow emission is powered by extended energy injection or delayed reverse shock emission followed by the emergence of forward shock emission.

The prompt optical emission that arrives with gamma-rays from a cosmic gamma-ray burst (GRB) is a signature of the engine powering the burst, the properties of the ultra-relativistic ejecta of the explosion, and the ejecta's interactions with the surroundings. Until now, only GRB 990123 had been detected at optical wavelengths during the burst phase. Its prompt optical emission was variable and uncorrelated with the prompt gamma-ray emission, suggesting that the optical emission was generated by a reverse shock arising from the ejecta's collision with the surrounding material. Here we report prompt optical emission from GRB 041219a. It is variable and correlated with the prompt gamma-rays, indicating a common origin for the optical light and the gamma-rays. Within the context of the standard fireball model of GRBs, we attribute this new optical component to internal shocks driven into the burst ejecta by variations of the inner engine. The correlated optical emission is a direct probe of the jet isolated from the medium. The timing of the uncorrelated optical emission is strongly dependent on the nature of the medium.

We present 97 QSO candidates in 48 Galactic Bulge fields of the Optical Gravitational Lensing Experiment II (OGLE-II) covering about 11 square degrees, which are selected via their variability. We extend light curves of variable objects which were detected in 3-year baseline in the OGLE-II variable star catalog to 4th year. We search for objects which are faint (16<I_0<18.5) and slowly variable during 4 years in this catalog by using the variogram/structure function. Finding the QSOs in such stellar-crowded and high extinction fields is challenging, but should be useful for the astrometric reference frame. Spectroscopic follow-up observations are required to confirm these QSO candidates. Follow-up observations are being prepared for four of these fields (BUL_SC1, 2, 32 and 45). Follow-up observations for other fields are strongly encouraged. The complete list and light curves of all 97 candidates are available in electronic format at http://www.astro.princeton.edu/~sumi/QSO-OGLEII/.

The internet has brought about great change in the astronomical community, but this interconnectivity is just starting to be exploited for use in instrumentation. Utilizing the internet for communicating between distributed astronomical systems is still in its infancy, but it already shows great potential. Here we present an example of a distributed network of telescopes that performs more efficiently in synchronous operation than as individual instruments. RAPid Telescopes for Optical Response (RAPTOR) is a system of telescopes at LANL that has intelligent intercommunication, combined with wide-field optics, temporal monitoring software, and deep-field follow-up capability all working in closed-loop real-time operation. The Telescope ALert Operations Network (TALON) is a network server that allows intercommunication of alert triggers from external and internal resources and controls the distribution of these to each of the telescopes on the network. TALON is designed to grow, allowing any number of telescopes to be linked together and communicate. Coupled with an intelligent alert client at each telescope, it can analyze and respond to each distributed TALON alert based on the telescopes needs and schedule.

We present new version of the OGLE-II catalog of eclipsing binary stars detected in the Small Magellanic Cloud, based on Difference Image Analysis catalog of variable stars in the Magellanic Clouds containing data collected from 1997 to 2000. We found 1351 eclipsing binary stars in the central 2.4 square degree area of the SMC. 455 stars are newly discovered objects, not found in the previous release of the catalog. The eclipsing objects were selected with the automatic search algorithm based on the artificial neural network. The full catalog is accessible from the OGLE Internet archive.

We search for microvariability in a sample of 485 Mira variables with high quality I-band light curves from the second generation Optical Gravitational Lensing Experiment (OGLE-II). Rapid variations with amplitudes in the ~0.2-1.1 mag range lasting hours to days were discovered in Hipparcos data by de Laverny et al. (1998). Our search is primarily sensitive to events with time-scales of about 1 day, but retains a few percent efficiency (per object) for detecting unresolved microvariability events as short as 2 hours. We do not detect any candidate events. Assuming that the distribution of the event time profiles is identical to that from the Hipparcos light curves we derive the 95% confidence level upper limit of 0.038 per year per star for the rate of such events (1 per 26 years per average object of the ensemble). The high event rates of the order of 1 per year per star implied by the Hipparcos study in the H_P band are excluded with high confidence by the OGLE-II data in the I band. Our non-detection could still be explained by much redder spectral response of the I filter compared to the H_P band or by population differences between the bulge and the solar neighborhood. In any case, the OGLE-II I-band data provide the first limit on the rate of the postulated microvariability events in Mira stars and offer new quantitative constraints on their properties. Similar limits are obtained for other pulse shapes and a range of the assumed time-scales and size-frequency distributions.

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

We present a proper motion (\mu) catalogue of 5,080,236 stars in 49 Optical Gravitational Lensing Experiment II (OGLE-II) Galactic bulge (GB) fields, covering a range of -11 deg. <l< 11 deg. and -6 deg. <b<3 deg., the total area close to 11 square degrees. The proper motion measurements are based on 138 - 555 I-band images taken during four observing seasons: 1997-2000. The catalogue stars are in the magnitude range 11 < I < 18 mag. In particular, the catalogue includes Red Clump Giants (RCGs) and Red Giants in the GB, and main sequence stars in the Galactic disc. The proper motions up to \mu = 500 mas/yr were measured with the mean accuracy of 0.8-3.5 mas/yr, depending on the brightness of a star. This catalogue may be useful for studying the kinematic of stars in the GB and the Galactic disk.

This is a progress report of the study of pulsating main-sequence stars in the LMC. Using the OGLE-II photometry supplemented by the MACHO photometry, we find 64 Beta Cephei stars in the LMC. Their periods are generally much longer than observed in stars of this type in the Galaxy (the median value is 0.27 d compared with 0.17 d in the Galaxy). In 20 stars with short periods attributable to the Beta Cephei-type instability, we also find modes with periods longer than 0.4 d. They are likely low-order g modes, which means that in these stars both kinds of variability, Beta Cephei and SPB, are observed. We also show examples of the multiperiodic SPB stars in the LMC, the first beyond our Galaxy.

We present the catalog of RR Lyrae stars discovered in a 4.5 square degrees area in the central parts of the Large Magellanic Cloud (LMC). Presented sample contains 7612 objects, including 5455 fundamental mode pulsators (RRab), 1655 first-overtone (RRc), 272 second-overtone (RRe) and 230 double-mode RR Lyr stars (RRd). Additionally we attach a list of several dozen other short-period pulsating variables. The catalog data include astrometry, periods, BVI photometry, amplitudes, and parameters of the Fourier decomposition of the I-band light curve of each object. We present density map of RR Lyr stars in the observed fields which shows that the variables are strongly concentrated toward the LMC center. The modal values of the period distribution for RRab, RRc and RRe stars are 0.573, 0.339 and 0.276 days, respectively. The period-luminosity diagrams for BVI magnitudes and for extinction insensitive index W_I are constructed. We provide the logP-I, logP-V and logP-W_I relations for RRab, RRc and RRe stars. The mean observed V-band magnitudes of RR Lyr stars in the LMC are 19.36 mag and 19.31 mag for ab and c types, respectively, while the extinction free values are 18.91 mag and 18.89 mag. We also found a large number of RR Lyr stars pulsating in two modes closely spaced in the power spectrum. We discovered three stars which simultaneously reveal RR Lyr-type and eclipsing-type variability. If any of these objects were an eclipsing binary system containing RR Lyr star, then for the first time the direct determination of the mass of RR Lyr variable would be possible. We also provide a list of six LMC star clusters which contain RR Lyr stars. The richest cluster, NGC 1835, hosts 84 RR Lyr variables.

We present the catalog of 2580 eclipsing binary stars detected in 4.6 square degree area of the central parts of the Large Magellanic Cloud. The photometric data were collected during the second phase of the OGLE microlensing search from 1997 to 2000. The eclipsing objects were selected with the automatic search algorithm based on an artificial neural network. Basic statistics of eclipsing stars are presented. Also, the list of 36 candidates of detached eclipsing binaries for spectroscopic study and for precise LMC distance determination is provided. The full catalog is accessible from the OGLE Internet archive.

We describe a unique mass determination for a microlensing event from the second phase of the Optical Gravitational Microlensing Experiment (OGLE-II). The event, sc26_2218, which is very bright (baseline magnitude I=15.10), appears to exhibit both parallax and finite source effects. The parallax effect allows us to determine the projected Einstein radius on the observer plane (approx. 3.8 AU), while the finite source effect allows us to determine the ratio of the angular source size and the angular Einstein radius. As the angular size of the star can be estimated using its color and magnitude, we can hence determine the angular Einstein radius (approx. 0.1 mas). By combining both the projected and angular Einstein radius we can determine the lens mass M = 0.050^{+0.016}_{-0.011} solar masses, independent of the source distance. The lens is therefore a brown dwarf candidate. However, the `parallax' signature is weak and so we cannot completely discount the possibility that these signatures originate from binary rotation of the source (which would prevent any estimate of the lens mass), rather than parallax. However, this can be tested by future spectroscopic observations. This event highlights the scientific returns for intense monitoring of bright microlensing events, since the parallax and finite source effects can be more easily identified due to their high signal-to-noise ratios.

We present the catalog of RR Lyrae stars from 2.4 square degrees of central parts of the Small Magellanic Cloud (SMC). The photometric data were collected during four years of the OGLE-II microlensing survey. Photometry of each star was obtained using the Difference Image Analysis (DIA) method. The catalog contains 571 objects, including 458 RRab, 56 RRc variables, and 57 double mode RR Lyr stars (RRd). Additionally we attach a list of variables with periods between 0.18-0.26 days, which are probably delta Sct stars. Period, BVI photometry, astrometry, amplitude, and parameters of the Fourier decomposition of the I-band light curve are provided for each object. We also present the Petersen diagram for double mode pulsators. We found that the SMC RR Lyr stars are fairly uniformly distributed over the studied area, with no clear correlation to any location. The most preferred periods for RRab and RRc stars are 0.589 and 0.357 days, respectively. We noticed significant excess of stars with periods of about 0.29 days, which might be second-overtone RR Lyr stars (RRe). The mean extinction free magnitudes derived for RRab stars are 18.97, 19.45 and 19.73 mag for the I, V and B-band, respectively. All presented data, including individual BVI observations, are available from the OGLE Internet archive.

The mining of Virtual Observatories (VOs) is becoming a powerful new method for discovery in astronomy. Here we report on the development of SkyDOT (Sky Database for Objects in the Time domain), a new Virtual Observatory, which is dedicated to the study of sky variability. The site will confederate a number of massive variability surveys and enable exploration of the time domain in astronomy. We discuss the architecture of the database and the functionality of the user interface. An important aspect of SkyDOT is that it is continuously updated in near real time so that users can access new observations in a timely manner. The site will also utilize high level machine learning tools that will allow sophisticated mining of the archive. Another key feature is the real time data stream provided by RAPTOR (RAPid Telescopes for Optical Response), a new sky monitoring experiment under construction at Los Alamos National Laboratory (LANL).

We present 3 years of photometry of the ``Double Hamburger'' lensed quasar, HE1104-1805, obtained on 102 separate nights using the OGLE 1.3-m telescope. Both the A and B images show variations, but with substantial differences in the lighcurves at all time delays. At the 310 day delay reported by Wisotzki and collaborators the difference lightcurve has an rms amplitude of 0.060 mag. The structure functions for the A and B images are quite different, with image A more than twice as variable as image B (a factor of 4 in structure function) on timescales of less than a month. Adopting microlensing as a working hypothesis for the uncorrelated variability, the short timescale argues for the relativistic motion of one or more components of the source. We argue that the small amplitude of the fluctuations is due to the finite size of the source with respect to the microlenses.

We report a measurement of the streaming motion of the stars in the Galactic bar with the Red Clump Giants (RCGs) using the data of the Optical Gravitational Lensing Experiment II (OGLE-II). We measure the proper motion of 46,961 stars and divide RCGs into bright and faint sub-samples which on average will be closer to the near and far side of the bar, respectively. We find that the far-side RCGs (4,979 stars) have a proper motion of \Delta <\mu> ~ 1.5 +- 0.11 mas yr^{-1} toward the negative l relative to the near-side RCGs (3,610 stars). This result can be explained by stars in the bar rotating around the Galactic center in the same direction as the Sun with v_b ~ 100 km s^{-1}. In the Disc Star (DS) and Red Giant (RG) samples, we do not find significant difference between bright and faint sub-samples. For those samples \Delta <\mu> \~ 0.3 +- 0.14 mas yr^{-1} and ~ 0.03 +- 0.14 mas yr^{-1}, respectively. It is likely that the average proper motion of RG stars is the same as that of the Galactic center. The proper motion of DSs with respect to RGs is ~ 3.3 mas yr^{-1} toward positive l. This value is consistent with the expectations for a flat rotation curve and Solar motion with respect to local standard of rest. RGs have proper motion approzimately equal to the average of bright and faint RCGs, which implies that they are on average near the center of the bar. This pilot project demonstrates that OGLE-II data may be used to study streaming motions of stars in the Galactic bar. We intend to extend this work to all 49 OGLE-II fields in the Galactic bulge region.

We present the VI photometric maps of the Galactic bulge. They contain VI photometry and astrometry of about 30 million stars from 49 fields of 0.225 square degree each in the Galactic center region. The data were collected during the second phase of the OGLE microlensing project. We discuss the accuracy of data and present color-magnitude diagrams of selected fields observed by OGLE in the Galactic bulge. The VI maps of the Galactic bulge are accessible electronically for the astronomical community from the OGLE Internet archive.

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

Fast variability of optical objects is an interesting though poorly explored subject in modern astronomy. Real-time data processing and identification of transient celestial events in the images is very important for such study as it allows rapid follow-up with more sensitive instruments. We discuss an approach which we have developed for the RAPTOR project, a pioneering closed-loop system combining real-time transient detection with rapid follow-up. RAPTOR's data processing pipeline is able to identify and localize an optical transient within seconds after the observation. The testing we performed so far have been confirming the effectiveness of our method for the optical transient detection. The software pipeline we have developed for RAPTOR can easily be applied to the data from other experiments.

The Rapid Telescope for Optical Response (RAPTOR) program consists of a network of robotic telescopes dedicated to the search for fast optical transients. The pilot project is composed of three observatories separated by approximately 38 kilometers located near Los Alamos, New Mexico. Each of these observatories is composed of a telescope, mount, enclosure, and weather station, all operating robotically to perform individual or coordinated transient searches. The telescopes employ rapidly slewing mounts capable of slewing a 250 pound load 180 degrees in under 2 seconds with arcsecond precision. Each telescope consists of wide-field cameras for transient detection and a narrow-field camera with greater resolution and sensitivity. The telescopes work together by employing a closed-loop system for transient detection and follow-up. Using the combined data from simultaneous observations, transient alerts are generated and distributed via the Internet. Each RAPTOR telescope also has the capability of rapidly responding to external transient alerts received over the Internet from a variety of ground-based and satellite sources. Each observatory may be controlled directly, remotely, or robotically while providing state-of-health and observational results to the client and the other RAPTOR observatories. We discuss the design and implementation of the spatially distributed RAPTOR system.

The Rapid Telescopes for Optical Response (RAPTOR) experiment is a spatially distributed system of autonomous robotic telescopes that is designed to monitor the sky for optical transients. The core of the system is composed of two telescope arrays, separated by 38 kilometers, that stereoscopically view the same 1500 square-degree field with a wide-field imaging array and a central 4 square-degree field with a more sensitive narrow-field "fovea" imager. Coupled to each telescope array is a real-time data analysis pipeline that is designed to identify interesting transients on timescales of seconds and, when a celestial transient is identified, to command the rapidly slewing robotic mounts to point the narrow-field ``fovea'' imagers at the transient. The two narrow-field telescopes then image the transient with higher spatial resolution and at a faster cadence to gather light curve information. Each "fovea" camera also images the transient through a different filter to provide color information. This stereoscopic monitoring array is supplemented by a rapidly slewing telescope with a low resolution spectrograph for follow-up observations of transients and a sky patrol telescope that nightly monitors about 10,000 square-degrees for variations, with timescales of a day or longer, to a depth about 100 times fainter. In addition to searching for fast transients, we will use the data stream from RAPTOR as a real-time sentinel for recognizing important variations in known sources. Altogether, the RAPTOR project aims to construct a new type of system for discovery in optical astronomy--one that explores the time domain by "mining the sky in real time".

We describe a highly unusual microlensing event, OGLE-1999-BUL-19, which exhibits multiple peaks in its light curve. The Einstein radius crossing time for this event is approximately one year, which is unusually long. We show that the motion of the Earth induces these multiple peaks in the light curve, since the relative transverse velocity of the lens projected into the observer plane is very small (v = 12.5 km/s). This is the lowest velocity so far published and we believe that this is the first multiple-peak parallax event ever observed. We also believe that this event may be exhibiting slight binary-source signatures in addition to these parallax-induced multiple peaks. With spectroscopic observations it is possible to test this `parallax plus binary-source' hypothesis and (if this hypothesis turns out to be correct) to simultaneously fit both models and obtain a measurement of the lens mass. Furthermore, spectroscopic observations could also supply information regarding the lens properties, possibly providing another avenue for determining the lens mass. We found that most of the I-band blending is probably caused by light from the lens or a binary companion to the source. However, in the V-band, there appears to be a second blended source 0.35" away from the lensed source. HST observations will be very useful for understanding the nature of the blends. We also suggest that a radial velocity survey of all parallax events will be very useful for further constraining the lensing kinematics and understanding the origins of these events and the excess of long events toward the bulge.

We present a Catalog of high proper motion (HPM) stars detected in the foreground of central parts of the Magellanic Clouds. The Catalog contains 2161 objects in the 4.5 square degree area towards the LMC, and 892 HPM stars in the 2.4 square degree area towards the SMC. The Catalog is based on observations collected during four years of the OGLE-II microlensing survey. The Difference Image Analysis (DIA) of the images provided candidate HPM stars with proper motion as small as 4 mas/yr. These appeared as pseudo-variables, and were all measured astrometrically on all CCD images, providing typically about 400 data points per star. The reference frame was defined by the majority of background stars, most of them members of the Magellanic Clouds. The reflex motion due to solar velocity with respect to the local standards of rest is clearly seen. The largest proper motion in our sample is 363 mas/yr. Parallaxes were measured with errors smaller than 20% for several stars.

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

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

We present the first edition of a catalog of variable stars from OGLE-II Galactic Bulge data covering 3 years: 1997-1999. Typically 200-300 I band data points are available in 49 fields between -11 and 11 degrees in galactic longitude, totaling roughly 11 square degrees in sky coverage. Photometry was obtained using the Difference Image Analysis (DIA) software and tied to the OGLE data base with the DoPhot package. The present version of the catalog comprises 221,801 light curves. In this preliminary work the level of contamination by spurious detections is still about 10%. Parts of the catalog have only crude calibration, insufficient for distance determinations. The next, fully calibrated, edition will include the data collected in year 2000. The data is accessible via FTP. Due to the data volume, we also distribute DAT tapes upon request.

We present the first edition of a catalog of variable stars found in the Magellanic Clouds using OGLE-II data obtained during four years: 1997-2000. The catalog covers about 7 square degrees of the sky - 21 fields in the Large Magellanic Cloud and 11 fields in the Small Magellanic Cloud. All variables were found with the Difference Image Analysis (DIA) software. The catalog is divided into two sections. The DC section contains FITS reference images (obtained by co-adding 20 best frames for each field) and profile photometry (DoPHOT) of all variable stars on those images. The AC section contains flux variations and magnitudes of detected variable stars obtained with DIA as well as with DoPHOT. The results are available in electronic form via FTP and through WWW interface from the OGLE internet archive. Typically, there are about 400 I-band data points and about 30 V and B-band data points for over 68,000 variables. A detailed analysis and classification of variable stars will be presented elsewhere. The FTP catalog contains approximately 2Gbytes of data.

We describe the Difference Image Analysis (DIA) algorithms and software used to analyze four years (1997-2000) of OGLE-II photometric monitoring of the Magellanic Clouds, the calibration, the photometric error analysis, and the search for variable stars. A preliminary analysis of photometric errors is based on the field LMC_SC2. A full catalog of about 70,000 variable stars will be presented in a separate publication.

We describe the discovery of the longest microlensing event ever observed, OGLE-1999-BUL-32, also independently identified by the MACHO collaboration as MACHO-99-BLG-22. This unique event has an Einstein radius crossing time of 641 days. The high quality data obtained with difference image analysis shows a small but significant parallax signature. This parallax effect allows one to determine the Einstein radius projected onto the observer plane as rE^hat ~ 29.2AU. The transverse velocity projected onto the observer plane is about 79km/s. We argue that the lens is likely to be have a mass of at least a few solar masses, i.e., it could be a stellar black hole. The black hole hypothesis can be tested using the astrometric microlensing signature with the soon-to-be installed Advanced Camera for Surveys on board the Hubble Space Telescope. Deep X-ray and radio images may also be useful for revealing the nature of the object.

We present a systematic search for parallax microlensing events among a total of 512 microlensing candidates in the OGLE II database for the 1997-1999 seasons. We fit each microlensing candidate with both the standard microlensing model and also a parallax model that accounts for the Earth's motion around the Sun. We then search for the parallax signature by comparing the chi^2 of the standard and parallax models. For the events which show a significant improvement, we further use the `duration' of the event and the signal-to-noise ratio as criteria to separate true parallax events from other noisy microlensing events. We have discovered one convincing new candidate, sc33_4505, and seven other marginal cases. The convincing candidate (sc33_4505) is caused by a slow-moving, and likely low-mass, object, similar to other known parallax events. We found that irregular sampling and gaps between observing seasons hamper the recovery of parallax events. We have also searched for long-duration events that do not show parallax signatures. The lack of parallax effects in a microlensing event puts a lower-limit on the Einstein radius projected onto the observer plane, which in turn imposes a lower limit on the lens mass divided by the relative lens-source parallax. Most of the constraints are however quite weak.

We present a sample of microlensing events discovered in the Difference Image Analysis (DIA) of the OGLE-II images collected during 3 observing seasons, 1997--1999. 4424 light curves pass our criteria on the presence of a brightening episode on top of a constant baseline. Among those, 512 candidate microlensing events were selected visually. We designed an automated procedure, which unambiguously selects up to 237 best events. Including 8 candidate events recovered by other means, a total of 520 light curves are presented in this work. In addition to microlensing events, the larger sample contains certain types of transients, but is also strongly contaminated by artifacts. All 4424 light curves in the weakly filtered group are available electronically, with the intent of showing the gray zone between microlensing events and variable stars, as well as artifacts, to some extent inevitable in massive data reductions. We welcome suggestions for improving the selection process before the full analysis of complete 4 seasons of the OGLE-II bulge data. Selection criteria for binary events can also be investigated with our extended sample.

The Difference Image Analysis (DIA) of the images obtained by the Optical Gravitational Lensing Experiment (OGLE-II) revealed a peculiar artifact in the sample of stars proposed as variable by Wozniak (2000) in one of the Galactic bulge fields: the occurrence of pairs of candidate variables showing anti-correlated light curves monotonic over a period of 3 years. This effect can be understood, quantified and related to the stellar proper motions. DIA photometry supplemented with a simple model offers an effective and easy way to detect high proper motion stars (HPM stars) in very dense stellar fields, where conventional astrometric searches are extremely inefficient.

We present a sample of 76 double-mode Cepheids detected in the 4.5 square degree area in the central part of the LMC. 19 stars from the sample pulsate in the fundamental mode and the first overtone while 57 objects are the first and second overtone pulsators. We analyze the period ratio of double-mode Cepheids and Fourier parameters of decomposition of the light curves of these objects. We also present location of different type Cepheids from the LMC in the color-magnitude diagram and show the distribution of their ${V-I}$ color indices.

We present the photometry and theoretical models for a Galactic bulge microlensing event OGLE-2000-BUL-43. The event is very bright with I=13.54 mag, and has a very long time scale, t_E=156 days. The long time scale and its light curve deviation from the standard shape strongly suggest that it may be affected by the parallax effect. We show that OGLE-2000-BUL-43 is the first discovered microlensing event, in which the parallax distortion is observed over a period of 2 years. Difference Image Analysis (DIA) using the PSF matching algorithm of Alard & Lupton enabled photometry accurate to 0.5%. All photometry obtained with DIA is available electronically. Our analysis indicates that the viewing condition from a location near Jupiter will be optimal and can lead to magnifications ~ 50 around January 31, 2001. These features offer a great promise for resolving the source (a K giant) and breaking the degeneracy between the lens parameters including the mass of the lens, if the event is observed with the imaging camera on the Cassini space probe.

I present an implementation of the difference image photometry based on the Alard & Lupton optimal PSF matching algorithm. The most important feature distinguishing this method from the ones using Fourier divisions is that equations are solved in real space and the knowledge of each PSF is not required for determination of the convolution kernel. I evaluate the method and software on 380 GB of OGLE-II bulge microlensing data obtained in 1997--1999 observing seasons. The error distribution is Gaussian to better than 99% with the amplitude only 17% above the photon noise limit for faint stars. Over the entire range of the observed magnitudes the resulting scatter is improved by a factor of 2--3 compared to DoPhot photometry, currently a standard tool in microlensing searches. For testing purposes the photometry of ~4600 candidate variable stars and sample difference image data are provided for BUL_SC1 field. In the candidate selection process, very few assumptions have been made about the specific types of flux variations, which makes this data set well suited for general variability studies, including the development of the classification schemes.

We present the BVI photometric maps of the Large Magellanic Cloud. They contain BVI photometry and astrometry of more than 7 million stars from the central parts of the LMC. The data were collected during the second phase of the OGLE microlensing project. We discuss the accuracy of the data and present color-magnitude diagrams of all 26 fields observed by OGLE in the LMC. The BVI maps of the LMC are accessible electronically for the astronomical community from the OGLE Internet archive.

To measure interstellar atomic deuterium abundances, we used the Interstellar Medium Absorption Profile Spectrograph (IMAPS) to obtain spectra of gamma2 Vel and zeta Pup over the wavelength interval 930-1150 A at a resolving power of 80,000. The interstellar D I features are resolved and cleanly separated from interstellar H I in the Ly-delta and Ly-epsilon profiles of both sight lines, and also in the Ly-gamma profile of zeta Pup. The D I profiles were modeled using a velocity template derived from several N I lines in the IMAPS spectra recorded at higher S/N. To find the best D I column density, we minimized chi-squares for model D I profiles that included not only the N(D I) as a free parameter, but also the effects of several potential sources of systematic error which could also be varied. For both stars, H I column densities were measured by analyzing Ly-alpha absorption profiles in a large number of IUE high dispersion spectra. Ultimately we found that D/H = 2.18(+0.36,-0.31)e-5 for gamma2 Vel and 1.42(+0.25,-0.23)e-5 for zeta Pup, values that contrast markedly with D/H derived in Paper I for delta Ori (the stated errors are 90% confidence limits). Evidently, the atomic D/H ratio in the ISM, averaged over path lengths of 250 to 500 pc, exhibits significant spatial variability. Furthermore, variations in D/H do not appear to be anticorrelated with N/H. Within the framework of standard Big Bang Nucleosynthesis, the large value of D/H found toward gamma2 Vel is equivalent to a cosmic baryon density of Omega_B h^2 = 0.023 (+-0.002), which we regard as an upper limit since there is no correction for the destruction of deuterium in stars.

In 1998 The Optical Gravitational Lensing Experiment (OGLE) successfully implemented automated data reductions for QSO 2237+0305. Using a new image subtraction method we achieved a differential photometry scatter of 1-5 % for images A-D respectively. Combined with a time sampling of 1-2 times a week this is sufficient for early detection of caustic crossings. Nearly real time photometry of QSO 2237+0305 is available from the OGLE web site http://www.astro.princeton.edu/~ogle/ogle2/huchra.html . During the 1999 observing season, the apparent V magnitude of the A, B, C and D images changed by 0.50, 0.15, 0.65 and 0.35 mag, respectively. Most likely however, none of the microlensing events involved a caustic crossing. The most rapid variation was 0.25 mag in 30 days, observed for image C. The alert system will continue to be active in the next observing season from late April until September 2000, when OGLE suspends operation for an upgrade. Observations will resume for season of 2001.

Absorption features from the Lyman and Werner bands of interstellar molecular hydrogen were recorded by the Interstellar Medium Absorption Profile Spectrograph (IMAPS) at a wavelength resolving power R=80,000 in the spectra of delta Ori and epsilon Ori. The objective was to find and study more examples of an unusual phenomenon found for one of the velocity components of H2 in the spectrum of zeta Ori by Jenkins & Peimbert (1997). Specifically, they detected a gradual shift in velocity and broadening for features arising from progressively higher rotational excitations J. This effect appears to be absent in the spectra of both delta and epsilon Ori, which are only a few degrees away in the sky from zeta Ori. The absence of atomic material at a large negative velocity in the spectra of delta and epsilon Ori (and its presence in zeta Ori) supports a proposal by Jenkins & Peimbert that the line of sight to zeta intercepts a bow shock facing away from us, perhaps created by the collision of wind-like material with some foreground obstruction. For both stars, the H2 absorption features are separated into two velocity components. Total H2 column densities toward delta and epsilon Ori are 5.5e14 and 1.9e16 cm^-2, respectively. The rotation temperatures of the molecules with J > 2 toward epsilon Ori indicate that the gas is in the general vicinity of the stars that emit UV fluxes capable of rotationally pumping the molecules. For the strongest component of H2 toward delta Ori, the pumping rate is lower and consistent with a general UV flux level in the plane of the Galaxy.

We present the Catalog of microlensing events detected toward the Galactic bulge in three observing seasons, 1997--1999, during the OGLE-II microlensing survey. The search for microlensing events was performed using a database of about 4*10^9 photometric measurements of about 20.5 million stars from the Galactic bulge. The Catalog comprises 214 cases of microlensing events found in the fields covering about 11 square degrees on the sky and distributed in different parts of the Galactic bulge. The sample includes 20 binary microlensing events, 14 of them are caustic crossing. In one case a double star is likely lensed. We present distribution of the basic parameters of microlensing events and show preliminary rate of microlensing in different regions of the Galactic bulge. The latter reveals clear dependence on the Galactic coordinates. The dependence on 'l' indicates that the majority of lenses toward the Galactic bulge are located in the Galactic bar. Models of the Galactic bar seem to reasonably predict the observed spatial distribution of microlensing events in the Galactic bulge. All data presented in the Catalog and photometry of all events are available from the OGLE Internet archive.

We present the catalog of star clusters found in the area of about 5.8 square degree in the central regions of the Large Magellanic Cloud. It contains data for 745 clusters. 126 of them are new objects. For each cluster equatorial coordinates, radius, approximate number of members and cross-identification are provided. Photometric data for all clusters presented in the catalog and Atlas consisting of finding charts and color-magnitude diagrams are available electronically from the OGLE Internet archive.

We present the Catalog of Cepheids from the SMC which contains data for 2049 objects detected in the 2.4 square degree area of central parts of the SMC. For each object period, BVI photometry, astrometry, and R_21, \phi_21 parameters of the Fourier decomposition of I-band light curve are provided. The Catalog is based on observations collected during the OGLE-II microlensing survey. Tests of completeness performed in overlapping parts of adjacent fields indicate that completeness of the Catalog is very high: about 92%. Statistics and distributions of basic parameters of Cepheids are also presented. All presented data, including individual BVI observations (about 4.7*10^5 BVI measurements), are available from the OGLE Internet archive.

We present the Catalog of Cepheids from the LMC. The Catalog contains 1333 objects detected in the 4.5 square degree area of central parts of the LMC. About 3.4*10^5 BVI measurements of these stars were collected during the OGLE-II microlensing survey. The Catalog data include period, BVI photometry, astrometry, and R_21, \phi_21 parameters of the Fourier decomposition of I-band light curve. The vast majority of objects from the Catalog are the classical Cepheids pulsating in the fundamental or first overtone mode. The remaining objects include Population II Cepheids and red giants with pulsation-like light curves. Tests of completeness performed in overlapping parts of adjacent fields indicate that completeness of the Catalog is very high: >96%. Statistics and distributions of basic parameters of Cepheids are also presented. Finally, we show the light curves of three eclipsing systems containing Cepheid detected among objects of the Catalog. All presented data, including individual BVI observations are available from the OGLE Internet archive.

We present Period-Luminosity-Color and Period-Luminosity relations of classical Cepheids constructed for about 1280 Cepheids from the LMC and 2140 from the SMC. High quality BVI observations (120-360 epochs in the I-band and 15-40 in the BV-bands) were collected during the OGLE-II microlensing experiment. The I-band diagrams of the LMC show very small scatter, \sigma=0.074 mag, indicating that Cepheid variables can potentially be a very good standard candle. We compare relations of fundamental mode Cepheids from the LMC and SMC and we do not find significant differences of slopes of the Period-Luminosity- Color and Period-Luminosity relations in these galaxies. For the first overtone Cepheids a small change of the slope of Period-Luminosity relation is possible. We determine the difference of distance moduli between the SMC and LMC with Cepheid relations and compare the result with difference obtained with other standard candles: RR Lyr and red clump stars. Results are very consistent and indicate that the values of zero points of the fundamental mode Cepheid relations are similar in these galaxies. The mean difference of distance moduli between the SMC and LMC is equal to \mu_SMC-\mu_LMC=0.51\pm0.03 mag. We calibrate the Period-Luminosity-Color and Period-Luminosity relations for classical, fundamental mode Cepheids using the observed LMC relations and adopting the short LMC distance modulus, \mu_LMC=18.22\pm0.05 mag, resulting from the recent determination with eclipsing system HV2274, RR Lyr and red clump stars and consistent with observations of Cepheids in NGC4258 galaxy, to which precise geometric distance is known. Finally, we determine a constraint on the absolute magnitude of Cepheids by comparison of their mean V-band magnitude with that of RR Lyr stars in both Magellanic Clouds.

We present UBVI photometry for 8530 stars in Baade's Window obtained during the OGLE-II microlensing survey. Among these are over one thousand red clump giants. 1391 of them have photometry with errors smaller than 0.04, 0.06, 0.12, and 0.20 mag in the I, V, B and U-band, respectively. We constructed a map of interstellar reddening. The corrected colors of the red clump giants: (U-B)_0, (B-V)_0, and (V-I)_0 are very well correlated, indicating that a single parameter determines the observed spread of their values, reaching almost 2 mag in the (U-B)_0. It seems most likely that heavy element content is the dominant parameter, but it is possible that another parameter: the age (or mass) of a star moves it along the same trajectory in the color--color diagram as the metallicity. The current ambiguity can be resolved with spectral analysis, and our catalog may be useful as a finding list of red clump giants. We point out that these K giants are more suitable for a fair determination of the distribution of metallicity than brighter M giants. We also present a compilation of UBVI data for 308 red clump giants near the Sun, for which Hipparcos parallaxes are more accurate than 10%. Spectral analysis of their metallicity may provide information about the local metallicity distribution as well as the extent to which mass (age) of these stars affects their colors. It is remarkable that in spite of a number of problems, stellar models agree with observations at the 0.1-0.2 mag level, making red clump giants not only the best calibrated but also the best understood standard candle.

We present $V$-band light curves of 61 variables from the core of the globular cluster M5 obtained using a newly developed image subtraction method (ISM). Four of these variables were previously unknown. Only 26 variables were found in the same field using photometry obtained with DoPHOT software. Fourier parameters of the ISM light curves have relative errors up to 20 times smaller than parameters measured from DoPHOT photometry. We conclude that the new method is very promising for searching for variable stars in the cores of the globular clusters and gives very accurate relative photometry with quality comparable to photometry obtained by HST. We also show that the variable V104 is not an eclipsing star as has been suggested, but is an RRc star showing non-radial pulsations.

We present results from 2 years of monitoring of Huchra's lens (QSO 2237+0305) with the 1.3 m Warsaw telescope on Las Campanas, Chile. Photometry in the V band was done using a newly developed method for image subtraction. Reliable subtraction without Fourier division removes all complexities associated with the presence of a bright lensing galaxy. With positions of lensed images adopted from HST measurements it is relatively easy to fit the variable part of the flux in this system, as opposed to modeling of the underlying galaxy. For the first time we observed smooth light variation over a period of a few months, which can be naturally attributed to microlensing. We also describe automated software capable of real time analysis of the images of QSO 2237+0305. It is expected that starting from the next observing season in 1999 an alert system will be implemented for high amplification events (HAE) in this object. Time sampling and photometric accuracy achieved should be sufficient for early detection of caustic crossings.

We present a sample of 93 double-mode Cepheids detected in the 2.4 square degree area in the central part of the SMC. 23 stars from the sample pulsate in the fundamental mode and the first overtone while 70 objects are the first and second overtone pulsators. This is the largest sample of such type Cepheids detected in one environment so far. We analyze period ratio of double-mode Cepheids and Fourier parameters of decomposition of the light curves of these objects. We also present location of different type Cepheids from the SMC in the color-magnitude diagram and show their distribution of V-I color indices. We find one object which is probably a blend, either physical or optical, of two Cepheids pulsating in the fundamental mode.

We present a sample of 13 candidates for single-mode second overtone Cepheids discovered in the 2.4 square degree area of the SMC bar during the OGLE microlensing search. All candidates passed four photometric tests: analysis of Fourier parameters of the light curve, R_21 and \phi_21, position in the period-luminosity diagram and location in the color-magnitude diagram. The stars form a very homogeneous group of objects with very similar light curves of low amplitude and almost sinusoidal shape. Their color indices B-V and V-I are also very similar and the objects determine the observational blue edge of the instability strip of the first overtone Cepheids. This is the first firm sample of the single-mode second overtone Cepheids detected so far.

Studies of the abundances of deuterium in different astrophysical sites are of fundamental importance to answering the question about how much deuterium was produced during big bang nucleosynthesis and what fraction of it was destroyed later. With this in mind, we used the Interstellar Medium Absorption Profile Spectrograph (IMAPS) on the ORFEUS-SPAS II mission to observe at a wavelength resolution of 4 km/s (FWHM) the L-delta and L-epsilon absorption features produced by interstellar atomic deuterium in the spectrum of delta Ori A. A chi-square analysis indicated that 0.96 < N(D I)< 1.45e15 cm^{-2} at a 90% level of confidence, and the gas is at a temperature of about 6000K. To obtain an accurate value of N(H I) needed for a determination of the atomic ratio of D to H, we measured the L-alpha absorption features in 57 spectra of delta Ori in the IUE archive. From our measurement of N(H I)= 1.56e20 cm^{-2}, we found that N(D I)/N(H I)= 7.4(+1.9,-1.3)e-6 (90% confidence). Our result for D/H contrasts with the more general finding along other lines of sight that D/H is approximately 1.5e-5. The underabundance of D toward delta Ori A is not accompanied by an overabundance of N or O relative to H, as one might expect if the gas were subjected to more stellar processing than usual.

We present the catalog of 1459 eclipsing binary stars detected in the central 2.4 square degree area of the Small Magellanic Cloud during the OGLE-II microlensing search. The sample includes objects brighter than I=20 mag with periods ranging from about 0.3 to 250 days. The average completeness of the catalog is about 80%. Statistics of the entire sample and well detached systems, suitable for distance determination, are also presented. The catalog, finding charts and BVI photometry data for all objects are available from the OGLE Internet archive.

We present {\it UBVI} photometry of the eclipsing binary HV2274 - the system which has been recently used for distance determination to the LMC by Guinan et al. (1998). We determine the interstellar reddening to the star, E(B-V)=0.149+/-0.015 mag, based on observed colors of the star. This value is in excellent agreement with the mean reddening towards HV2274 obtained from photometry of the red clump stars in the surrounding field. The reddening is almost twice as large as determined by Guinan et al. (1998). We discuss the consequences of reddening underestimate. Most likely HV2274 is located much closer with the distance modulus to the star and the LMC: m-M = 18.22+/-0.13 mag supporting the short distance scale to the LMC. Such a distance modulus is in excellent agreement with the recent distance determinations with RR Lyr and red clump stars.

We present photometry of the unique binary microlensing event MACHO-SMC-98-1 collected by the OGLE group. Particularly interesting observation was made close to the first caustic crossing which was not covered by observations of other groups. It allows to test proposed models of which Model~1 proposed by PLANET group seems to be in the best agreement with the OGLE observations.

We present three color, BVI maps of the Small Magellanic Cloud. The maps contain precise photometric and astrometric data for about 2.2 million stars from the central regions of the SMC bar covering ~2.4 square degrees on the sky. Mean brightness of stars is derived from observations collected in the course of the OGLE-II microlensing search from about 130, 30 and 15 measurements in the I, V and B-bands, respectively. Accuracy of the zero points of photometry is about 0.01 mag, and astrometry 0.15 arcsec (with possible systematic error up to 0.7 arcsec). Color-magnitude diagrams of observed fields are also presented. The maps of the SMC are the first from the series of similar maps covering other OGLE fields: LMC, Galactic bulge and Galactic disk. The data are very well suited for many projects, particularly for the SMC which has been neglected photometrically for years. Because of potentially great impact on many astrophysical fields we decided to make the SMC data available to the astronomical community from the OGLE Internet archive.

We present the catalog of clusters found in the area of ~2.4 square degree in the central regions of the Small Magellanic Cloud. The catalog contains data for 238 clusters, 72 of them are new objects. For each cluster equatorial coordinates, radii, approximate number of members, cross-identification, finding chart and color magnitude diagrams: V-(B-V) and V-(V-I) are provided. Photometric data for all clusters presented in the catalog are available from the OGLE Internet archive.

We study X-ray and soft gamma-ray spectral properties of BLRGs using data from Ginga, ASCA, OSSE and EXOSAT. The X-ray spectra are well fitted by an intrinsic power-law continuum with an energy index of alpha ~ 0.7, moderately absorbed by a cold medium. In addition, the Ginga spectra show Fe K-alpha lines with an average equivalent width of ~ 100 eV, and, in some cases, Compton reflection humps. However, the latter components are significantly weaker than both those seen in radio-quiet Seyferts and those expected if the Fe K-alpha lines were due to reflection. Some ASCA and EXOSAT spectra show soft X-ray excesses below ~ 3 keV. When that component is taken into account, the Fe K-alpha lines in the ASCA data become unresolved with equivalent widths < 200 eV. Multiple observations of 3C 390.3 and 3C 382 show the Fe K-alpha line approximately constant in flux but accompanied by strong continuum variations. This indicates the bulk of the line is formed by matter at a distance much larger than an accretion-disk scale, consistent with the ASCA line width measurements. The column density of the matter required to account for the observed line fluxes is N_H > 10^23 cm^-2. Such a medium is in the line-of-sight in 3C 445 but it has to be out of it in other objects, in which the observed N_H are substantially lower. Thus, a cold medium with that N_H and covering a large solid angle is out of the line-of-sight in most objects. The spectra of BLRGs break and become softer above ~ 100 keV, as shown by a simultaneous ASCA/OSSE observation of 3C 120 and by the OSSE spectra being on average much softer than the X-ray spectra. Finally, we find the X-ray and gamma-ray spectral properties of Cen A intrinsically very similar to BLRGs studied here.

We search OGLE-I photometric database for stars, which, as defined by formal criteria adopted by OGLE-I microlensing search, showed variability during only one out of 3 or 4 observing seasons. The results include 17 previously reported microlensing events, 2 newly discovered candidate events and 15 intrinsically variable stars that have a potential of contaminating samples of microlensing events. Based on photometry obtained in 1992 and 1993 OGLE #10 was tentatively included in the list of microlensing candidates, however its light curve in 1994 and 1995 shows many characteristics of the variable stars found in our search, and most likely it is not a microlensing event. For all stars which passed our tests, we provide 44 x 44 arcsec (101 x 101 pixels) centered subframes from each OGLE-I frame in I band. It is the first time when images used to derive photometry of microlensing events are available in convenient format to astronomical community.

We present a new distance determination to the Large and Small Magellanic Clouds using the newly developed red clump stars method (Paczynski and Stanek 1998). This new, single-step, Hipparcos calibrated method seems to be one of the most precise techniques of distance determination with very small statistical error due to large number of red clump stars usually available. The distances were determined independently along four lines-of-sight located at opposite sides of each Magellanic Cloud. The results for each line-of-sight are very consistent. For the SMC we obtain the distance modulus: m-M=18.56+/-0.03+/-0.06 mag (statistical and systematic errors, respectively) and for the LMC: m-M=18.08+/-0.03+/-0.12} mag where systematic errors are mostly due to uncertainty in reddening estimates. Both distances will be refined and systematic errors reduced when accurate reddening maps for our fields are available. Distance moduli to both Magellanic Clouds are ~0.4 mag smaller than generally accepted values. The modulus to the LMC is in good agreement with the recent determinations from RR Lyrae type stars and upper limit resulting from the SN1987A echo. We suspect that the distance to the LMC and SMC is shorter by about 15% than previously assumed: 42 kpc and 52 kpc, respectively. We also present our color-magnitude diagrams around the red clump for the LMC and SMC. We identify vertical red clump, first noted by Zaritsky and Lin (1997), in the color-magnitude diagram of both Magellanic Clouds and we interpret it as an evolutionary feature rather than unknown stellar population between the LMC and our Galaxy.

We measure the astrometric shifts of the light centroids for the microlensing events in the OGLE-1 database. Of the 15 consistently detected events, 7 had significant shifts which we were able to measure with a fair degree of confidence. Those events with large shifts are also expected to be unresolved with a background ``blend'', and thus, we suggest that we have identified events in the OGLE-1 catalog which are strongly blended. Though we concentrate on the OGLE-1 database, and use the DoPHOT package in order to perform our analysis, we suggest that this shift is a generic effect, and should be observable in any crowded field.

We present photometric observations of the MACHO-SMC-1 microlensing candidate collected by the OGLE-2 project. We show light curves of both components of the 1.6 arcsec blend: microlensed star and its optical companion. We find the contribution of the optical companion to the total flux to be 24% and confirm presence of the small amplitude periodic oscillations in the light curve of the lensed star with the period of 5.096 days and amplitude 0.05 mag. The lensed star is probably an ellipsoidal binary system.

The current modelling of single microlensing light curves neglects the possibility that only a fraction of the light is due to the lensed star, the remaining being due to a close, unresolved blend, which may be related or unrelated to the lens. Unfortunately, the effects of blending are significant as all microlensing experiments choose very crowded fields as their targets. In this paper we point out a strong degeneracy of the fitting procedure which makes it practically impossible to detect the presence of a blend by purely photometric means, except in a small part of the parameter space. Some blends may be detected by astrometric means, but the majority have to be corrected for statistically. The luminosity function reaching well below the ground based detection limit (with the HST) would be very also important and this could be determined with the repeating microlensing events. If no correction is made then the event time scales, the lens masses, and the optical depth are all systematically underestimated.

We propose a new method to investigate the coefficient of the selective extinction, based on two band photometry. This method uses red clump stars as a means to construct the reddening curve. We apply this method to the OGLE color-magnitude diagrams to investigate the variations of the selective extinction towards various parts of the Galactic bulge. We find that $A_{_V}/E_{_{V-I}}$ coefficient is within the errors the same for $l=\pm 5\deg$ OGLE fields. Therefore, the difference of $\sim 0.37\;mag$ in the extinction adjusted apparent magnitude of the red clump stars in these fields (Stanek et al.~1994, 1995) cannot be assigned to a large-scale gradient of the selective extinction coefficient. This strengthens the implication of this difference as indicator of the presence of the bar in our Galaxy. However using present data we cannot entirely exclude the possibility of $\sim 0.2\;mag/mag$ variations of the selective extinction coefficient on the large scales across the bulge.