We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.

We explore a unique electromagnetic signature of stellar-mass compact-object binaries long before they are detectable in gravitational waves. We show that gravitational lensing of light emitting components of a compact-object binary, by the other binary component, could be detectable in the nearby universe. This periodic lensing signature could be detected from present and future X-ray observations, identifying the progenitors of binaries that merge in the LIGO band, and also unveiling populations that do not merge, thus providing a tracer of the compact-object binary population in an enigmatic portion of its life. We argue that periodically repeating lensing flares could be observed for ~100 ksec orbital-period binaries with the future Lynx X-ray mission, possibly concurrent with gravitational wave emission in the LISA band. Binaries with longer orbital periods could be more common and be detectable as single lensing flares, though with reliance on a model for the flare that can be tested by observations of succeeding flares. Non-detection of such events, even with existing X-ray observations, will help to constrain the population of EM bright compact-object binaries.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

We investigated the time-variability and spectral properties of the eclipsing X-ray source Circinus Galaxy X-1 (GG X-1), using Chandra, XMM-Newton and ROSAT. We phase-connected the lightcurves observed over 20 years, and obtained a best-fitting period $P = (25,970.0 \pm 0.1)$ s $\approx$7.2 hr, and a period derivative $\dot{P}/P = ( 10.2\pm4.6) \times 10^{-7}$ yr$^{-1}$. The X-ray lightcurve shows asymmetric eclipses, with sharp ingresses and slow, irregular egresses. The eclipse profile and duration vary substantially from cycle to cycle. We show that the X-ray spectra are consistent with a power-law-like component, absorbed by neutral and ionized Compton-thin material, and by a Compton-thick, partial-covering medium, responsible for the irregular dips. The high X-ray/optical flux ratio rules out the possibility that CG X-1 is a foreground Cataclysmic Variable; in agreement with previous studies, we conclude that it is the first example of a compact ultraluminous X-ray source fed by a Wolf-Rayet star or stripped Helium star. Its unocculted luminosity varies between $\approx$4 $\times 10^{39}$ erg s$^{-1}$ and $\approx$3 $\times 10^{40}$ erg s$^{-1}$. Both the donor star and the super-Eddington compact object drive powerful outflows: we suggest that the occulting clouds are produced in the wind-wind collision region and in the bow shock in front of the compact object. Among the rare sample of Wolf-Rayet X-ray binaries, CG X-1 is an exceptional target for studies of super-critical accretion and close binary evolution; it is also a likely progenitor of gravitational wave events.

The Roche lobe formalism describes mass transfer from one star to another. We develop an extension to hierarchical triples, considering the case in which a star donates mass to a companion which is itself a binary. The L1 point moves as the inner binary rotates, and the Roche lobe pulsates with the period of the inner binary. Signatures of mass transfer may therefore be imprinted with the orbital period of the inner binary. For some system parameters, the pulsing Roche lobe can drive mass transfer at high rates. Systems undergoing this type of mass transfer include those with inner binaries consisting of compact objects that will eventually merge, as well as progenitors of Type Ia supernovae.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans.

NASA's WFIRST mission will perform a wide-field, NIR survey of the Galactic Bulge to search for exoplanets via the microlensing techniques. As the mission is due to launch in the mid-2020s, around half-way through the LSST Main Survey, we have a unique opportunity to explore synergistic science from two landmark programs. LSST can survey the entire footprint of the WFIRST microlensing survey in a single Deep Drilling Field. Here we explore the great scientific potential of this proposal and recommend the most effective observing strategies.

Close-orbit binaries consisting of two compact objects are a center of attention because of the detection of gravitational-radiation-induced mergers. The creation of close, compact-object binaries involves physical processes that are not yet well understood; there are open questions about the manner in which two compact objects come to be close enough to merge within a Hubble time. Here we explore an important, and likely common physical process: mass transfer from a third star in a wider, hierarchical orbit. Mass added to the close binary's components can reduce the time to merger and can even change the nature of an accretor, transforming a white dwarf to a neutron star and/or a neutron star to a black hole. Some accreting white dwarfs in close binaries may even explode as Type Ia supernovae. Given the ubiquity of higher-order multiples, the evolutionary channels we lay out may be important pathways to gravitational mergers, including Type Ia supernovae. Fortunately, these pathways also lead to testable predictions.

We recommend that the search for exoplanets around binary stars be extended to include X-ray binaries in which the accretor is a white dwarf, neutron star, or black hole. We present a novel idea for detecting planets bound to such mass transfer binaries: we propose that the X-ray light curves of these binaries be inspected for signatures of transiting planets. X-ray transits may be the only way to detect planets around some systems, while providing a complementary approach to optical and/or radio observations in others. Any planets associated with X-ray binaries must be in stable orbits. We consider the range of allowable separations and find that orbital periods can be hours or longer, while transit durations extend upward from about a minute for Earth-radius planets in very close orbits, to hours for Jupiter-radius planets in wider orbits. The search for planets around mass transfer binaries could begin at once with existing X-ray observations of these systems. If and when a planet is detected around an X-ray binary, the size and mass of the planet may be readily measured, and it may also be possible to study the transmission and absorption of X-rays through its atmosphere. Finally, a noteworthy application of our proposal is that the same technique could be used to search for signals from extraterrestrial intelligence. If an advanced exocivilization placed a Dyson sphere or similar structure in orbit around the accretor of an X-ray binary in order to capture energy, such an artificial structure might cause detectable transits in the X-ray light curve.

Nearly 150 massive black hole binary (MBHB) candidates at sub-pc orbital separations have been reported in recent literature. Nevertheless, the definitive detection of even a single such object remains elusive. If at least one of the black holes is accreting, the light emitted from its accretion disc will be lensed by the other black hole for binary orbital inclinations near to the line of sight. This binary self-lensing could provide a unique signature of compact MBHB systems. We show that, for MBHBs with masses in the range 10^6 - 10^10 solar masses and with orbital periods less than ~10 yr, strong lensing events should occur in one to 10s of percent of MBHB systems that are monitored for an entire orbit. Lensing events will last from days for the less massive, shorter period MBHBs up to a year for the most massive ~10 year orbital period MBHBs. At small inclinations of the binary orbit to the line of sight, lensing must occur and will be accompanied by periodicity due to the relativistic Doppler boost. Flares at the same phase as the otherwise average flux of the Doppler modulation would be a smoking gun signature of self-lensing and can be used to constrain binary parameters. For MBHBs with separation >100 Schwarzschild radii, we show that finite-sized source effects could serve as a probe of MBH accretion disc structure. Finally, we stress that our lensing probability estimate implies that ~10 of the known MBHB candidates identified through quasar periodicity should exhibit strong lensing flares.

Gravitational lensing provides a means to measure mass that does not rely on detecting and analysing light from the lens itself. Compact objects are ideal gravitational lenses, because they have relatively large masses and are dim. In this paper we describe the prospects for predicting lensing events generated by the local population of compact objects, consisting of 250 neutron stars, 5 black holes, and approximately 35,000 white dwarfs. By focusing on a population of nearby compact objects with measured proper motions and known distances from us, we can measure their masses by studying the characteristics of any lensing event they generate. Here we concentrate on shifts in the position of a background source due to lensing by a foreground compact object. With HST, JWST, and Gaia, measurable centroid shifts caused by lensing are relatively frequent occurrences. We find that 30-50 detectable events per decade are expected for white dwarfs. Because relatively few neutron stars and black holes have measured distances and proper motions, it is more difficult to compute realistic rates for them. However, we show that at least one isolated neutron star has likely produced detectable events during the past several decades. This work is particularly relevant to the upcoming data releases by the Gaia mission and also to data that will be collected by JWST. Monitoring predicted microlensing events will not only help to determine the masses of compact objects, but will also potentially discover dim companions to these stellar remnants, including orbiting exoplanets.

We compile observations of early-type binaries identified via spectroscopy, eclipses, long-baseline interferometry, adaptive optics, common proper motion, etc. Each observational technique is sensitive to companions across a narrow parameter space of orbital periods P and mass ratios q = M_comp/M_1. After combining the samples from the various surveys and correcting for their respective selection effects, we find the properties of companions to O-type and B-type main-sequence (MS) stars differ among three regimes. First, at short orbital periods P < 20 days (separations a < 0.4 AU), the binaries have small eccentricities e < 0.4, favor modest mass ratios <q> = 0.5, and exhibit a small excess of twins q > 0.95. Second, the companion frequency peaks at intermediate periods log P (days) = 3.5 (a = 10 AU), where the binaries have mass ratios weighted toward small values q = 0.2-0.3 and follow a Maxwellian "thermal" eccentricity distribution. Finally, companions with long orbital periods log P (days) = 5.5-7.5 (a = 200-5,000 AU) are outer tertiary components in hierarchical triples, and have a mass ratio distribution across q = 0.1-1.0 that is nearly consistent with random pairings drawn from the initial mass function. We discuss these companion distributions and properties in the context of binary star formation and evolution. We also reanalyze the binary statistics of solar-type MS primaries, taking into account that (30+/-10)% of single-lined spectroscopic binaries likely contain white dwarf companions instead of low-mass stellar secondaries. The mean frequency of stellar companions with q > 0.1 and log P (days) < 8.0 per primary increases from 0.50+/-0.04 for solar-type MS primaries to 2.1+/-0.3 for O-type MS primaries. We fit joint probability density functions f(M_1,q,P,e) to the corrected distributions, which can be incorporated into binary population synthesis studies.

Globular clusters are ancient stellar populations with no star formation or core-collapse supernovae. Several lines of evidence suggest that globular clusters are rich in planets. If so, and if advanced civilizations can develop there, then the distances between these civilizations and other stars would be far smaller than typical distances between stars in the Galactic disk. The relative proximity would facilitate interstellar communication and travel. However, the very proximity that promotes interstellar travel also brings danger, since stellar interactions can destroy planetary systems. However, by modeling globular clusters and their stellar populations, we find that large regions of many globular clusters can be thought of as "sweet spots" where habitable-zone planetary orbits can be stable for long times. We also compute the ambient densities and fluxes in the regions within which habitable-zone planets can survive. Globular clusters are among the best targets for searches for extraterrestrial intelligence (SETI). We use the Drake equation to compare globular clusters to the Galactic disk, in terms of the likelihood of housing advanced communicating civilizations. We also consider free-floating planets, since wide-orbit planets can be ejected and travel freely through the cluster. A civilization spawned in a globular cluster may have opportunities to establish self-sustaining outposts, thereby reducing the probability that a single catastrophic event will destroy the civilization or its descendants. Although individual civilizations within a cluster may follow different evolutionary paths, or even be destroyed, the cluster may always host some advanced civilization, once a small number of them have managed to jump across interstellar space.

The formation of relativistic jets by an accreting compact object is one of the fundamental mysteries of astrophysics. While the theory is poorly understood, observations of relativistic jets from systems known as microquasars have led to a well-established phenomenology. Relativistic jets are not expected from sources with soft or supersoft X-ray spectra, although two such systems are known to produce relatively low-velocity bipolar outflows. Here we report optical spectra of an ultraluminous supersoft X-ray source (ULS) in the nearby galaxy M81 (M81 ULS-1) showing blueshifted broad H\alpha\ emission lines, characteristic of baryonic jets with relativistic speeds. The time variable jets have projected velocities ~17 per cent of the speed of light, and seem similar to those in the prototype microquasar SS 433. Such relativistic jets are not expected to be launched from white dwarfs, but an origin from a black hole or neutron star in M81 ULS-1 is hard to reconcile with its constant soft X-rays. The completely unexpected presence of relativistic jets in a ULS challenges the canonical theories for jet formation, but may possibly be explained by a long speculated super-critically accreting black hole with optically thick outflows

Gravitational microlensing events are powerful tools for the study of stellar populations. In particular, they can be used to discover and study a variety of binary systems. A large number of binary lenses have already been found through microlensing surveys and a few of these systems show strong evidence of orbital motion on the timescale of the lensing event. We expect that more binary lenses of this kind will be detected in the future. For binaries whose orbital period is comparable to the event duration, the orbital motion can cause the lensing signal to deviate drastically from that of a static binary lens. The most striking property of such light curves is the presence of quasi- periodic features, which are produced as the source traverses the same regions in the rotating lens plane. These repeating features contain information about the orbital period of the lens. If this period can be extracted, then much can be learned about the lensing system even without performing time-consuming, detailed light curve modeling. However, the relative transverse motion between the source and the lens significantly complicates the problem of period extraction. To resolve this difficulty, we present a modification of the standard Lomb-Scargle periodogram analysis. We test our method for four representative binary lens systems and demonstrate its efficiency in correctly extracting binary orbital periods.

We analyze 221 eclipsing binaries (EBs) in the Large Magellanic Cloud with B-type main-sequence (MS) primaries ($M_1$ $\approx$ 4 - 14 M$_{\odot}$) and orbital periods $P$ = 20 - 50 days that were photometrically monitored by the Optical Gravitational Lensing Experiment. We utilize our three-stage automated pipeline to (1) classify all 221 EBs, (2) fit physical models to the light curves of 130 detached well-defined EBs from which unique parameters can be determined, and (3) recover the intrinsic binary statistics by correcting for selection effects. We uncover two statistically significant trends with age. First, younger EBs tend to reside in dustier environments with larger photometric extinctions, an empirical relation that can be implemented when modeling stellar populations. Second, younger EBs generally have large eccentricities. This demonstrates that massive binaries at moderate orbital periods are born with a Maxwellian "thermal" orbital velocity distribution, which indicates they formed via dynamical interactions. In addition, the age-eccentricity anticorrelation provides a direct constraint for tidal evolution in highly eccentric binaries containing hot MS stars with radiative envelopes. The intrinsic fraction of B-type MS stars with stellar companions $q$ = $M_2$/$M_1$ $>$ 0.2 and orbital periods $P$ = 20 - 50 days is (7 $\pm$ 2)%. We find early-type binaries at $P$ = 20 - 50 days are weighted significantly toward small mass ratios $q$ $\approx$ 0.2 - 0.3, which is different than the results from previous observations of closer binaries with $P$ $<$ 20 days. This indicates that early-type binaries at slightly wider orbital separations have experienced substantially less coevolution and competitive accretion during their formation in the circumbinary disk.

We report on the discovery of a new X-ray pulsator, Swift J201424.9+152930 (Sw J2014). Owing to its X-ray modulation at 491 s, it was discovered in a systematic search for coherent signals in the archival data of the Swift X-ray Telescope. To investigate the nature of Sw J2014, we performed multi-wavelength follow-up observations with space-borne (Swift and XMM-Newton) and ground-based (the 1.5-m Loiano Telescope and the 3.6-m Telescopio Nazionale Galileo) instruments. The X-ray spectrum of Sw J2014 can be described by a hard and highly absorbed power law. The optical observations made it possible to single out the optical counterpart to this source, which displays several variable emission lines and total eclipses lasting ~20 min. Total eclipses of similar length were observed also in X-rays. The study of the eclipses, allowed us to infer a second periodicity of 3.44 h, which we interpret as the orbital period of a close binary system. We also found that the period has not significantly changed over a ~7 yr timespan. Based on the timing signatures of Sw J2014, and its optical and X-ray spectral properties, we suggest that it is a close binary hosting an accreting magnetic white dwarf. The system is therefore a cataclysmic variable of the intermediate polar type and one of the very few showing deep eclipses.

At intervals as short as ten thousand years, each white dwarf (WD) passes within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravitational tidal forces tear the planetoid apart; its metal-rich debris falls onto the WD, enriching the atmosphere. A third of WDs exhibit atmospheric "pollution". For roughly every hundred planetoid disruptions, a planetoid collides with a WD. We simulate a small number of collisions, in which "death-by-dynamics" refers to the fate of the planetoid. We also compute the energies and likely durations of a broad sample of collision events, and identify detection strategies at optical and X-ray wavelengths. Collisions with the most massive planetoids can be detected in external galaxies. Some may trigger nuclear burning. If one in $\sim 10^7-10^8$ of WD-planetoid collisions creates the conditions needed for a Type Ia supernova (SN~Ia), "death-by-dynamics" would also refer to the fate of the WD, and could provide a novel channel for the production of SN~Ia. We consider the circumstances under which the rate of SNe~Ia can be increased by interactions with planetoids.

Early B-type main-sequence (MS) stars (M$_1$ = 5-16 M$_{\odot}$) with closely orbiting low-mass stellar companions (q = M$_2$/M$_1$ < 0.25) can evolve to produce Type Ia supernovae, low-mass X-ray binaries, and millisecond pulsars. However, the formation mechanism and intrinsic frequency of such close extreme mass-ratio binaries have been debated, especially considering none have hitherto been detected. Utilizing observations of the Large Magellanic Cloud galaxy conducted by the Optical Gravitational Lensing Experiment, we have discovered a new class of eclipsing binaries in which a luminous B-type MS star irradiates a closely orbiting low-mass pre-MS companion that has not yet fully formed. The primordial pre-MS companions have large radii and discernibly reflect much of the light they intercept from the B-type MS primaries ($\Delta$I$_{\rm refl}$ = 0.02-0.14 mag). For the 18 definitive MS + pre-MS eclipsing binaries in our sample with good model fits to the observed light curves, we measure short orbital periods P = 3.0-8.5 days, young ages $\tau$ = 0.6-8 Myr, and small secondary masses M$_2$ = 0.8-2.4 M$_{\odot}$ (q = 0.07-0.36). The majority of these nascent eclipsing binaries are still associated with stellar nurseries, e.g. the system with the deepest eclipse $\Delta$I$_1$ = 2.8 mag and youngest age $\tau$ = 0.6$\pm$0.4 Myr is embedded in the bright H II region 30 Doradus. After correcting for selection effects, we find that (2.0$\pm$0.6)% of B-type MS stars have companions with short orbital periods P = 3.0-8.5 days and extreme mass ratios q = 0.06-0.25. This is $\approx$10 times greater than that observed for solar-type MS primaries. We discuss how these new eclipsing binaries provide invaluable insights, diagnostics, and challenges for the formation and evolution of stars, binaries, and H II regions.

Gravitational lensing events provide unique opportunities to discover and study planetary systems and binaries. Here we build on previous work to explore the role that orbital motion can play in both identifying and learning more about multiple-mass systems that serve as gravitational lenses. We find that a significant fraction of planet-lens and binary-lens light curves are influenced by orbital motion. Furthermore, the effects of orbital motion extend the range of binaries for which lens multiplicity can be discovered and studied. Orbital motion will play an increasingly important role as observations with sensitive photometry, such as those made by the space missions Kepler, Transiting Exoplanet Survey Satellite, (TESS), and WFIRST discover gravitational lensing events. Similarly, the excellent astrometric measurements made possible by GAIA will allow it to study the effects of orbital motion. Frequent observations, such as those made possible with the Korean Microlensing Telescope Network, KMTNet, will also facilitate the study of orbital motion during gravitational lensing events. Finally, orbital motion will typically play a significant role in the characteristics of lensing events in which the passage of a specific nearby star in front of a background star can be predicted.

This article offers a short guide to the steps scientists can take to ensure that their data and associated analyses continue to be of value and to be recognized. In just the past few years, hundreds of scholarly papers and reports have been written on questions of data sharing, data provenance, research reproducibility, licensing, attribution, privacy, and more, but our goal here is not to review that literature. Instead, we present a short guide intended for researchers who want to know why it is important to "care for and feed" data, with some practical advice on how to do that.

In order to understand the rates and properties of Type Ia and Type Ib/c supernovae, X-ray binaries, gravitational wave sources, and gamma ray bursts as a function of galactic environment and cosmic age, it is imperative that we measure how the close binary properties of O and B-type stars vary with metallicity. We have studied eclipsing binaries with early-B main-sequence primaries in three galaxies with different metallicities: the Large and Small Magellanic Clouds (LMC and SMC, respectively) as well as the Milky Way (MW). The observed fractions of early-B stars which exhibit deep eclipses 0.25 < Delta(m) (mag) < 0.65 and orbital periods 2 < P (days) < 20 in the MW, LMC, and SMC span a narrow range of (0.7-1.0)%, which is a model independent result. After correcting for geometrical selection effects and incompleteness toward low-mass companions, we find for early-B stars in all three environments: (1) a close binary fraction of (22+/-5)% across orbital periods 2 < P (days) < 20 and mass ratios q = M_2/M_1 > 0.1, (2) an intrinsic orbital period distribution slightly skewed toward shorter periods relative to a distribution that is uniform in log P, (3) a mass-ratio distribution weighted toward low-mass companions, and (4) a small, nearly negligible excess fraction of twins with q > 0.9. Our fitted parameters derived for the MW eclipsing binaries match the properties inferred from nearby, early-type spectroscopic binaries, which further validates our results. There are no statistically significant trends with metallicity, demonstrating that the close binary properties of massive stars do not vary across metallicities -0.7 < log(Z/Z_sun) < 0.0 beyond the measured uncertainties.

We report on observations of a luminous supersoft X-ray source (SSS) in M31, r1-25, that has exhibited spectral changes to harder X-ray states. We document these spectral changes. In addition, we show that they have important implications for modeling the source. Quasisoft states in a source that has been observed as an SSS represent a newly- discovered phenomenon. We show how such state changers could prove to be examples of unusual black hole or neutron star accretors. Future observations of this and other state changers can provide the information needed to determine the nature(s) of these intriguing new sources.

The Einstein rings and proper motions of nearby stars tend to be large. Thus, every year some foreground stars within a few hundred parsecs of Earth induce gravitational lensing events in background stars. In some of these cases, the events may exhibit evidence of planets orbiting the nearby star. In fact, planets can even be discovered during relatively distant passages. Here, we study the lensing signatures associated with planets orbiting nearby high-proper-motion stars. We find the following. (1) Wide-orbit planets can be detected for all distances of closest approach between the foreground and background stars, potentially producing independent events long before and/or after the closest approach. (2) Close-orbit planets can be detected for intermediate distances of closest approach, producing quasiperiodic signatures that may occur days or weeks before and after the stellar-lens event. (3) Planets in the so-called 'zone for resonant lensing' can significantly increase the magnification when the distance of closest approach is small, making the stellar-lens event easier to detect, while simultaneously providing evidence for planets. Because approaches close enough to allow planets to be detected can be predicted, we can plan observing strategies to take advantage of the theoretical framework built in this paper, which describes the sequence of expected effects in terms of a sequence of detection regimes.

Type Ia supernovae (SNe Ia) play important roles in our study of the expansion and acceleration of the Universe, but because we do not know the exact nature or natures of the progenitors, there is a systematic uncertainty that must be resolved if SNe Ia are to become more precise cosmic probes. No progenitor system has ever been identified either in the pre- or post-explosion images of a Ia event. There have been recent claims for and against the detection of ex-companion stars in several SNe Ia remnants. These studies, however, usually ignore the angular momentum gain of the progenitor WD, which leads to a spin-up phase and a subsequent spin-down phase before explosion. For spin-down timescales greater than 10^5 years, the donor star could be too dim to detect by the time of explosion. Here we revisit the current limits on ex-companion stars to SNR 0509-67.5, a 400 year old remnant in the Large Magellanic Cloud (LMC). If the effects of possible angular momentum gain on the WD are included, a wide range of single-degenerate progenitor models are allowed for this remnant. We demonstrate that the current absence of evidence for ex-companion stars in this remnant, as well as other SNe Ia remnants, does not necessarily provide the evidence of absence for ex-companions. We discuss potential ways to identify such ex-companion stars through deep imaging observations.

The detections of both X-ray and radio emission from the cluster G1 in M31 have provided strong support for existing dynamical evidence for an intermediate mass black hole (IMBH) of mass 1.8 +/- 0.5 x 10^4 solar masses at the cluster center. However, given the relatively low significance and astrometric accuracy of the radio detection, and the non-simultaneity of the X-ray and radio measurements, this identification required further confirmation. Here we present deep, high angular resolution, strictly simultaneous X-ray and radio observations of G1. While the X-ray emission (L_X = 1.74^{+0.53}_{-0.44} x 10^{36} (d/750 kpc)^2 erg/s in the 0.5-10 keV band) remained fully consistent with previous observations, we detected no radio emission from the cluster center down to a 3-sigma upper limit of 4.7 microJy/beam. Our favored explanation for the previous radio detection is flaring activity from a black hole low mass X-ray binary (LMXB). We performed a new regression of the Fundamental Plane of black hole activity, valid for determining black hole mass from radio and X-ray observations of sub-Eddington black holes, finding log M_{BH} = (1.638 +/- 0.070)log L_R - (1.136 +/- 0.077)log L_X - (6.863 +/- 0.790), with an empirically-determined uncertainty of 0.44 dex. This constrains the mass of the X-ray source in G1, if a black hole, to be <9700 solar masses at 95% confidence, suggesting that it is a persistent LMXB. This annuls what was previously the most convincing evidence from radiation for an IMBH in the Local Group, though the evidence for an IMBH in G1 from velocity dispersion measurements remains unaffected by these results.

We report the discovery of an X-ray/UV stellar flare from the source LMC 335, captured by XMM-Newton in the field of the Large Magellanic Cloud. The flare event was recorded continuously in X-ray for its first 10 hours from the precursor to the late decay phases. The observed fluxes increased by more than two orders of magnitude at its peak in X-ray and at least one in the UV as compared to quiescence. The peak 0.1-7.0 keV X-ray flux is derived from the two-temperature APEC model to be ~(8.4 +/- 0.6) x 10^-12 erg cm-2 s-1. Combining astrometric information from multiple X-ray observations in the quiescent and flare states, we identify the NIR counterpart of LMC 335 as the 2MASS source J05414534-6921512. The NIR color relations and spectroscopic parallax characterize the source as a Galactic K7-M4 dwarf at a foreground distance of (100 - 264) pc, implying a total energy output of the entire event of ~(0.4 - 2.9) x 10^35 erg. This report comprises detailed analyses of this late-K / early-M dwarf flare event that has the longest time coverage yet reported in the literature. The flare decay can be modeled with two exponential components with timescales of ~28 min and ~4 hours, with a single component decay firmly ruled out. The X-ray spectra during flare can be described by two components, a dominant high temperature component of ~40-60MK and a low temperature component of ~10MK, with a flare loop length of about 1.1-1.3 stellar radius.

The successful prediction of lensing events is a new and exciting enterprise that provides opportunities to discover and study planetary systems. The companion paper investigates the underlying theory. This paper is devoted to outlining the components of observing programs that can discover planets orbiting stars predicted to make a close approach to a background star. If the time and distance of closest approach can be well predicted, then the system can be targeted for individual study. In most cases, however, the predictions will be imprecise, yielding only a set of probable paths of approach and event times. We must monitor an ensemble of such systems to ensure discovery, a strategy possible with observing programs similar to a number of current surveys, including PTF and Pan-STARRS; nova searches, including those conducted by amateurs; ongoing lensing programs such as MOA and OGLE; as well as MEarth, Kepler and other transit studies. If well designed, the monitoring programs will be guaranteed to either discover planets in orbits with semi-major axes smaller than about two Einstein radii, or else to rule out their presence. Planets on wider orbits may not all be discovered, but if they are common, will be found among the events generated by ensembles of potential lenses. We consider the implications for VB 10, the first star to make a predicted approach to a background star that is close enough to allow planets to be discovered. VB 10 is not an ideal case, but it is well worth studying. A more concise summary of this work, and information for observers can be found at https://www.cfa.harvard.edu/~jmatthews/vb10.html.

Searches for planets via gravitational lensing have focused on cases in which the projected separation, a, between planet and star is comparable to the Einstein radius, R_E. This paper considers smaller orbital separations and demonstrates that evidence of close-orbit planets can be found in the low-magnification portion of the light curves generated by the central star. We develop a protocol to discover hot Jupiters as well as Neptune-mass and Earth-mass planets in the stellar habitable zone. When planets are not discovered, our method can be used to quantify the probability that the lens star does not have planets within specified ranges of the orbital separation and mass ratio. Nearby close-orbit planets discovered by lensing can be subject to follow-up observations to study the newly-discovered planets or to discover other planets orbiting the same star. Careful study of the low-magnification portions of lensing light curves should produce, in addition to the discoveries of close-orbit planets, definite detections of wide-orbit planets through the discovery of "repeating" lensing events. We show that events exhibiting extremely high magnification can effectively be probed for planets in close, intermediate, and wide distance regimes simply by adding several-time-per-night monitoring in the low-magnification wings, possibly leading to gravitational lensing discoveries of multiple planets occupying a broad range of orbits, from close to wide, in a single planetary system.

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

The explosion of a Type Ia supernova, SN 2011ef, in the nearby Pinwheel galaxy (M101 at 6.4 Mpc) provides an opportunity to study pre-explosion images and search for the progenitor, which should consist of a white dwarf (WD), possibly surrounded by an accretion disk, in orbit with another star. We report on our use of deep Chandra observations to limit the luminosity and temperature of the pre-explosion white dwarf (WD). It is found that if the spectrum was a blackbody, then WDs of highest possible temperatures and luminosities are excluded but, even if the WD was emitting at the Eddington luminosity, values of kT less than roughly 60 eV are permitted. This allows the progenitor to be an accreting nuclear-burning WD with an expanded photosphere. Pre-SN HST observations were used to derive a lower limit of about 10 eV for the expanded photosphere. Li et al.\, (2011) have already ruled out the possibility of a giant donor. We consider the combined emission from the WD, disk, and donor, and find that even the combined emission from a bright subgiant, WD and disk would not likely have been observed prior to explosion, and neither would some local candidates for the nuclear-burning WD model.

In the single degenerate scenario for Type Ia supernova (SNeIa), a white dwarf (WD) must gain a significant amount of matter from a companion star. Because the accreted mass carries angular momentum, the WD is likely to achieve fast spin periods, which can increase the critical mass, $M_{crit}$, needed for explosion. When $M_{crit}$ is higher than the maximum mass achieved by the WD, the WD must spin down before it can explode. This introduces a delay between the time at which the WD has completed its epoch of mass gain and the time of the explosion. Matter ejected from the binary during mass transfer therefore has a chance to become diffuse, and the explosion occurs in a medium with a density similar to that of typical regions of the interstellar medium. Also, either by the end of the WD's mass increase or else by the time of explosion, the donor may exhaust its stellar envelope and become a WD. This alters, generally diminishing, explosion signatures related to the donor star. Nevertheless, the spin-up/spin-down model is highly predictive. Prior to explosion, progenitors can be super-$M_{Ch}$ WDs in either wide binaries with WD companions, or else in cataclysmic variables. These systems can be discovered and studied through wide-field surveys. Post explosion, the spin-up/spin-down model predicts a population of fast-moving WDs, low-mass stars, and even brown dwarfs. In addition, the spin-up/spin-down model provides a paradigm which may be able to explain both the similarities and the diversity observed among SNeIa.

Based on new and archival Chandra observations of the Sombrero galaxy (M 104), we study the diffuse X-ray emission in and around its massive stellar bulge. The 2-6 keV unresolved emission from the bulge region closely follows the K-band star light and most likely arises from unresolved stellar sources. At lower energies, however, the unresolved emission reaches a galactocentric radius of at least 23 kpc, significantly beyond the extent of the starlight, clearly indicating the presence of diffuse hot gas. We isolate the emission of the gas by properly accounting for the emission from unresolved stellar sources, predominantly cataclysmic variables and coronally active binaries, whose quasi-universal X-ray emissivity is recently established. We find a gas temperature of ~0.6 keV with little variation across the field of view, except for a lower temperature of ~0.3 keV along the stellar disk. We measure a total intrinsic 0.3-2 keV luminosity of ~2e39 erg/s, which is comparable to the prediction by the latest galaxy formation models for disk galaxies as massive as Sombrero. However, such numerical models do not fully account for internal feedback processes, such as nuclear feedback and stellar feedback, against accretion from the intergalactic medium. On the other hand, we find no evidence for either the nucleus or the very modest star-forming activities in the disk to be a dominant heating source for the diffuse gas. We also show that neither the expected energy released by Type Ia supernovae nor the expected mass returned by evolved stars is recovered by observations. We argue that in Sombrero a galactic-scale subsonic outflow of hot gas continuously removes much of the "missing" energy and mass input from the bulge region. The observed density and temperature distributions of such an outflow, however, continues to pose challenges to theoretical studies.

We present a study of discrete X-ray sources in and around the bulge-dominated, massive Sa galaxy, Sombrero (M104), based on new and archival Chandra observations with a total exposure of ~200 ks. With a detection limit of L_X = 1E37 erg/s and a field of view covering a galactocentric radius of ~30 kpc (11.5 arcminute), 383 sources are detected. Cross-correlation with Spitler et al.'s catalogue of Sombrero globular clusters (GCs) identified from HST/ACS observations reveals 41 X-rays sources in GCs, presumably low-mass X-ray binaries (LMXBs). We quantify the differential luminosity functions (LFs) for both the detected GC and field LMXBs, whose power-low indices (~1.1 for the GC-LF and ~1.6 for field-LF) are consistent with previous studies for elliptical galaxies. With precise sky positions of the GCs without a detected X-ray source, we further quantify, through a fluctuation analysis, the GC LF at fainter luminosities down to 1E35 erg/s. The derived index rules out a faint-end slope flatter than 1.1 at a 2 sigma significance, contrary to recent findings in several elliptical galaxies and the bulge of M31. On the other hand, the 2-6 keV unresolved emission places a tight constraint on the field LF, implying a flattened index of ~1.0 below 1E37 erg/s. We also detect 101 sources in the halo of Sombrero. The presence of these sources cannot be interpreted as galactic LMXBs whose spatial distribution empirically follows the starlight. Their number is also higher than the expected number of cosmic AGNs (52+/-11 [1 sigma]) whose surface density is constrained by deep X-ray surveys. We suggest that either the cosmic X-ray background is unusually high in the direction of Sombrero, or a distinct population of X-ray sources is present in the halo of Sombrero.

We report the most accurate X-ray position of the X-ray source in the giant globular cluster G1 in M31 by using the Chandra X-ray Observatory, Hubble Space Telescope (HST), and Canada-France-Hawaii Telescope (CFHT). G1 is clearly detected with Chandra and by cross-registering with HST and CFHT images, we derive a 1-sigma error radius of 0.15 arcsec, significantly smaller than the previous measurement by XMM-Newton. We conclude that the X-ray emission of G1 is likely to come from within the core radius of the cluster. We have considered a number of possibilities for the origin of the X-ray emission but can rule all but two scenarios out: it could be due to either accretion onto a central intermediate-mass black hole (IMBH), or an ordinary low-mass X-ray binary (LMXB). Based on the X-ray luminosity and the Bondi accretion rate, an IMBH accreting from the cluster gas seems unlikely and we suggest that the X-rays are due to accretion from a companion. Alternatively, the probability that a 1.5 solar masses cluster LMXB lies within the 95 per cent X-ray error circle is about 0.7. Therefore we cannot rule out a single LMXB as the origin of the X-ray emission. While we cannot distinguish between different models with current observations, future high-resolution and high-sensitivity radio imaging observations will reveal whether there is an IMBH at the centre of G1.

In order for a white dwarf (WD) to achieve the Chandrasekhar mass, M_C, and explode as a Type Ia supernova (SNIa), it must interact with another star, either accreting matter from or merging with it. The failure to identify the types of binaries which produce SNeIa is the "progenitor problem". Its solution is required if we are to utilize the full potential of SNeIa to elucidate basic cosmological and physical principles. In single-degenerate models, a WD accretes and burns matter at high rates. Nuclear-burning WDs (NBWDs) with mass close to M_C are hot and luminous, potentially detectable as supersoft x-ray sources (SSSs). In previous work we showed that > 90-99% of the required number of progenitors do not appear as SSSs during most of the crucial phase of mass increase. The obvious implication is that double-degenerate (DD) binaries form the main class of progenitors. We show in this paper, however, that many binaries that later become DDs must pass through a long-lived NBWD phase during which they are potentially detectable as SSSs. The paucity of SSSs is therefore not a strong argument in favor of DD models. Those NBWDs that are the progenitors of DD binaries are likely to appear as symbiotic binaries for intervals > 10^6 years. In fact, symbiotic pre-DDs should be common, whether or not the WDs eventually produce SNeIa. The key to solving the progenitor problem lies in understanding the appearance of NBWDs. Most do not appear as SSSs most of the time. We therefore consider the evolution of NBWDs to address the question of what their appearance may be and how we can hope to detect them.

Kepler's first major discoveries are two hot objects orbiting stars in its field. These may be the cores of stars that have each been eroded or disrupted by a companion star. The companion, which is the star monitored today, is likely to have gained mass from its now-defunct partner, and can be considered to be a blue straggler. KOI-81 is almost certainly the product of stable mass transfer; KOI-74 may be as well, or it may be the first clear example of a blue straggler created throughthree-body interactions. We show that mass transfer binaries are common enough that Kepler should discover ~1000 white dwarfs orbiting main sequence stars. Most, like KOI-74 and KOI-81, will be discovered through transits, but many will be discovered through a combination of gravitational lensing and transits, while lensing will dominate for a subset. In fact, some events caused by white dwarfs will have the appearance of "anti-transits" --i.e., short-lived enhancements in the amount of light received from the monitored star. Lensing and other mass measurements methods provide a way to distinguish white dwarf binaries from planetary systems. This is important for the success of Kepler's primary mission, in light of the fact that white dwarf radii are similar to the radii of terrestrial planets, and that some white dwarfs will have orbital periods that place them in the habitable zones of their stellar companions. By identifying transiting and/or lensing white dwarfs, Kepler will conduct pioneering studies of white dwarfs and of the end states of mass transfer. It may also identify orbiting neutron stars or black holes. The calculations inspired by the discovery of KOI-74 and KOI-81 have implications for ground-based wide-field surveys as well as for future space-based surveys.

The population of stellar black holes (SBHs) in the Galaxy and galaxies generally is poorly known in both number and distribution. SBHs are the fossil record of the massive stars in galaxy evolution and may have produced some (if not all) of the intermediate mass (\gsim100\Msun) black holes (IMBHs) and, in turn, the central supermassive black holes (SMBHs) in galactic nuclei. For the first time, a Galaxy-wide census of accreting black holes, and their more readily recognizable tracer population, accreting neutron stars (NSs), could be measured with a wide-field hard X-ray imaging survey and soft X-ray and optical/IR prompt followup -- as proposed for the EXIST mission.

Do white dwarfs host asteroid systems? Although several lines of argument suggest that white dwarfs may be orbited by large populations of asteroids, transits would provide the most direct evidence. We demonstrate that the Kepler mission has the capability to detect transits of white dwarfs by asteroids. Because white-dwarf asteroid systems, if they exist, are likely to contain many asteroids orbiting in a spatially extended distribution, discoveries of asteroid transits can be made by monitoring only a small number of white dwarfs, compatible with Kepler's primary mission, which is to monitor stars with potentially habitable planets. Possible future missions that survey ten times as many stars with similar sensitivity and minute-cadence monitoring can establish the characteristics of asteroid systems around white dwarfs, such as the distribution of asteroid sizes and semimajor axes. Transits by planets would be more dramatic, but the probability that they will occur is lower. Ensembles of planetary moons and/or the presence of rings around planets can also produce transits detectable by Kepler. The presence of moons and rings can significantly increase the probability that Kepler will discover planets orbiting white dwarfs, even while monitoring only a small number of them.

Ongoing microlensing observations by OGLE and MOA regularly identify and conduct high-cadence sampling of lensing events with Einstein diameter crossing time, tau_E, of 16 or fewer days. Events with estimated values of tau_E of one to two days have been detected. Short duration events tend to be generated by low-mass lenses or by lenses with high transverse velocities. We compute the expected rates, demonstrate the expected ranges of parameters for lenses of different mass, and develop a protocol for observing and modeling short-duration events. Relatively minor additions to the procedures presently used will increase the rate of planet discovery, and also discover or place limits on the population of high-speed dim stars and stellar remnants in the vicinity of the Sun.

Short duration lensing events tend to be generated by low-mass lenses or by lenses with high transverse velocities. Furthermore, for any given lens mass and speed, events of short duration are preferentially caused by nearby lenses (mesolenses) that can be studied in detail, or else by lenses so close to the source star that finite-source-size effects may be detected, yielding information about both the Einstein ring radius and the surface of the lensed star. Planets causing short-duration events may be in orbits with any orientation, and may have semimajor axes smaller than an AU, or they may reach the outer limits of their planetary systems, in the region corresponding to the Solar System's Oort Cloud. They can have masses larger than Jupiter's or smaller than Pluto's. Lensing therefore has a unique potential to expand our understanding of planetary systems. A particular advantage of lensing is that it can provide precision measurements of system parameters, including the masses of and projected separation between star and planet. We demonstrate how the parameters can be extracted and show that a great deal can be learned. For example, it is remarkable that the gravitational mass of nearby free-floating planet-mass lenses can be measured by complementing observations of a photometric event with deep images that detect the planet itself. A fraction of short events may be caused by high-velocity stars located within a kpc. Many high-velocity lenses are likely to be neutron stars that received large natal kicks. Other high-speed stars may be members of the halo population. Still others may be hypervelocity stars that have been ejected from the Galactic Center, or runaway stars escaped from close binaries, possibly including the progenitor binaries of Type Ia supernovae.

In a canonical model, the progenitors of Type Ia supernovae (SNe Ia) are accreting, nuclear-burning white dwarfs (NBWDs), which explode when the white dwarf reaches the Chandrasekhar mass, M_C. Such massive NBWDs are hot (kT ~100 eV), luminous (L ~ 10^{38} erg/s), and are potentially observable as luminous supersoft X-ray sources (SSSs). During the past several years, surveys for soft X-ray sources in external galaxies have been conducted. This paper shows that the results falsify the hypothesis that a large fraction of progenitors are NBWDs which are presently observable as SSSs. The data also place limits on sub-M_C models. While Type Ia supernova progenitors may pass through one or more phases of SSS activity, these phases are far shorter than the time needed to accrete most of the matter that brings them close to M_C.

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.

Observations of hundreds of supersoft x-ray sources (SSSs) in external galaxies have shed light on the diversity of the class and on the natures of the sources. SSSs are linked to the physics of Type Ia supernovae and accretion-induced collapse, ultraluminous x-ray sources and black holes, the ionization of the interstellar medium, and tidal disruption by supermassive black holes. The class of SSSs has an extension to higher luminosities: ultraluminous SSSs have luminosities above 10^39 erg/s. There is also an extension to higher energies: quasisoft x-ray sources (QSSs) emit photons with energies above 1 eV, but few or none with energies above 2 keV. Finally, a significant fraction of the SSSs found in external galaxies switch states between observations, becoming either quasisoft or hard. For many systems ``supersoft'' refers to a temporary state; SSSs are sources, possibly including a variety of fundamentally different system types, that pass through such a state. We review those results derived from extragalactic data and related theoretical work that are most surprising and that suggest directions for future research.

Over the next decade, we can expect time domain astronomy to flourish at optical and radio wavelengths. In parallel with these efforts, a dedicated transient "machine" operating at higher energies (X-ray band through soft gamma-rays) is required to reveal the unique subset of events with variable emission predominantly visible above 100 eV. Here we focus on the transient phase space never yet sampled due to the lack of a sensitive, wide-field and triggering facility dedicated exclusively to catching high energy transients and enabling rapid coordinated multi-wavelength follow-up. We first describe the advancements in our understanding of known X-ray transients that can only be enabled through such a facility and then focus on the classes of transients theoretically predicted to be out of reach of current detection capabilities. Finally there is the exciting opportunity of revealing new classes of X-ray transients and unveiling their nature through coordinated follow-up observations at longer wavelengths.

Ultraluminous supersoft X-ray sources (ULSSS) exhibit supersoft spectra with blackbody temperatures of 50-100 eV and bolometric luminosities above $10^{39}$ erg s$^{-1}$, and are possibly intermediate mass black holes (IMBHs) of $\ge10^3 M_\odot$ or massive white dwarfs that are progenitors of type Ia supernovae. In this letter we report our optical studies of such a source in M81, M81-ULS1, with HST archive observations. M81-ULS1 is identified with a point-like object, the spectral energy distribution of which reveals a blue component in addition to the companion of an AGB star. The blue component is consistent with the power-law as expected from the geometrically-thin accretion disk around an IMBH accretor, but inconsistent with the power-law as expected from the X-ray irradiated flared accretion disk around a white dwarf accretor. This result is strong evidence that M81-ULS1 is an IMBH instead of a white dwarf.

This paper is devoted to exploring how we can discover and study nearby (< 1-2 kpc) planetary and binary systems by observing their action as gravitational lenses. Lensing can extend the realm of nearby binaries and planets that can be systematically studied to include dark and dim binaries, and face-on systems. As with more traditional studies, which use light from the system, orbital parameters (including the total mass, mass ratio, and orbital separation) can be extracted from lensing data, Also in common with these traditional methods, individual systems can be targeted for study. We discuss the specific observing strategies needed in order to optimize the discovery and study of nearby planetary and binary systems by observing their actions as lenses.

We demonstrate that gravitational lensing can be used to discover and study planets in the habitable zones of nearby dwarf stars. If appropriate software is developed, a new generation of monitoring programs will automatically conduct a census of nearby planets in the habitable zones of dwarf stars. In addition, individual nearby dwarf stars can produce lensing events at predictable times; careful monitoring of these events can discover any planets located in the zone of habitability. Because lensing can discover planets (1) in face-on orbits, and (2) in orbit around the dimmest stars, lensing techniques will provide complementary information to that gleaned through Doppler and/or transit investigations. The ultimate result will be a comprehensive understanding of the variety of systems with conditions similar to those that gave rise to life on Earth.

Nearby masses can have a high probability of lensing stars in a distant background field. High-probability lensing, or mesolensing, can therefore be used to dramatically increase our knowledge of dark and dim objects in the solar neighborhood, where it can discover and study members of the local dark population (free-floating planets, low-mass dwarfs, white dwarfs, neutron stars, and stellar mass black holes). We can measure the mass and transverse velocity of those objects discovered (or already known), and determine whether or not they are in binaries with dim companions. We explore these and other applications of mesolensing, including the study of forms of matter that have been hypothesized but not discovered, such as intermediate-mass black holes, dark matter objects free-streaming through the Galactic disk, and planets in the outermost regions of the solar system. In each case we discuss the feasibility of deriving results based on present-day monitoring systems, and also consider the vistas that will open with the advent of all-sky monitoring in the era of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), and the Large Synoptic Survey Telescope (LSST).

Binary microlensing light curves have a variety of morphologies. Many are indistinguishable from point lens light curves. Of those that deviate from the point lens form, caustic crossing light curves have tended to dominate identified binary lens events. Other distinctive signatures of binary lenses include significant asymmetry, multiple peaks, and repeating events. We have quantified, using high resolution simulations, the theoretically expected relative numbers of each type of binary lens event, based on its measurable characteristics. We find that a microlensing survey with current levels of photometric uncertainty and sampling should find at least as many non-caustic crossing binary lens events as caustic crossing events; in future surveys with more sensitive photometry, the contribution of distinctive non-caustic crossing events will be even greater. To try to explain why caustic crossing light curves appear to be so dominant among the published binary lensing events, we investigate the influence of several physical effects, including blending, sampling rate, and various binary populations.

We report the discovery of an eclipsing X-ray binary in M101, the first such system to be discovered outside the Local Group. Based on a sequence of 25 Chandra observations that sample a wide range of orbital phases, we find a period of 32.688 $\pm$ 0.002 hours, which we interpret as an orbital period. The folded light curve exhibits an eclipse lasting about 8 hours, suggesting a compact orbit in an nearly edge-on configuration. The X-ray binary has an average luminosity of $L_X(0.3-8 {\rm keV}) \approx 1.3\times10^{38}$ erg/sec, with only one out of the 25 observations significantly lower in flux than the average light curve.The presence of the eclipse and the $\sim$ 1.4-day orbital period suggests that this source is an analog of the well studied eclipsing X-ray binary pulsars Her X-1 or LMC X-4. Combining the Chandra data and the HST ACS/WFC images, we have identified several possible optical counterparts, including an O5-O3 star with V = 25.0. Follow-up optical monitoring observations should be able to identify the donor and further constrain the orbital properties.

First discovered in the Magellanic Clouds and in the Milky Way, the largest pools of luminous supersoft X-ray sources (SSSs) now known lie in M31 and in more distant galaxies. Hundreds of newly-discovered SSSs are helping us to test models for Type Ia supernovae and to identify SSSs that may represent a wider range of physical systems, including accreting intermediate-mass black holes. In this short report we list ten intriguing facts about distant SSSs.

XMM-Newton and Chandra have ushered in a new era for the study of dwarf galaxies in the Local Group. We provide an overview of the opportunities, challenges, and some early results. The large number of background sources relative to galaxy sources is a major theme. Despite this challenge, the identification of counterparts has been possible, providing hints that the same mechanisms producing X-ray sources in larger galaxies are active in dwarf galaxies. A supersoft X-ray source within 2" of the supermassive black hole in M32 may be a remnant of the tidal disruption of a giant, although other explanations cannot be ruled out.

In a variety of astronomical situations, there is a relatively high probability that a single isolated lens will produce a detectable event. The high probability is caused by some combination of a large Einstein angle, fast angular motion, and a dense background field. We refer to high-probability lenses as mesolenses. The mesolensing regime is complementary to the regime in which the optical depth is high, in that it applies to isolated lenses instead of to dense lens fields. Mesolensing is complementary to microlensing, because it is well suited to the study of different lens populations, and also because different observing and analysis strategies are required to optimize the discovery and study of mesolenses. Planetary and stellar masses located within 1-2 kpc are examples of mesolenses. We show that their presence can be detected against a wide variety of background fields, using distinctive signatures of both time variability and spatial effects. Time signatures can be identical to those of microlensing, but can also include baseline jitter, extreme apparent chromaticity, events well fit by lens models in which several independent sources are simultaneously lensed, and sequences of events. Spatial signatures include astrometric lensing of surface brightness fluctuations, as well as patterns of time variability that sweep across the background field as the lens moves. Wide-field monitoring programs, such as Pan-STARRS and LSST, are well-suited to the study of mesolensing. In addition, high resolution observations of the region behind a known mass traveling across a background field can use lensing effects to measure the lens mass and study its multiplicity.

We report the results of Chandra and XMM-Newton observations of a new outburst of an ultra-luminous X-ray source in M101. During a Chandra monitoring observation of M101, M101 ULX-1 was found to be in outburst in 2004 December, the second outburst in 2004. The peak bolometric luminosity is about 3e40 erg/s (7e39 erg/s in 0.3-7 keV). The outburst spectra are very soft and can generally be fitted with a blackbody model with temperatures of 40-80 eV, similar to supersoft X-ray sources in the Milky Way and in the Magellanic Clouds. In one Chandra observation, the source spectrum appears to be harder with a temperature of 150 eV. Such a spectral state is rarely seen in M101 ULX-1 and no X-ray source in the Milky Way shows this kind of spectrum. However, such an unusual spectral state very likely belongs to a new class of X-ray sources, quasisoft X-ray sources, recently discovered in nearby galaxies. M101 ULX-1 returned to supersoft state in a subsequent XMM-Newton observation. Based on the two outbursts in 2004, the extremely high luminosity (L_bol=1e40-1e41 erg/s), very soft X-ray spectra (kT=40-150 eV), transient behavior, and state transition provide strong evidence that M101 ULX-1 harbors an intermediate-mass black hole.

We report the results of Chandra and XMM-Newton observations of an ultra-luminous supersoft X-ray source in M101. M101 ULX-1 underwent 2 outbursts in 2004 during which the peak bolometric luminosities reached 1e41 erg/s. The outburst spectra were very soft and can generally be fitted with a blackbody model with temperatures of 50-160 eV. In two of the observations, absorption edges at 0.33 keV, 0.56 keV, 0.66 keV, and 0.88 keV were found. A cool accretion disk was also found in the 2004 December outburst. During the low luminosity state, a power-law tail was seen up to 7 keV. It is clear the source changed from a low/hard state to a high/soft state. In addition, it showed at least 5 outbursts between 1996 and 2004. This is the first ultra-luminous X-ray source for which recurrent outbursts with state transitions similar to Galactic X-ray binaries have been observed. From the Hubble Space Telescope data, we found an optical counterpart to the source. During the 2004 outbursts, we also performed radio and ground-based optical observations. All the results strongly suggest that the accreting object is a > 2800 solar mass black hole.

We report the results of Chandra and XMM-Newton observations of a new outburst of an ultraluminous supersoft X-ray source in M101. CXOU J140332.3+542103 was observed in a low luminosity state (L_x~1e37 erg/s) between 2004 January and May. The low-state X-ray spectra were relatively hard; the combined low state spectrum can be fitted with a combination of a power-law with photon index of 1.4 and a blackbody of 63 eV. During 2004 July, the source underwent a strong outburst and the peak 0.3-7 keV luminosity reached 3e40 erg/s, with a bolometric luminosity of about 1e41 erg/s. The outburst spectra were very soft and can generally be fitted with a blackbody model with temperatures of 50-100 eV. In two of the observations, absorption edges at 0.33 keV, 0.56 keV, 0.66 keV, and 0.88 keV were found. An XMM-Newton observation was also performed during the decay of the outburst, and a power-law tail was seen in addition to the supersoft spectrum. We consider different accretion models; one involving an intermediate-mass black hole can explain the observations.

Chandra observations of nearby galaxies have revealed a number of X-ray sources characterized by high luminosities ($L_X > \sim10^{36}$ erg s$^{-1}$) and $k T$ in the range $100 - 350$ eV. These ``quasisoft X-ray sources'' (QSSs) are harder than luminous supersoft X-ray sources (SSSs), whose characteristic temperatures are tens of eV, but, with little or no emission above 2 keV, they are significantly softer than most canonical X-ray sources. They are likely to include a range of physical systems; some may be common systems in unusual states: neutron stars or stellar-mass black holes (BHs), or SNRs. Others may be accreting BHs of intermediate mass. We have analyzed {\it Chandra} data from $19 $ near-by galaxies to identify QSSs with sufficiently high count rates to allow spectral fits. Six of these galaxies have been studied in great detail; in these we find 89 SSSs and 122 QSSs. In this paper we present spectra for QSSs those with more than 50 counts. We also use data from six of these galaxies to study the broadband spectral distributions of QSSs and to compare them with those of other X-ray sources. Since QSSs in the Milky Way are likely to be missed by our selection procedures, we discuss other signatures by which they can be found. If some of the QSSs for which we have spectra are accreting BHs, the lower bounds for the BH mass range from roughly $30 M_\odot$ to $1000 M_\odot$.

We report on very soft X-ray sources (VSSs) in M31. In a survey which was most sensitive to soft sources in four 8'x8' regions covered by Chandra's ACIS-S S3 CCD, we find 33 VSSs that appear to belong to M31. Fifteen VSSs have spectral characteristics mirroring the supersoft X-ray sources studied in the Magellanic Cloud and Milky Way ($k T_{eff} \leq 100$ eV); we therefore call these ``classical'' supersoft sources, or simply supersoft sources (SSSs). Eighteen VSSs may either have small (< 10%) hard components, or slightly higher effective temperatures (but still < 350 eV). We refer to these VSSs as quasisoft sources (QSSs). While hot white dwarf models may apply to SSSs, the effective temperatures of QSSs are too high, unless, e.g., the radiation emanates from only a small portion of surface. Two of the SSSs were first detected and identified as such through ROSAT observations. One SSS and one QSS may be identified with symbiotics, and 2 SSSs with supernova remnants. Both SSSs and QSSs in the disk are found near star-forming regions, possibly indicating that they are young. VSSs in the outer disk and halo are likely to be old systems; in these regions, there are more QSSs than SSSs, which is opposite to what is found in fields closer to the galaxy center. The largest density of bright VSSs is in the bulge; some of the bulge sources are close enough to the nucleus to be remnants of the tidal disruption of a giant by the massive central black hole. By using Chandra data in combination with ROSAT and XMM observations, we find most VSSs to be highly variable, fading from or brightening toward detectability on time scales of months. There is evidence for VSSs with low luminosities ($\sim 10^{36}$ erg s$^{-1}$).

We apply a uniform procedure to select very soft sources from point sources observed by Chandra in 4 galaxies. This sample includes one elliptical galaxy (NGC 4967), 2 face-on spirals (M101 and M83), and an interacting galaxy (M51). We have found very soft X-ray sources (VSSs) in every galaxy. Some of these fit the criteria for canonical supersoft sources (SSSs), while others are somewhat harder. These latter have characteristic values of kT < 300 eV; we refer to them as quasisoft sources (QSSs). We found a combined total of 149 VSSs in the 4 galaxies we considered; 77 were SSSs and 72 were QSSs. (See the paper for the original long abstract)

We have obtained 17 epochs of Chandra High Resolution Camera (HRC) snapshot images, each covering most of the M31 disk. The data cover a total baseline of 2.5 years and contain a mean effective exposure of 17 ks. We have measured the mean fluxes and long-term lightcurves for 166 objects detected in these data. At least 25% of the sources show significant variability. The cumulative luminosity function (CLF) of the disk sources is well-fit by a power-law with a slope comparable to those observed in typical elliptical galaxies. The CLF of the bulge is a broken power law similar to measurements made by previous surveys. We note several sources in the southwestern disk with L_X > 10^{37} erg/s . We cross-correlate all of our sources with published optical and radio catalogs, as well as new optical data, finding counterpart candidates for 55 sources. In addition, 17 sources are likely X-ray transients. We analyze follow-up HST WFPC2 data of two X-ray transients, finding F336W (U-band equivalent) counterparts. In both cases, the counterparts are variable. In one case, the optical counterpart is transient with F336W = 22.3 +/- 0.1 mag. The X-ray and optical properties of this object are consistent with a ~10 solar mass black hole X-ray nova with an orbital period of ~20 days. In the other case, the optical counterpart varies between F336W = 20.82 +/- 0.06 mag and F336W = 21.11 +/- 0.02 mag. Ground-based and HST observations show this object is bright (V = 18.8 +/- 0.1) and slightly extended. Finally, the frequency of bright X-ray transients in the M31 bulge suggests that the ratio of neutron star to black hole primaries in low-mass X-ray binaries (NS/BH) is ~1.

We report on possible associations between 6 globular clusters in the Virgo elliptical galaxy NGC 4472 (M49) and bright ($L_x > 10^{38}$ erg s$^{-1}$) very soft X-ray sources (VSSs). Two of the VSSs have broad-band spectral properties consistent with those of luminous supersoft X-ray sources (SSSs). The other VSSs are somewhat harder harder, possibly with values of kT between roughly 150 eV and 250 eV. These sources may be too hot to be explained by the white dwarf models so promising for SSSs; they are members of the newly-established class of quasisoft sources (QSSs). We examine white dwarf, neutron star, and black hole models for the VSSs. One of the SSSs is hot and bright enough to be a possible progenitor of a Type Ia supernova, while the most natural model for the other VSSs is one in which the accretors are intermediate-mass black holes. Whatever their physical natures, these sources are unlike any X-ray sources in Galactic globular clusters. No Galactic globular cluster houses a quasisoft source and only one, M3, contains a dim ($L_x \sim 10^{36}$ erg s$^{-1}$) transient SSS.

We introduce a procedure to identify very soft X-ray sources (VSSs) in external galaxies. Our immediate goal was to formulate a systematic procedure to identify luminous supersoft X-ray sources (SSSs), so as to allow comparisons among galaxies and to study environmental effects. The focus of this paper is on the design of the selection algorithm and on its application to simulated data. In the companion paper we test it by applying it to sources discovered through Chandra observations of 4 galaxies. We find that, in its application to both simulated and real data, our procedure also selects somewhat harder sources, which we call quasisoft. Whereas values of kT for SSSs are typically tens of eV, some quasisoft sources (QSSs) may have kT as high as ~250-300 eV. The dominant spectral component of other QSSs may be as soft as SSS spectra, but the spectra may also include a low-luminosity harder component. We sketch physical models for both supersoft and quasisoft sources. Some SSSs are likely to be accreting white dwarfs; some of these may be progenitors of Type Ia supernovae. Most QSSs may be too hot to be white dwarfs. They, together with a subset of SSSs, may be neutron stars or, perhaps most likely, accreting intermediate-mass black holes.

We present new optical and {\it Chandra} observations of the field containing the ultraluminous X-ray source NGC 1313 X-2. On an ESO 3.6 m image, the {\it Chandra} error box embraces a $R=21.6$ point-like object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of NGC 1313 X-2 is $\sim 500$. The value of $f_X/f_{opt}$, the X-ray variability history and the spectral distribution derived from a re-analysis of the {\it ROSAT}, {\it ASCA} and {\it XMM} data indicate a luminous X-ray binary in NGC 1313 as a likely explanation for NGC 1313 X-2. If the X-ray soft component observed in the {\it XMM} EPIC spectrum originates from an accretion disk, the inferred mass of the compact remnant is $\approx 100 M_\odot$, making it an intermediate mass black hole. The derived optical luminosity ($L\approx 10^5 L_\odot$) is consistent with that of a $\approx 15-20 M_\odot$ companion. The properties of the environment of NGC 1313 X-2 are briefly discussed.

We present new optical and {\it Chandra} observations of the field containing the ultraluminous X-ray source NGC1313 X-2. On an ESO 3.6 m image, the {\it Chandra} error box embraces a $R=21.6$ mag stellar-like object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of NGC 1313 X-2 is $\sim 500$. The value of $f_X/f_{opt}$, the X-ray variability history and spectral distribution indicate a luminous X-ray binary in NGC 1313 as a likely explanation for NGC 1313 X-2. The inferred optical luminosity ($L\approx 10^5 L_\odot$) is consistent with that of a $\approx 10 M_\odot$ companion.

We report on the population of point sources discovered during an 18.5 ksec Chandra ACIS-S observation of the Sombrero Galaxy. We present the luminosity function, the spectra of the 6 brightest sources, consider correlations with globular clusters (GCs) and with planetary nebulae (PNe), and study the galaxy's population of SSSs. We detected 122 sources, 22 of them are identified as luminous supersoft X-ray sources (SSSs). There is an over density of SSSs within 1.5 kpc of the nucleus, which is itself the brightest X-ray source. SSSs are also found in the disk and halo, with 1 SSS in a globular cluster (GC). Several sources in Sombrero's halo are good candidates for SSS models in which the accretor is a nuclear-burning white dwarf. In total, 32 X-ray sources are associated with GCs. The majority of sources with luminosity > 1e38 erg/s are in GCs. These results for M104, an Sa galaxy, are similar to what has been found for elliptical galaxies and for the late-type spiral M31. We find that those optically bright GCs with X-ray sources house only the brightest X-ray sources. We find that, in common with other galaxies, there appears to be a positive connection between young (metal-rich) GCs and X-ray sources, but that the brightest X-ray sources are equally likely to be in metal-poor GCs. We propose a model which can explain the trends seen in the data sets from the Sombrero and other galaxies. Thermal-time scale mass transfer can occur in some of the the younger clusters in which the turn-off mass is slightly greater than $0.8 M_\odot$; multiplicity may play a role in some of the most massive clusters; accretion from giant stars may be the dominant mechanism in some older, less massive and less centrally concentrated clusters.

We report the results of XMM-Newton observations for an especially luminous supersoft X-ray source (SSS) with bolometric luminosity of 10^39 erg/s in the spiral galaxy NGC 300. The source was detected as a SSS in 1992 and disappeared in subsequent X-ray observations. The source was active again during recent XMM-Newton observations. It appeared to be very soft (kT~60 eV) and very luminous (~10^38 - 10^39 erg/s). The two XMM-Newton observations also reveal that the source went from a ``high'' state to a ``low'' state in 6 days. We also found a 5.4-hr periodicity during the ``low'' state. We consider white dwarf, black hole, and neutron star models to explain the nature of the source.

We use a set of conservative criteria to identify luminous supersoft X-ray sources(SSSs) in external galaxies. We test this approach on blackbody models, and then apply it to {\it Chandra} data from 4 galaxies: an elliptical (NGC 4967), 2 face-on spirals (M101 and M83), and an interacting galaxy (M51). We find SSSs in every galaxy, and estimate that the total SSS population of sources with $L > 10^{37}$ erg s$^{-1}$ in each galaxy is at least several hundred and could be significantly larger. In addition, we discover two intriguing features of galactic populations of SSSs. First, there are significant subpopulations of high-luminosity sources; the bolometric luminosities can exceed $10^{39}$ erg s$^{-1}.$ Second, in the spiral galaxies M101, M83 and M51, SSSs appear to be associated with the spiral arms. This may indicate that some SSSs are young systems, possibly younger than $10^8$ years.

Three different M31 disk fields, spanning a range of stellar populations, were observed by Chandra. We report the X-ray point source luminosity function (LF) of each region, and the LF of M31's globular clusters, and compare these with each other and with the LF of the galaxy's bulge. To interpret the results we also consider tracers of the stellar population, such as OB associations and supernova remnants. We find differences in the LFs among the fields, but cannot definitively relate them to the stellar content of the fields. We find that stellar population information, average and maximum source luminosities, X-ray source densities, and slopes of the LF are useful in combination.

We report on \chandra observations of the central region of M31. By combining eight \chandra ACIS-I observations taken between 1999 and 2001, we have identified 204 X-ray sources within the central $\sim 17'\times17'$ region of M31, with a detection limit of $\sim 2\times10^{35}$\lum. Of these 204 sources, 22 are identified with globular clusters, 2 with supernova remnants, 9 with planetary nebula, and 9 as supersoft sources. By comparing individual images, about 50% of the sources are variable on time scales of months. We also found 13 transients, with light curves showing a variety of shapes. We also extracted the energy spectra of the 20 brightest sources; they can be well fit by a single power-law with a mean photon index of 1.8. The spectral shapes of 12 sources are shown to be variable, suggesting that they went through state changes. The luminosity function of all the point sources is consistent with previous observations (a broken power-law with a luminosity break at $1.7\times10^{37}$\lum). However, when the X-ray sources in different regions are considered separately, different luminosity functions are obtained. This indicates that the star-formation history might be different in different regions.

We have conducted Chandra observations of ~2560 square arcmin (~131 kpc^2) of of M31, and find that the most luminous X-ray sources in our fields are in globular clusters (GCs). Of the 28 GC X-ray sources in our fields, 17 are newly discovered. Approximately 1/3 of all the sources have L_x([0.5-7] keV) > 10^{37} ergs/s; approximately 1/10 of all the sources have L_x([0.5-7] keV) close to or above 10^{38} erg/s. The most luminous source, in the globular cluster Bo 375, is consistently observed to have L_x greater than 2 x 10^{38} ergs/s. (1) We present data on the spectra and/or light curves of the 5 most luminous M31 GC sources. (2)We explore possible explanations for the high X-ray luminosities of the brightest sources. These include that the X-ray sources may be composites, the radiation we receive may be beamed, metallicity effects could be at work, or the sources may be accreting black holes. We weigh each of these possibilities against the data. In addition, we introduce a neutron star model in which mass transfer proceeds on the thermal time scale of the donor star. Our model can produce luminosities of several times 10^{38} ergs/s, and leads to a set of well-defined predictions. (3) We compute the X-ray luminosity function and the distribution of counts in wavebands that span the range of energies to which Chandra is sensitive. We find the peak X-ray luminosity is higher and that systems with L_x > 10^{37} ergs/s constitute a larger fraction of all GC sources than in the MW. (4) We study the possible reasons for this difference between M31 and Galactic GC X-ray sources and identify three promising explanations.

We have been monitoring M31 with HRC and ACIS onboard \chandra regularly in the past two years. By combining eight \chandra ACIS-I observations taken between 1999 to 2001, we have identified 204 X-ray sources within the central $\sim 17\arcmin\times17\arcmin$ region of M31, with a detection limit of $\sim 1.6\times10^{35}$\lum. Of these 204 sources, 21 are identified with globular clusters, 2 with supernova remnants (one of them is spatially resolved with \chandra), and 8 with planetary nebula. By comparing individual images, about 50% of the sources are variable in time scales of months. We also found 14 transients. Combining all the available transients found in literatures, there are 25 transients in M31 and M32 detected by \chandra and \xmm; we present some of the long-term lightcurves by using the HRC, ACIS and \xmm data. The spectral shape of 12 sources is shown to be variable, suggesting that they went through state changes. The luminosity function of all the point sources is consistent with previous observations (a broken power-law with a luminosity break at $1.7\times10^{37}$\lum). However, when the X-ray sources in different regions are considered separately, different luminosity functions are obtained. This indicates that the star formation history might be different in different regions.

We point out that the tidal disruption of a giant may leave a luminous (10^35-10^39 ergs/s), hot (10-100 eV) stellar core. The ``supersoft'' source detected by Chandra at the center of M31 may be such a core; whether or not it is, the observations have shown that such a core is detectable, even in the center of a galaxy. We therefore explore the range of expected observational signatures and how they may be used to (1) test the hypothesis that the M31 source is a remnant of tidal stripping and (2) discover evidence of black holes and disruption events in other galaxies.

The problem of formulating a fully consistent quantum gravity theory (QGT) has not yet been solved. Even before we are able to work out the details of a complete theory, however, we do know some important qualitative features to be expected in any quantum theory. Fluctuations of the metric, for example, are expected and are associated with fluctuations of the lightcone. Lightcone fluctuations affect the arrival time of signals from distant sources in potentially measurable ways, broadening the spectra. In this paper, we start with a thermal spectrum and derive the form of spectral changes expected in a wide class of QGTs. We apply these results, valid for any thermal spectrum, to the cosmic microwave background (CMB). The CMB offers two advantages: (1) deviations from a thermal spectrum are well constrained, and (2) the radiation emanates from the most distant source of light, the surface of last scattering. We use existing CMB data to derive an upper bound on the value of $\Delta t,$ the mean spread in arrival times due to metric fluctuations: ${\Delta t} < 2.1 \times 10^{-14}$ s at the 95% confidence limit. This limit applies to a wide range of QGTs, and falsifies those predicting a larger spread in arrival times. We find this limit rules out at least one QGT, the 5-dimensional theory in which the "extra" dimension is flat. Tests of other models may also be possible, depending on the results of calculations to predict values of \Delta t,$ and also a second time scale, $\tau_c,$ the correlation time, which is the characteristic time scale of the metric fluctuations. We show that stronger limits on the value of $\Delta t,$ hence on lightcone fluctuations, can likely be derived through observations of of higher-T sources, e.g., in the X-ray and gamma-ray regimes.

This paper uses the caustic crossing events in the microlens data sets to explore the nature and location of the lenses. We conclude that the large majority of lenses, whether they are luminous or dark, are likely to be binaries. Further, we demonstrate that blending is an important feature of all the data sets. An additional interpretation suggested by the data, that the caustic crossing events along the directions to the Magellanic Clouds are due to lenses located in the Clouds, implies that most of the LMC/SMC events to date are due to lenses in the Magellanic Clouds. All of these conclusions can be tested. If they are correct, a large fraction of lenses along the direction to the LMC may be ordinary stellar binary systems, just as are the majority of the lenses along the direction to the Bulge. Thus, a better understanding of the larger-than-anticipated value derived for the Bulge optical depth may allow us to better interpret the large value derived for the optical depth to the LMC. Indeed, binarity and blending in the data sets may illuminate connections among several other puzzles: the dearth of binary-source light curves, the dearth of non-caustic-crossing perturbed binary-lens events, and the dearth of obviously blended point-lens events.

In the companion paper we began the task of systematically studying the detection of planets in wide orbits ($a > 1.5 R_E$) via microlensing surveys. In this paper we continue, focusing on repeating events. We find that, if all planetary systems are similar to our own Solar System, reasonable extensions of the present observing strategies would allow us to detect 3-6 repeating events per year along the direction to the Bulge. Indeed, if planetary systems with multiple planets are common, then future monitoring programs which lead to the discovery of thousands of stellar-lens events will likely discover events in which several different planets within a single system serve as lenses, with light curves exhibiting multiple repetitions. In this paper we discuss observing strategies to maximize the discovery of all wide-orbit planet-lens events. We also compare the likely detection rates of planets in wide orbits to those of planets located in the zone for resonant lensing. We find that, depending on the values of the planet masses and stellar radii of the lensed sources (which determine whether or not finite source size is important), and also on the sensitivity of the photometry used by observers, the detection of planets in wide orbits may be the primary route to the discovery of planets via microlensing. We also discuss how the combination of resonant and wide-orbit events can help us to learn about the distribution of planetary system properties (S 6.1). In addition, by determining the fraction of short-duration events due to planets, we indirectly derive information about the fraction of all short-duration events that may be due to low-mass MACHOs (S 6.2).

Are microlensing searches likely to discover planets that harbor life? Given our present state of knowledge, this is a difficult question to answer. We therefore begin by asking a more narrowly focused question: are conditions on planets discovered via microlensing likely to be similar to those we experience on Earth? In this paper I link the microlensing observations to the well-known "Goldilocks Problem" (conditions on the Earth-like planets need to be "just right"), to find that Earth-like planets discovered via microlensing are likely to be orbiting stars more luminous than the sun. This means that light from the planetary system's central star may contribute a significant fraction of the baseline flux relative to the star that is lensed. Such blending of light from the lens with light from the lensed source can, in principle, limit our ability to detect these events. This turns out not to be a significant problem, however. A second consequence of blending is the opportunity to determine the spectral type of the lensed spectral type of the lensed star. This circumstance, plus the possibility that finite-source-size effects are important, implies that some meaningful follow-up observations are likely to be possible for a subset Earth-like planets discovered via microlensing. In addition, calculations indicate that reasonable requirements on the planet's density and surface gravity imply that the mass of Earth-like planets is likely to be within a factor of $\sim 15$ of an Earth mass.

We propose and evaluate the feasibility of a new strategy to search for planets via microlensing. This new strategy is designed to detect planets in "wide" orbits, i.e., with orbital separation, $a$ greater than $\sim 1.5 R_E$. Planets in wide orbits may provide the dominant channel for the microlensing discovery of planets, particularly low-mass (e.g., Earth-mass) planets. This paper concentrates on events in which a single planet serves as a lens, leading to an isolated event of short duration. We point out that a distribution of events due to lensing by stars with wide-orbit planets is necessarily accompanied by a distribution of shorter- duration events. The fraction of events in the latter distribution is proportional to the average value of $\sqrt{q}$, where $q$ is the ratio between \pl and stellar masses. The position of the peak or peaks also provides a measure of the mass ratios typical of planetary systems. We study detection strategies that can optimize our ability to discover isolated short-duration events due to lensing by planets, and find that monitoring employing sensitive photometry is particularly useful. If planetary systems similar to our own are common, even modest changes in detection strategy should lead to the discovery of a few isolated events of short duration every year. We therefore also address the issue of the contamination due to stellar populations of any microlensing signal due to low-mass MACHOs. We describe how, even for isolated events of short duration, it will be possible to test the hypothesis that the lens was a planet instead of a low-mass MACHO, if the central star of the planetary system contributes a measurable fraction of the baseline flux.

We propose and evaluate the feasibility of a new strategy to search for planets via microlensing observations. This new strategy is designed to detect planets in "wide" orbits, i.e., with orbital separation, a, greater than ~1.5 R_E. Planets in wide orbits may provide the dominant channel for the discovery of planets via microlensing, particularly low-mass (e.g., Earth-mass) planets. Because the ongoing microlensing observations and extensions of them should be able to discover planets in wide orbits, we provide a foundation for the search through detailed calculations and simulations that quantify the expected results and compare the relative benefits of various search strategies. If planetary systems similar to our own or to some of the known extra-solar systems are common, then the predicted detection rates of wide-orbit events are high, generally in the range 2-10% of the present detection rate for apparently single events by stars. The expected high rates should allow the microlensing observing teams to either place significant limits on the presence of planetary systems in the Galactic Bulge, or begin to probe the population in detail within the next few years. We also address the issues of (1) whether planets discovered via microlensing are likely to harbor life, (2) the feasibility of follow-up observations to learn more about planet microlenses, and (3) the contamination due to stellar populations of any microlensing signal due to low-mass MACHOs.

We have obtained spectra of the yellow giant S1040 in the open cluster M67 using the Goddard High-Resolution Spectrograph (GHRS) and the Faint Object Spectrograph on the Hubble Space Telescope. S1040 is a single-lined spectroscopic binary with a 42.8d period that occupies a "red straggler" position in the M67 color-magnitude diagram (CMD), 0.2 mag blueward of the giant branch. A detection of S1040 at 1620 A with the Ultraviolet Imaging Telescope provided evidence that the secondary is a hot white dwarf, and thus that the anomalous location of S1040 in the CMD is likely due to a prior episode of mass-transfer. Our GHRS spectrum shows a broad Ly alpha absorption profile that confirms the white dwarf identification of the S1040 secondary. A model atmosphere fit to the GHRS spectrum yields Teff = 16,160 K, log g = 6.7, and a mass of 0.22 Msun, for an assumed cluster distance of 820 pc and reddening of E(B-V) = 0.02. The unusually low mass derived for the white dwarf implies that it must have a helium core, and that a mass-transfer episode must have begun while the progenitor was on the lower giant branch. We construct a plausible mass-transfer history for S1040 in which it originated as a short (~2d) period binary, and evolved through a blue straggler phase to reach its current state.

The microlensing monitoring programs have studied large numbers of standard light curves which seem to be due to lensing by a dark point mass. Theory predicts that many microlensing events should display significant deviations from the standard form. Lens binarity in particular is expected to be common. So far, however, only a handful of light curves exhibit evidence that the lens is a binary; all of these display dramatic deviations from the standard light curve, exhibiting pronounced multiple peaks and/or caustic crossings. Binary-lens events in which the light curve is less dramatically perturbed should also exist in the data set. Why, then, have we not detected them? The answer may lie in the fact that the perturbations, though often significant, tend to be less distinctive than those associated with caustic crossings. We present a method to determine whether a light curve is due to lensing by a binary. The method works for both gently and dramatically perturbed binary-lens light curves. Our method identifies all degenerate solutions-- i.e., all possible lensing events that might have given rise to the observed light curve. It also enables us to eliminate from consideration large ranges of possible false positive identifications associated with light curves that might mimic microlensing by a binary. This method, or a generalization of it, can also be applied to the analysis of light curves that deviate from the standard point-mass lens form because of astronomical effects other than lens binarity.

The ability to perform detailed evolutionary calculations is essential to the development of a well-defined and testable binary model. Unfortunately, traditional evolutionary calculations cannot be used to follow a significant fraction of possible close-binary supersoft sources (CBSSs). It is therefore important to examine the input physics carefully, to be sure that all relevant and potentially important physical processes are included. In this paper we continue a line of research begun last year, and explore the role that winds are expected to play in the evolution of CBSSs. We find that at least a subset of the systems that seemed to be candidates for common envelope evolution may survive, if radiation emitted by white dwarf drives winds from the system. We study the effects of winds on the binary evolution of CBSSs, and compute the number and characteristics of CBSSs expected to be presently active in galaxies such as our own or M31.

In some luminous supersoft X-ray sources (SSSs), hydrogen accretes onto the surface of a white dwarf at rates more-or-less compatible with steady nuclear burning. The white dwarfs in these systems therefore have a good chance to grow in mass. Here we review what is known about the rate of Type Ia supernovae that may be associated with SSSs. Observable consequences of the conjecture that SSSs can be progenitors of Type Ia supernovae are also discussed.

We present the light curve of an unusual variable object, OGLE #7, detected during the OGLE search for microlensing events. After one season of being in a low, normal state, the star brightened by more than 2~mag with a characteristic double-maximum shape, and returned to normal brightness after 60 days. We consider possible explanations of the photometric behavior of OGLE #7. The binary microlens model seems to be the most likely explanation -- it reproduces well the observed light curve and explains the observed colors of OGLE #7. The characteristic time scale of the OGLE #7 event, $t_E$, is equal to 80 days, the longest observed to date. The binary microlens model predicts that the spectrum of the star should be composite, with roughly 50% of its light in the $I$-band coming from a non-lensed source.