This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical MHD, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust \it{quantitative} measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to HI observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent MHD simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.

The interstellar medium (ISM) is a magnetised system in which transonic or supersonic turbulence is driven by supernova explosions. This leads to the production of intermittent, filamentary structures in the ISM gas density, whilst the associated dynamo action also produces intermittent magnetic fields. The traditional theory of random functions, restricted to second-order statistical moments (or power spectra), does not adequately describe such systems. We apply topological data analysis (TDA), sensitive to all statistical moments and independent of the assumption of Gaussian statistics, to the gas density fluctuations in a magnetohydrodynamic (MHD) simulation of the multi-phase ISM. This simulation admits dynamo action, so produces physically realistic magnetic fields. The topology of the gas distribution, with and without magnetic fields, is quantified in terms of Betti numbers and persistence diagrams. Like the more standard correlation analysis, TDA shows that the ISM gas density is sensitive to the presence of magnetic fields. However, TDA gives us important additional information that cannot be obtained from correlation functions. In particular, the Betti numbers per correlation cell are shown to be physically informative. Magnetic fields make the ISM more homogeneous, reducing the abundance of both isolated gas clouds and cavities, with a stronger effect on the cavities. Remarkably, the modification of the gas distribution by magnetic fields is captured by the Betti numbers even in regions more than 300 pc from the midplane, where the magnetic field is weaker and correlation analysis fails to detect any signatures of magnetic effects.

Current astrophysical models of the interstellar medium assume that small scale variation and noise can be modelled as Gaussian random fields or simple transformations thereof, such as lognormal. We use topological methods to investigate this assumption for three regions of the southern sky. We consider Gaussian random fields on two-dimensional lattices and investigate the expected distribution of topological structures quantified through Betti numbers. We demonstrate that there are circumstances where differences in topology can identify differences in distributions when conventional marginal or correlation analyses may not. We propose a non-parametric method for comparing two fields based on the counts of topological features and the geometry of the associated persistence diagrams. When we apply the methods to the astrophysical data, we find strong evidence against a Gaussian random field model for each of the three regions of the interstellar medium that we consider. Further, we show that there are topological differences at a local scale between these different regions.

The invariance of the speed of light implies a series of consequences related to our perception of simultaneity and of time itself. Whilst these consequences are experimentally well studied for subluminal speeds, the kinematics of superluminal motion lack direct evidence. Using high temporal resolution imaging techniques, we demonstrate that if a source approaches an observer at superluminal speeds, the temporal ordering of events is inverted and its image appears to propagate backwards. If the source changes its speed, crossing the interface between sub- and super-luminal propagation, we observe image pair annihilation and creation. These results show that it is not possible to unambiguously determine the kinematics of an event from imaging and time-resolved measurements alone.

The Cosmology Data Management system (CosmoDM) is an automated and flexible data management system for the processing and calibration of data from optical photometric surveys. It is designed to run on supercomputers and to minimize disk I/O to enable scaling to very high throughput during periods of reprocessing. It serves as an early prototype for one element of the ground-based processing required by the Euclid mission and will also be employed in the preparation of ground based data needed in the eROSITA X-ray all sky survey mission. CosmoDM consists of two main pipelines. The first is the single-epoch or detrending pipeline, which is used to carry out the photometric and astrometric calibration of raw exposures. The second is the co- addition pipeline, which combines the data from individual exposures into deeper coadd images and science ready catalogs. A novel feature of CosmoDM is that it uses a modified stack of As- tromatic software which can read and write tile compressed images. Since 2011, CosmoDM has been used to process data from the DECam, the CFHT MegaCam and the Pan-STARRS cameras. In this paper we shall describe how processed Pan-STARRS data from CosmoDM has been used to optically confirm and measure photometric redshifts of Planck-based Sunyaev-Zeldovich effect selected cluster candidates.

This report exists to provide high-level guidance for the strategic and engineering development of Data Management and Preservation plans for 'Big Science' data. Although the report's nominal audience is therefore rather narrow, we intend the document to be of use to other planners and data architects who wish to implement good practice in this area. For the purposes of this report, we presume that the reader is broadly persuaded (by external fiat if nothing else) of the need to preserve research data appropriately, and that they have both sophisticated technical support and the budget to support developments. The goal of the document is not to provide mechanically applicable recipes, but to allow the user to develop and lead a high-level plan which is appropriate to their organisation. Throughout, the report is informed where appropriate by the OAIS reference model.

We present the global group properties of two samples of galaxy groups containing 39 high quality X-ray selected systems and 38 optically (spectroscopically) selected systems in coincident spatial regions at 0.12<z<0.79. Only nine optical systems are associable with X-ray systems. We discuss the confusion inherent in the matching of both galaxies to extended X-ray emission and of X-ray emission to already identified optical systems. Extensive spectroscopy has been obtained and the resultant redshift catalog and group membership are provided here. X-ray, dynamical, and total stellar masses of the groups are also derived and presented. We explore the effects of applying three different kinds of radial cut to our systems: a constant cut of 1 Mpc and two r200 cuts, one based on the velocity dispersion of the system and the other on the X-ray emission. We find that an X-ray based r200 results in less scatter in scaling relations and less dynamical complexity as evidenced by results of the Anderson-Darling and Dressler-Schectman tests, indicating that this radius tends to isolate the virialized part of the system. The constant and velocity dispersion based cuts can overestimate membership and can work to inflate velocity dispersion and dynamical and stellar mass. We find Lx-sigma and Mstellar-Lx scaling relations for X-ray and optically selected systems are not dissimilar. The mean fraction of mass found in stars for our systems is approximately 0.014 with a logarithmic standard deviation of 0.398 dex. We also define and investigate a sample of groups which are X-ray underluminous given the total group stellar mass. For these systems the fraction of stellar mass contributed by the most massive galaxy is typically lower than that found for the total population of groups implying that there may be less IGM contributed from the most massive member in these systems. (Abridged)

We present quantitative morphology measurements of a sample of optically selected group galaxies at 0.3 < z < 0.55 using the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and the GIM2D surface brightness--fitting software package. The group sample is derived from the Canadian Network for Observational Cosmology Field Redshift survey (CNOC2) and follow-up Magellan spectroscopy. We compare these measurements to a similarly selected group sample from the Millennium Galaxy Catalogue (MGC) at 0.05 < z < 0.12. We find that, at both epochs, the group and field fractional bulge luminosity (B/T) distributions differ significantly, with the dominant difference being a deficit of disk--dominated (B/T < 0.2) galaxies in the group samples. At fixed luminosity, z=0.4 groups have ~ 5.5 +/- 2 % fewer disk--dominated galaxies than the field, while by z=0.1 this difference has increased to ~ 19 +/- 6 %. Despite the morphological evolution we see no evidence that the group environment is actively perturbing or otherwise affecting the entire existing disk population. At both redshifts, the disks of group galaxies have similar scaling relations and show similar median asymmetries as the disks of field galaxies. We do find evidence that the fraction of highly asymmetric, bulge--dominated galaxies is 6 +/- 3 % higher in groups than in the field, suggesting there may be enhanced merging in group environments. We replicate our group samples at z=0.4 and z=0 using the semi-analytic galaxy catalogues of Bower et al (2006). This model accurately reproduces the B/T distributions of the group and field at z=0.1. However, the model does not reproduce our finding that the deficit of disks in groups has increased significantly since z=0.4.

We have obtained near-infrared imaging of 58 galaxy groups, in the redshift range 0.1<z<0.6, from the William Herschel Telescope and from the Spitzer IRAC data archive. The groups are selected from the CNOC2 redshift survey, with additional spectroscopy from the Baade telescope (Magellan). Our group samples are statistically complete to K(Vega}=17.7 (INGRID) and [3.6](AB)=19.9 (IRAC). From these data we construct near-infrared luminosity functions, for groups in bins of velocity dispersion, up to 800 km/s, and redshift. The total amount of near-infrared luminosity per group is compared with the dynamical mass, estimated from the velocity dispersion, to compute the mass-to-light ratio, M/Lk. We find that the M/Lk values in these groups are in good agreement with those of their statistical descendants at z=0, with no evidence for evolution beyond that expected for a passively evolving population. There is a trend of M/Lk with group mass, which increases from M/Lk~10 for groups with sigma<250 km/s to M/Lk~100 for sigma=425-800 km/s. This trend is weaker, but still present, if we estimate the total mass from weak lensing measurements. In terms of stellar mass, stars make up ~2 per cent of the mass in the smallest groups, and ~1 per cent in the most massive groups. We also use the near-infrared data to consider the correlations between stellar populations and stellar masses, for group and field galaxies at 0.1<z<0.6. We find that fewer group galaxies show strong [OII] emission compared with field galaxies of the same stellar mass and at the same redshift. We conclude that most of the stellar mass in these groups was already in place by z~0.4, with little environment-driven evolution to the present day.