The intermediate polar FO Aquarii (FO Aqr) experienced its first-reported low-accretion states in 2016, 2017, and 2018, and using newly available photographic plates, we identify pre-discovery low states in 1965, 1966, and 1974. The primary focus of our analysis, however, is an extensive set of time-series photometry obtained between 2002 and 2018, with particularly intensive coverage of the 2016-2018 low states. After computing an updated spin ephemeris for the white dwarf (WD), we show that its spin period began to increase in 2014 after having spent 27 years decreasing; no other intermediate polar has experienced a sign change of its period derivative, but FO Aqr has now done so twice. Our central finding is that the recent low states all occurred shortly after the WD began to spin down, even though no low states were reported in the preceding quarter-century, when it was spinning up. Additionally, the system's mode of accretion is extremely sensitive to the mass-transfer rate, with accretion being almost exclusively disk-fed when FO Aqr is brighter than V~14 and substantially stream-fed when it is not. Even in the low states, a grazing eclipse remains detectable, confirming the presence of a disk-like structure (but not necessarily a Keplerian accretion disk). We relate these various observations to theoretical predictions that during the low state, the system's accretion disk dissipates into a non-Keplerian ring of diamagnetic blobs. Finally, a new XMM-Newton observation from a high state in 2017 reveals an anomalously soft X-ray spectrum and diminished X-ray luminosity compared to pre-2016 observations.

Binary stars provide a valuable test of stellar structure and evolution, because the masses of the individual stellar components can be derived with high accuracy and in a model-independent way. In this work, we study Spica, an eccentric double-lined spectroscopic binary system with a beta Cep type variable primary component. We use state-of-the-art modelling tools to determine accurate orbital elements of the binary system and atmospheric parameters of both stellar components. We interpret the short-period variability intrinsic to the primary component, detected on top of the orbital motion both in the photometric and spectroscopic data. The non-LTE based spectrum analysis reveals two stars of similar atmospheric chemical composition consistent with the present day cosmic abundance standard defined by Nieva&Przybilla (2012). The masses and radii of the stars are found to be 11.43+/-1.15 M_sun and 7.21+/-0.75 M_sun, and 7.47+/-0.54 R_sun and 3.74+/-0.53 R_sun for the primary and secondary, respectively. We find the primary component to pulsate in three independent modes, of which one is identified as a radial mode, while the two others are found to be non-radial, low degree l modes. The frequency of one of these modes is an exact multiple of the orbital frequency, and the l=m=2 mode identification suggests a tidal nature for this particular mode. We find a very good agreement between the derived dynamical and evolutionary masses for the Spica system to within the observational errors of the measured masses. The age of the system is estimated to be 12.5+/-1 Myr.

Based on gas-phase laboratory spectra at 6 K, Campbell et al. (2015) confirmed that the diffuse interstellar bands (DIBs) at 9632.7 and 9577.5A are due to absorption by the fullerene ion C60+. They also reported the detection of two other, weaker bands at 9428.5 and 9365.9A. These lie in spectral regions heavily contaminated by telluric water vapour lines. We acquired CFHT ESPaDOnS spectra of HD183143 close to the zenith and chopped with a nearby standard to correct for the telluric line absorption which enabled us to detect a DIB at 9365.9A of relative width and strength comparable to the laboratory absorption. There is a DIB of similar strength and FWHM at 9362.5A. A stellar emission feature at 9429A prevented detection of the 9428.5A band. However, a CFHT archival spectrum of HD169454, where emission is absent at 9429A, clearly shows the 9428.5A DIB with the expected strength and width. These results further confirm C60+ as a DIB carrier.

We investigate the nature of the long-period radial velocity variations in Alpha Tau first reported over 20 years ago. We analyzed precise stellar radial velocity measurements for Alpha Tau spanning over 30 years. An examination of the Halpha and Ca II 8662 spectral lines, and Hipparcos photometry was also done to help discern the nature of the long-period radial velocity variations. Our radial velocity data show that the long-period, low amplitude radial velocity variations are long-lived and coherent. Furthermore, Halpha equivalent width measurements and Hipparcos photometry show no significant variations with this period. Another investigation of this star established that there was no variability in the spectral line shapes with the radial velocity period. An orbital solution results in a period of P = 628.96 +/- 0.90 d, eccentricity, e = 0.10 +/- 0.05, and a radial velocity amplitude, K = 142.1 +/- 7.2 m/s. Evolutionary tracks yield a stellar mass of 1.13 +/- 0.11 M_sun, which corresponds to a minimum companion mass of 6.47 +/- 0.53 M_Jup with an orbital semi-major axis of a = 1.46 +/- 0.27 AU. After removing the orbital motion of the companion, an additional period of ~ 520 d is found in the radial velocity data, but only in some time spans. A similar period is found in the variations in the equivalent width of Halpha and Ca II. Variations at one-third of this period are also found in the spectral line bisector measurements. The 520 d period is interpreted as the rotation modulation by stellar surface structure. Its presence, however, may not be long-lived, and it only appears in epochs of the radial velocity data separated by $\sim$ 10 years. This might be due to an activity cycle. The data presented here provide further evidence of a planetary companion to Alpha Tau, as well as activity-related radial velocity variations.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

Context: The central star of the HR 8799 system is a gamma Doradus-type pulsator. The system harbours four planetary-mass companions detected by direct imaging, and is a good solar system analogue. The masses of the companions are not known accurately, because the estimation depends strongly on the age of the system, which is also not known with sufficient accuracy. Asteroseismic studies of the star might help to better constrain the age of HR 8799. We organized an extensive photometric and multi-site spectroscopic observing campaign for studying the pulsations of the central star. Aims: The aim of the present study is to investigate the pulsation properties of HR 8799 in detail via the ultra-precise 47-d-long nearly continuous photometry obtained with the MOST space telescope, and to find as many independent pulsation modes as possible, which is the prerequisite of an asteroseismic age determination. Methods: We carried out Fourier analysis of the wide-band photometric time series. Results: We find that resonance and sudden amplitude changes characterize the pulsation of HR 8799. The dominant frequency is always at f1 = 1.978 c/d. Many multiples of one ninth of the dominant frequency appear in the Fourier spectrum of the MOST data: n/9 f1, where n={1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 13, 14, 17, 18}. Our analysis also reveals that many of these peaks show strong amplitude decrease and phase variations even on the 47-d time-scale. The dependencies between the pulsation frequencies of HR 8799 make the planned subsequent asteroseismic analysis rather difficult. We point out some resemblance between the light curve of HR 8799 and the modulated pulsation light curves of Blazhko RR Lyrae stars.

This paper reports the results from three targeted searches of Milagro TeV sky maps: two extragalactic point source lists and one pulsar source list. The first extragalactic candidate list consists of 709 candidates selected from the Fermi-LAT 2FGL catalog. The second extragalactic candidate list contains 31 candidates selected from the TeVCat source catalog that have been detected by imaging atmospheric Cherenkov telescopes (IACTs). In both extragalactic candidate lists Mkn 421 was the only source detected by Milagro. This paper presents the Milagro TeV flux for Mkn 421 and flux limits for the brighter Fermi-LAT extragalactic sources and for all TeVCat candidates. The pulsar list extends a previously published Milagro targeted search for Galactic sources. With the 32 new gamma-ray pulsars identified in 2FGL, the number of pulsars that are studied by both Fermi-LAT and Milagro is increased to 52. In this sample, we find that the probability of Milagro detecting a TeV emission coincident with a pulsar increases with the GeV flux observed by the Fermi-LAT in the energy range from 0.1 GeV to 100 GeV.

TeV flaring activity with time scales as short as tens of minutes and an orphan TeV flare have been observed from the blazar Markarian 421 (Mrk 421). The TeV emission from Mrk 421 is believed to be produced by leptonic synchrotron self-Compton (SSC) emission. In this scenario, correlations between the X-ray and the TeV fluxes are expected, TeV orphan flares are hardly explained and the activity (measured as duty cycle) of the source at TeV energies is expected to be equal or less than that observed in X-rays if only SSC is considered. To estimate the TeV duty cycle of Mrk 421 and to establish limits on its variability at different time scales, we continuously observed Mrk 421 with the Milagro observatory. Mrk 421 was detected by Milagro with a statistical significance of 7.1 standard deviations between 2005 September 21 and 2008 March 15. The observed spectrum is consistent with previous observations by VERITAS. We estimate the duty cycle of Mrk 421 for energies above 1 TeV for different hypothesis of the baseline flux and for different flare selections and we compare our results with the X-ray duty cycle estimated by Resconi et al. 2009. The robustness of the results is discussed.

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

We use precise radial velocity measurements and photometric data to derive the frequency spacing of the p-mode oscillation spectrum of the planet-hosting star Beta Gem. This spacing along with the interferometric radius for this star is used to derive an accurate stellar mass. A long time series of over 60 hours of precise stellar radial velocity measurements of Beta Gem were taken with an iodine absorption cell and the echelle spectrograph mounted on the 2m Alfred Jensch Telescope. Complementary photometric data for this star were also taken with the MOST microsatellite spanning 3.6 d. A Fourier analysis of the radial velocity data reveals the presence of up to 17 significant pulsation modes in the frequency interval 10-250 micro-Hz. Most of these fall on a grid of equally-spaced frequencies having a separation of 7.14 +/- 0.12 micro-Hz. An analysis of 3.6 days of high precision photometry taken with the MOST space telescope shows the presence of up to 16 modes, six of which are consistent with modes found in the spectral (radial velocity) data. This frequency spacing is consistent with high overtone radial pulsations; however, until the pulsation modes are identified we cannot be sure if some of these are nonradial modes or even mixed modes. The radial velocity frequency spacing along with angular diameter measurements of Beta Gem via interferometry results in a stellar mass of M = 1.91 +/- 0.09 solar masses. This value confirms the intermediate mass of the star determined using stellar evolutionary tracks. Beta Gem is confirmed to be an intermediate mass star. Stellar pulsations in giant stars along with interferometric radius measurements can provide accurate determinations of the stellar mass of planet hosting giant stars. These can also be used to calibrate stellar evolutionary tracks.

We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of anisotropy in the motions and no radial variation in the proper motions with respect to distance from the cluster center. We also set an upper limit on any central black hole to be ~150 Msun at 90% confidence level.

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. MGRO J2019+37 is detected with a significance of 12.3 standard deviations ($\sigma$), and is better fit by a power law with an exponential cutoff than by a simple power law, with a probability $>98$% (F-test). The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of $7^{+5}_{-2}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$, a spectral index of $2.0^{+0.5}_{-1.0}$ and a cutoff energy of $29^{+50}_{-16}$ TeV. MGRO J2031+41 is detected with a significance of 7.3$\sigma$, with no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of $2.4^{+0.6}_{-0.5}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$ and a spectral index of $3.08^{+0.19}_{-0.17}$. The overall flux is subject to an $\sim$30% systematic uncertainty. The systematic uncertainty on the power law indices is $\sim$0.1. A comparison with previous results from TeV J2032+4130, MGRO J2031+41 and MGRO J2019+37 is also presented.

We present a modern ephemeris and modern light curve of the first-discovered Cepheid variable in M31, Edwin Hubble's M31-V1. Observers of the American Association of Variable Star Observers undertook these observations during the latter half of 2010. The observations were in support of an outreach program by the Space Telescope Science Institute's Hubble Heritage project, but the resulting data are the first concentrated observations of M31-V1 made in modern times. AAVSO observers obtained 214 V-band, Rc-band, and unfiltered observations from which a current ephemeris was derived. The ephemeris derived from these observations is JD(Max) = 2455430.5(+/-0.5) + 31.4 (+/-0.1) E. The period derived from the 2010 data are in agreement with the historic values of the period, but the single season of data precludes a more precise determination of the period or measurement of the period change using these data alone. However, using an ephemeris based upon the period derived by Baade and Swope we are able to fit all of the observed data acceptably well. Continued observations in the modern era will be very valuable in linking these modern data with data from the 1920s-30s and 1950s, and will enable us to measure period change in this historic Cepheid. In particular, we strongly encourage intensive observations of this star around predicted times of maximum to constrain the date of maximum to better than 0.5 days.

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV - 100 TeV gamma ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between ~1 and ~100 TeV. The TeV emission is believed to be caused by inverse-Compton up-scattering scattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fit. The fit is consistent with values measured by IACTs between 1 and 20 TeV. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4 to 2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20 to 40 TeV.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

We present strong evidence that the broad, diffuse interstellar bands (DIBs) at 4881 and 5450\,\AA are caused by the $B\,^1$B$_1$\,$\leftarrow$\,$X\,^1$A$_1$ transition of H$_2$CCC (l-C$_3$H$_2$). The large widths of the bands are due to the short lifetime of the $B\,^1$B$_1$ electronic state. The bands are predicted from absorption measurements in a neon matrix and observed by cavity ring-down in the gas phase and show exact matches to the profiles and wavelengths of the two broad DIBs. The strength of the 5450\,\AA DIB leads to a l-C$_3$H$_2$ column density of $\sim5\times10^{14}$ cm$^{-2}$ towards HD\,183143 and $\sim2\times10^{14}$\,cm$^{-2}$ to HD\,206267. Despite similar values of $E$($B-V$), the 4881 and 5450\,\AA DIBs in HD\,204827 are less than one third their strength in HD\,183143, while the column density of interstellar C$_3$ is unusually high for HD\,204827 but undetectable for HD\,183143. This can be understood if C$_3$ has been depleted by hydrogenation to species such as l-C$_3$H$_2$ towards HD\,183143. There are also three rotationally resolved sets of triplets of l-C$_3$H$_2$ in the 6150$-$6330\,\AA region. Simulations, based on the derived spectroscopic constants and convolved with the expected instrumental and interstellar line broadening, show credible coincidences with sharp, weak DIBs for the two observable sets of triplets. The region of the third set is too obscured by the $\alpha$-band of telluric O$_2$.

The aim of this work is to link the broad lambda 5450 diffuse interstellar band (DIB) to a laboratory spectrum recorded through an expanding acetylene plasma. Cavity ring-down direct absorption spectra and astronomical observations of HD 183143 with the HERMES spectrograph on the Mercator Telescope in La Palma and the McKellar spectrograph on the DAO 1.2 m Telescope are compared. In the 543-547 nm region a broad band is measured with a band maximum at 545 nm and FWHM of 1.03(0.1) nm coinciding with a well-known diffuse interstellar band at lambda 5450 with FWHM of 0.953 nm. A coincidence is found between the laboratory and the two independent observational studies obtained at higher spectral resolution. This result is important, as a match between a laboratory spectrum and a - potentially lifetime broadened - DIB is found. A series of additional experiments has been performed in order to unambiguously identify the laboratory carrier of this band. This has not been possible. The laboratory results, however, restrict the carrier to a molecular transient, consisting of carbon and hydrogen.

We present the result of a search of the Milagro sky map for spatial correlations with sources from a subset of the recent Fermi Bright Source List (BSL). The BSL consists of the 205 most significant sources detected above 100 MeV by the Fermi Large Area Telescope. We select sources based on their categorization in the BSL, taking all confirmed or possible Galactic sources in the field of view of Milagro. Of the 34 Fermi sources selected, 14 are observed by Milagro at a significance of 3 standard deviations or more. We conduct this search with a new analysis which employs newly-optimized gamma-hadron separation and utilizes the full 8-year Milagro dataset. Milagro is sensitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10-50 TeV depending on the source spectrum and declination. These results extend the observation of these sources far above the Fermi energy band. With the new analysis and additional data, multi-TeV emission is definitively observed associated with the Fermi pulsar, J2229.0+6114, in the Boomerang Pulsar Wind Nebula (PWN). Furthermore, an extended region of multi-TeV emission is associated with the Fermi pulsar, J0634.0+1745, the Geminga pulsar.

Context: Asteroseismology has great potential for the study of metal-poor stars due to its sensitivity to determine stellar ages. Aims: Our goal was to detect p-mode oscillations in the metal-poor sub-dwarf 85 Peg A and to search for other variability on longer timescales. Methods: We have obtained continuous high-precision photometry of the binary system 85 Pegasi with the MOST space telescope in two seasons (2005 & 2007). Furthermore, we redetermined vsini for 85 Peg A using high resolution spectra obtained through the ESO archive, and used photometric spot modeling to interpret long periodic variations. Results: Our frequency analysis yields no convincing evidence for p-modes significantly above a noise level of 4 ppm. Using simulated p-mode patterns we provide upper RMS amplitude limits for 85 Peg A. The light curve shows evidence for variability with a period of about 11 d and this periodicity is also seen in the follow up run in 2007; however, as different methods to remove instrumental trends in the 2005 run yield vastly different results, the exact shape and periodicity of the 2005 variability remain uncertain. Our re-determined vsini value for 85 Peg A is comparable to previous studies and we provide realistic uncertainties for this parameter. Using these values in combination with simple photometric spot models we are able to reconstruct the observed variations. Conclusions: The null-detection of p-modes in 85 Peg A is consistent with theoretical values for pulsation amplitudes in this star. The detected long-periodic variation must await confirmation by further observations with similar or better precision and long-term stability. If the 11 d periodicity is real, rotational modulation of surface features on one of the components is the most likely explanation.

In this paper, we present new spectropolarimetric observations of the planet-hosting star Tau Bootis, using ESPaDOnS and Narval spectropolarimeters at Canada-France-Hawaii Telescope (CFHT) and Telescope Bernard Lyot (TBL), respectively. We detected the magnetic field of the star at three epochs in 2008. It is a weak magnetic field of only a few Gauss, oscillating between a predominant toroidal component in January and a dominant poloidal component in June and July. A magnetic polarity reversal was observed relative to the magnetic topology in June 2007. This is the second such reversal observed in two years on this star, suggesting that Tau Boo has a magnetic cycle of about 2 years. This is the first detection of a magnetic cycle for a star other than the Sun. The role of the close-in massive planet in the short activity cycle of the star is questioned. Tau Boo has strong differential rotation, a common trend for stars with shallow convective envelope. At latitude 40 deg., the surface layer of the star rotates in 3.31 d, equal to the orbital period. Synchronization suggests that the tidal effects induced by the planet may be strong enough to force at least the thin convective envelope into corotation. Tau Boo shows variability in the Ca H & K and Halpha throughout the night and on a night to night time scale. We do not detect enhancement in the activity of the star that may be related to the conjunction of the planet. Further data is needed to conclude about the activity enhancement due to the planet.

Results are presented of a harmonic analysis of the large scale cosmic-ray anisotropy as observed by the Milagro observatory. We show a two-dimensional display of the sidereal anisotropy pro jections in right ascension generated by the fitting of three harmonics to 18 separate declination bands. The Milagro observatory is a water Cherenkov detector located in the Jemez mountains near Los Alamos, New Mexico. With a high duty cycle and large field-of-view, Milagro is an excellent instrument for measuring this anisotropy with high sensitivity at TeV energies. The analysis is conducted using a seven year data sample consisting of more than 95 billion events, the largest such data set in existence. We observe an anisotropy with a magnitude around 0.1% for cosmic rays with a median energy of 6 TeV. The dominant feature is a deficit region of depth (2.49 +/- 0.02 stat. +/- 0.09 sys.)x10^(-3) in the direction of the Galactic North Pole centered at 189 degrees right ascension. We observe a steady increase in the magnitude of the signal over seven years.

In the late 1970s and early 1980s, the introduction of solid-state, signal-generating detectors and absorption cells to impose wavelength fiducials directly on the starlight, the errors in stellar radial velocity (RV) measurements were reduced to the point where Doppler searches for planets became feasible. In 1980 we began to use a hydrogen fluoride gas cell with the CFHT coud\'{e} spectrograph and, for 12 years, monitored RVs of some 29 solar-type stars. Since extra-solar planets were expected to resemble Jupiter in both mass and orbit, we were awarded only three or four two-night observing runs each year. In 1988 we highlighted a potential planetary companion to $\gamma$ Cep (K1 IV), in 1993 one to $\beta$ Gem (K0 III), and another to $\epsilon$ Eri (K2 V) in 1992. The putative planets all resembled Jovian systems with periods and masses of 2.5 yr and 1.4 $M_{J}$, 1.6 yr and 2.6 $M_{J}$, and 6.9 yr and 0.9 $M_{J}$, respectively. All three were subsequently confirmed from more extensive data by the Texas group led by Cochran and Hatzes who derived the currently accepted orbital elements. None of the systems is simple and some still question $\epsilon$ Eri b.

We present precise photometry of the pulsating Herbig Ae star HD 142666 obtained in two consecutive years with the MOST (Microvariability & Oscilations of STars) satellite. Previously, only a single pulsation period was known for HD 142666. The MOST photometry reveals that HD 142666 is multi-periodic. However, the unique identification of pulsation frequencies is complicated by the presence of irregular variability caused by the star's circumstellar dust disk. The two light curves obtained with MOST in 2006 and 2007 provided data of unprecedented quality to study the pulsations in HD 142666 and also to monitor the circumstellar variability. We attribute 12 frequencies to pulsation. Model fits to the three frequencies with the highest amplitudes lie well outside the uncertainty box for the star's position in the HR diagram based on published values. The models suggest that either (1) the published estimate of the luminosity of HD 142666, based on a relation between circumstellar disk radius and stellar luminosity, is too high and/or (2) additional physics such as mass accretion may be needed in our models to accurately fit both the observed frequencies and HD 142666's position in the HR diagram.

A field star, HD 61199 (V ~ 8), simultaneously observed with Procyon by the MOST (Microvariability & Oscillations of STars) satellite in continuous runs of 34, 17, and 34 days in 2004, 2005, and 2007, was found to pulsate in 11 frequencies in the delta Scuti range with amplitudes from 1.7 down to 0.09 mmag. The photometry also showed variations with a period of about four days. To investigate the nature of the longer period, 45 days of time-resolved spectroscopy was obtained at the Thueringer Landessternwarte Tautenburg in 2004. The radial velocity measurements indicate that HD 61199 is a triple system. A delta Scuti pulsator with a rich eigenspectrum in a multiple system is promising for asteroseismology. Our objectives were to identify which of the stars in the system is the delta Scuti variable and to obtain the orbital elements of the system and the fundamental parameters of the individual components, which are constrained by the pulsation frequencies of the delta Scuti star. Classical Fourier techniques and least-squares multi-sinusoidal fits were applied to the MOST photometry to identify the pulsation frequencies. The groundbased spectroscopy was analysed with least-squares-deconvolution (LSD) techniques, and the orbital elements derived with the KOREL and ORBITX routines. Asteroseismic models were also generated. The photometric and spectroscopic data are compatible with a triple system consisting of a close binary with an orbital period of 3.57 days and a delta Scuti companion (HD 61199,A) as the most luminous component. The delta Scuti star is a rapid rotator with about vsin i = 130 km/s and an upper mass limit of about 2.1 Msun. For the close binary components, we find they are of nearly equal mass, with lower mass limits of about 0.7 Msun.

An analysis of 7 years of Milagro data performed on a 10-degree angular scale has found two localized regions of excess of unknown origin with greater than 12 sigma significance. Both regions are inconsistent with gamma-ray emission with high confidence. One of the regions has a different energy spectrum than the isotropic cosmic-ray flux at a level of 4.6 sigma, and it is consistent with hard spectrum protons with an exponential cutoff, with the most significant excess at ~10 TeV. Potential causes of these excesses are explored, but no compelling explanations are found.

MOST (Microvariability & Oscillations of STars) and ASAS (All Sky Automated Survey) observations have been used to characterize photometric variability of TW Hya on time scales from a fraction of a day to 7.5 weeks and from a few days to 8 years, respectively. The two data sets have very different uncertainties and temporal coverage properties and cannot be directly combined, nevertheless, they suggests a global variability spectrum with "flicker noise" properties, i.e. with amplitudes a ~ 1/sqrt(f), over >4 decades in frequency, in the range f = 0.0003 to 10 cycles per day (c/d). A 3.7 d period is clearly present in the continuous 11 day, 0.07 d time resolution, observations by MOST in 2007. Brightness extrema coincide with zero-velocity crossings in periodic (3.56 d) radial velocity variability detected in contemporaneous spectroscopic observations of Setiawan et al. (2008) and interpreted as caused by a planet. The 3.56/3.7 d periodicity was entirely absent in the second, four times longer MOST run in 2008, casting doubt on the planetary explanation. Instead, a spectrum of unstable single periods within the range of 2 - 9 days was observed; the tendency of the periods to progressively shorten was well traced using the wavelet analysis. The evolving periodicities and the overall flicker-noise characteristics of the TW Hya variability suggest a combination of several mechanisms, with the dominant ones probably related to the accretion processes from the disk around the star.

We present two series of MOST (Microvariability & Oscillations of STars) space-based photometry, covering nearly continuously 10 days in 2004 and 30 days in 2007, of selected variable stars in the upper Main Sequence of the old open cluster M67. New high-precision light curves were obtained for the blue-straggler binary/triple systems AH Cnc, ES Cnc and EV Cnc. The precision and phase coverage of ES Cnc and EV Cnc is by far superior to any previous observations. The light curve of ES Cnc is modelled in detail, assuming two dark photospheric spots and Roche geometry. An analysis of the light curve of AH Cnc indicates a low mass ratio (q about 0.13) and a high inclination angle for this system. Two new long-period eclipsing binaries, GSC 814-323 and HD75638 (non-members of M67) were discovered. We also present ground-based DDO spectroscopy of ES Cnc and of the newly found eclipsing binaries. Especially interesting is HD75638, a member of a visual binary, which must itself be a triple or a higher-multiplicity system. New light curves of two delta Scuti pulsators, EX Cnc and EW Cnc, have been analyzed leading to detection of 26 and 8 pulsation frequencies of high temporal stability.

Diffuse $\gamma$-ray emission produced by the interaction of cosmic-ray particles with matter and radiation in the Galaxy can be used to probe the distribution of cosmic rays and their sources in different regions of the Galaxy. With its large field of view and long observation time, the Milagro Gamma Ray Observatory is an ideal instrument for surveying large regions of the Northern Hemisphere sky and for detecting diffuse $\gamma$-ray emission at very high energies. Here, the spatial distribution and the flux of the diffuse $\gamma$-ray emission in the TeV energy range with a median energy of 15 TeV for Galactic longitudes between 30$^\circ$ and 110$^\circ$ and between 136$^\circ$ and 216$^\circ$ and for Galactic latitudes between -10$^\circ$ and 10$^\circ$ are determined. The measured fluxes are consistent with predictions of the GALPROP model everywhere except for the Cygnus region ($l\in[65^\circ,85^\circ]$). For the Cygnus region, the flux is twice the predicted value. This excess can be explained by the presence of active cosmic ray sources accelerating hadrons which interact with the local dense interstellar medium and produce gamma rays through pion decay.

Measuring the albedo of an extrasolar planet provides insights into its atmospheric composition and its global thermal properties, including heat dissipation and weather patterns. Such a measurement requires very precise photometry of a transiting system sampling fully many phases of the secondary eclipse. Spacebased optical photometry of the transiting system HD 209458 from the MOST (Microvariablity and Oscillations of STars) satellite, spanning 14 and 44 days in 2004 and 2005 respectively, allows us to set a sensitive limit on the optical eclipse of the hot exosolar giant planet in this system. Our best fit to the observations yields a flux ratio of the planet and star of 7 $\pm$ 9 ppm (parts per million), which corresponds to a geometric albedo through the MOST bandpass (400-700 nm) of $A_g$ = 0.038 $\pm$ 0.045. This gives a 1$\sigma$ upper limit of 0.08 for the geometric albedo and a 3$\sigma$ upper limit of 0.17. HD 209458b is significantly less reflective than Jupiter (for which $A_g$ would be about 0.5). This low geometric albedo rules out the presence of bright reflective clouds in this exoplanet's atmosphere. We determine refined parameters for the star and exoplanet in the HD 209458 system based on a model fit to the MOST light curve.

The Canadian MOST satellite is a unique platform for observations of bright transiting exoplanetary systems. Providing nearly continuous photometric observations for up to 8 weeks, MOST can produce important observational data to help us learn about the properties of exosolar planets. We review our current observations of HD 209458, HD 189733 with implications towards the albedo and our progress towards detecting reflected light from an exoplanet.

The MOST (Microvariability and Oscillations of Stars) satellite has discovered SPBe (Slowly Pulsating Be) oscillations in the stars HD 127756 (B1/B2 Vne) and HD 217543 (B3 Vpe). For HD 127756, 30 significant frequencies are identified from 31 days of nearly continuous photometry; for HD 217543, up to 40 significant frequencies from 26 days of data. In both cases, the oscillations fall into three distinct frequency ranges, consistent with models of the stars. The variations are caused by nonradial g-modes (and possibly r-modes) distorted by rapid rotation and excited by the opacity mechanism near the iron opacity bump. A comparison of pulsation models and observed frequency groups yields a rotation frequency for each star, independently of vsini. The rotation rates of these stars, as well as those of the SPBe stars previously discovered by MOST, HD 163868 and $\beta$ CMi, are all close to their critical values.

Context: We present 31.2 days of nearly continuous MOST photometry of the roAp star 10Aql. Aims:The goal was to provide an unambiguous frequency identification for this little studied star, as well as to discuss the detected frequencies in the context of magnetic models and analyze the influence of the magnetic field on the pulsation. Methods: Using traditional Fourier analysis techniques on three independent data reductions, intrinsic frequencies for the star are identified. Theoretical non-adiabatic axisymmetric modes influenced by a magnetic field having polar field strengths Bp = 0-5kG were computed to compare the observations to theory. Results: The high-precision data allow us to identify three definite intrinsic pulsation frequencies and two other candidate frequencies with low S/N. Considering the observed spacings, only one (50.95microHz) is consistent with the main sequence nature of roAp stars. The comparison with theoretical models yields a best fit for a 1.95Msun model having solar metallicity, suppressed envelope convection, and homogenous helium abundance. Furthermore, our analysis confirms the suspected slow rotation of the star and sets new lower limits to the rotation period (Prot>1 month) and inclination (i>30\pm10deg.). Conclusions:The observed frequency spectrum is not rich enough to unambiguously identify a model. On the other hand, the models hardly represent roAp stars in detail due to the approximations needed to describe the interactions of the magnetic field with stellar structure and pulsation. Consequently, errors in the model frequencies needed for the fitting procedure can only be estimated. Nevertheless, it is encouraging that models which suppress convection and include solar metallicity, in agreement with current concepts of roAp stars, fit the observations best.

Aims. Despite photometry and spectroscopy of its oscillations obtained over the past 25 years, the pulsation frequency spectrum of the rapidly oscillating Ap (roAp) star gamma Equ has remained poorly understood. Better time-series photometry, combined with recent advances to incorporate interior magnetic field geometry into pulsational models, enable us to perform improved asteroseismology of this roAp star. Methods. We obtained 19 days of continuous high-precision photometry of gamma Equ with the MOST (Microvariability & Oscillations of STars) satellite. The data were reduced with two different reduction techniques and significant frequencies were identified. Those frequencies were fitted by interpolating a grid of pulsation models that include dipole magnetic fields of various polar strengths. Results. We identify 7 frequencies in gamma Equ that we associate with 5 high-overtone p-modes and 1st and 2nd harmonics of the dominant p-mode. One of the modes and both harmonics are new discoveries for this star. Our best model solution (1.8 M_sun, log T_eff ~ 3.882; polar field strength ~8.1 kG) leads to unique mode identifications for these frequencies (ell = 0, 1, 2 and 4). This is the first purely asteroseismic fit to a grid of magnetic models. We measure amplitude and phase modulation of the primary frequency due to beating with a closely spaced frequency which had never been resolved. This casts doubts on theories that such modulation - unrelated to the rotation of the star - is due to a stochastic excitation mechanism.

We have measured transit times for HD 189733b passing in front of its bright (V = 7.67) chromospherically active and spotted parent star. Nearly continuous broadband optical photometry of this system was obtained with the MOST (Microvariability & Oscillations of STars) space telescope during 21 days in August 2006, monitoring 10 consecutive transits. We have used these data to search for deviations from a constant orbital period which can indicate the presence of additional planets in the system that are as yet undetected by Doppler searches. There are no transit timing variations above the level of ${\pm}45$ s, ruling out super-Earths (of masses $1 - 4 M_{\earth}$) in the 1:2 and 2:3 inner resonances and planets of 20 $M_{\earth}$ in the 2:1 outer resonance of the known planet. We also discuss complications in measuring transit times for a planet that transits an active star with large star spots, and how the transits can help constrain and test spot models. This has implications for the large number of such systems expected to be discovered by the CoRoT and Kepler missions.

(abridged) There is considerable interest in the possible interaction between parent stars and giant planetary companions in 51 Peg-type systems. We demonstrate from MOST satellite photometry and Ca II K line emission that there has been a persistent, variable region on the surface of tau Boo A which tracked its giant planetary companion for some 440 planetary revolutions and lies ~68deg (phi=0.8) in advance of the sub-planetary point. The light curves are folded on a range of periods centered on the planetary orbital period and phase dependent variability is quantified by Fourier methods and by the mean absolute deviation (MAD) of the folded data for both the photometry and the Ca II K line reversals. The region varies in brightness on the time scale of a rotation by ~1 mmag. In 2004 it resembled a dark spot of variable depth, while in 2005 it varied between bright and dark. Over the 123 planetary orbits spanned by the photometry the variable region detected in 2004 and in 2005 are synchronised to the planetary orbital period within 0.0015 d. The Ca II K line in 2001, 2002 and 2003 also shows enhanced K-line variability centered on phi=0.8, extending coverage to some 440 planetary revolutions. The apparently constant rotation period of the variable region and its rapid variation make an explanation in terms of conventional star spots unlikely. The lack of complementary variability at phi=0.3 and the detection of the variable region so far in advance of the sub-planetary point excludes tidal excitation, but the combined photometric and Ca II K line reversal results make a good case for an active region induced magnetically on the surface of tau Boo A by its planetary companion.

We have obtained new spectropolarimetric observations of the planet-hosting star tauBootis, using the ESPaDOnS and NARVAL spectropolarimeters at the Canada-France-Hawaii Telescope and Telescope Bernard-Lyot. With this data set, we are able to confirm the presence of a magnetic field at the surface of tauBoo and map its large-scale structure over the whole star. The overall polarity of the magnetic field has reversed with respect to our previous observation (obtained a year before), strongly suggesting that tauBoo is undergoing magnetic cycles similar to those of the Sun. This is the first time that a global magnetic polarity switch is observed in a star other than the Sun; we speculate that the magnetic cycle period of tauBoo is much shorter than that of the Sun. Our new data also allow us to confirm the presence of differential rotation from the latitudinal shearing that the magnetic structure is undergoing. The differential rotation surface shear that tauBoo experiences is found to be 6 to 10 times larger than that of the Sun. We propose that the short magnetic cycle period is due to the strong level of differential rotation. With a rotation period of 3.0 and 3.9 d at the equator and pole respectively, tauBoo appears as the first planet-hosting star whose rotation (at intermediate latitudes) is synchronised with the orbital motion of its giant planet (period 3.3 d). Assuming that this synchronisation is not coincidental, it suggests that the tidal effects induced by the giant planet can be strong enough to force the thin convective enveloppe (though not the whole star) into corotation and thus to play a role in the activity cycle of tauBoo.

We report on the measurement of transit times for the HD 209458 planetary system from photometry obtained with the MOST (Microvariability & Oscillations of STars) space telescope. Deviations from a constant orbital period can indicate the presence of additional planets in the system that are yet undetected, potentially with masses approaching an Earth mass. The MOST data sets of HD 209458 from 2004 and 2005 represent unprecedented time coverage with nearly continuous observations spanning 14 and 43 days and monitoring 3 transits and 12 consecutive transits, respectively. The transit times we obtain show no variations on three scales: (a) no long-term change in P since before 2004 at the 25 ms level, (b) no trend in transit timings during the 2005 run, and (c) no individual transit timing deviations above 80 sec level. Together with previously published transit times from Agol & Steffen (2007), this allows us to place limits on the presence of additional close-in planets in the system, in some cases down to below an Earth mass. This result, along with previous radial velocity work, now eliminates the possibility that a perturbing planet could be responsible for the additional heat source needed to explain HD 209458b's anomalous low density.

Evidence suggesting an observable magnetic interaction between a star and its hot Jupiter appears as a cyclic variation of stellar activity synchronized to the planet's orbit. In this study, we monitored the chromospheric activity of 7 stars with hot Jupiters using new high-resolution echelle spectra collected with ESPaDOnS over a few nights in 2005 and 2006 from the CFHT. We searched for variability in several stellar activity indicators (Ca II H, K, the Ca II infrared triplet, Halpha, and He I). HD 179949 has been observed almost every year since 2001. Synchronicity of the Ca II H & K emission with the orbit is clearly seen in four out of six epochs, while rotational modulation with P_rot=7 days is apparent in the other two seasons. We observe a similar phenomenon on upsilon And, which displays rotational modulation (P_rot=12 days) in September 2005, in 2002 and 2003 variations appear to correlate with the planet's orbital period. This on/off nature of star-planet interaction (SPI) in the two systems is likely a function of the changing stellar magnetic field structure throughout its activity cycle. Variability in the transiting system HD 189733 is likely associated with an active region rotating with the star, however, the flaring in excess of the rotational modulation may be associated with its hot Jupiter. As for HD 179949, the peak variability as measured by the mean absolute deviation for both HD 189733 and tau Boo leads the sub-planetary longitude by 70 degrees. The tentative correlation between this activity and the ratio of Mpsini to the planet's rotation period, a quantity proportional to the hot Jupiter's magnetic moment, first presented in Shkolnik et al. 2005 remains viable. This work furthers the characterization of SPI, improving its potential as a probe of extrasolar planetary magnetic fields.

The Microvariability and Oscillations of STars (MOST) photometric satellite has already undertaken more than 64 primary campaigns which include some clusters and has obtained observations of >850 secondary stars of which ~180 are variable. More than half of the variables pulsate, with the majority being of B-type. Since 2006 January, MOST has operated with only a single CCD for both guiding and science. The resulting increase in read-out cadence has improved precision for the brightest stars. The 2007 light curve for Procyon confirms the lack of predicted p-modes with photometric amplitudes exceeding 8 ppm as we found in 2004 and 2005. p-modes have been detected in other solar-type stars as well as pre-main sequence objects, roAp and delta Scuti variables. g-modes have been detected in a range of slowly pulsating B stars, Be stars and beta Cephei variables. Differential rotation has been defined for several spotted solar-type stars and limits set to the albedo of certain transiting planets and the presence of other perturbing planets. The mission is expected to continue as long as the experiment operates.

We report null results on a two year photometric search for outburst predictors in SS Cyg. Observations in Johnson V and Cousins I were obtained almost daily for multiple hours per night for two observing seasons. The accumulated data are put through various statistical and visual analysis techniques but fails to detect any outburst predictors. However, analysis of 102 years of AAVSO archival visual data led to the detection of a correlation between a long term quasi-periodic feature at around 1,000-2,000 days in length and an increase in outburst rate.

We have made a comprehensive transit search for exoplanets down to ~1.5 - 2 Earth radii in the HD 189733 system, based on 21-days of nearly uninterrupted broadband optical photometry obtained with the MOST (Microvariability & Oscillations of STars) satellite in 2006. We have searched these data for realistic limb-darkened transits from exoplanets other than the known hot Jupiter, HD 189733b, with periods ranging from about 0.4 days to one week. Monte Carlo statistical tests of the data with synthetic transits inserted into the data-set allow us to rule out additional close-in exoplanets with sizes ranging from about 0.15 - 0.31 RJ (Jupiter radii), or 1.7 - 3.5 RE (Earth radii) on orbits whose planes are near that of HD 189733b. These null results constrain theories that invoke lower-mass hot Super-Earth and hot Neptune planets in orbits similar to HD 189733b due to the inward migration of this hot Jupiter. This work also illustrates the feasibility of discovering smaller transiting planets around chromospherically active stars.

Highly precise and nearly uninterrupted optical photometry of the RR Lyrae star AQ Leo was obtained with the MOST (Microvariability & Oscillations of STars) satellite over 34.4 days in February-March 2005. AQ Leo was the first known double-mode RR Lyrae pulsator (RRd star). Three decades after its discovery, MOST observations have revealed that AQ Leo oscillates with at least 42 frequencies, of which 32 are linear combinations (up to the sixth order) of the radial fundamental mode and its first overtone. Evidence for period changes of these modes is found in the data. The other intrinsic frequencies may represent an additional nonradial pulsation mode and its harmonics (plus linear combinations) which warrant theoretical modeling. The unprecedented number of frequencies detected with amplitudes down to millimag precision also presents an opportunity to test nonlinear theories of mode growth and saturation in RR Lyrae pulsators.

The Hyades K2V+WD system 471 Tau is a prototype post-common envelope system and a likely cataclysmic binary progenitor. We present 10 days of nearly continuous optical photometry by the MOST (Microvariability & Oscillations of STars) satellite and partly simultaneous optical spectroscopy from DDO (David Dunlap Observatory) of the binary. The photometric data indicate that the spot coverage of the K dwarf component was less than observed in the past, suggesting that we monitored the star close to a minimum in its activity cycle. Despite the low spot activity, we still detected seven flare-like events whose estimated energies are among the highest ever observed in V471 Tau and whose times of occurrence do not correlate with the binary orbital phase. A detailed O-C analysis of the times of eclipse over the last ~35 years reveals timing variations which could be explained in several ways, including perturbations by an as-yet-undetected third body in the system or by a small orbital eccentricity inducing slow apsidal motion. The DDO spectra result in improved determinations of the K dwarf projected rotation velocity, V_K sin i = 92 km s^-1, and the orbital amplitude, K_K=150.5 km s^-1. The spectra also allow us to measure changes in H_alpha emission strength and radial velocity (RV) variations. We measure a larger H_alpha velocity amplitude than found previously suggesting that the source of the emission in V471 Tau was less concentrated around the sub-white-dwarf point on the K star than had been observed in previous studies.

Using MOST (Microvariability and Oscillations of STars) satellite guide star photometry, we have discovered a metallic A star showing hybrid p- and g-mode pulsations. HD 114839 was observed nearly continuously for 10 days in March, 2005. We identify frequencies in three groups: the first centered near 2 cycles/day, in the gamma Dor pulsation range, and two others near 8 and 20, both in the delta Scuti range. This is only the fourth known such hybrid pulsator, including another MOST discovery (Rowe et al. 2006, this issue).

The MOST photometric space mission discovered an eclipsing binary among its guide stars in June 2006 which combines a relatively large eccentricity e = 0.20 with an orbital period of only 2.27 days. HD 313926 appears to consist of two early-type stars of spectral type B3 - B7. It has a largest eccentricity among known early-type binaries with periods less than 3.5 d. Despite the large components indicated by its spectral type and light curve model, and its short period, the orbit of HD 313926 has not yet circularised so it is probably very young, even compared to other young B stars.

Recent rapid localizations of short, hard gamma-ray bursts (GRBs) by the Swift and HETE satellites have led to the observation of the first afterglows and the measurement of the first redshifts from this type of burst. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Seventeen short duration (< 5 s) GRBs detected by satellites occurred within the field of view of the Milagro gamma-ray observatory between 2000 January and 2006 December. We have searched the Milagro data for >100 GeV counterparts to these GRBs and find no significant emission correlated with these bursts. Due to the absorption of high-energy gamma rays by the extragalactic background light (EBL), detections are only expected for redshifts less than ~0.5. While most long duration GRBs occur at redshifts higher than 0.5, the opposite is thought to be true of short GRBs. Lack of a detected VHE signal thus allows setting meaningful fluence limits. One GRB in the sample (050509b) has a likely association with a galaxy at a redshift of 0.225, while another (051103) has been tentatively linked to the nearby galaxy M81. Fluence limits are corrected for EBL absorption, either using the known measured redshift, or computing the corresponding absorption for a redshift of 0.1 and 0.5, as well as for the case of z=0.

A survey of Galactic gamma-ray sources at a median energy of ~20 TeV has been performed using the Milagro Gamma Ray Observatory. Eight candidate sources of TeV emission are detected with pre-trials significance $>4.5\sigma$ in the region of Galactic longitude $l\in[30^\circ,220^\circ]$ and latitude $b\in[-10^\circ,10^\circ]$. Four of these sources, including the Crab nebula and the recently published MGRO J2019+37, are observed with significances $>4\sigma$ after accounting for the trials involved in searching the 3800 square degree region. All four of these sources are also coincident with EGRET sources. Two of the lower significance sources are coincident with EGRET sources and one of these sources is Geminga. The other two candidates are in the Cygnus region of the Galaxy. Several of the sources appear to be spatially extended. The fluxes of the sources at 20 TeV range from ~25% of the Crab flux to nearly as bright as the Crab.

We first reported evidence for differential rotation of Kappa1 Ceti in Paper I. In this paper we demonstrate that the differential rotation pattern closely matches that for the Sun. This result is based on additional MOST (Microvariability & Oscillations of STars) observations in 2004 and 2005, to complement the 2003 observations discussed in Paper I. Using StarSpotz, a program developed specifically to analyze MOST photometry, we have solved for k, the differential rotation coefficient, and P_{EQ}, the equatorial rotation period using the light curves from all three years. The spots range in latitude from 10 to 75 degrees and k = 0.090^{+0.006}_{-0.005} -- less than the solar value but consistent with the younger age of the star. k is also well constrained by the independent spectroscopic estimate of vsini. We demonstrate independently that the pattern of differential rotation with latitude in fact conforms to solar. Details are given of the parallel tempering formalism used in finding the most robust solution which gives P_{EQ} = 8.77^{+0.03}_{-0.04} days -- smaller than that usually adopted, implying an age < 750 My. Our values of P_{EQ} and k can explain the range of rotation periods determined by others by spots or activity at a variety of latitudes. Historically, Ca II activity seems to occur consistently between latitudes 50 and 60 degrees which might indicate a permanent magnetic feature. Knowledge of k and P_{EQ} are key to understanding the dynamo mechanism and rotation structure in the convective zone as well assessing age for solar-type stars. We recently published values of k and P_{EQ} for epsilon Eri based on MOST photometry and expect to analyze MOST light curves for several more spotted, solar-type stars.

We have made a comprehensive transit search for exoplanets down to about 2 Earth radii in the HD 209458 system, based on nearly uninterrupted broadband optical photometry obtained with the MOST (Microvariability and Oscillations of Stars) satellite, spanning 14 days in 2004 and 44 days in 2005. We have searched these data for limb-darkened transits at periods other than that of the known giant planet, from about 0.5 days to 2 weeks. Monte Carlo statistical tests of the data with synthetic transits inserted allow us to rule out additional close-in exoplanets with sizes ranging from about 0.20-0.36 RJ (Jupiter radii), or 2.2-4.0 RE (Earth radii) on orbits whose planes are near that of HD 209458b. These null results constrain theories that invoke lower mass planets in orbits similar to HD 209458b to explain its anomalously large radius, and those that predict "hot Earths" due to the inward migration of HD 209458b.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.

We present 10.27$ \rm{hrs}$ of photometry of the roAp star HD 134214 obtained by the MOST satellite. The star is shown to be monoperiodic and oscillating at a frequency of 2948.97$ \pm 0.55 \mu\rm{Hz}$. This is consistent with earlier ground based photometric campaigns (e.g. Kreidl et al. 1994). We do not detect any of the additional frequencies identified in the recent spectroscopic study by Kurtz et al. (2006) down to an amplitude limit of 0.36$ \rm{mmag}$ ($2\sigma$ significance limit).

We present a total of 57 days of contiguous, high-cadence photometry (14 days in 2004 and 43 in 2005) of the star BD+18 4914 obtained with the MOST satellite. We detect 16 frequencies down to a signal-to-noise of 3.6 (amplitude \~ 0.5 mmag). Six of these are less than 3 cycles/day, and the other ten are between 7 and 16 cycles/day. We intrepret the low frequencies as g-mode $\gamma$ Doradus-type pulsations and the others as $\delta$ Scuti-type p-modes, making BD+18 4914 one of the few known hybrid pulsators of its class. If the g-mode pulsations are high-overtone non-radial modes with identical low degree l, we can assign a unique mode classification of n={12, 20, 21, 22, 31, 38} based on the frequency ratio method.

Hubble Space Telescope observations of the nearby (3.22 pc), K2 V star epsilon Eridani have been combined with ground-based astrometric and radial velocity data to determine the mass of its known companion. We model the astrometric and radial velocity measurements simultaneously to obtain the parallax, proper motion, perturbation period, perturbation inclination, and perturbation size. Because of the long period of the companion, \eps b, we extend our astrometric coverage to a total of 14.94 years (including the three year span of the \HST data) by including lower-precision ground-based astrometry from the Allegheny Multichannel Astrometric Photometer. Radial velocities now span 1980.8 -- 2006.3. We obtain a perturbation period, P = 6.85 +/- 0.03 yr, semi-major axis, alpha =1.88 +/- 0.20 mas, and inclination i = 30.1 +/- 3.8 degrees. This inclination is consistent with a previously measured dust disk inclination, suggesting coplanarity. Assuming a primary mass M_* = 0.83 M_{\sun}, we obtain a companion mass M = 1.55 +/- 0.24 M_{Jup}. Given the relatively young age of epsilon Eri (~800 Myr), this accurate exoplanet mass and orbit can usefully inform future direct imaging attempts. We predict the next periastron at 2007.3 with a total separation, rho = 0.3 arcsec at position angle, p.a. = -27 degrees. Orbit orientation and geometry dictate that epsilon Eri b will appear brightest in reflected light very nearly at periastron. Radial velocities spanning over 25 years indicate an acceleration consistent with a Jupiter-mass object with a period in excess of 50 years, possibly responsible for one feature of the dust morphology, the inner cavity.

The Microvariability and Oscillations of stars (MOST) satellite has detected low-amplitude light variations ($\Delta m\sim$1 mmag) in the Be star $\beta$ CMi (B8Ve). The observations lasted 41 days and the variations have typical periods $\sim 0.3$ days. We demonstrate that the dominant frequencies are consistent with prograde high-order g-modes of $m=-1$ excited by the Fe-bump of opacity in an intermediate-mass ($\approx 3.5 M_\odot$) star with a nearly critical rotation period of 0.38 days. This is the first detection of nonradial g-mode pulsations in a Be star later than B6 leading to the possibility that pulsations are excited in all classical Be stars.

The {\it Microvariability and Oscillations of Stars (MOST)} satellite observed the B supergiant HD 163899 (B2 Ib/II) for 37 days as a guide star and detected 48 frequencies $\la$ 2.8 c d$^{-1}$ with amplitudes of a few milli-magnitudes (mmag) and less. The frequency range embraces g- and p-mode pulsations. It was generally thought that no g-modes are excited in less luminous B supergiants because strong radiative damping is expected in the core. Our theoretical models, however, show that such g-modes are excited in massive post-main-sequence stars, in accordance with these observations. The nonradial pulsations excited in models between $20M_\odot$ at $\log T_{\rm eff} \approx 4.41$ and $15M_\odot$ at $\log T_{\rm eff} \approx 4.36$ are roughly consistent with the observed frequency range. Excitation by the Fe-bump in opacity is possible because g-modes can be partially reflected at a convective zone associated with the hydrogen-burning shell, which significantly reduces radiative damping in the core. The {\it MOST} light curve of HD 163899 shows that such a reflection of g-modes actually occurs, and reveals the existence of a previously unrecognized type of variable, slowly pulsating B supergiants (SPBsg) distinct from $\alpha$ Cyg variables. Such g-modes have great potential for asteroseismology.

We present precise stellar radial velocity measurements for the K giant star Beta Gem spanning over 25 years. These data show that the long period low amplitude radial velocity variations found by Hatzes & Cochran (1993) are long-lived and coherent. An examination of the Ca II K emission, spectral line shapes from high resolution data (R = 210,000), and Hipparcos photometry show no significant variations of these quantities with the RV period. These data confirm the planetary companion hypothesis suggested by Hatzes & Cochran (1993). An orbital solution assuming a stellar mass of 1.7 M_sun yields a period, P = 589.6 days, a minimum mass of 2.3 M_Jupiter, and a semi-major axis, and a = 1.6 AU. The orbit is nearly circular (e = 0.02). Beta Gem is the seventh intermediate mass star shown to host a sub-stellar companion and suggests that planet-formation around stars much more massive than the sun may common.

Based on the long-term radial-velocity surveys carried out with the McDonald Observatory 2.7m Harlan J. Smith Telescope from 1988 to the present, we derive upper limits to long-period giant planet companions for 31 nearby stars. Data from three phases of the McDonald Observatory 2.7m planet-search program have been merged together and for 17 objects, data from the pioneering Canada-France-Hawaii Telescope (CFHT) radial-velocity program have also been included in the companion-limits determination. For those 17 objects, the baseline of observations is in excess of 23 years, enabling the detection or exclusion of giant planets in orbits beyond 8 AU. We also consider the possibility of eccentric orbits in our computations. At an orbital separation of 5.2 AU, we can exclude on average planets of M sin i > 2.0+/-1.1 Mjup (e=0) and M sin i > 4.0+/-2.8 Mjup (e=0.6) for 25 of the 31 stars in this survey. However, we are not yet able to rule out "true Jupiters," i.e. planets of M sin i ~ 1 Mjup in 5.2 AU orbits. These limits are of interest for the Space Interferometry Mission (SIM), Terrestrial Planet Finder (TPF), and Darwin missions which will search for terrestrial planets orbiting nearby stars, many of which are included in this work.

We report the discovery of a new slowly pulsating B (SPB) star, with the largest number of detected frequencies to date by more than a factor of three, based on 37 days of MOST (Microvariability & Oscillations of STars) satellite guide star photometry. The star HD 163830 (V = 9.3, B5II/III) varies in twenty detected frequencies in the range $0.035 - 1.06$ d$^{-1}$ ($0.4 - 12.3 ~ \mu$Hz) with amplitudes from 0.7 to 7.6 mmag ($S/N$ from {4 to 43}). Eighteen of these frequencies are consistent with low-degree, high-order nonradial g-modes {of seismic models of an evolved 4.5$M_{\odot}$ star. We are unable to identify one unique model due to lack of mode identifications. The lowest two frequencies may be associated with the rotation of HD 163830, but firm proof of this must await future spectroscopic data.

We present space-based photometry of the transiting exoplanetary system HD 209458 obtained with the MOST (Microvariablity and Oscillations of STars) satellite, spanning 14 days and covering 4 transits and 4 secondary eclipses. The HD 209458 photometry was obtained in MOST's lower-precision Direct Imaging mode, which is used for targets in the brightness range $6.5 < V < 13$. We describe the photometric reduction techniques for this mode of observing, in particular the corrections for stray Earthshine. We do not detect the secondary eclipse in the MOST data, to a limit in depth of 0.053 mmag (1 \sigma). We set a 1 \sigma upper limit on the planet-star flux ratio of 4.88 x 10^-5 corresponding to a geometric albedo upper limit in the MOST bandpass (400 to 700 nm) of 0.25. The corresponding numbers at the 3 \sigma level are 1.34 x 10^-4 and 0.68 respectively. HD 209458b is half as bright as Jupiter in the MOST bandpass. This low geometric albedo value is an important constraint for theoretical models of the HD209458b atmosphere, in particular ruling out the presence of reflective clouds. A second MOST campaign on HD 209458 is expected to be sensitive to an exoplanet albedo as low as 0.13 (1 sigma), if the star does not become more intrinsically variable in the meantime.

The beta Cephei star delta Ceti was considered one of the few monoperiodic variables in the class. Despite (or perhaps because of) its apparently simple oscillation spectrum, it has been challenging and controversial to identify this star's pulsation mode and constrain its physical parameters seismically. Broadband time-resolved photometry of delta Ceti spanning 18.7 days with a duty cycle of about 65% obtained by the MOST (Microvariability & Oscillations of STars) satellite -- the first scientific observations ever obtained by MOST -- reveals that the star is actually multiperiodic. Besides the well-known dominant frequency of f1 = 6.205886/d, we have discovered in the MOST data its first harmonic 2f1 and three other frequencies (f2 = 3.737/d, f3 = 3.673/d and f4 = 0.318/d), all detected with S/N > 4. In retrospect, f2 was also present in archival spectral line profile data but at lower S/N. We present seismic models whose modes match exactly the frequencies f1 and f2. Only one model falls within the common part of the error boxes of the star's observed surface gravity and effective temperature from photometry and spectroscopy. In this model, f1 is the radial (l = 0) first overtone and f2 is the g2 (l = 2, m = 0) mode. This model has a mass of 10.2+/-0.2 Msun and an age of 17.9+/-0.3 million years, making delta Ceti an evolved beta Cephei star. If f2 and f3 are rotationally split components of the same g2 mode, then the star's equatorial rotation velocity is either 27.6 km/s or half this value. Given its vsini of about 1 km/s, this implies we are seeing delta Ceti nearly pole-on.

The cataclysmic variable ASAS J002511+1217.2 was discovered in outburst by the All-Sky Automated Survey in September 2004, and intensively monitored by AAVSO observers through the following two months. Both photometry and spectroscopy indicate that this is a very short-period system. Clearly defined superhumps with a period of 0.05687 +/- 0.00001 days (1-sigma) are present during the superoutburst, 5 to 18 days following the ASAS detection. We observe a change in superhump profile similar to the transition to ``late superhumps'' observed in other short-period systems; the superhump period appears to increase slightly for a time before returning to the original value, with the resulting superhump phase offset by approximately half a period. We detect variations with a period of 0.05666 +/- 0.00003 days (1-sigma) during the four-day quiescent phase between the end of the main outburst and the single echo outburst. Weak variations having the original superhump period reappear during the echo and its rapid decline. Time-resolved spectroscopy conducted nearly 30 days after detection and well into the decline yields an orbital period measurement of 82 +/- 5 minutes. Both narrow and broad components are present in the emission line spectra, indicating the presence of multiple emission regions. The weight of the observational evidence suggests that ASAS J002511+1217.2 is a WZ Sge-type dwarf nova, and we discuss how this system fits into the WZ classification scheme.

Eight consecutive low-frequency radial p-modes are identified in the G0 IV star eta Bootis based on 27 days of ultraprecise rapid photometry obtained by the MOST (Microvariability & Oscillations of Stars) satellite. The MOST data extend smoothly to lower overtones the sequence of radial p-modes reported in earlier groundbased spectroscopy by other groups. The lower-overtone modes from the MOST data constrain the interior structure of the model of eta Boo. With the interior fit anchored by the lower-overtone modes seen by MOST, standard models are not able to fit the higher-overtone modes with the same level of accuracy. The discrepancy is similar to the discrepancy that exists between the Sun's observed p-mode frequencies and the p-mode frequencies of the standard solar model. This discrepancy promises to be a powerful constraint on models of 3D convection.

The synchronously rotating G stars in the detached, short-period (0.7 d), partially eclipsing binary, ER Vul, are the most chromospherically active solar-type stars known. We have monitored activity in the Ca II H & K reversals for almost an entire orbit. Rucinski's Broadening Function Formalism allows the photospheric contribution to be objectively subtracted from the highly blended spectra. The power of the BF technique is also demonstrated by the good agreement of radial velocities with those measured by others from less crowded spectral regions. In addition to strong Ca II emission from the primary and secondary, there appears to be a high-velocity stream flowing onto the secondary where it stimulates a large active region on the surface 30 - 40 degrees in advance of the sub-binary longitude. A model light curve with a spot centered on the same longitude also gives the best fit to the observed light curve. A flare with approximately 13% more power than at other phases was detected in one spectrum. We suggest ER Vul may offer a magnified view of the more subtle chromospheric effects synchronized to planetary revolution seen in certain `51 Peg'-type systems.

We monitored the chromospheric activity in the Ca II H & K lines of 13 solar-type stars (including the Sun); 8 of them over three years at the CFHT and 5 in a single run at the VLT. Ten of the 13 targets have close planetary companions. All of the stars observed at the CFHT show long-term (months to years) changes in H & K intensity levels. Four stars display short-term (days) cyclical activity. For two, HD 73256 and kappa^1 Ceti, the activity is likely associated with an active region rotating with the star, however, the flaring in excess of the rotational modulation may be associated with a hot jupiter. A planetary companion remains a possibility for kappa^1 Ceti. For the other two, HD 179949 and upsilon And, the cyclic variation is synchronized to the hot jupiter's orbit. For both stars this synchronicity with the orbit is clearly seen in two out of three epochs. The effect is only marginal in the third epoch at which the seasonal level of chromospheric activity had changed for both stars. Short-term chromospheric activity appears weakly dependent on the mean K-line reversal intensities for the sample of 13 stars. Also, a suggestive correlation exists between this activity and the M_p sin(i) of the star's hot jupiter. Because of their small separation (<= 0.1 AU), many of the hot jupiters lie within the Alfv\'en radius of their host stars which allows a direct magnetic interaction with the stellar surface. We discuss the conditions under which a planet's magnetic field might induce activity on the stellar surface and why no such effect was seen for the prime candidate, tau Boo. This work opens up the possibility of characterizing planet-star interactions, with implications for extrasolar planet magnetic fields and the energy contribution to stellar atmospheres.

We present time-series observations, spectra and archival outburst data of a newly-discovered variable star in Hercules, Var Her 04. Its orbital period, mass ratio, and outburst amplitude resemble those of the UGWZ-type subclass of UGSU dwarf novae. However, its supercycle and outburst light curve defy classification as a clear UGWZ. Var Her 04 is most similar to the small group of possible hydrogen-burning ``period bouncers'', dwarf novae that have passed beyond the period minimum and returned.

About 30.5 days of nearly uninterrupted broadband photometry of the solar-type star Kappa1 Ceti, obtained with the MOST (Microvariability & Oscillations of STars) satellite, shows evidence for two large starspots with different rotation periods of 8.9 and approximately 9.3 days (DeltaOmega/Omega~4%). Ground based measurements in 2002 and 2003 of Ca II H & K emission reveal variations in chromospheric activity with a period of about 9.3 days. The data were obtained during the MOST commissioning phase. When the data are combined with historical observations, they indicate that the 9.3-day spot has been stable in its period for over 30 years. The photometry, with a sampling rate of approximately once per minute, was also used to search for acoustic (p-mode) oscillations in the star. We detect no clear evidence for p-modes in the Kappa1 Ceti photometry, with a noise level around 7-9 mu_mag at frequencies in the range 0.5-4 mHz (3-sigma detection limit of 21 - 27 mu_mag). There were no flares or planetary transits during the 30.5 days of MOST monitoring with light amplitudes greater than 2 mmag (durations greater than 200 minutes) and 3 mmag (2-200 min durations). While this rules out any close-in planets of >=0.5 Jupiter diameters with an orbital inclination close to 90 deg, the scatter in differential radial velocities permit a close giant planet in a more highly inclined orbit.

We have examined the directional cross-correlation of statistical `hot-spots' between a Northern Sky TeV Gamma Ray Survey by the Milagro Observatory and a similar survey by the Tibet Array. We find the directions of these hot-spots are angularly uncorrelated between the two surveys for large angular separations (Delta theta > 4 degrees), but there appears to be a statistically significant correlation between hot-spot directions for Delta theta < 1.5 degrees. Independent simulations indicate the chance probability for the occurrence of this correlation is approximately 10^-4, implying the existence of one or more previously unobserved TeV gamma-ray sources in these directions. The data sets are consistent with both point-like sources or diffuse sources with extent of 1 - 2 degrees.

We report on the detection of a planetary companion in orbit around the primary star of the binary system $\gamma$ Cephei. High precision radial velocity measurements using 4 independent data sets spanning the time interval 1981--2002 reveal long-lived residual radial velocity variations superimposed on the binary orbit that are coherent in phase and amplitude with a period or 2.48 years (906 days) and a semi-amplitude of 27.5 m s$^{-1}$. We performed a careful analysis of our Ca II H & K S-index measurements, spectral line bisectors, and {\it Hipparcos} photometry. We found no significant variations in these quantities with the 906-d period. We also re-analyzed the Ca II $\lambda$8662 {\AA} measurements of Walker et al. (1992) which showed possible periodic variations with the ``planet'' period when first published. This analysis shows that periodic Ca II equivalent width variations were only present during 1986.5 -- 1992 and absent during 1981--1986.5. Furthermore, a refined period for the Ca II $\lambda$8662 {\AA} variations is 2.14 yrs, significantly less than residual radial velocity period. The most likely explanation of the residual radial velocity variations is a planetary mass companion with $M$ sin $i$ = 1.7 $M_{Jupiter}$ and an orbital semi-major axis of $a_2$ $=$ 2.13 AU. This supports the planet hypothesis for the residual radial velocity variations for $\gamma$ Cep first suggested by Walker et al. (1992). With an estimated binary orbital period of 57 years $\gamma$ Cep is the shortest period binary system in which an extrasolar planet has been found. This system may provide insights into the relationship between planetary and binary star formation.

An all-sky survey performed with the Tibet Air Shower Array (Tibet AS) has found a number of potential point sources of TeV gamma rays. If they are steady sources, the implied Tibet AS fluxes should be visible with strong significance to the Whipple 10-m gamma-ray telescope (E > 400 GeV) with only a short (5 hour) exposure. We have observed four candidate directions from the Tibet-II HD dataset for ~5 hours each with the Whipple telescope. In addition, we observed a new candidate direction from the Tibet-III Phase 1 dataset for 7.5 hours. We have found no corresponding excesses at the flux levels implied, and we have set upper limits for each candidate.

We have detected the synchronous enhancement of Ca II H & K emission with the short-period planetary orbit in HD 179949. High-resolution spectra taken on three observing runs extending more than a year show the enhancement coincides with phi ~ 0 (the sub-planetary point) of the 3.093-day orbit with the effect persisting for more than 100 orbits. The synchronous enhancement is consistent with planet-induced chromospheric heating by magnetic rather than tidal interaction. Something which can only be confirmed by further observations. Independent observations are needed to determine whether the stellar rotation is sychronous with the planet's orbit. Of the five 51 Peg-type systems monitored, HD 179949 shows the greatest chromospheric H & K activity. Three others show significant nightly variations but the lack of any phase coherence prevents us saying whether the activity is induced by the planet. Our two standards, tau Ceti and the Sun, show no such nightly variations.

We have measured the radial velocities of five 51 Peg-type stars and one star with constant velocity. Our measurements, on 20 \AA centered at 3947 \AA, were conventional using Th/Ar comparison spectra taken every 20 or 40 minutes between the stellar exposures. Existing IRAF routines were used for the reduction. We find $\sigma_{RV}$ $\leq$ 20 m s$^{-1}$, provided 4 measurements (out of 72) with residuals $>5\sigma_{RV}$ are neglected. The observations were made with the CFHT Gecko spectrograph, fiber-fed with the CAFE system (R$\sim$110,000). $\sigma_{RV}$ $\leq$10 m s$^{-1}$ seems possible with additional care. This study was incidental to the main program and so not exhaustive but the small value of $\sigma_{RV}$ implies that the fiber feed/image slicer system on Gecko + CAFE, essentially eliminates the long standing problem of guiding errors in radial velocity measurements. We are not suggesting this conventional approach for serious Doppler planet searches (especially with Gecko which has such a small multiplex gain), but the precision is valuable for observations made in spectral regions remote from telluric lines or captive-gas fiducials. Instrument builders might consider the advantages of the CAFE optics which incorporate agitation and invert the object and pupil for slit and grating illumination in future spectrograph designs.

We present the first results obtained at CFHT with the TRIDENT infrared camera, dedicated to the detection of faint companions close to bright nearby stars. The camera's main feature is the acquisition of three simultaneous images in three wavelengths (simultaneous differential imaging) across the methane absorption bandhead at 1.6 micron, that enables a precise subtraction of the primary star PSF while keeping the companion signal. The main limitation is non-common path aberrations between the three optical paths that slightly decorrelate the PSFs. Two types of PSF calibrations are combined with the differential simultaneous imaging technique to further attenuate the PSF: reference star subtraction and instrument rotation to smooth aberrations. It is shown that a faint companion with a DeltaH of 10 magnitudes would be detected at 0.5 arcsec from the primary.

We present the first results obtained at CFHT with the TRIDENT infrared camera, dedicated to the detection of faint companions close to bright nearby stars. Its main feature is the acquisition of three simultaneous images in three wavelengths (simultaneous differential imaging) across the methane absorption bandhead at 1.6 micron, that enables a precise subtraction of the primary star PSF while keeping the companion signal. Gl229 and 55Cnc observations are presented to demonstrate TRIDENT subtraction performances. It is shown that a faint companion with a DeltaH of 10 magnitudes would be detected at 0.5 arcsec from the primary.

We present the first results obtained at CFHT with the TRIDENT infrared camera, dedicated to the detection of faint companions close to bright nearby stars. Its main feature is the acquisition of three simultaneous images in three wavelengths (simultaneous differential imaging) across the methane absorption bandhead at 1.6 micron, that enables a precise subtraction of the primary star PSF while keeping the companion signal. Thirty-five stars have been observed in two observing missions, with no detection so far. It is shown that a faint companion with a DeltaH of 10 magnitudes would be detected at 0.5 arcsec from the primary.

We have detected modulation of the Ca II H&K reversal structure in four out of five 51 Peg-type stars whose planets have orbital periods between 3 and 4 days. We observe two effects in the K-core: (1) a broad 3-A variation at ~1% level and (2) changes on a scale of 0.5 A (~1-3%) in each of the three reversal features. The nightly variations are coherent in both H and K. From differential radial velocities measured to better than 20 m/s, up-to-date phases were extracted. The enhancements in the reversals tend to be greatest at the sub-planetary points which may imply that there is a magnetic interaction between the star's outer layers and the magnetosphere of the planet. These high-S/N (500 per pixel in the continuum) and high-resolution (R = 110,000) data are too few to confirm orbital synchronization.

The BL Lac object H1426+428 was recently detected as a high energy gamma-ray source by the VERITAS collaboration (Horan et al. 2002). We have reanalyzed the 2001 portion of the data used in the detection in order to examine the spectrum of H1426+428 above 250 GeV. We find that the time-averaged spectrum agrees with a power law of the shape dF/dE = 10^(-7.31 +- 0.15(stat) +- 0.16(syst)) x E^(-3.50 +- 0.35(stat) +- 0.05(syst)) m^(-2)s^(-1)TeV^(-1) The statistical evidence from our data for emission above 2.5 TeV is 2.6 sigma. With 95% c.l., the integral flux of H1426+428 above 2.5 TeV is larger than 3% of the corresponding flux from the Crab Nebula. The spectrum is consistent with the (non-contemporaneous) measurement by Aharonian et al. (2002) both in shape and in normalization. Below 800 GeV, the data clearly favours a spectrum steeper than that of any other TeV Blazar observed so far indicating a difference in the processes involved either at the source or in the intervening space.

We have made a sensitive search for the origin bands in the known electronic transitions of the linear carbon chains C_4 and C_5 at 3789 and 5109 A towards zeta Oph (A_V <= 1). The incentive was a recent detection of C_3 in this interstellar cloud with a column density of 1.6 x 10^12 cm^-2 plus the availability of laboratory gas phase spectra of C_4 and C_5. Further, some models of diffuse interstellar clouds predict that the abundance of these latter species should be within an order of magnitude of C_3. Despite achieving S/N of 2300 to 2600 per pixel at a resolution of ~110,000, the searches were negative, leading to 3 sigma upper limits to the column density of N(C_5) = 2 x 10^11 cm^-2 and N(C_4) = 4 x 10^12-13 cm^-2 where these values rely on theoretically calculated oscillator strengths. The implication of these limits are discussed on the choice of molecules for study in future attempts to identify the carriers of the stronger diffuse interstellar bands.

The smallest polyatomic carbon chain, C$_{3}$, has been identified in interstellar clouds (A$_{v}\sim$1 mag) towards $\zeta$ Ophiuchi, 20 Aquilae, and $\zeta$ Persei by detection of the origin band in its $A^{1}\Pi_{u}-X^{1}\Sigma^{+}_{g}$ electronic transition, near 4052\AA. Individual rotational lines were resolved up to $J$=30 enabling the rotational level column densities and temperature distributions to be determined. The inferred limits for the total column densities ($\sim$1 to 2$\times10^{12}$ cm$^{-2}$) offer a strong incentive to laboratory and astrophysical searches for the longer carbon chains. Concurrent searches for C$_2^{+}$, C$_2^{-}$ and C$_3^{-}$ were negative but provide sensitive estimates for their maximum column densities.

High precision radial velocity (RV) measurements spanning the years 1980.8--2000.0 are presented for the nearby (3.22 pc) K2 V star $\epsilon$ Eri. These data, which represent a combination of six independent data sets taken with four different telescopes, show convincing variations with a period of $\approx$ 7 yrs. A least squares orbital solution using robust estimation yields orbital parameters of period, $P$ = 6.9 yrs, velocity $K$-amplitude $=$ 19 {\ms}, eccentricity $e$ $=$ 0.6, projected companion mass $M$ sin $i$ = 0.86 $M_{Jupiter}$, and semi-major axis $a_2$ $=$ 3.3 AU. Ca II H&K S-index measurements spanning the same time interval show significant variations with periods of 3 and 20 yrs, yet none at the RV period. If magnetic activity were responsible for the RV variations then it produces a significantly different period than is seen in the Ca II data. Given the lack of Ca II variation with the same period as that found in the RV measurements, the long-lived and coherent nature of these variations, and the high eccentricity of the implied orbit, Keplerian motion due to a planetary companion seems to be the most likely explanation for the observed RV variations. The wide angular separation of the planet from the star (approximately 1 arc-second) and the long orbital period make this planet a prime candidate for both direct imaging and space-based astrometric measurements.

We present high spatial resolution, near-IR images in J, H, K of the nucleus of NGC3227, obtained with the Adaptive Optics bonnette on CFHT. The ~0.15 arcsec (17pc) resolution allows structures to be probed in the core region. Dust obscuration becomes significantly less pronounced at longer wavelengths, revealing the true geometry of the core region. We are able to identify two main features in our maps: (i) a spiraling association of knots with a counterpart in an HST F606W image; (ii) a smaller scale annulus, orthogonal to the spiral of knots. These features may provide a means to transport material inwards to fuel the active nucleus.

We have searched, so far in vain, for brown dwarfs or giant planets in the vicinity of several nearby stars using cameras built at Universite de Montreal attached to the CFHT Adaptive Optics system. Here, we show how pairs of images taken simultaneously through separate filters on, and off, the 1.6 micron methane absorption band, can be used to completely defeat speckles, the dominant noise source, while revealing cool brown dwarfs and giant planets in high contrast. We achieve the photon shot-noise limit for data taken without a focal plane mask but, so far, have had less success with data where the image is occulted by a mask.

As part of a project to map nearby Seyfert galaxies in the near-IR with adaptive optics, we present high spatial resolution, near-IR images in J,H, and K of the nuclei of NGC 3227 and NGC 2992, obtained with the Adaptive Optics Bonnette (AOB) on CFHT. The ~0.15" resolution allows us to probe structures in the core region at unprecedented scales. With NGC 3227, we are able to identify an inwards spiraling starburst in all three near-IR bands. NGC 2992 shows evidence for emission along a radio loop. Compared with HST optical images, dust obscuration becomes significantly less pronounced at longer wavelengths, revealing the true geometry of the core regions. The observed structures may help to elucidate fueling mechanisms for the central engine, as well as providing insight into the unification paradigms. The results are tempered with the discovery of AOB-related artifacts in the central arcsec of observed AGN galaxies which take on the appearance of spiral/disk structures.

First imaging results are obtained with a new CCD mosaic prototype (3K x 3k, 15 \mu m pixels). The CCDs are aligned using an etched socket alignment technique. Three different measurements of the alignment are made using star images, test p attern images, and microscope analysis. The CCDs have an angular misalignement o f less than 30 ppm. The composite device is flat to within \pm 3 \mu m, with rows/columns oriented to within 20 ppm. The use of an existing technology with built in precision reduces many of the difficulties and expenses typically encountered with mosaic detector construction. A new camera being built for the UBC liquid mirror telescope is al so described.

We report the successful development and operation of a large astronomical liquid-mirror telescope. Employing a rotating 2.7-meter diameter mirror with a surface of liquid mercury, the telescope images a half-degree diameter field centered at the zenith. Located near Vancouver, British Columbia, it is equipped with a low-noise 2048x2048-pixel CCD detector, operating in TDI mode, which produces continuous imaging of a 20 arcminute-widestrip of sky with 2 minute integration time. Images with FWHM of 2 arcseconds or less are regularly obtained. This image quality is limited only by atmospheric seeing and star-trail curvature. The telescope is equipped with a series of narrow-band filters, designed to produce 40-point spectral energy distributions from 0.4 to 1.0 microns of all detected objects. These will allow classification and redshift estimation of approximately 25,000 galaxies and 1,000 quasars to a limiting magnitude of R~21.