Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

We calculate the cooling times at constant density for halos with virial temperatures from 100 K to 10^5 K that originate from a 3-sigma fluctuation of a CDM power spectrum in three different cosmologies. Our intention is to determine the first objects that can cool to low temperatures, but not to follow their dynamical evolution. We identify two generations of halos: those with low virial temperatures, Tvir < 9000 K that remain largely neutral, and those with larger virial temperatures that become ionized. The lower-temperature, lower-mass halos are the first to cool to 75 percent of their virial temperature. The precise temperature and mass of the first objects are dependent upon the molecular hydrogen (H2) cooling function and the cosmological model. The higher-mass halos collapse later but, in this paradigm, cool much more efficiently once they have done so, first via electronic transitions and then via molecular cooling: in fact, a greater residual ionization once the halos cool below 9000 K results in an enhanced H2 production and hence a higher cooling rate at low temperatures than for the lower-mass halos, so that within our constant-density model it is the former that are the first to cool to really low temperatures. We discuss the possible significance of this result in the context of CDM models in which the shallow slope of the initial fluctuation spectrum on small scales leads to a wide range of halo masses (of differing overdensities) collapsing over a small redshift interval. This ``crosstalk'' is sufficiently important that both high- and low-mass halos collapse during the lifetimes of the massive stars which may be formed at these epochs. Further investigation is thus required to determine which generation of halos plays the dominant role in early structure formation.

We model a one-dimensional shock-tube using smoothed particle hydrodynamics and investigate the consequences of having finite shock-width in numerical simulations. We investigate the cooling of gas during passage through the shock for different cooling regimes. For a shock temperature of 10^5K, the maximum temperature of the gas is much reduced and the cooling time was reduced by a factor of 2. At lower temperatures, we are especially interested in the production of molecular Hydrogen and so we follow the ionization level and H_2 abundance across the shock. This regime is particularly relevent to simulations of primordial galaxy formation for halos in which the virial temperature of the galaxy is sufficiently high to partially re-ionize the gas. The effect of in-shock cooling is substantial: the maximum temperature the gas reaches compared to the theoretical temperature was found to vary between 0.15 and 0.81 for the simulations performed. The downstream ionization level is reduced from the theoretical level by a factor of between 2.4 and 12.5, and the resulting H_2 abundance was found to be reduced to a fraction of 0.45 to 0.74 of its theoretical value. At temperatures above 10^5K, radiative shocks are unstable and will oscillate. We reproduce these oscillations and find good agreement with the previous work of Chevalier and Imamura (1982), and Imamura, Wolff and Durisen (1984). The effect of in-shock cooling in such shocks is difficult to quantify, but is undoubtedly present.

We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).

We report on work done by the Virgo consortium, an international collaboration set up in order to study the formation and evolution of Large Scale Structure using N-body simulations on the latest generation of parallel supercomputers. We show results of 256^3 particle simulations of the formation of clusters in four Dark Matter models with different cosmological parameters. Normalizing the models such that one obtains the correct abundance of rich clusters yields an interesting result: The peculiar velocities of the clusters are almost independent of \Omega, and depend only weakly on \Gamma, the shape parameter of the power spectrum. Thus, it is nearly impossible to distinguish between high and low \Omega models on the basis of the peculiar velocities.

We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $\Omega_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $\Omega_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.