We investigate various phenomenological schemes for the rapid generation of 3D mock galaxy catalogues with a given power spectrum and bispectrum. We apply the fast bispectrum estimator \MODALLSS{} to these mock galaxy catalogues and compare to $N$-body simulation data analysed with the halo-finder \texttt{ROCKSTAR} (our benchmark data). We propose an assembly bias model for populating parent halos with subhalos by using a joint lognormal-Gaussian probability distribution for the subhalo occupation number and the halo concentration. This prescription enabled us to recover the benchmark power spectrum from $N$-body simulations to within 1\% and the bispectrum to within 4\% across the entire range of scales of the simulation. A small further boost adding an extra galaxy to all parent halos above the mass threshold $M>2\times10^{14}\,h^{-1} M_\odot$ obtained a better than 1\% fit to both power spectrum and bispectrum in the range $K/3<1.1\,h\,\text{Mpc}^{-1}$, where $K=k_1+k_2+k_3$. This statistical model should be applicable to fast dark matter codes, allowing rapid generation of mock catalogues which simultaneously reproduce the halo power spectrum and bispectrum obtained from $N$-body simulations. We also investigate alternative schemes using the Halo Occupation Distribution (HOD) which depend only on halo mass, but these yield results deficient in both the power spectrum (2\%) and the bispectrum (>4\%) at $k,K/3 \approx 0.2\,h\,\text{Mpc}^{-1}$, with poor scaling for the latter. Efforts to match the power spectrum by modifying the standard four-parameter HOD model result in overboosting the bispectrum (with a 10\% excess). We also characterise the effect of changing the halo profile on the power spectrum and bispectrum.

Cosmological information from forthcoming galaxy surveys, such as LSST and Euclid, will soon exceed that available from the CMB. Higher order correlation functions, like the bispectrum, will be indispensable for realising this potential. The interpretation of this data faces many challenges because gravitational collapse of matter is a complex non-linear process, typically modelled by computationally expensive N-body simulations. Proposed alternatives using fast dark matter codes (e.g. 2LPT or particle-mesh) are primarily evaluated on their ability to reproduce clustering statistics linked to the matter power spectrum. The accuracy of these codes can be tested in more detail by looking at higher-order statistics, and in this paper we will present an efficient and optimal methodology (MODAL-LSS) to reconstruct the full bispectrum of any 3D density field. We make quantitative comparisons between a number of fast dark matter codes and Gadget at redshift $z=0.5$. This will serve as an important diagnostic tool for dark matter/halo mock catalogues and lays the foundation for realistic high precision analysis with the galaxy bispectrum. In particular, we show that the lack of small-scale power in the bispectrum of fast codes can be ameliorated by a simple `boosting' technique for the power spectrum. We also investigate the covariance of the MODAL-LSS bispectrum estimator, demonstrating the plateauing of non-Gaussian errors in contrast to simple Gaussian extrapolations. This has important consequences for the extraction of information from the bispectrum and hence parameter estimation. Finally we make quantitative comparisons of simulation bispectra with theoretical models, discussing the initial parameters required to create mock catalogues with accurate bispectra.

We present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy. Since the images of gravitationally lensed quasars are the superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and $gi$ multicolour photometric observations from the Dark Energy Survey (DES), near IR $JK$ photometry from the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have identified a candidate system with two catalogue components with $i_{AB}=18.61$ and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of $z=2.739\pm0.003$ and a foreground early type galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass $M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days. The relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.