{'fraction_non_alphanumeric': 0.029441624365482234, 'fraction_numerical': 0.013127139652933538, 'mean_word_length': 4.609323268441266, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Models of gravitational lens candidates from Space Warps CFHTLS

We report modelling follow-up of recently-discovered gravitational-lens candidates in the Canada France Hawaii Telescope Legacy Survey. Lens modelling was done by a small group of specially-interested volunteers from the SpaceWarps citizen-science community who originally found the candidate lenses. Models are categorised according to seven diagnostics indicating (a) the image morphology and how clear or indistinct it is, (b) whether the mass map and synthetic lensed image appear to be plausible, and (c) how the lens-model mass compares with the stellar mass and the abundance-matched halo mass. The lensing masses range from ~10^11 Msun to>10^13 Msun. Preliminary estimates of the stellar masses show a smaller spread in stellar mass (except for two lenses): a factor of a few below or above ~10^11 Msun. Therefore, we expect the stellar-to-total mass fraction to decline sharply as lensing mass increases. The most massive system with a convincing model is J1434+522 (SW05). The two low-mass outliers are J0206-095 (SW19) and J2217+015 (SW42); if these two are indeed lenses, they probe an interesting regime of very low star-formation efficiency. Some improvements to the modelling software (SpaghettiLens), and discussion of strategies regarding scaling to future surveys with more and frequent discoveries, are included.


INTRODUCTION
By a curious coincidence, the typical escape velocity of massive galaxies -of order a few hundred km/s -is such that v 2 esc /c 2 , expressed as an angular distance, is comparable to the apparent sizes of galaxies at cosmological distances. This coincidence is fortunate, because it makes the gravitational lensing deflection angle of distant galaxies (α ∼ 2v 2 esc /c 2 ) comparable to their size on the sky. As a result, strong lensing by galaxies produces images that probe their host dark matter halos, providing a useful tool to understand galaxy formation. While there is a general consensus about the basic mechanisms at play, involving gravitational collapse, fragmentation, and mergers of dark-matter clumps, into which gas fell, cooling through radiative processes to form dense clouds and eventually stars, there is much debate about the details (for a summary, see Silk & Mamon 2012). In particular, the nature of dark matter remains mysterious: most researchers take it to be a collisionless non-relativistic fluid (cold dark matter or CDM) readily studied by simulations (for example, the influential Millennium simulation, Springel et al. 2005). However, alternative scenarios, where dark matter has exotic dynamical properties (Saxton & Ferreras 2010;Schive et al. 2016), or is not really matter at all, but a modification of gravity (McGaugh et al. 2016), have also been considered.
All this motivates the use of strong gravitational lensing over galaxy scales to study the mutual dynamics of dark matter, gas and stars. Several studies in recent years have done so (see, e.g., Koopmans et al. 2009;Leier et al. 2011Leier et al. , 2012Leier et al. , 2016Bruderer et al. 2016) but it is desirable to enlarge the samples from tens of lensing galaxies to thousands. Doing so requires both finding more lenses and also modelling their mass distribution. Recent searches through the CFHT Lens Survey (CFHTLS, Heymans et al. 2012) using arc-finders (e.g., More et al. 2012;Maturi et al. 2014;Gavazzi et al. 2014;Sonnenfeld et al. 2017) by either machine learning methods (e.g., Paraficz et al. 2016;Lanusse et al. 2017) or visual inspection by citizen-science volunteers through the Space Warps project (More et al. 2016) have, between them, discovered an average of four lenses per square degree, so one can be optimistic about finding many thousands of lenses in the next generation of wide-field surveys, from ground-based surveys such as LSST in the optical window and SKA in radio, to space-based missions such as Euclid and WFIRST (Oguri & Marshall 2010;Collett 2015). The availability of the recently published year 1 result of the Dark Energy Survey (Drlica-Wagner et al. 2017), which covered 1800 deg 2 , suggests that the number of lenses could grow by at least an order of magnitude within a very short time period.
The expected flood of new lens discoveries will need a similarly large modelling effort to reconstruct their mass distributions. Lens-modelling robots have started to be developed (Nightingale et al. 2017;Hezaveh et al. 2017) but so far are able to handle only very clean systems. For typical lenses, human input is still needed. With that in mind Küng et al. (2015) developed a new modelling strategy, implemented as the SpaghettiLens 1 system. The idea is to collaborate with experienced members of the citizen-science community, who have already participated in lens discovery through Space Warps, as well as several other projects involving astronomical data. The system was tested on a sample of simulated lenses, which were part of the training and testing set in Space Warps.
This paper follows up that study by applying Spaghet-tiLens to candidates discovered through Space Warps. We present results from the modelling of 58 of the 59 lens candidates reported by More et al. (2016). Each lens candidate was modelled following a collaborative refinement process, where anyone interested could improve the analysis by modifying an existing model or creating a new one. Note the difference with respect to the main Space Warps project, where volunteers from a group of 10 4 people make independent contributions. Each person is presented with a random selection of survey-patches and invited to (in effect) vote on each. The system estimates the skill level of each volunteer according to test-patches interspersed with the real data, and weights their votes accordingly (Marshall et al. 2016). There is an active forum for volunteers, but since everyone is seeing different data samples with minimal overlap, the forum has little if any influence on the votes. In SpaghettiLens, the number of volunteers is significantly lower, but the level of interaction is higher. The resulting model represents a consensus among contributors, as to the best that could be achieved with the available data and software.
We emphasize that the interpretation of the results presented here is tentative, because the systems are lens candidates at this stage, not secure lenses. Moreover the 1 http://labs.spacewarps.org/spaghetti/ Figure 1. Number of models generated over time in days since the launch of SpaghettiLens; the black line shows all models, the coloured lines show the five most active candidates. A total of 377 models were generated for 58 lens candidates, eight users contributed more than five models per person.
candidate-lens redshifts have large uncertainties, while the candidate-source redshifts can only be guessed at present. Nevertheless it is interesting to explore the trends observed with the already available data, based on some 300 models created by 8 volunteers for 58 candidate lenses found in the 150 deg 2 CFHTLS.
This paper is organized as follows: Section 2 introduces the candidate lenses, their models and the diagnostics applied to them. The following sections elaborate on the diagnostics. Image morphology diagnostics are explained in Section 3 and diagnostics based on the mass models are discussed in Section 4. Stellar masses are presented in Section 5, to compare stellar and lensing masses. Finally, section 6 summarises and tabulates the diagnostics in Table 1.
We include three appendices devoted to various technical issues relating to the modelling. The online supplement gives the results from all the modelled systems generated for all the lensing candidates.

THE CANDIDATES AND MODELS
Space Warps is a citizen science gravitational lens search (Marshall et al. 2016). Its first run searched the CFHT Legacy Survey, a survey carried out in five optical bands (u * ,g ,r ,i ,z ) covering a total area of 160 deg 2 divided up into four fields W1 -W4. The MegaPrime camera, used for the survey, has a field of view of 1 deg 2 , with 0.186 pixels. The flux limit of the survey is g < 25.6 AB (5σ) with a typical seeing FWHM ∼ 0.7 (Erben et al. 2013). The cutouts shown to the volunteers were colour composites of the g , r and i bands from randomly chosen regions. The programme resulted in 59 new lens candidates, of which 29 are considered high quality. Moreover, 82 previously known lens candidates were "rediscovered". The candidates have broadly similar redshifts (0.2 < z < 1.0), Einstein radii between 0.7 and 5 , and arc fluxes in g-band AB magnitudes between 22 and 26. These properties are similar to those found by previous robotic searches like ArcFinder and RingFinder (More et al. 2016).
A first version of SpaghettiLens was launched in May 2013, while the first run of Space Warps was ongoing. Fig. 1 shows the total number of models generated since launch for our lens candidates. A small number of volunteers started creating models of possible lens candidates from Space Warps that were debated in the discussion forum. In the beginning of August 2013 (around day 70 since release) a modelling challenge for the candidate later to be identified as SW 01 was launched. A group of five volunteers presented their results in an online letter converging on a set of 30 models. In April 2015 (day 700), the first major version of SpaghettiLens was released. At the same time, the list of identified candidates was made available (as a preprint of More et al. 2016). At this stage, the volunteers were asked to create models for all the lens candidates that were still missing. The work done for candidate SW 29 show that the volunteers converged on a favourable model much faster, after 15 models, generated in 30 days. The total modelling effort resulted in a set of 377 SpaghettiLens models for all but one of the 59 Space Warps lens candidates; candidate SW 39 was not modelled by any volunteer. Over its runtime, the system logged eight major contributors (defined as volunteers that generated, at least, five models).
In parallel, Küng et al. (2015) tested the system on simulated lenses and identified some areas for improvement. In § A1 we introduce fitting of the surface brightness profile of the source. This feature has not been included yet in Spaghet-tiLens, but we carried the analysis during post-processing for a few interesting cases. In § A2 we show that fine-grained mass maps within the central region relieves a tendency in the earlier work for mass profiles to be too shallow. In § A3 we consider the possibility of fitting a parametric lens model to the ensemble. We characterise each model with seven diagnostics, grouped into three categories, whose purpose is to help identify the most likely cases of a gravitational lens, as well as flag the most interesting candidates for future follow-up observations. The diagnostics are as follows.
• Firstly we have diagnostics based on morphology. Section 3 and Fig. 2 explain this diagnostic in more detail. We consider: (i) Whether the images are unblended: Distinct, unblended images are an advantage in modelling, although not essential.
(ii) Whether all images are discernible. The topography of an arrival-time surface, as encoded by a spaghetti diagram, may require more images than are visible, in which case the modeller has to insert conjectural image positions. (iii) Whether the lens is isolated. (iv) The image morphology concisely described: double or quads, further sub-categorised to indicate the elongation direction of the lensing mass.
• Secondly we have mass models, covered in Section 4. We assess the following points: (v) Whether the mass map is reasonable (see Fig. 5). (vi) Whether the arrival-time surface and synthetic image are plausible. In particular, an unsatisfactory model is flagged when additional images are implied in regions where they are not observed (see Figs. 3, 4 6).
• Thirdly, whether the total implied lensing mass determined from the lens model is plausible, given the photometry of the lensing galaxy. In this work we define the lensing mass as the sum of all mass tiles in the model. These mass tiles reach out to typically twice the radius of the outermost image. Section 5 explains how we compare the lensing mass with the stellar mass. Galaxies may have halos extending well beyond the mass maps in the models, so the calculated lensing mass is only a lower bound on the total galaxy mass. Hence, the lensing mass should be somewhere between the stellar mass (lower bound) and the total halo mass (upper bound). Extrapolation of the lens model to estimate the full virial mass may be possible (cf. Leier et al. 2012), but is not attempted in this work.
(vii) We then introduce a so-called halo index (H), that quantifies how the lensing mass compares with these two constraints (see Fig. 7).

IMAGE MORPHOLOGY
The main input of a modeller to the process is a markup of the candidate lens system, which we call a spaghetti diagram. This is a sketch of the arrival-time surface from a point-like source, with proposed locations of maxima, minima and saddle points, and an implied time-ordering of the images. Such information encodes a starting proposal for the mass distribution. A spaghetti diagram is thus a completely abstract construction, and moreover it refers to a simplified system with a point source. However, spaghetti diagrams are intuitive because they tend to resemble the form of lensed arcs (see Fig.3 of Ferreras et al. 2008), and of course they are simple to draw, and easy to vary and refine in an open collaborative environment. This makes them very practical for non-professional lens enthusiasts in the citizen-science community. Details and tests are given in Küng et al. (2015).
We now discuss the diagnostics that can be taken from the process of drawing the spaghetti diagram, even before detailed mass-modelling takes place. Fig. 2 shows nine examples, each consisting of a cutout of the Space Warps image of a lens candidate, marked up with a spaghetti diagram. The volunteers are initially presented a square image with side 82 (440 pixels) in full size, but they have the ability to zoom in individually. The cutouts presented in Figs. 2 to 5 are rescaled during post-processing, relative to the outermost image.
All the examples shown in Fig. 2 identify five locations: the centre of the main lensing galaxy, which is also a maximum of the arrival time surface (red dots); two minima (blue dots); and two saddle points (green dots and also selfintersections of the curves). Only two (SW 05 and SW 42) of the nine systems, however, have five distinct features in plausible locations. In the other seven examples, an arc is interpreted as a blend of three or four images. This defines the 'unblended images' diagnostic. Note that this characterisation could be different if the spaghetti diagram were different. For example, SW 28 has also been modelled with the arc on the right interpreted as a single image, rather than as three images as shown in the figure; for such a model, the 'unblended images' diagnostic would be true.
The second diagnostic checks whether all images are visible. For example, we see in SW 58 at the top left of Fig. 2 that an image at the second minimum is conjectural and does not correspond to any visible feature.
The third diagnostic tests whether the lens is isolated, or whether other galaxies could contribute to the lensing, based on visual inspection of the field. For this, we do not consider stars or other clearly foregound objects. Additional galaxies contributing to the lensing mass can be marked by the volunteer alongside the spaghetti diagrams. An example can be seen as the grey dot and circle in the cutout with SW 57 at the top right of Fig. 2. Objects marked in this way are modelled as point masses. Other possible contributors to lensing are galaxies or groups that are not in the immediate vicinity of the lensed images, yet massive enough to exert an influence. These are accounted for by allowing a constant but adjustable external shear.
We remark that the lack of expected lensed images, or the presence of blended images or perturbing galaxies do not imply that a given candidate is unlikely to be a lens. This means, rather, that the models are more uncertain and perhaps could be more easily improved by trying further variations in the markup.
The fourth diagnostic is based on the fact that the arrangement of lensed images of a pointlike source through a non-circular gravitational lens depends on the location of the source relative to the long and short axes of the lens. This dependence is quite robust and independent of many other details of the lensing mass distribution. Sources close to being dead-centre behind a lens tend to produce quads; sources at larger transverse distance tend to produce doubles. We denote these as Q and D, respectively, and add a prefix to the Q systems, as follows: we write LQ if the source is inferred as displaced along the long axis of the lens, SQ if displaced along the short axis, IQ if inclined to both axes, CQ if only very little or no displacement is evident. Although the unlensed source and its location are obviously not seen, the LQ, SQ, IQ and CQ cases correspond to easily-seen image morphologies (see, e.g., Saha & Williams 2003).
• The simplest is the LQ case: this creates a saddle point and two minima in an arc, with the second saddle point on the other side of the lensing galaxy, closer to the galaxy than the arc. SW 58 and SW 28 in the upper row of Fig. 2 are typical examples of LQ. If the source were moved outwards along the short axis, the minimum-saddle-minimum set would merge into a single minimum, leaving a D system; the transition is known as a cusp catastrophe.
• The middle row of Fig. 2 shows three IQ systems: SW 05, SW 42 and SW 19. This type has a characteristic asymmetry, often with two images close together. If the source were moved outwards, the minimum-saddle pair would merge and mutually cancel, leaving again a D system. This transition is known as a fold catastrophe.
• The lower row of Fig. 2 shows three SQ systems: SW 09, SW 29 and SW 02. The failed model, SW 57, at top right may also belong to this category. Here the images have a fairly symmetric arrangement with an arc and a counter-image on the other side, but the spaghetti diagram is completely asymmetric. If the source were moved outwards along the long axis, the saddle-minimum-saddle set would merge into a single saddle -another form of cusp catastrophe. SQ can be visually distinguished from LQ by the relative distances of the arc and the counter-image. For SQ the arc is closer, for LQ the counter-image is closer.
CQ systems are often called 'cross' or 'Einstein-cross' systems; IQ systems are sometimes called 'folds', with 'cusp' commonly used for both SQ and LQ. The labels 'short-axis quad' and so on are not standard in the literature, but the morphological classification they express is familiar to experienced modellers. Hence they can be useful to researchers wishing to apply other modelling methods to the same systems.

MASS MODELS
Once a spaghetti diagram has been drawn on a web browser, it is forwarded to a server-side numerical framework, which searches for mass maps consistent with the image locations, parities and time ordering, given by the modeller. The mass maps are made up of mass tiles and are free-form, but are required to be concentrated around the identified lens centre (see Coles et al. 2014, for the precise formulation of the search problem). The modeller can also set various options for the search, such as the number of mass tiles and the extent of the mass map; all options have defaults. Typically, there are ∼ 500 mass tiles arranged in a disc, centred on the lensing galaxy and extending to twice the radius of the outermost image. By default, the mass distributed is required to have 180 • rotation symmetry, but this option can be unset. Assuming mass distributions can be found -which in practice is usually the case -a statistical ensemble of 200 two-dimensional mass maps is returned. From each mass map, further quantities such as lens potentials or enclosed masses can be derived. Thus, there are ensembles of 200 values for the surface density at any point, for the enclosed mass within a given radius and so on. In this paper we present the averages and ranges of different quantities, but other quantities such as 90% confidence ranges could also be computed. The whole ensemble of mass maps, along with derived quantities and uncertainties, makes up one SpaghettiLens model.
The mass will naturally depend on the lens and source redshifts, which are unknown when a lens candidate is first identified. However, this is not a problem, because a model can be trivially rescaled to use better redshift values, as they become available. The mass normalisation of the models is proportional, through the angular critical density, to the factor dL dS d −1 LS , where dL, dS and dLS are the usual angulardiameter distances. In the redshift ranges typical of galaxy lensing, the normalisation factor is roughly proportional to the lens redshift, and weakly dependent on the source redshift. This work applies photometric redshifts from the CFHTLS pipeline (Coupon et al. 2009) to the candidate lensing galaxies; the values range from z = 0.2 to 1 (see Table 1). The source redshifts are assumed as z = 2, unless an unambiguous photometric redshift is available. The lens redshifts entail rather big uncertainties, up to a few tens of percent (see Fig. 5 in Coupon et al. 2009). The lensing masses would be uncertain at the same level. On the other hand, the lensing masses in the sample range from 10 11 M to 10 13 M ; in comparison, the redshift uncertainty is not so important in this preliminary analysis.
The ensemble of mass maps can be post-processed in many different ways. Four different graphical quantities are particularly useful. Fig. 3 shows the arrival-time surfaces corresponding to the spaghetti diagrams of Fig. 2. The arrival-time contours look like machine-made elaborations of the input spaghetti diagrams. If the saddle-point contours in the arrival time are qualitatively the same as the curves in the spaghetti diagram (the detailed shape of the spaghetti curves is unimportant), it immediately suggests a successful model. On the other hand, if the arrival-time surface has unexpected minima or saddle points, and especially if the unexpected features are far from identified lens images, that signals an improbable model. Fig. 4 shows what we call synthetic images, meaning reconstructions of the extended lensed features by fitting for a source. These were generated by a new method, explained in Appendix A1, implemented during the offline post-processing after the modelling process was complete. The synthetic images provided by SpaghettiLens during the collaborative modelling and discussion were more crude; those are included in the online supplement.
The arrival-time surface and synthetic image are summarised by one diagnostic, the most important of all: are the lensed features satisfactorily reproduced? This diagnostic remains a judgment call by modellers. A useful quantitative criterion for whether the synthetic image is consistent with the data would need to allow for PSF dependence and unmodelled substructure -otherwise all models would be summarily rejected, something left for a future implementation in SpaghettiLens. Fig. 5 shows the projected mass maps of the sample. In fact, this figure only shows the ensemble-average mass maps, and not the variation within the ensemble, from which uncertainties can be inferred. The same applies to the arrivaltime surfaces and synthetic images in the previous two figures. The uncertainties will be shown in a concise form in Fig. 6. Fig. 5 makes evident the tiled nature of the mass model. The tiles can be smoothed over by interpolation, and this was done in the mass maps during the modelling process, available in the online version. It is interesting, however, to note the tiling artifacts, if only as a reminder that the substructure in the mass distribution is very uncertain, even if some integrated quantities are well constrained. How the free-form mass maps relate to parameterised lens models is discussed in Appendix A3. Note that although the mass distribution can have discontinuous jumps, the lens equation and arrival-time surface are continuous. Fig. 6 shows the enclosed-mass profiles, expressed as the average convergence, κ, within circles of given projected radius. Uncertainties are included (see the figure caption for details). Appendix A2 describes the improvements made since our earlier work (Küng et al. 2015), to allow for steeper profiles in the inner regions. The enclosed mass is typically best constrained at the notional Einstein radius, becoming more uncertain at larger and smaller radii.
The mass maps and mass profiles are the basis of a further diagnostic: are the mass distributions plausible? This is also a judgment call made by modellers, but it showed to be a powerful diagnostic, summarizing three aspects. The overall shape is forced to be 180 • -rotation symmetric, usually a plausible assumption, but volunteers can deactivate this constraint. The profile slope turned out to be a good indicator of plausible models, as can seen by contrasting the model for SW 57 with the rest of the sample: The missing core in Fig. 5 and the flat profile in Fig. 6 disqualify this model. The clumpiness of the mass map is another useful indicator. Flat profile slopes can often be identified directly in the mass map, where the mass tiles form a checkerboard pattern.
More experienced volunteers applied these diagnostics already during the process of creating models as a criteria of a successful model, the evaluation presented in this work however was generated by the authors during post processing.

STELLAR AND HALO MASS ESTIMATES
The stellar masses of the lens galaxies are derived by comparison of the photometric data with stellar M/L estimates from population synthesis models. In principle, a detailed analysis of the spectral energy distribution is needed to derive accurate stellar masses (e.g. Gallazzi & Bell 2009;Taylor et al. 2011). However, at the redshifts probed, the photometric data mainly constrain the information-rich 4000 Å break region, whose strength depends sensitively on age and metallicity, thereby providing a strong constraint on the stellar M/L. Hence, estimates to within ∼0.3 dex in log(M stel /M ) can be derived with a single colour, preferably tracing a rest-frame colour similar to U − V (see Fig. 1 of Ferreras et al. 2008), assuming a universal initial mass function (IMF). There is evidence from detailed absorption line strength analysis that massive galaxies can feature a non-standard IMF (e.g. Ferreras et al. 2013). However, these variations -towards a bottom-heavy distribution -are typically found in the cores of massive early-type galaxies (La ). The effect of these variations on the stellar M/L of lensing systems is still rather controversial (Smith et al. 2015;Leier et al. 2016).
In this paper we further simplify the analysis by assuming a relationship between the apparent total magnitude and stellar mass, at the redshift of the lens. For typical stellarpopulation parameters, the variation of this relation is at most ∆ logMs ∼1 dex. A further potential systematic can arise from contamination of the light of the lensing galaxy by the lensed background source. Reducing or eliminating the latter would require detailed fitting of light distributions for each candidate (see Leier et al. 2011), which we have not yet attempted. Nonetheless, because the lensing masses range over two orders of magnitude, it is still interesting to compare them with rough estimates of stellar mass.
We make use of the Bruzual & Charlot (2003) models to derive two functional forms of the stellar mass with respect to the i -band magnitudes. The models have solar metallicity, with a Chabrier IMF, and assume two different age trends: a "young" model, with a constant 500 Myr age at all redshifts, and an "old" model where the age is the oldest one possible at each redshift, adopting a standard ΛCDM model with H0 = 70 km s −1 Mpc −1 and Ωm = 0.3. Fig. 7 shows a comparison of stellar and lensing mass in our sample. The comparatively large span of the error bars in stellar mass (horizontal axis) shows the range between the masses derived using the two age trends respectively, and lies between 0.4 and 0.8 dex. It will be improved in future work by the use of available optical and NIR magnitudes to derive more accurate constraints on the stellar populations. In addition, we also derive halo masses for the lenses using an abundance-matching formula. This technique matches the distribution function of observed stellar mass in galaxies with that of dark-halo masses from N -body simulations, to define a simple relation between stellar mass and halo mass. We emphasize that a halo mass from abundance matching should be considered an "average" estimate, and a significant scatter can be expected as galaxies with the same stellar mass can be found in different environments. We refer the reader to Leier et al. (2012) for an assessment of the effect of abundance matching on the derivation of dark matter halo properties in lensing galaxies. We follow the prescription of Moster et al. (2010), namely: that relates the observed lensing to stellar mass, with the global ratio expected if the host halo corresponds to the average value derived by abundance matching. Several cases for H can be considered: • H < 0 is unphysical because M < M stel .
• H = 0 means the stellar mass exactly accounts for the lensing mass (i.e. no dark matter affects the lensing model).
• 0 < H < 1 is the typical situation, where the lens includes stars and dark matter, but not the full halo.
• H = 1 means that the lens consists of the entire halo.
• H > 1 is in tension with abundance-matching, because the lensing mass exceeds the expected halo mass.
The halo-matching index expresses whether the lensing mass is plausible given the flux received from the candidate lensing galaxy. Fig. 7 and Table 1 show that most of the candidates have stellar and lensing masses typical of massive ellipticals 2 . SW 05 is one of the most massive of all the candidates, corresponding to a galaxy-group mass scale. It is a particularly attractive system for follow-up observations at higher resolution, as it is a large system with clear multiple-image features. Modelling leaves little doubt that it is a lens. SW 04 seems to be even more massive, but the diagnostics leave some doubts about the validity of this model. The two lowest-mass systems, SW 19 and SW 42, are important if they are indeed lenses, as they would be low-mass lenses dominated by dark matter. All the modelled systems have reasonable stellarmass fractions, except for two cases where the stellar-mass fraction is too low (H > 1): these are SW 42 and SW 57. In the case of SW 57, the model has poor diagnostics and should be discarded. The model for SW 42, on the other hand, is quite convincing -except for the high halo-matching index. If SW 42 turned out not to be a lens, that would support the halo-matching index as an effective criterion to discriminate models. 2 The mass values themselves are given in the online version of Table 1

SUMMARY AND CONCLUSIONS
We report on a first set of mass distributions and followup diagnostics for the Space Warps lens candidates created with a novel approach that aims to be scalable by orders of magnitude to prepare for the many thousands of lenses the next generation of wide field surveys will yield (e.g. Euclid, WFIRST).
Over the past few years, the way of discovering lenses has changed with the introduction of machine learning and citizen science methods, combined with the coverage of large areas of the sky by modern surveys. The way lensing mass models are constructed also needs to change, in order to be prepared for the increasing influx of lens candidates. The work in this paper represents a hybrid approach between the classical style -where a small team of experts invest many hours into the creation of a single model -and a citizen science project -where a crowd of amateur volunteers make independent contributions. The authors of this paper are a collaboration of professionals and experienced citizen-science volunteers, aiming to create early-stage lens models as soon as a lens candidate is found.
To assist volunteers in constraining lensing models, we introduce a set of diagnostics that help asses the validity of the models. At a later stage, we encourage modellers to apply those diagnostics as feedback on the plausibility of their assumptions, and to suggest additional diagnostics.
The diagnostics (i) -(iv) (see Section. 2) turned out to be useful measures of the difficulty in modelling a system, but they did not constitute necessary conditions for a promising model. They can help select systems to introduce novice volunteers to the modelling process. In contrast, diagnostics (v) and (vi) can be considered as necessary criteria for a good model. Volunteers employed those ones to evaluate their models, and turned out to be easy enough to grasp by new volunteers. The halo-matching index diagnostic (vii, H), is an interesting criterion that might be useful for the modellers, but needs further investigation. Table 1 is a summary of our results. It characterises each modelled system with seven diagnostics, indicating (a) the image morphology and how clear or indistinct it is, (b) whether the mass map and synthetic lensed image appear to be plausible, and (c) how the model mass compares with the estimated stellar and full-halo masses. Missing entries are due unavailable imaging data, whereas missing rows are due to models that were not created for this particular candidate.
The trend in Fig. 7, where higher-mass galaxies get progressively more dark-matter dominated, is expected (see, e.g. Ferreras et al. 2005), as is the span of about one order of magnitude for the stellar mass and the two orders of magnitude in lensing mass. With future data, it will be interesting to compare the enclosed stellar and lensing mass as a function of radius, going from the star-dominated inner regions to the outer dark halos. Leier et al. (2011) illustrate this behaviour in their Fig. 5, but the present sample could go an order of magnitude higher in mass.
The quick creation of many models for the Space Warps candidates successfully showed that a subset of citizen scientists are interested in being involved in more challenging tasks that take some time to learn. The next step involves recruiting more lensing enthusiasts, as soon as the next round of Space Warps is started. In the meantime, the improvements shown in the Appendix will be integrated in the standard version of SpaghettiLens. Photometric fitting could also be integrated into SpaghettiLens. This would allow experienced citizen scientists to generate photometric redshifts and stellar masses, and thus generate preliminary dark-matter maps as soon as a lens-candidate is identified. Figure 2. Nine of the lens candidates marked up with spaghetti diagrams. Red, blue and green dots are proposed locations for maxima, minima and saddle points of the arrival time respectively. The curves help guide the placement of the dots, but their precise appearance has no significance. This images are screenshots from the SpaghettiLens user interface, which applies interpolation to background images. The scaling is adjusted to fit the other images. This selection includes the best-modelled systems, but also one case (SW 57 at upper right) of unsuccessful modelling. Since the modelling process is collaborative among the volunteers, with anyone welcome to contribute new models or modify existing ones, there are variant spaghetti diagrams for all the modelled systems. The online supplement displays all the models presented for discussion during this work.     The upper-left shaded region is unphysical according to abundance matching (see Section 5) because it gives M > M halo . That is to say, the unshaded region is 0 < H < 1.

A1 Improved synthetic images
The mass maps produced by the current implementation of SpaghettiLens are based on images of point-like features. No information about extended images is used, except in so far as they help the user identify images of point-like features. The synthetic images offered to users are rudimentary, corresponding to conical light profiles (i.e. circular light profiles with brightness decreasing linearly with radius).
We have now developed a prototype to improve the generation of synthetic images, as illustrated in Fig. A1. Areas containing lensed images are selected (green frames in the figure). The selected areas should be as free as possible from other sources, including the lensing galaxy. Pixels within the selected areas are mapped to a grid on the source plane, using bending angles given by the mass model. The mapping from lens-plane pixels to source-plane grid cells is many-toone, because of image multiplicity and magnification. The brightness of each source-plane pixel is set to the mean of all the lens-plane pixels mapped to it. Finally, the mapping is run back to the lens plane. The result is a synthetic image. In effect, one is reconstructing a source-plane brightness map by least-squares.
The procedure is not yet implemented in SpaghettiLens but can be applied during post-processing. The new synthetic images could be used to improve the mass reconstruction, by weighing the ensemble of maps according to how good the synthetic images are.

A2 Sub-sampling of the central region
The models of simulated lenses in Küng et al. (2015) showed a tendency to produce density profiles which were too shallow, resulting typically in overestimates of the Einstein radius.
Allowing some extra mass tiles in central region, thus allowing the mass profile to rise more steeply near the centre, was suggested as a possible cure. Fig. A2 shows an experiment with smaller mass tiles in the inner region. Replacing the very central mass tile with 9 smaller tiles allows for steeper central profiles. Doing the same for the 25 innermost mass tiles allows for even steeper central profiles, eliminating the systematic shallow profiles. However, it does not provide a completely satisfactory solution, because (a) it increases the number of mass tiles by 40% and significantly increases the computational time, and (b) the square boundary between areas with different tile sizes is rather undesirable. The main modelling work in this paper was, however, done before the experiments with smaller mass tiles was complete. Some of the models presented in this paper apply the intermediate option (corresponding to the middle panel in Fig. A2) while others use the old system. We note this paper mainly concerns the enclosed mass in the outer regions, so shallowness in the central region should not be an issue.

A3 Parametrization of pixel models
In order to fit the set of pixelated models to a single parametrized model, a programme was written that took a parametrized function and subtracted from it the mean and the principal components (PCs) of the data. This created the residual function. The number of principal components in the analysis was varied, to test how this affected the output. It was found that 5 PCs gave a reasonable approximation. A masking function was added, selecting only the data points that fell inside the image of the lens, and the PCs were clipped in order to keep the values inside the region of the ensemble. Any value higher than the clip was set to be the clip value. This was chosen to 2.5 because, assuming that the data follows a Gaussian error distribution, almost all the values for the variance should lie between 2 and 3 standard deviations from the mean. Minimising the residuals function produces the set of parameters that fit the parametrized function to the original pixelated ensemble most closely. A least squares fit was used to perform this minimisation. The parametrized model function was obtained from the gravitational potential of an isothermal ellipsoid mass distribution (Keeton 2001). This model is frequently used to describe gravitational lenses, giving good fits to the observations. The isothermal ellipsoid model outputs three useful parameters: the radius of the Einstein ring, the ellipticity of the model and the angle of the ellipticity from the vertical, giving the orientation of the galaxy. By applying this model to simulated lenses for which the values of these parameters were already known, it was possible to assess the accuracy of the methodology, before applying the model to the candidate lensing galaxies. Preliminary results on the recovery of Einstein radii are shown in Fig. A3. This paper has been typeset from a T E X/L A T E X file prepared by the author. Figure A2. Model improvement resulting from the use of smaller mass tiles in the inner region of the mass model. Shown here are the average enclosed κ within a given projected radius, for three different reconstructions of a simulated lens (sim) from Space Warps. In each panel, the dashed blue curve is the correct answer. The orange band represents the statistical ensemble from SpaghettiLens, the orange line shows the ensemble mean. Locations of images (maximum, saddle point, minimum) are marked with vertical arrows. The radial value at κ = 1 is the effective Einstein radius, r E . The upper panel is taken from Küng et al. (2015), see Fig. 3 of that paper. The middle panel is the result when the innermost mass tile is replaced by 9 smaller tiles. The lower panel results from replacing the innermost 5×5 tiles with 9 smaller tiles each. Figure A3. Comparison of Einstein radii r E obtained from mass tiles to the results from a parameterized model to test the performance of the algorithm. The parameterized model was generated using principal component analysis on the ensemble of models. The blue dashed line represents a perfect recovery of r E .