A long-held vision has been to realize diffraction-limited optical aperture synthesis over kilometer baselines. This will enable imaging of stellar surfaces and their environments, and reveal interacting gas flows in binary systems. An opportunity is now opening up with the large telescope arrays primarily erected for measuring Cherenkov light in air induced by gamma rays. With suitable software, such telescopes could be electronically connected and also used for intensity interferometry. Second-order spatial coherence of light is obtained by cross correlating intensity fluctuations measured in different pairs of telescopes. With no optical links between them, the error budget is set by the electronic time resolution of a few nanoseconds. Corresponding light-travel distances are approximately one meter, making the method practically immune to atmospheric turbulence or optical imperfections, permitting both very long baselines and observing at short optical wavelengths. Previous theoretical modeling has shown that full images should be possible to retrieve from observations with such telescope arrays. This project aims at verifying diffraction-limited imaging experimentally with groups of detached and independent optical telescopes. In a large optics laboratory, artificial stars were observed by an array of small telescopes. Using high-speed photon-counting solid-state detectors, intensity fluctuations were cross-correlated over up to 180 baselines between pairs of telescopes, producing coherence maps across the interferometric Fourier-transform plane. These measurements were used to extract parameters about the simulated stars, and to reconstruct their two-dimensional images. As far as we are aware, these are the first diffraction-limited images obtained from an optical array only linked by electronic software, with no optical connections between the telescopes.

Highest resolution imaging in astronomy is achieved by interferometry, connecting telescopes over increasingly longer distances, and at successively shorter wavelengths. Here, we present the first diffraction-limited images in visual light, produced by an array of independent optical telescopes, connected electronically only, with no optical links between them. With an array of small telescopes, second-order optical coherence of the sources is measured through intensity interferometry over 180 baselines between pairs of telescopes, and two-dimensional images reconstructed. The technique aims at diffraction-limited optical aperture synthesis over kilometre-long baselines to reach resolutions showing details on stellar surfaces and perhaps even the silhouettes of transiting exoplanets. Intensity interferometry circumvents problems of atmospheric turbulence that constrain ordinary interferometry. Since the electronic signal can be copied, many baselines can be built up between dispersed telescopes, and over long distances. Using arrays of air Cherenkov telescopes, this should enable the optical equivalent of interferometric arrays currently operating at radio wavelengths.

Optical interferometry has been successful at achieving milliarcsecond resolution on bright stars. Imaging performance can improve greatly by increasing the number of baselines, which has motivated proposals to build large (~ 100 m) optical interferometers with tens to hundreds of telescopes. It is also desirable to adaptively correct atmospheric turbulence to obtain direct phased images of astrophysical sources. When a natural guide star is not available, we investigate the feasibility of using a modified laser-guide-star technique that is suitable for large diluted apertures. The method consists of using sub-sets of apertures to create an array of artificial stars in the sodium layer and collecting back-scattered light with the same sub-apertures. We present some numerical and laboratory simulations that quantify the requirements and sensitivity of the technique.

Stellar intensity interferometers will achieve stellar imaging with a tenth of a milli- arcsecond resolution in the optical band by taking advantage of the large light collect- ing area and broad range of inter-telescope distances offered by future gamma-ray Air Cherenkov Telescope (ACT) arrays. Up to now, studies characterizing the capabilities of intensity interferometers using ACTs have not accounted for realistic effects such as telescope mirror extension, detailed photodetector time response, excess noise, and night sky contamination. In this paper, we present the semi-classical quantum optics Monte-Carlo simulation we developed in order to investigate these experimental limi- tations. In order to validate the simulation algorithm, we compare our first results to models for sensitivity and signal degradation resulting from mirror extension, pulse shape, detector excess noise, and night sky contamination.

Using kilometric arrays of air Cherenkov telescopes, intensity interferometry may increase the spatial resolution in optical astronomy by an order of magnitude, enabling images of rapidly rotating stars with structures in their circumstellar disks and winds, or mapping out patterns of nonradial pulsations across stellar surfaces. Intensity interferometry (pioneered by Hanbury Brown and Twiss) connects telescopes only electronically, and is practically insensitive to atmospheric turbulence and optical imperfections, permitting observations over long baselines and through large airmasses, also at short optical wavelengths. The required large telescopes with very fast detectors are becoming available as arrays of air Cherenkov telescopes, distributed over a few square km. Digital signal handling enables very many baselines to be synthesized, while stars are tracked with electronic time delays, thus synthesizing an optical interferometer in software. Simulated observations indicate limiting magnitudes around m(v)=8, reaching resolutions ~30 microarcsec in the violet. The signal-to-noise ratio favors high-temperature sources and emission-line structures, and is independent of the optical passband, be it a single spectral line or the broad spectral continuum. Intensity interferometry provides the modulus (but not phase) of any spatial frequency component of the source image; for this reason image reconstruction requires phase retrieval techniques, feasible if sufficient coverage of the interferometric (u,v)-plane is available. Experiments are in progress; test telescopes have been erected, and trials in connecting large Cherenkov telescopes have been carried out. This paper reviews this interferometric method in view of the new possibilities offered by arrays of air Cherenkov telescopes, and outlines observational programs that should become realistic already in the rather near future.

We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.

With its unprecedented light-collecting area for night-sky observations, the Cherenkov Telescope Array (CTA) holds great potential for also optical stellar astronomy, in particular as a multi-element intensity interferometer for realizing imaging with sub-milliarcsecond angular resolution. Such an order-of-magnitude increase of the spatial resolution achieved in optical astronomy will reveal the surfaces of rotationally flattened stars with structures in their circumstellar disks and winds, or the gas flows between close binaries. Image reconstruction is feasible from the second-order coherence of light, measured as the temporal correlations of arrival times between photons recorded in different telescopes. This technique (once pioneered by Hanbury Brown and Twiss) connects telescopes only with electronic signals and is practically insensitive to atmospheric turbulence and to imperfections in telescope optics. Detector and telescope requirements are very similar to those for imaging air Cherenkov observatories, the main difference being the signal processing (calculating cross correlations between single camera pixels in pairs of telescopes). Observations of brighter stars are not limited by sky brightness, permitting efficient CTA use during also bright-Moon periods. While other concepts have been proposed to realize kilometer-scale optical interferometers of conventional amplitude (phase-) type, both in space and on the ground, their complexity places them much further into the future than CTA, which thus could become the first kilometer-scale optical imager in astronomy.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.

It has been suggested that the ejection to interplanetary space of terrestrial crustal material, accelerated in a large impact, may result in the interchange of biological material between Earth and other Solar System bodies. In this paper, we analyze the fate of debris ejected from Earth by means of numerical simulations of the dynamics of a large collection of test particles. This allows us to determine the probability and conditions for the collision of ejecta with other planets of the Solar System. We also estimate the amount of particles falling-back to Earth as a function of time after being ejected. We find that, in general, the collision rates of Earth ejecta with Venus and the Moon, as well as the fall-back rates, are consistent with results reported in the literature. By considering a larger number of particles than in all previous calculations we have also determined directly the collision probability with Mars and, for the first time, computed collision probabilities with Jupiter. We find that the collision probability with Mars is greater than values determined from collision cross section estimations previously reported.

Intensity interferometry permits very long optical baselines and the observation of sub-milliarcsecond structures. Using planned kilometric arrays of air Cherenkov telescopes at short wavelengths, intensity interferometry may increase the spatial resolution achieved in optical astronomy by an order of magnitude, inviting detailed studies of the shapes of rapidly rotating hot stars with structures in their circumstellar disks and winds, or mapping out patterns of nonradial pulsations across stellar surfaces. Signal-to-noise in intensity interferometry favors high-temperature sources and emission-line structures, and is independent of the optical passband, be it a single spectral line or the broad spectral continuum. Prime candidate sources have been identified among classes of bright and hot stars. Observations are simulated for telescope configurations envisioned for large Cherenkov facilities, synthesizing numerous optical baselines in software, confirming that resolutions of tens of microarcseconds are feasible for numerous astrophysical targets.

Kilometric-scale optical imagers seem feasible to realize by intensity interferometry, using telescopes primarily erected for measuring Cherenkov light induced by gamma rays. Planned arrays envision 50--100 telescopes, distributed over some 1--4 km$^2$. Although array layouts and telescope sizes will primarily be chosen for gamma-ray observations, also their interferometric performance may be optimized. Observations of stellar objects were numerically simulated for different array geometries, yielding signal-to-noise ratios for different Fourier components of the source images in the interferometric $(u,v)$-plane. Simulations were made for layouts actually proposed for future Cherenkov telescope arrays, and for subsets with only a fraction of the telescopes. All large arrays provide dense sampling of the $(u,v)$-plane due to the sheer number of telescopes, irrespective of their geographic orientation or stellar coordinates. However, for improved coverage of the $(u,v)$-plane and a wider variety of baselines (enabling better image reconstruction), an exact east-west grid should be avoided for the numerous smaller telescopes, and repetitive geometric patterns avoided for the few large ones. Sparse arrays become severely limited by a lack of short baselines, and to cover astrophysically relevant dimensions between 0.1--3 milliarcseconds in visible wavelengths, baselines between pairs of telescopes should cover the whole interval 30--2000 m.