We discuss the estimation of the interferometric visibility (fringe contrast) for the exozodi survey conducted at the CHARA array with the JouFLU beam combiner. We investigate the use of the statistical median to estimate the uncalibrated visibility from an ensemble of fringe exposures. Under a broad range of operating conditions, numerical simulations indicate that this estimator has a smaller bias compared to other estimators. We also propose an improved method for calibrating visibilities, which not only takes into account the time-interval between observations of calibrators and science targets, but also the uncertainties of the calibrators' raw visibilities. We test our methods with data corresponding to stars that do not display the exozodi phenomenon. The results of our tests show that the proposed method yields smaller biases and errors. The relative reduction in bias and error is generally modest, but can be as high as $\sim 20-40\%$ for the brightest stars of the CHARA data, and statistically significant at the $95\%$ confidence level (CL).

We report the results of high-angular-resolution observations that search for exozodiacal light in a sample of main sequence stars and sub-giants. Using the "jouvence" of the fiber linked unit for optical recombination (JouFLU) at the center for high angular resolution astronomy (CHARA) telescope array, we have observed a total of 44 stars. Out of the 44 stars, 33 are new stars added to the initial, previously published survey of 42 stars performed at CHARA with the fiber linked unit for optical recombiation (FLUOR). Since the start of the survey extension, we have detected a K-band circumstellar excess for six new stars at the ~ 1\% level or higher, four of which are known or candidate binaries, and two for which the excess could be attributed to exozodiacal dust. We have also performed follow-up observations of 11 of the stars observed in the previously published survey and found generally consistent results. We do however detect a significantly larger excess on three of these follow-up targets: Altair, $\upsilon$ And and $\kappa$ CrB. Interestingly, the last two are known exoplanet host stars. We perform a statistical analysis of the JouFLU and FLUOR samples combined, which yields an overall exozodi detection rate of $21.7^{+5.7}_{-4.1}\%$. We also find that the K-band excess in FGK-type stars correlates with the existence of an outer reservoir of cold ($\lesssim 100\,$K) dust at the $99\%$ confidence level, while the same cannot be said for A-type stars.

Future large arrays of telescopes, used as intensity interferometers, can be used to image the surfaces of stars with unprecedented angular resolution. Fast-rotating, hot stars are particularly attractive targets for intensity interferometry since shorter (blue) wavelength observations do not pose additional challenges. Starting from realistic surface brightness simulations of fast-rotating stars, we discuss the capabilities of future intensity interferometers for imaging effects such as gravity darkening and rotational deformation. We find that two-telescope intensity correlation data allow reasonably good imaging of these phenomena, but can be improved with additional higher order (e.g. three-telescope) correlation data, which contain some Fourier phase information.

Most neighboring stars are still detected as point sources and are beyond the angular resolution reach of current observatories. Methods to improve our understanding of stars at high angular resolution are investigated. Air Cherenkov telescopes (ACTs), primarily used for Gamma-ray astronomy, enable us to increase our understanding of the circumstellar environment of a particular system. When used as optical intensity interferometers, future ACT arrays will allow us to detect stars as extended objects and image their surfaces at high angular resolution. Optical stellar intensity interferometry (SII) with ACT arrays, composed of nearly 100 telescopes, will provide means to measure fundamental stellar parameters and also open the possibility of model-independent imaging. A data analysis algorithm is developed and permits the reconstruction of high angular resolution images from simulated SII data. The capabilities and limitations of future ACT arrays used for high angular resolution imaging are investigated via Monte-Carlo simulations. Simple stellar objects as well as stellar surfaces with localized hot or cool regions can be accurately imaged. Finally, experimental efforts to measure intensity correlations are expounded. The functionality of analog and digital correlators is demonstrated. Intensity correlations have been measured for a simulated star emitting pseudo-thermal light, resulting in angular diameter measurements. The StarBase observatory, consisting of a pair of 3 m telescopes separated by 23 m, is described.

Recent proposals have been advanced to apply imaging air Cherenkov telescope arrays to stellar intensity interferometry (SII). Of particular interest is the possibility of model-independent image recovery afforded by the good (u, v)-plane coverage of these arrays, as well as recent developments in phase retrieval techniques. The capabilities of these instruments used as SII receivers have already been explored for simple stellar objects, and here the focus is on reconstructing stellar images with non-uniform radiance distributions. We find that hot stars (T > 6000 K) containing hot and/or cool localized regions (T \sim 500 K) as small as \sim 0.1 mas can be imaged at short wavelengths ({\lambda} = 400 nm).

VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.

Optical stellar intensity interferometry with air Cherenkov telescope arrays, composed of nearly 100 telescopes, will provide means to measure fundamental stellar parameters and also open the possibility of model-independent imaging. In addition to sensitivity issues, a main limitation of image recovery in intensity interferometry is the loss of phase of the complex degree of coherence during the measurement process. Nevertheless, several model-independent phase reconstruction techniques have been developed. Here we implement a Cauchy-Riemann based algorithm to recover images from simulated data. For bright stars (m_v~6) and exposure times of a few hours, we find that scale features such as diameters, oblateness and overall shapes are reconstructed with uncertainties of a few percent. More complex images are also well reconstructed with high degrees of correlation with the pristine image. Results are further improved by using a forward algorithm.

The X-ray binary LSI +61 303, consisting of a main sequence Be star and a compact object has been detected in the TeV range with MAGIC and VERITAS, and showed a clear intensity modulation as a function of the orbital phase. We describe a gamma-ray attenuation model and apply it to this system. Our first result is that interaction of high energy photons with the background radiation produced by the main sequence star alone does not account for the observed modulation. We then include interactions between very high energy radiation and matter, and are able to constrain fundamental parameters of the system such as the mass of the compact object and the density of circumstellar matter around the Be star. In our analysis of the TeV data, we find that the compact object has mass $M_2>2.5M_{\odot}$ at the 99% confidence level, implying it is most likely a black hole. However, we find a column density which conflicts with results from X-ray observations, suggesting that attenuation may not play an important role in the modulation.

Sub milli-arcsecond imaging in the visible band will provide a new perspective in stellar astrophysics. Even though stellar intensity interferometry was abandoned more than 40 years ago, it is capable of imaging and thus accomplishing more than the measurement of stellar diameters as was previously thought. Various phase retrieval techniques can be used to reconstruct actual images provided a sufficient coverage of the interferometric plane is available. Planned large arrays of Air Cherenkov telescopes will provide thousands of simultaneously available baselines ranging from a few tens of meters to over a kilometer, thus making imaging possible with unprecedented angular resolution. Here we investigate the imaging capabilities of arrays such as CTA or AGIS used as Stellar Intensity Interferometry receivers. The study makes use of simulated data as could realistically be obtained from these arrays. A Cauchy-Riemann based phase recovery allows the reconstruction of images which can be compared to the pristine image for which the data were simulated. This is first done for uniform disk stars with different radii and corresponding to various exposure times, and we find that the uncertainty in reconstructing radii is a few percent after a few hours of exposure time. Finally, more complex images are considered, showing that imaging at the sub-milli-arc-second scale is possible.

Experiments are in progress to prepare for intensity interferometry with arrays of air Cherenkov telescopes. At the Bonneville Seabase site, near Salt Lake City, a testbed observatory has been set up with two 3-m air Cherenkov telescopes on a 23-m baseline. Cameras are being constructed, with control electronics for either off- or online analysis of the data. At the Lund Observatory (Sweden), in Technion (Israel) and at the University of Utah (USA), laboratory intensity interferometers simulating stellar observations have been set up and experiments are in progress, using various analog and digital correlators, reaching 1.4 ns time resolution, to analyze signals from pairs of laboratory telescopes.

The imminent availability of large arrays of large light collectors deployed to exploit atmospheric Cherenkov radiation for gamma-ray astronomy at more than 100GeV, motivates the growing interest in application of intensity interferometry in astronomy. Indeed, planned arrays numbering up to one hundred telescopes will offer close to 5,000 baselines, ranging from less than 50m to more than 1000m. Recent and continuing signal processing technology developments reinforce this interest. Revisiting Stellar Intensity Interferometry for imaging is well motivated scientifically. It will fill the short wavelength (B/V bands) and high angular resolution (< 0.1mas) gap left open by amplitude interferometers. It would also constitute a first and important step toward exploiting quantum optics for astronomical observations, thus leading the way for future observatories. In this paper we outline science cases, technical approaches and schedule for an intensity interferometer to be constructed and operated in the visible using gamma-ray astronomy Air Cherenkov Telescopes as receivers.

The space-time correlations of streams of photons can provide fundamentally new channels of information about the Universe. Today's astronomical observations essentially measure certain amplitude coherence functions produced by a source. The spatial correlations of wave fields has traditionally been exploited in Michelson-style amplitude interferometry. However the technology of the past was largely incapable of fine timing resolution and recording multiple beams. When time and space correlations are combined it is possible to achieve spectacular measurements that are impossible by any other means. Stellar intensity interferometry is ripe for development and is one of the few unexploited mechanisms to obtain potentially revolutionary new information in astronomy. As we discuss below, the modern use of stellar intensity interferometry can yield unprecedented measures of stellar diameters, binary stars, distance measures including Cepheids, rapidly rotating stars, pulsating stars, and short-time scale fluctuations that have never been measured before.

When rotation is not taken into account, the measurement of the Gravitational Redshift can provide unique information about the compactness ($M/R$) of the star. Rotation alters the gravitational redshift rendering thereby a unique determination of the compactness parameter impossible. Nevertheless, it can be shown that by using some theoretical input, useful information on, say, the radii of compact rotating objects can still be extracted. Moreover, by measuring the gravitational redshift one can infer the maximum angular velocity of the object. As it is well known, the minimum observed periods of rotation are found in millisecond pulsars. Here we show that millisecond periods are actually a semi-theoretical limit that can be found by General Relativistic arguments corresponding to the maximum angular velocity. We apply our method to compact objects such as pulsars, white dwarfs and neutron stars.