The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$-$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more strongly correlated with galaxy structure than either stellar populations or star formation rate. These results lead to an approximate lower bound on the dispersion in halo masses among $ \log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are very tightly coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.

The measurement of an annual modulation in the event rate of direct dark matter detection experiments is a powerful tool for dark matter discovery. Indeed, several experiments have already claimed such a discovery in the past decade. While most of them have later revoked their conclusions, and others have found potentially contradictory results, one still stands today. This paper explains the potential as well as the challenges of annual modulation measurements, and gives an overview of past, present and future direct detection experiments.

Muons and other ionizing radiation produced by cosmic rays and radiative decays affect CMOS/CCD sensor. When particles colliding with sensors atoms cause specific kind of noise on images recorded by cameras. We present a concept and preliminary implementation of method for recognizing those events and algorithms for image processing and their classification by machine learning. Our method consists of analyzing the shape of traces present in images recorded by a camera sensor and metadata related to an image like camera model, GPS location of camera, vertical and horizontal orientation of a camera sensor, timestamp of image acquisition, and other events recognized near-by sensors. The so created feature vectors are classified as either a muon-like event, an electron-like event or the other event, possibly noise. For muon-like events our method estimates azimuth of a muon track. Source of the data is database of CREDO (Cosmic-Ray Extremely Distributed Observatory) project and ESO (European Southern Observatory) archives. The telescope dark frames from ESO are analysed. CREDO project collected so far over 2 millions images of events from many kinds of cameralike: smartphones camera, laptop webcams and Internet of Things cameras localised around the globe.

The search for Cosmic-Ray Ensembles (CRE), groups of correlated cosmic rays that might be distributed over very large areas, even of the size of the planet, requires a globally spread and dense network of detectors, as proposed by the Cosmic-Ray Extremely Distributed Observatory (CREDO) Collaboration. This proposal motivates an effort towards exploring the potential of using even very much diversified detection technologies within one system, with detection units located even in hard-to-reach places, where, nevertheless, the sensors could work independently - without human intervention. For these reasons we have developed a dedicated communication solution enabling the connection of many different types of detectors, in a range of environments. The proposed data transmission system uses radio waves as an information carrier on the 169MHz frequency band in a contrast to the typical commercially used frequencies in a IoT systems (868MHz). The connectivity within the system is based on the star topology, which ensures the least energy consumption. The solution is now being prepared to being implemented using the prototype detection system based on the CosmicWatch open hardware design: a portable, pocket size, and economy particle detector using the scintillation technique. Our prototype detector is equipped with a dedicated software that integrates it with the already operational CREDO server system.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) uses the hunt for particle cascades from deep space as a vehicle for a unique "bottom-up" approach to scientific research. By engaging the non-specialist public of all ages as "citizen scientists" we create opportunities for lifelong learning for individuals as well as for cooperation and the sharing of common educational tools amongst institutions. The discoveries of these citizen scientists will feed directly into a pioneering new area of scientific research oriented on Cosmic Ray Ensembles (CRE). The detection (or non-detection) of such particle groups promises to open up a new method for exploring our universe, and a new channel on the multi-messenger stage, oriented on both astro- and geo-investigations. The opportunities this would create for cross-disciplinary research are significant and beneficial for individuals, networks of institutions and the global communities of both professional scientists and science enthusiasts.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may spread over a significant fraction of the Earth, is a scientific terra incognita. The key idea of CREDO is to combine existing cosmic-ray detectors (large professional arrays, educational instruments, individual detectors, such as smartphones, etc.) into a worldwide network, thus enabling a global analysis. The second goal of CREDO involves a large number of participants (citizen science!), assuring the geographical spread of the detectors and managing manpower necessary to deal with vast amount of data to search for evidence for cosmic-ray ensembles. In this paper the status and perspectives of the project are presented.

Propagation of ultra-high energy photons in the galactic and intergalactic space gives rise to cascades comprising thousands of photons. Using Monte Carlo simulations, we investigate the development of such cascades in the solar magnetosphere, and find that the photons in the cascades are distributed over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We also perform similar study for cascades starting as far as 10 Mpc away from us using relevant magnetic field models. A few photons correlated in time are expected to arrive at the Earth from the latter type of cascade. We present our simulation results and discuss the prospects for detection of these cascades with the Cosmic-Ray Extremely Distributed Observatory.

Propagation of ultra-high energy photons in the solar magnetosphere gives rise to cascades comprising thousands of photons. We study the cascade development using Monte Carlo simulations and find that the photons in the cascades are spatially extended over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We compare results from simulations which use two models of the solar magnetic field, and show that although signatures of such cascades are different for the models used, for practical detection purpose in the ground-based detectors, they are similar.

We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte Carlo (MCMC) methods by restricting plausible regions and exploring parameter space efficiently. However, PRISM can additionally be used as a standalone alternative to MCMC for model analysis, providing insight into the behavior of complex scientific models. With PRISM, the time spent on evaluating a model is minimized, providing developers with an advanced model analysis for a fraction of the time required by more traditional methods. This paper provides an overview of the different techniques and algorithms that are used within PRISM. We demonstrate the advantage of using the Bayes linear approach over a full Bayesian analysis when analyzing complex models. Our results show how much information can be captured by PRISM and how one can combine it with MCMC methods to significantly speed up calibration processes (>15 times faster). PRISM is an open-source Python package that is available under the BSD 3-Clause License (BSD-3) at https://github.com/1313e/PRISM and hosted at https://prism-tool.readthedocs.io. PRISM has also been reviewed by "The Journal of Open Source Software" (https://doi.org/10.21105/joss.01229).

One of the main objectives of cosmic-ray studies are precise measurements of energy and chemical composition of particles with extreme energies. Large and sophisticated detectors are used to find events seen as showers starting in the Earth's atmosphere with recorded energies larger than 100 EeV. However, a Cosmic-Ray Ensemble (CRE) developing before reaching the Earth as a bunch of correlated particles may spread over larger areas and requires an extended set of detectors to be discovered. The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a solution to find such phenomena. Even simple detectors measuring the particle arrival time only are useful in this approach, as they are sufficient both to provide candidate CRE events and to determine the direction from which they are arriving.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may overspread a significant fraction of the Earth, is a scientific terra incognita. Instead of developing and commissioning a completely new global detector infrastructure, CREDO proposes approaching the global cosmic-ray analysis objectives with all types of available detectors, from professional to pocket size, merged into a worldwide network. With such a network it is possible to search for evidences of correlated cosmic-ray ensembles. One of the observables that can be investigated in CREDO is a number of spatially isolated events collected in a small time window which could shed light on fundamental physics issues. The CREDO mission and strategy requires active engagement of a large number of participants, also non-experts, who will contribute to the project by using common electronic devices (e.g. smartphones). In this note the status and perspectives of the project is presented.

We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to modulate gas cooling and star formation efficiencies, the SAM can reproduce the cosmic evolution of galaxy properties predicted by the hydrodynamic simulation. These include the galaxy stellar mass function, total baryonic mass, star-forming gas mass and star formation rate at $z\sim5-11$. However, this agreement is only possible by reducing the star formation threshold relative to that suggested by local observations. Otherwise, too much star-forming gas is trapped in quenched dwarf galaxies. We further find that dwarf galaxies rapidly build up their star-forming reservoirs in the early universe ($z>10$), with the relevant time-scale becoming significantly longer towards lower redshifts. This indicates efficient accretion in cold mode in these low-mass objects at high redshift. Note that the improved SAM, which has been calibrated against hydrodynamic simulations, can provide more accurate predictions of high-redshift dwarf galaxy properties that are essential for reionization study.

We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.

We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and increases again for higher-mass haloes. However, the galaxy gas accretion does not flatten at intermediate halo masses when AGN feedback is switched off. To better understand these trends, we develop a physically motivated semi-analytic model of galaxy gas accretion. We show that the flattening is produced by the rate of gas cooling from the hot halo. The ratio of the cooling radius and the virial radius does not decrease continuously with increasing halo mass as generally thought. While it decreases up to ${\sim}10^{13}M_{\odot}$ haloes, it increases for higher halo masses, causing an upturn in the galaxy gas accretion rate. This may indicate that in high-mass haloes AGN feedback is not sufficiently efficient. When there is no AGN feedback, the density of the hot halo is higher, the ratio of the cooling and virial radii does not decrease as much and the cooling rate is higher. Changes in the efficiency of stellar feedback can also increase or decrease the accretion rates onto galaxies. The trends can plausibly be explained by the re-accretion of gas ejected by progenitor galaxies and by the suppression of black hole growth, and hence AGN feedback, by stellar feedback.

Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $z\sim0.6-7$, the black hole mass function at $z\lesssim0.5$, the global ionizing emissivity at $z\sim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${\sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $z\sim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $z\sim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $z\sim6$. We conclude that quasars are unlikely to make a significant contribution to reionization.

We investigate the clustering properties of Lyman-break galaxies (LBGs) at $z\sim6$ - $8$. Using the semi-analytical model {\scshape Meraxes} constructed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation (DRAGONS) project, we predict the angular correlation function (ACF) of LBGs at $z\sim6$ - $8$. Overall, we find that the predicted ACFs are in good agreement with recent measurements at $z\sim 6$ and $z\sim 7.2$ from observations consisting of the Hubble eXtreme Deep Field (XDF), the Hubble Ultra-Deep Field (HUDF) and Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. We confirm the dependence of clustering on luminosity, with more massive dark matter haloes hosting brighter galaxies, remains valid at high redshift. The predicted galaxy bias at fixed luminosity is found to increase with redshift, in agreement with observations. We find that LBGs of magnitude $M_{{\rm AB(1600)}} < -19.4$ at $6\lesssim z \lesssim 8$ reside in dark matter haloes of mean mass $\sim 10^{11.0}$- $10^{11.5} M_{\rm \odot}$, and this dark matter halo mass does not evolve significantly during reionisation.

We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters.

Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.

Direct detection of regions of ionized hydrogen (HII) has been suggested as a promising probe of cosmic reionization. Observing the redshifted 21-cm signal of hydrogen from the epoch of reionization (EoR) is a key scientific driver behind new-generation, low-frequency radio interferometers. We investigate the feasibility of combining low-frequency observations with the Square Kilometre Array and near infra-red survey data of the Wide-Field Infrared Survey Telescope to detect cosmic reionization by imaging HII bubbles surrounding massive galaxies during the cosmic dawn. While individual bubbles will be too small to be detected, we find that by stacking redshifted 21-cm spectra centred on known galaxies, it will be possible to directly detect the EoR at $z \sim 9-12$, and to place qualitative constraints on the evolution of the spin temperature of the intergalactic medium (IGM) at $z \geq 9$. In particular, given a detection of ionized bubbles using this technique, it is possible to determine if the IGM surrounding them is typically in absorption or emission. Determining the globally-averaged neutral fraction of the IGM using this method will prove more difficult due to degeneracy with the average size of HII regions.

Using a series of high-resolution hydrodynamical simulations we show that during the rapid growth of high-redshift (z > 5) galaxies, reserves of molecular gas are consumed over a time-scale of 300Myr, almost independent of feedback scheme. We find that there exists no such simple relation for the total gas fractions of these galaxies, with little correlation between gas fractions and specific star formation rates. The bottleneck or limiting factor in the growth of early galaxies is in converting infalling gas to cold star-forming gas. Thus, we find that the majority of high redshift dwarf galaxies are effectively in recession, with demand (of star formation) never rising to meet supply (of gas), irrespective of the baryonic feedback physics modelled. We conclude that the basic assumption of self-regulation in galaxies - that they can adjust total gas consumption within a Hubble time - does not apply for the dwarf galaxies thought to be responsible for providing most UV photons to reionize the high redshift Universe. We demonstrate how this rapid molecular time-scale improves agreement between semi-analytic model predictions of the early Universe and observed stellar mass functions.

The SKA and its pathfinders will enable studies of HI emission at higher redshifts than ever before. In moving beyond the local Universe, this will require the use of cosmologically appropriate formulae that have traditionally been simplified to their low-redshift approximations. In this paper, we summarise some of the most important relations for tracing HI emission in the SKA era, and present an online calculator to assist in the planning and analysis of observations (hifi.icrar.org).

We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${\sim}6$ and ${\sim}11$), halos that host dwarf galaxies ($\lesssim 10^{9} \mathrm{M_\odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${\sim}90\%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${\sim}30\%$ to 60\%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations.

Using Hubble data, including new grism spectra, Oesch et al. recently identified GN-z11, an $M_\textrm{UV}$=-21.1 galaxy at $z$=11.1 (just 400Myr after the big bang). With an estimated stellar mass of $\sim$10$^9$M$_{\odot}$, this galaxy is surprisingly bright and massive, raising questions as to how such an extreme object could form so early in the Universe. Using \Meraxes{}, a semi-analytic galaxy-formation model developed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) programme, we investigate the potential formation mechanisms and eventual fate of GN-z11. The volume of our simulation is comparable to that of the discovery observations and possesses two analogue galaxies of similar luminosity to this remarkably bright system. Existing in the two most massive subhaloes at $z$=11.1 ($M_\textrm{vir}$=1.4$\times 10^{11}$M$_{\odot}$ and 6.7$\times 10^{10}$M$_{\odot}$), our model analogues show excellent agreement with all available observationally derived properties of GN-z11. Although they are relatively rare outliers from the full galaxy population at high-$z$, they are no longer the most massive or brightest systems by $z$=5. Furthermore, we find that both objects possess relatively smooth, but extremely rapid mass growth histories with consistently high star formation rates and UV luminosities at $z{>}11$, indicating that their brightness is not a transient, merger-driven feature. Our model results suggest that future wide-field surveys with the \textit{James Webb Space Telescope} may be able to detect the progenitors of GN-z11 analogues out to $z{\sim}$13--14, pushing the frontiers of galaxy-formation observations to the early phases of cosmic reionization and providing a valuable glimpse of the first galaxies to reionize the Universe on large scales.

We investigate high-redshift galaxy sizes using a semi-analytic model constructed for the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation project. Our fiducial model, including strong feedback from supernovae and photoionization background, accurately reproduces the evolution of the stellar mass function and UV luminosity function. Using this model, we study the size--luminosity relation of galaxies and find that the effective radius scales with UV luminosity as $R_\mathrm{e}\propto L^{0.25}$ at $z{\sim}5$--$9$. We show that recently discovered very luminous galaxies at $z{\sim}7$ (Bowler et al. 2016) and $z{\sim}11$ (Oesch et al. 2016) lie on our predicted size--luminosity relations. We find that a significant fraction of galaxies at $z>8$ will not be resolved by JWST, but GMT will have the ability to resolve all galaxies in haloes above the atomic cooling limit. We show that our fiducial model successfully reproduces the redshift evolution of average galaxy sizes at $z>5$. We also explore galaxy sizes in models without supernova feedback. The no-supernova feedback models produce galaxy sizes that are smaller than observations. We therefore confirm that supernova feedback plays an important role in determining the size--luminosity relation of galaxies and its redshift evolution during reionization.

We introduce Meraxes, a new, purpose-built semi-analytic galaxy formation model designed for studying galaxy growth during reionization. Meraxes is the first model of its type to include a temporally and spatially coupled treatment of reionization and is built upon a custom (100 Mpc)$^3$ N-body simulation with high temporal and mass resolution, allowing us to resolve the galaxy and star formation physics relevant to early galaxy formation. Our fiducial model with supernova feedback reproduces the observed optical depth to electron scattering and evolution of the galaxy stellar mass function between $z$=5 and 7, predicting that a broad range of halo masses contribute to reionization. Using a constant escape fraction and global recombination rate, our model is unable to simultaneously match the observed ionizing emissivity at $z{\lesssim}6$. However, the use of an evolving escape fraction of 0.05-0.1 at $z{\sim}6$, increasing towards higher redshift, is able to satisfy these three constraints. We also demonstrate that photoionization suppression of low mass galaxy formation during reionization has only a small effect on the ionization history of the inter-galactic medium. This lack of `self-regulation' arises due to the already efficient quenching of star formation by supernova feedback. It is only in models with gas supply-limited star formation that reionization feedback is effective at regulating galaxy growth. We similarly find that reionization has only a small effect on the stellar mass function, with no observationally detectable imprint at $M_{\rm *}{>}10^{7.5}\,{\rm M_{{\odot}}}$. However, patchy reionization has significant effects on individual galaxy masses, with variations of factors of 2-3 at $z$=5 that correlate with environment.

We use the Dark-ages, Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) framework to investigate the effect of galaxy-formation physics on the morphology and statistics of ionized hydrogen (HII) regions during the Epoch of Reioinization (EoR). DRAGONS self-consistently couples a semi-analytic galaxy-formation model with the inhomogeneous ionizing UV background, and can therefore be used to study the dependence of morphology and statistics of reionization on feedback phenomena of the ionizing source galaxy population. Changes in galaxy-formation physics modify the sizes of HII regions and the amplitude and shape of 21-cm power spectra. Of the galaxy physics investigated, we find that supernova feedback plays the most important role in reionization, with HII regions up to $\approx 20$ per cent smaller and a fractional difference in the amplitude of power spectra of up to $\approx 17$ per cent at fixed ionized fraction in the absence of this feedback. We compare our galaxy-formation-based reionization models with past calculations that assume constant stellar-to-halo mass ratios and find that with the correct choice of minimum halo mass, such models can mimic the predicted reionization morphology. Reionization morphology at fixed neutral fraction is therefore not uniquely determined by the details of galaxy formation, but is sensitive to the mass of the haloes hosting the bulk of the ionizing sources. Simple EoR parametrizations are therefore accurate predictors of reionization statistics. However, a complete understanding of reionization using future 21-cm observations will require interpretation with realistic galaxy-formation models, in combination with other observations.

In this paper we present calculations of the UV luminosity function from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) project, which combines N-body, semi-analytic and semi-numerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from $z{\sim}5$ through to $z{\sim}10$. We investigate the luminosity--star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of $0.1$--$0.3$ dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the luminosity function below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at $M_\mathrm{UV}{\gtrsim}{-14}$. We find that $48$ ($17$) per cent of the total UV flux at $z{\sim}6$ ($10$) has been detected above an observational limit of $M_\mathrm{UV}{\sim}{-17}$, and that galaxies fainter than $M_\mathrm{UV}{\sim}{-17}$ are the main source of ionizing photons for reionization. We investigate the luminosity--stellar mass relation, and find a correlation for galaxies with $M_\mathrm{UV}{<}{-14}$ that has the form $M_*{\propto}10^{-0.47M_\mathrm{UV}}$, in good agreement with observations, but which flattens for fainter galaxies. We determine the luminosity--halo mass relation to be $M_\mathrm{vir}{\propto}10^{-0.35M_\mathrm{UV}}$, finding that galaxies with $M_\mathrm{UV}{=}{-20}$ reside in host dark matter haloes of $10^{11.0\pm 0.1}\mathrm{M_\odot}$ at $z{\sim}6$, and that this mass decreases towards high redshift.

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $\Delta c_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $\alpha$, depends on peak height, $\nu$, in a manner that is accurately described by $\alpha {=}0.0070\nu^2{+}0.1839$. The distribution of the spin parameter, $\lambda$, has a weak dependence on equilibrium state; $\lambda$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulations. Because the semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. The resulting concentration-mass (c-M) relations are found to agree with the simulations, and because the model applies only to relaxed halos, they do not exhibit the upturn at high masses or high redshifts found by some recent works. We predict a change of slope in the z=0 c-M relation at a mass scale of $10^{11}\rm{M}_{\odot}$. We find that this is due to the change in the functional form of the halo MAH, which goes from being dominated by an exponential (for high-mass halos) to a power-law (for low-mass halos). During the latter phase, the core radius remains approximately constant, and the concentration grows due to the drop of the background density. We also analyse how the c-M relation predicted by this work affects the power produced by dark matter annihilation, finding that at z = 0 the power is two orders of magnitude lower than that obtained from extrapolating best-fitting c-M relations. We provide fits to the c-M relations as well as numerical routines to compute concentrations and MAHs.

We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^{\alpha}e^{\beta z}$, where the parameters $\alpha$ and $\beta$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $\alpha$, $\beta$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.

Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime. However, in the low-redshift regime the mass history follows a power law because the growth of density perturbations is halted in the dark energy dominated era due to the accelerated expansion of the Universe. We provide an analytic model that follows the expression $M(z) = M_{0}(1+z)^{af(M_{0})}e^{-f(M_{0})z}$, where $M_{0} = M(z=0)$, $a$ depends on cosmology and $f(M_{0})$ depends only on the linear matter power spectrum. The analytic model does not rely on calibration against numerical simulations and is suitable for any cosmology. We compare our model with the latest empirical models for the mass accretion history in the literature and find very good agreement. We provide numerical routines for the model online.

We use high-resolution simulations of cosmological volumes to model galaxy formation at high-redshift, with the goal of studying the photon budget for reionization. We demonstrate that galaxy formation models that include a strong, thermally coupled supernovae scheme reproduce current observations of star formation rates and specific star formation rates, both during and after the reionization era. These models produce enough UV photons to sustain reionization at z<8 (z<6) through a significant population of faint, unobserved, galaxies for an assumed escape fraction of 20% (5%). This predicted population is consistent with extrapolation of the faint end of observed UV luminosity functions. We find that heating from a global UV/X-ray background after reionization causes a dip in the total global star formation rate density in galaxies below the current observational threshold. Finally, while the currently observed specific star formation rates are incapable of differentiating between supernovae feedback models, sufficiently deep observations will be able to use this diagnostic in the future to investigate galaxy formation at high redshift.

Galaxy redshift surveys using optical telescopes have, in combination with other cosmological probes, enabled precision measurements of the nature of dark energy. We show that radio telescopes are rapidly becoming competitive with optical facilities in spectroscopic surveys of large numbers of galaxies. Two breakthroughs are driving this change. Firstly, individual radio telescopes are more efficient at mapping the sky thanks to the large field-of-view of new phased-array feeds. Secondly, ever more dishes can be correlated in a cost-effective manner with rapid increases in computing power. The next decade will see the coming of age of the 21cm radio wavelength as a cosmological probe as first the Pathfinders then, ultimately, the Square Kilometre Array is constructed. The latter will determine precise 3D positions for a billion galaxies, mapping the distribution of matter in the Universe over the last 12 billion years. This radio telescope will be able to constrain the equation of state of dark energy, and its potential evolution, to a precision rivalling that of future optical facilities such as DESI and Euclid.

We present new wideband radio observations with the Australia Telescope Compact Array of a sample of 12 giant radio galaxies. The radio observations are part of a larger radio-optical study aimed at relating the radio structures with the ambient medium on large scales. With projected linear sizes larger than 0.7 Mpc, these objects are ideal candidates for the study of the Warm-Hot Intergalactic Medium (WHIM). The sample includes sources with sizes spanning 0.8 to 3.2 Mpc and total powers of 1.2*10^24 to 4.0*10^26 W Hz^-1 at 2.1 GHz. Redshifts were limited to z<0.15 to permit spectroscopic observations of the hosts and neighbouring galaxies, which were obtained using the AAOmega spectrograph on the Anglo-Australian Telescope. We derive lobe energy densities from the radio observations via equipartition arguments. The inferred pressures in the lobes of the giant radio sources, which range from 1.1*10^-15 to 2.0*10^-14 Pa (80 to 1500 cm^-3 K), are lower than previously inferred from X-ray observations of dense filaments. Comparison with the OverWhelmingly Large Simulations (OWLS) suggests that the WHIM in pressure balance with the radio lobes has a temperature in excess of ~10^6.5 K or a particle overdensity in the range 50 to 500. This study highlights the capability of next generation surveys, such as the Evolutionary Map of the Universe (EMU) survey with the Australian Square Kilometre Array Pathfinder (ASKAP), to study populations of giant radio sources at lower surface brightness and thereby discriminate between models for the cosmological evolution of the intergalactic medium and examine the validity of cosmological hydrodynamical simulations.

We use numerical simulations to investigate how the statistical properties of dark matter (DM) haloes are affected by the baryonic processes associated with galaxy formation. We focus on how these processes influence the spin and shape of a large number of DM haloes covering a wide range of mass scales, from galaxies to clusters at redshifts zero and one, extending to dwarf galaxies at redshift two. The haloes are extracted from the OverWhelmingly Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations run with a range of feedback prescriptions. We find that the median spin parameter in DM-only simulations is independent of mass, redshift and cosmology. At z=0 baryons increase the spin of the DM in the central region (<=0.25 r_200) by up to 30 per cent when feedback is weak or absent. This increase can be attributed to the transfer of angular momentum from baryons to the DM, but is no longer present at z=2. We also present fits to the mass dependence of the DM halo shape at both low and high redshift. At z=0 the sphericity (triaxiality) is negatively (positively) correlated with halo mass and both results are independent of cosmology. Interestingly, these mass-dependent trends are markedly weaker at z=2. While the cooling of baryons acts to make the overall DM halo more spherical, stronger feedback prescriptions tend to reduce the impact of baryons by reducing the central halo mass concentration. More generally, we demonstrate a strongly positive (negative) correlation between halo sphericity (triaxiality) and galaxy formation efficiency, with the latter measured using the central halo baryon fraction. In conclusion, our results suggest that the effects of baryons on the DM halo spin and shape are minor when the effects of cooling are mitigated, as required by realistic models of galaxy formation, although they remain significant for the inner halo.

The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) which will probe the mass and dynamics of over 600,000 galaxies. A much deeper small-area HI survey, called DINGO, aims to trace the evolution of HI from z = 0 - 0.43, a cosmological volume of 40 million Mpc^3, detecting potentially 100,000 galaxies. The high-sensitivity 30 antenna ASKAP core (diameter ~2 km) will provide an angular resolution of 30 arcsec (at z=0). Our simulations show that the majority of galaxies detected in WALLABY (87.5%) will be resolved. About 5000 galaxies will be well resolved, i.e. more than five beams (2.5 arcmin) across the major axis, enabling kinematic studies of their gaseous disks. This number would rise to 160,000 galaxies if all 36 ASKAP antennas could be used; the additional six antennas provide baselines up to 6 km, resulting in an angular resolution of 10 arcsec. For DINGO this increased resolution is highly desirable to minimise source confusion; reducing confusion rates from a maximum of 10% of sources at the survey edge to 3%. We estimate that the sources detected by WALLABY and DINGO will span four orders of magnitude in total halo mass (from 10^{11} to 10^{15} Msol) and nearly seven orders of magnitude in stellar mass (from 10^{5} to 10^{12} Msol), allowing us to investigate the process of galaxy formation across the last four billion years.

This is a design study into the capabilities of the Australian Square Kilometre Array Pathfinder in performing a full-sky low redshift neutral hydrogen survey, termed WALLABY, and the potential cosmological constraints one can attain from measurement of the galaxy power spectrum. We find that the full sky survey will likely attain 0.6 million redshifts which, when combined with expected Planck CMB data, will constrain the Dark Energy equation of state to 20%, representing a coming of age for radio observations in creating cosmological constraints.

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.

We discuss the possibility of performing a substantial spectroscopic galaxy redshift survey selected via the 21cm emission from neutral hydrogen using the Five-hundred metre Aperture Spherical Telescope (FAST) to be built in China. We consider issues related to the estimation of the source counts and optimizations of the survey, and discuss the constraints on cosmological models that such a survey could provide. We find that a survey taking around two years could detect ~10^7 galaxies with an average redshift of ~0.15 making the survey complementary to those already carried out at optical wavelengths. These conservative estimates have used the z=0 HI mass function and have ignored the possibility of evolution. The results could be used to constrain Gamma = (Omega_m h) to 5 per cent and the spectral index, n_s, to 7 per cent independent of cosmic microwave background data. If we also use simulated power spectra from the Planck satellite, we can constrain w to be within 5 per cent of -1.

We present new photometric observations in the BVRI bands of the double-lined eclipsing binary BD+05 706 conducted over three observing seasons, as well as new X-ray observations obtained with the ROSAT satellite covering a full orbital cycle (P = 18.9 days). A detailed light-curve analysis of the optical data shows the system to be semidetached, confirming indications from an earlier analysis by Torres et al. (1998), with the less massive and cooler star filling its Roche lobe. The system is a member of the rare class of cool Algol systems, which are different from the "classical" Algol systems in that the mass-gaining component is also a late-type star rather than a B- or A-type star. By combining the new photometry with a reanalysis of the spectroscopic observations reported by Torres et al. (1998) we derive accurate absolute masses for the components of M1 = 2.633 +/- 0.028 Msun and M2 = 0.5412 +/- 0.0093 Msun, radii of R1 = 7.55 +/- 0.20 Rsun and R2 = 11.02 +/- 0.21 Rsun, as well as effective temperatures of 5000 +/- 100 K and 4640 +/- 150 K for the primary and secondary, respectively. There are obvious signs of activity (spottedness) in the optical light curve of the binary. Our X-ray light curve clearly shows the primary eclipse but not the secondary eclipse, suggesting that the primary star is the dominant source of the activity in the system. The depth and duration of the eclipse allow us to infer some of the properties of the X-ray emitting region around that star.

We conducted an extensive search for Wolf-Rayet stars (W-Rs) in the SMC, using the same interference filter imaging techniques that have proved successful in finding W-Rs in more distant members of the Local Group. Photometry of some 1.6 million stellar images resulted in some 20 good candidates, which we then examined spectroscopically. Two of these indeed proved to be newly found W-Rs, bringing the total known in the SMC from 9 to 11. Other finds included previously unknown Of-type stars (one as early as O5f?p)),the recovery of the Luminous Blue Variable S18, and the discovery of a previously unknown SMC symbiotic star. More important, however, is the fact that there does not exist a significant number of W-Rs waiting to be discovered in the SMC. The number of W-Rs in the SMC is a factor of 3 lower than in the LMC (per unit luminosity), and we argue this is the result of the SMC's low metallicity on the evolution of the most massive stars.