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.

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).