Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

The recent demonstration of a real-time direct imaging radio interferometry correlator represents a new capability in radio astronomy. However wide field imaging with this method is challenging since wide-field effects and array non-coplanarity degrade image quality if not compensated for. Here we present an alternative direct imaging correlation strategy using a Direct Fourier Transform (DFT), modelled as a linear operator facilitating a matrix multiplication between the DFT matrix and a vector of the electric fields from each antenna. This offers perfect correction for wide field and non-coplanarity effects. When implemented with data from the Long Wavelength Array (LWA), it offers comparable computational performance to previously demonstrated direct imaging techniques, despite having a theoretically higher floating point cost. It also has additional benefits, such as imaging sparse arrays and control over which sky co-ordinates are imaged, allowing variable pixel placement across an image. It is in practice a highly flexible and efficient method of direct radio imaging when implemented on suitable arrays. A functioning Electric Field Direct imaging architecture using the DFT is presented, alongside an exploration of techniques for wide-field imaging similar to those in visibility based imaging, and an explanation of why they do not fit well to imaging directly with the digitized electric field data. The DFT imaging method is demonstrated on real data from the LWA telescope, alongside a detailed performance analysis, as well as an exploration of its applicability to other arrays.

In radio interferometry imaging, the gridding procedure of convolving visibilities with a chosen gridding function is necessary to transform visibility values into uniformly sampled grid points. We propose here a parameterised family of "least-misfit gridding functions" which minimise an upper bound on the difference between the DFT and FFT dirty images for a given gridding support width and image cropping ratio. When compared with the widely used spheroidal function with similar parameters, these provide more than 100 times better alias suppression and RMS misfit reduction over the usable dirty map. We discuss how appropriate parameter selection and tabulation of these functions allow for a balance between accuracy, computational cost and storage size. Although it is possible to reduce the errors introduced in the gridding or degridding process to the level of machine precision, accuracy comparable to that achieved by CASA requires only a lookup table with 300 entries and a support width of 3, allowing for a greatly reduced computation cost for a given performance.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

Methods currently in use for locating and characterising sources in radio interferometry maps are designed for processing images, and require interferometric maps to be preprocessed so as to resemble conventional images. We demonstrate a Bayesian code - BaSC - that takes into account the interferometric visibility data despite working with more computationally manageable image domain data products. This method is better able to discriminate nearby sources than the commonly used SExtractor, and has potential even in more complicated cases. We also demonstrate the correctness of the Bayesian resolving formula for simulated data, and its implications for source discrimination at distances below the full width half maximum of the restoring beam.

We present a technique to automatically minimise the re-computation when a data analysis program is iteratively changed, or added to, as is often the case in exploratory data analysis in astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and referentially transparent tasks. We describe the implementation of this technique for the CASA radio astronomy data reduction package. We also propose a technique for optimising efficiency of storage of memoized intermediate data products using copy-on-write and block level de-duplication and measure their practical efficiency. We find the minimal recomputation technique improves the efficiency of data analysis while reducing the possibility for user error and improving the reproducibility of the final result. It also aids exploratory data analysis on batch-schedule cluster computer systems.

I show that a software framework intended primarily for training of neural networks, PyTorch, is easily applied to a general function minimisation problem in science. The qualities of PyTorch of ease-of-use and very high efficiency are found to be applicable in this domain and lead to two orders of magnitude improvement in time-to-solution with very small software engineering effort.

Detecting neutral Hydrogen (H\,{\sc i}) via the 21~cm line emission from the intergalactic medium at $z\gtrsim 6$ has been identified as one of the most promising probes of the epoch of cosmic reionization -- a major phase transition of the Universe. However, these studies face severe challenges imposed by the bright foreground emission from cosmic objects. Current techniques require precise instrumental calibration to separate the weak H\,{\sc i} line signal from the foreground continuum emission. We propose to mitigate this calibration requirement by using measurements of the interferometric bi-spectrum phase. Bi-spectrum phase is unaffected by antenna-based direction-independent calibration errors and hence for a compact array it depends on the sky brightness distribution only (subject to the usual thermal-like noise). We show that the bi-spectrum phase of foreground synchrotron continuum has a characteristically smooth spectrum relative to the cosmological line signal. The two can be separated effectively by exploiting this spectral difference using Fourier techniques, while eliminating the need for precise antenna-based calibration of phases introduced by the instrument, and the ionosphere, inherent in existing approaches. Using fiducial models for continuum foregrounds, and for the cosmological H\,{\sc i} signal, we show the latter should be detectable in bi-spectrum phase spectra, with reasonable significance at $|k_\parallel| \gtrsim 0.5\,h$~Mpc$^{-1}$, using existing instruments. Our approach will also benefit other H\,{\sc i} intensity mapping experiments that face similar challenges, such as those measuring Baryon Acoustic Oscillations (BAO).

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

We present a technique to automatically minimise the re-computation when a data processing program is iteratively changed, or added to, as is often the case in exploratory data analysis in radio astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and referentially transparent tasks. We describe a prototype implementation for the CASA data reduction package. This technique improves the efficiency of data analysis while reducing the possibility for user error and improving the reproducibility of the final result.

HERA is a highly-redundant transit interferometer with 14\,m-diameter parabolic dish elements. We exploit the fact that the Galactic centre transits through the main beam of the telescope to attempt a conventional self-calibration approach to imaging and calibration. The Galactic centre provides a bright source which, we show, can be approximated as a point source sufficiently well to initialise the self-calibration loop and derive initial delays and antenna frequency-independent phases. Subsequent iteration using a more complex sky model derived from the data itself then converges to a reasonable bandpass calibration. The calibration solutions have good stability properties. We show therefore that the conventional self-calibration is a feasible parallel approach in addition to the redundant calibration already planned for HERA. The conventional imaging and calibration is useful as a cross-check to the alternatives being pursued in the HERA project, as a way of quantifying the performance of the hardware on the ground (and potentially identifying problems) and as a route to imaging and removing brighter continuum sources before power spectrum analysis.

This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

Fluctuating properties of the atmosphere, and in particular its water vapour content, give rise to phase fluctuations of astronomical signals which, if uncorrected, lead to rapid deterioration of performance of (sub)-mm interferometers on long baselines. The Atacama Large Millimetre/submillimeter Array (ALMA) uses a 183 GHz Water Vapour Radiometer (WVR) system to help correct these fluctuations and provide much improved performance on long baselines and at high frequencies. Here we describe the design of the overall ALMA WVR system, the choice of design parameters and the data processing strategy. We also present results of initial tests that demonstrate both the large improvement in phase stability that can be achieved and the very low contribution to phase noise from the WVRs. Finally, we describe briefly the main limiting factors to the accuracy of phase correction seen in these initial tests; namely, the degrading influence of cloud and the residual phase fluctuations that are most likely to be due to variations in the density of the dry component of the air

Atacama Large Millimeter/submillimeter Array (ALMA) will be the world largest mm/submm interferometer, and currently the Early Science is ongoing, together with the commissioning and science verification (CSV). Here we present a study of the temporal phase stability of the entire ALMA system from antennas to the correlator. We verified the temporal phase stability of ALMA using data, taken during the last two years of CSV activities. The data consist of integrations on strong point sources (i.e., bright quasars) at various frequency bands, and at various baseline lengths (up to 600 m). From the observations of strong quasars for a long time (from a few tens of minutes, up to an hour), we derived the 2-point Allan Standard Deviation after the atmospheric phase correction using the 183 GHz Water Vapor Radiometer (WVR) installed in each 12 m antenna, and confirmed that the phase stability of all the baselines reached the ALMA specification. Since we applied the WVR phase correction to all the data mentioned above, we also studied the effectiveness of the WVR phase correction at various frequencies, baseline lengths, and weather conditions. The phase stability often improves a factor of 2 - 3 after the correction, and sometimes a factor of 7 improvement can be obtained. However, the corrected data still displays an increasing phase fluctuation as a function of baseline length, suggesting that the dry component (e.g., N2 and O2) in the atmosphere also contributes the phase fluctuation in the data, although the imperfection of the WVR phase correction cannot be ruled out at this moment.

This memo describes the software engineering and technical details of the design and implementation of the wvrgcal program and associated libraries. This program performs off-line correction of atmospheric phase fluctuations in ALMA observations, using the 183 GHz Water Vapour Radiometers (WVRs) installed on the ALMA 12 m dishes. The memo can be used as a guide for detailed study of the source code of the program for purposes of further development or maintenance.

Radio emission at around 90 GHz from star-forming galaxies is expected to be strongly dominated by the free-free component due to ionising radiation from massive, short-lived, stars. We present high surface-brightness sensitivity observations at 90 GHz of the nearby star-forming galaxy Messier 66 with resolution of about 9 arcsec (corresponding to a physical scale of about 500 pc) and analyse these observations in combination with archival lower frequency radio and mid-infrared measurements. For the four regions for which the observations support our models we find that the free-free component indeed dominates the emission at 90 GHz, making up 76--90 per cent of the luminosity at this frequency but with the data also consistent with all of the emission being due to free-free. The estimates of free-free luminosities are also consistent, within measurement and decomposition errors, with star-formation rates derived from lower radio frequencies and mid-infrared observations. In our analysis we consider both power-law and curved spectra for the synchrotron component but do not find evidence to support the curved model in preference to the power-law.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

The success of any ALMA phase-calibration strategy, which incorporates phase transfer, depends on a good understanding of how the atmospheric path delay changes with frequency (e.g. Holdaway & Pardo 2001). We explore how the wet dispersive path delay varies for realistic atmospheric conditions at the ALMA site using the ATM transmission code. We find the wet dispersive path delay becomes a significant fraction (>5 per cent) of the non-dispersive delay for the high-frequency ALMA bands (>160 GHz, Bands 5 to 10). Additionally, the variation in dispersive path delay across ALMA's 4-GHz contiguous bandwidth is not significant except in Bands 9 and 10. The ratio of dispersive path delay to total column of water vapour does not vary significantly for typical amounts of water vapour, water vapour scale heights and ground pressures above Chajnantor. However, the temperature profile and particularly the ground-level temperature are more important. Given the likely constraints from ALMA's ancillary calibration devices, the uncertainty on the dispersive-path scaling will be around 2 per cent in the worst case and should contribute about 1 per cent overall to the wet path fluctuations at the highest frequencies.

I describe an approach to fitting and comparison of radio spectra based on Bayesian analysis and realised using a new implementation of the nested sampling algorithm. Such an approach improves on the commonly used maximum-likelihood fitting of radio spectra by allowing objective model selection, calculation of the full probability distributions of the model parameters and provides a natural mechanism for including information other than the measured spectra through priors. In this paper I cover the theoretical background, the algorithms used and the implementation details of the computer code. I also briefly illustrate the method with some previously published data for three near-by galaxies. In forthcoming papers we will present the results of applying this analysis larger data sets, including some new observations, and the physical conclusions that can be made. The computer code as well as the overall approach described here may also be useful for analysis of other multi-chromatic broad-band observations and possibly also photometric redshift estimation. All of the code is publicly available, licensed under the GNU General Public License, at http://www.mrao.cam.ac.uk/~bn204/galevol/speca/index.html

By observing bright and compact astronomical sources while also taking data with the 183 GHz Water Vapour Radiometers, ALMA will be able to measure the `empirical' relationship between fluctuations in the phase of the astronomical signal and the fluctuations of sky brightness around 183 GHz. Simulations of such measurements assuming only thermal noise in the astronomical and WVR receivers are presented and it is shown that accurate determination of the empirical relationship should be possible in a relatively short time. It is then proposed that the best way of using these empirical coefficients is to include them as a constraint on a physical model of the atmosphere -- this allows them to be used for longer period of time, increasing the efficiency of observing. This approach fits naturally into the analysis framework presented in the previous memo, which has now been extended to implement it. The technique is illustrated via simulations and on a short data set collected at the SMA.

We present a Bayesian approach to calculating the coefficients that convert the outputs of ALMA 183 GHz water-vapour radiometers into estimates of path fluctuations which can then be used to correct the observed interferometric visibilities. The key features of the approach are a simple, thin-layer, three-parameter model of the atmosphere; using the absolute measurements from the radiometers to constrain the model; priors to incorporate physical constraints and ancillary information; and a Markov Chain Monte Carlo characterisation of the posterior distribution including full distributions for the phase correction coefficients. The outcomes of the procedure are therefore estimates of the coefficients and their confidence intervals. We illustrate the technique with simulations showing some degeneracies that can arise and the importance of priors in tackling them. We then apply the technique to an hour-long test observation at the Sub-Millimetre Array and find that the technique is stable and that, in this case, its performance is close to optimal. The modelling is described in detail in the appendices and all of the implementation source code is made publicly available under the GPL.

We define a volume limited sample of over 14,000 early-type galaxies (ETGs) selected from data release six of the Sloan Digital Sky Survey. The density of environment of each galaxy is robustly measured. By comparing narrow band spectral line indices with recent models of simple stellar populations (SSPs) we investigate trends in the star formation history as a function of galaxy mass (velocity dispersion), density of environment and galactic radius. We find that age, metallicity and alpha-enhancement all increase with galaxy mass and that field ETGs are younger than their cluster counterparts by ~2 Gyr. We find negative radial metallicity gradients for all masses and environments, and positive radial age gradients for ETGs with velocity dispersion over 180 km/s. Our results are qualitatively consistent with a relatively simple picture for ETG evolution in which the low-mass halos accreted by a proto-ETG contained not only gas but also a stellar population. This fossil population is preferentially found at large radii in massive ETGs because the stellar accretions were dissipationless. We estimate that the typical, massive ETG should have been assembled at z < 3.5. The process is similar in the cluster and the field but occurred earlier in dense environments.

We present a semi-empirical model for the infrared emission of dust around star-forming sites in galaxies. Our approach combines a simple model of radiative transfer in dust clouds with a state-of-the-art model of the microscopic optical properties of dust grains pioneered by Draine & Li. In combination with the Starburst 99 stellar spectral synthesis package, this framework is able to produce synthetic spectra for galaxies which extend from the Lyman limit through to the far-infrared. We use it to probe how model galaxy spectra depend upon the physical characteristics of their dust grain populations, and upon the energy sources which heat that dust. We compare the predictions of our model with the 8- and 24-micron luminosities of sources in the Spitzer First Look Survey, and conclude by using the models to analyse the relative merits of various colour diagnostics in distinguishing systems out to a redshift of 2 with ongoing star formation from those with only old stellar populations.

We present a technique for the accurate estimation of large-scale errors in an antenna surface using astronomical sources and detectors. The technique requires several out-of-focus images of a compact source and the signal-to-noise ratio needs to be good but not unreasonably high. For a given pattern of surface errors, the expected form of such images can be calculated directly. We show that it is possible to solve the inverse problem of finding the surface errors from the images in a stable manner using standard numerical techniques. To do this we describe the surface error as a linear combination of a suitable set of basis functions (we use Zernike polynomials). We present simulations illustrating the technique and in particular we investigate the effects of receiver noise and pointing errors. Measurements of the 15-m James Clerk Maxwell telescope made using this technique are presented as an example. The key result is that good measurements of errors on large spatial scales can be obtained if the input images have a signal-to-noise ratio of order 100 or more. The important advantage of this technique over transmitter-based holography is that it allows measurements at arbitrary elevation angles, so allowing one to characterise the large scale deformations in an antenna as a function of elevation.

We describe phase-retrieval holography measurements of the 100-m diameter Green Bank Telescope using astronomical sources and an astronomical receiver operating at a wavelength of 7 mm. We use the technique with parameterization of the aperture in terms of Zernike polynomials and employing a large defocus, as described by Nikolic, Hills & Richer (2006). Individual measurements take around 25 minutes and from the resulting beam maps (which have peak signal to noise ratios of 200:1) we show that it is possible to produce low-resolution maps of the wavefront errors with accuracy around a hundredth of a wavelength. Using such measurements over a wide range of elevations, we have calculated a model for the wavefront-errors due to the uncompensated gravitational deformation of the telescope. This model produces a significant improvement at low elevations, where these errors are expected to be the largest; after applying the model, the aperture efficiency is largely independent of elevation. We have also demonstrated that the technique can be used to measure and largely correct for thermal deformations of the antenna, which often exceed the uncompensated gravitational deformations during daytime observing. We conclude that the aberrations induced by gravity and thermal effects are large-scale and the technique used here is particularly suitable for measuring such deformations in large millimetre wave radio telescopes.

We present Giant Metrewave Radio Telescope 21-cm HI observations towards a sample of compact radio sources behind galaxy groups, to search for cool HI. The results -- from high dynamic range spectra for 8 lines-of-sight through 7 galaxy groups -- do not show any evidence for absorption by cool HI. At a resolution of 20 km/s, the optical depth upper limits obtained were between 0.0075 and 0.035 (3sigma); these correspond to upper limits of a few times 10**23 m**-2 for the column density of any cool HI along these lines of sight (assuming a spin temperature of 100 K).

Using the third data release of the Sloan Digital Sky Survey (SDSS) we have rigorously defined a volume limited sample of early-type galaxies in the redshift range z < 0.1. We have defined the density of the local environment for each galaxy using a method which takes account of the redshift bias introduced by survey boundaries if traditional methods are used. At luminosities greater than our absolute r-band magnitude cutoff of -20.45 the mean density of environment shows no trend with redshift. We calculate the Lick indices for the entire sample and correct for aperture effects and velocity dispersion in a model independent way. Although we find no dependence of redshift or luminosity with environment we do find that the mean velocity dispersion, sigma, of early-type galaxies in dense environments tends to be higher than in low density environments. Taking account of this effect we find that several indices show small but very significant trends with environment that are not the result of the correlation between indices and velocity dispersion. The statistical significance of the data is sufficiently high to reveal that models accounting only for alpha-enhancement struggle to produce a consistent picture of age and metallicity of the sample galaxies, whereas a model that also includes carbon enhancement fares much better. We find that early-type galaxies in the field are younger than those in environments typical of clusters but that neither metallicity, alpha-enhancement nor carbon enhancement are influenced by the environment. The youngest early-type galaxies in both field and cluster environments are those with the lowest sigma. However, there is some evidence that the objects with the largest sigma are slightly younger, especially in denser environments.

We have combined up-to-date stellar population synthesis models, a simple radiative transfer approach, and a fully comprehensive dust model with the aim of developing a simple but quantitative way of interpreting the mid-infrared spectra of galaxies. We apply these models to the observed correlations of mid-infrared luminosities (at 8 and 24 micron) with the star-formation rate of normal galaxies and find that the observations are naturally reproduced by our models. We further find that the observed 24 micron correlation places a weak constraint on relative distribution of dust and stars.

The effect of galaxy interactions on star formation has been investigated using Data Release 1 of the Sloan Digital Sky Survey (SDSS). Both the imaging and spectroscopy data products have been used to construct a catalogue of nearest companions to a volume limited (0.03 < z < 0.1) sample of galaxies drawn from the main galaxy sample of SDSS. Of the 13973 galaxies in the volume limited sample, we have identified 12492 systems with companions at projected separations less than 300 kpc. Star-formation rates for the volume-limited sample have been calculated from extinction and aperture corrected H-alpha luminosities and, where available, IRAS data. Specific star formation rates were calculated by estimating galaxy masses from z-band luminosities, and r-band concentration indices were used as an indicator of morphological class. The mean specific star-formation rate is found to be strongly enhanced for projected separations of less than 30 kpc. For late-type galaxies the correlation extends out to projected separations of 300 kpc and is most pronounced in actively star-forming systems. The specific star-formation rate is observed to decrease with increasing recessional velocity difference, but the magnitude of this effect is small compared to that associated with the projected separation. We also observe a tight relationship between the concentration index and pair separation; the mean concentration index is largest for pairs with separations of approximately 75 kpc and declines rapidly for separations smaller than this. This is interpreted as being due to the presence of tidally-triggered nuclear starbursts in close pairs. Further, we find no dependence of star formation enhancement on the morphological type or mass of the companion galaxy