In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales.

In this paper, we study the effects of a recent drop in star formation rate (SFR) on the spectra of epoch of reionization (EoR) galaxies, and the resulting degeneracy with the spectral features produced by extreme Lyman continuum leakage. In order to study these effects in the wavelength range relevant for the upcoming James Webb Space Telescope (JWST), we utilize synthetic spectra of simulated EoR galaxies from cosmological simulations together with synthetic spectra of partially quenched mock galaxies. We find that rapid declines in the SFR of EoR galaxies could seriously affect the applicability of methods that utilize the equivalent width of Balmer lines and the ultraviolet spectral slope to assess the escape fraction of EoR galaxies. In order to determine if the aforementioned degeneracy can be avoided by using the overall shape of the spectrum, we generate mock NIRCam observations and utilize a classification algorithm to identify galaxies that have undergone quenching. We find that while there are problematic cases, JWST/NIRCam or NIRSpec should be able to reliably identify galaxies with redshifts $z\sim 7$ that have experienced a significant decrease in the SFR (by a factor 10-100) in the past 50-100 Myr with a success rate $\gtrsim 85\%$. We also find that uncertainties in the dust-reddening effects on EoR galaxies significantly affect the performance of the results of the classification algorithm. We argue that studies that aim to characterize the dust extinction law most representative in the EoR would be extremely useful.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

Recent observations of galaxies at $z \gtrsim 7$, along with the low value of the electron scattering optical depth measured by the Planck mission, make galaxies plausible as dominant sources of ionizing photons during the epoch of reionization. However, scenarios of galaxy-driven reionization hinge on the assumption that the average escape fraction of ionizing photons is significantly higher for galaxies in the reionization epoch than in the local Universe. The NIRSpec instrument on the James Webb Space Telescope (JWST) will enable spectroscopic observations of large samples of reionization-epoch galaxies. While the leakage of ionizing photons will not be directly measurable from these spectra, the leakage is predicted to have an indirect effect on the spectral slope and the strength of nebular emission lines in the rest-frame ultraviolet and optical. Here, we apply a machine learning technique known as lasso regression on mock JWST/NIRSpec observations of simulated $z=7$ galaxies in order to obtain a model that can predict the escape fraction from JWST/NIRSpec data. Barring systematic biases in the simulated spectra, our method is able to retrieve the escape fraction with a mean absolute error of $\Delta f_{\mathrm{esc}} \approx 0.12$ for spectra with $S/N\approx 5$ at a rest-frame wavelength of 1500 {\AA} for our fiducial simulation. This prediction accuracy represents a significant improvement over previous similar approaches.

The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.

A proposed method for dealing with foreground emission in upcoming 21-cm observations from the epoch of reionization is to limit observations to an uncontaminated window in Fourier space. Foreground emission can be avoided in this way, since it is limited to a wedge-shaped region in $k_{\parallel}, k_{\perp}$ space. However, the power spectrum is anisotropic owing to redshift-space distortions from peculiar velocities. Consequently, the 21-cm power spectrum measured in the foreground avoidance window---which samples only a limited range of angles close to the line-of-sight direction---differs from the full spherically-averaged power spectrum which requires an average over \emph{all} angles. In this paper, we calculate the magnitude of this "wedge bias" for the first time. We find that the bias is strongest at high redshifts, where measurements using foreground avoidance will over-estimate the power spectrum by around 100 per cent, possibly obscuring the distinctive rise and fall signature that is anticipated for the spherically-averaged 21-cm power spectrum. In the later stages of reionization, the bias becomes negative, and smaller in magnitude ($\lesssim 20$ per cent). The effect shows only a weak dependence on spatial scale and reionization topology.

This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm rad~m^{-2}}$ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350~{\rm MHz}. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm \mu G}$ and its spatial variation across the 3C196 field is $0.1~{\rm \mu G}$. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density ($n_e B_\parallel>6.2~{\rm cm^{-3}\mu G}$) compared to its surroundings.

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$.

We study the prospects for constraining the ionized fraction of the intergalactic medium (IGM) at $z>6$ with the next generation of large Ly$\alpha$ emitter surveys. We make predictions for the upcoming Subaru Hyper Suprime-Cam (HSC) Ly$\alpha$ survey and a hypothetical spectroscopic survey performed with the James Webb Space Telescope (JWST). Considering various scenarios where the observed evolution of the Ly$\alpha$ luminosity function of Ly$\alpha$ emitters at $z>6$ is explained partly by an increasingly neutral IGM and partly by intrinsic galaxy evolution, we show how clustering measurements will be able to distinguish between these scenarios. We find that the HSC survey should be able to detect the additional clustering induced by a neutral IGM if the global IGM neutral fraction is greater than $\sim$20 per cent at $z=6.5$. If measurements of the Ly$\alpha$ equivalent widths (EWs) are also available, neutral fractions as small as 10 per cent may be detectable by looking for correlation between the EW and the local number density of objects. In this case, if it should turn out that the IGM is significantly neutral at $z=6.5$ and the intrinsic EW distribution is relatively narrow, the observed EWs can also be used to construct a map of the locations and approximate sizes of the largest ionized regions. For the JWST survey, the results appear a bit less optimistic. Since such surveys probe a large range of redshifts, the effects of the IGM will be mixed up with any intrinsic galaxy evolution that is present, making it difficult to disentangle the effects. However, we show that a survey with the JWST will have a possibility of observing a large group of galaxies at $z\sim7$, which would be a strong indication of a partially neutral IGM.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

We present a detailed comparison of three different simulations of the epoch of reionization (EoR). The radiative transfer simulation (${\rm C}^2$-RAY) among them is our benchmark. Radiative transfer codes can produce realistic results, but are computationally expensive. We compare it with two semi-numerical techniques: one using the same halos as ${\rm C}^2$-RAY as its sources (Sem-Num), and one using a conditional Press-Schechter scheme (CPS+GS). These are vastly more computationally efficient than ${\rm C}^2$-RAY, but use more simplistic physical assumptions. We evaluate these simulations in terms of their ability to reproduce the history and morphology of reionization. We find that both Sem-Num and CPS+GS can produce an ionization history and morphology that is very close to ${\rm C}^2$-RAY, with Sem-Num performing slightly better compared to CPS+GS. We also study different redshift space observables of the 21-cm signal from EoR: the variance, power spectrum and its various angular multipole moments. We find that both semi-numerical models perform reasonably well in predicting these observables at length scales relevant for present and future experiments. However, Sem-Num performs slightly better than CPS+GS in producing the reionization history, which is necessary for interpreting the future observations.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($\Delta z$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $\Delta z = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $\Delta z = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively.

The fraction of ionizing photons (fesc) that escape from z>6 galaxies is an important parameter when assessing the role of these objects in the reionization of the Universe, but the opacity of the intergalactic medium precludes a direct measurement of fesc for individual galaxies at these epochs. We argue, that since fesc regulates the impact of nebular emission on the spectra of galaxies, it should nonetheless be possible to indirectly probe fesc well into the reionization epoch. As a first step, we demonstrate that by combining measurements of the rest-frame UV slope beta with the equivalent width of the Hb emission line, galaxies with very high Lyman continuum escape fractions (fesc>0.5) should be identifiable up to z~9 through spectroscopy with the upcoming James Webb Space Telescope (JWST). By targeting strongly lensed galaxies behind low-redshift galaxy clusters, JWST spectra of sufficiently good quality can be obtained for M(1500)<-16.0 galaxies at z~7 and for M(1500)<-17.5 galaxies at z~9. Dust-obscured star formation may complicate the analysis, but supporting observations with ALMA or the planned SPICA mission may provide useful constraints on the dust properties of these galaxies.

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

Dark halo substructure may reveal itself through secondary, small-scale gravitational lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we explore the prospects of using Very Long Baseline Interferometry (VLBI) observations of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions produced by various forms of dark substructures in the 1e3-1e8 Msolar mass range. We present lensing simulations relevant for the angular resolutions attainable with the existing European VLBI Network (EVN), the global VLBI array, and an upcoming observing mode in which the Atacama Large Millimeter Array (ALMA) is connected to the global VLBI array. While observations of this type would not be sensitive to standard cold dark matter subhalos, they can be used to detect more compact forms of halo substructure predicted in alternative structure formation scenarios. By mapping ~5 strongly lensed systems, it should be possible to detect or robustly rule out primordial black holes in the 1e3-1e6 Msolar mass range if they constitute >1% percent of the dark matter in these lenses. Ultracompact minihalos are harder to detect using this technique, but 1e6-1e8 Msolar ultracompact minihalos could in principle be detected if they constitute >10% of the dark matter.

We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart.

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.

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.

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.

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.

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.

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.

We have made an asteroseismic analysis of the variable blue stragglers in the open cluster M67. The data set consists of photometric time series from eight sites using nine 0.6-2.1 meter telescopes with a time baseline of 43 days. In two stars, EW Cnc and EX Cnc, we detect the highest number of frequencies (41 and 26) detected in delta Scuti stars belonging to a stellar cluster, and EW Cnc has the second highest number of frequencies detected in any delta Scuti star. We have computed a grid of pulsation models that take the effects of rotation into account. The distribution of observed and theoretical frequencies show that in a wide frequency range a significant fraction of the radial and non-radial low-degree modes are excited to detectable amplitudes. Despite the large number of observed frequencies we cannot constrain the fundamental parameters of the stars. To make progress we need to identify the degrees of some of the modes either from multi-colour photometry or spectroscopy.

This is the first step in a project to combine studies of eclipsing binaries and oscillating stars to probe the interior of Blue Stragglers (BS). This may imply a way to discriminate observationally between different birth mechanisms of BS stars. We study the open cluster NGC2506 which contains oscillating BS stars and detached eclipsing binaries for which accurate parameters can be derived. This will tightly constrain the cluster isochrone and provide an absolute mass, radius and luminosity-scale for the cluster stars along with the cluster age, metallicity and distance. The present work focuses on obtaining the light curves of the binaries and determine their orbital periods, on obtaining power spectra of the oscillating BS stars to select targets for follow-up studies, and on searching for gamma Doradus type variables which are also expected to be present in the cluster. With a two-colour, dual-site photometric campaign we obtained 3120 CCD-images of NGC2506 spread over four months. We analysed the BI time-series of the oscillating stars and used simulations to derive statistical uncertainties of the resulting frequencies, amplitudes and phases. A preliminary mode-identification was performed using frequency ratios for the oscillating BS stars, and amplitude ratios and phase differences for a population of newly detected gamma Doradus stars. We quadrupled the number of known variables in NGC2506 by discovering 3 new oscillating BS stars, 15 gamma Doradus stars and four new eclipsing binaries. The orbital periods of 2 known, detached eclipsing binaries were derived. We discovered a BS star with both p-mode and g-mode variability and we confronted our gamma Doradus observations with state-of-the-art seismic models, but found significant discrepancy between theory and observations.

Asteroseismology of stellar clusters is potentially a powerful tool. The assumption of a common age, distance, and chemical composition provides constraints on each cluster member, which significantly improves the asteroseismic output. Driven by this great potential, we carried out multi-site observations aimed at detecting solar-like oscillations in the red giant stars in the open cluster M67 (NGC 2682) (Stello et al. 2006). Here we present the first analysis of our data, which show evidence of excess power in the Fourier spectra, shifting to lower frequencies for more luminous stars, consistent with expectations from oscillations. If the observed power excesses were due to stellar oscillations, this result would show great prospects for asteroseismology in stellar clusters.

Measuring solar-like oscillations in an ensemble of stars in a cluster, holds promise for testing stellar structure and evolution more stringently than just fitting parameters to single field stars. The most ambitious attempt to pursue these prospects was by Gilliland et al. (1993) who targeted 11 turn-off stars in the open cluster M67 (NGC 2682), but the oscillation amplitudes were too small (<20micromag) to obtain unambiguous detections. Like Gilliland et al. (1993) we also aim at detecting solar-like oscillations in M67, but we target red giant stars with expected amplitudes in the range 50-500micromag and periods of 1 to 8 hours. We analyse our recently published photometry measurements, obtained during a six-week multisite campaign using nine telescopes around the world. The observations are compared with simulations and with estimated properties of the stellar oscillations. Noise levels in the Fourier spectra as low as 27micromag are obtained for single sites, while the combined data reach 19micromag, making this the best photometric time series of an ensemble of red giant stars. These data enable us to make the first test of the scaling relations (used to estimate frequency and amplitude) with an homogeneous ensemble of stars. The detected excess power is consistent with the expected signal from stellar oscillations, both in terms of its frequency range and amplitude. However, our results are limited by apparent high levels of non-white noise, which cannot be clearly separated from the stellar signal.

We report on an ambitious multi-site campaign aimed at detecting stellar variability, particularly solar-like oscillations, in the red giant stars in the open cluster M67 (NGC 2682). During the six-week observing run, which comprised 164 telescope nights, we used nine 0.6-m to 2.1-m class telescopes located around the world to obtain uninterrupted time-series photometry. We outline here the data acquisition and reduction, with emphasis on the optimisation of the signal-to-noise of the low amplitude (50-500 micromag) solar-like oscillations. This includes a new and efficient method for obtaining the linearity profile of the CCD response at ultra high precision (~10 parts per million). The noise in the final time series is 0.50 mmag per minute integration for the best site, while the noise in the Fourier spectrum of all sites combined is 20 micromag. In addition to the red giant stars, this data set proves to be very valuable for studying high-amplitude variable stars such as eclipsing binaries, W UMa systems and delta Scuti stars.

Recent theoretical yields and chemical evolution models demonstrate that intermediate-mass AGB stars cannot reproduce the observed abundance distributions of O, Na, Mg, and Al. As a further observational test of this finding, we present elemental abundance ratios [X/Fe] for 20 elements in 38 bright giants of the globular cluster NGC 6752. Our mean abundance ratios [X/Fe] are in good agreement with previous studies of this cluster and are also consistent with other globular clusters and field stars at the same metallicity. The mean abundance ratios [Ba/Eu] and [La/Eu] exhibit values, in agreement with field stars at the same metallicity, that lie approximately midway between the pure r-process and the solar (s-process + r-process) mix, indicating that AGB stars have played a role in the chemical evolution of the proto-cluster gas. For the first time, we find possible evidence for an abundance variation for elements heavier than Al in this cluster. We find a correlation between [Si/Fe] and [Al/Fe] which is consistent with the abundance anomalies being synthesized via proton captures at high temperatures. Leakage from the Mg-Al chain into 28Si may explain the Si excess in stars with the highest [Al/Fe]. We identify correlations between [Y/Fe] and [Al/Fe], [Zr/Fe] and [Al/Fe], and [Ba/Fe] and [Al/Fe] suggesting that Y, Zr, and Ba abundances may increase by about 0.1 dex as Al increases by about 1.3 dex. While the correlations are statistically significant, the amplitudes of the variations are small. If the small variations in Y, Zr, and Ba are indeed real, then the synthesis of the Al anomalies must have taken place within an unknown class of stars that also ran the s-process. [Abridged]