We derive robust constraints on primordial non-Gaussianity (PNG) using the
clustering of 800,000 photometric quasars from the Sloan Digital Sky Survey in
the redshift range 0.5<z<3.5. These measurements rely on the novel technique
of {\it extended mode projection} to control the impact of spatially-varying
systematics in a robust fashion, making use of blind analysis techniques. This
allows the accurate measurement of quasar halo bias at the largest scales,
while discarding as little as possible of the data. The standard local-type PNG
parameters fNL and gNL both imprint a k−2
scale-dependent effect in the bias. Constraining these individually, we obtain
−49<fNL<31 and −2.7×105<gNL<1.9×105, while
their joint constraints lead to −105<fNL<72 and
−4.0×105<gNL<4.9×105 (all at 95% CL) . Introducing a
running parameter nfNL to constrain b(k)∝k−2+nfNL and a generalised PNG amplitude
f~NL, we obtain −45.5exp(3.7nfNL)<f~NL<34.4exp(3.3nfNL) at 95% CL. These
results incorporate uncertainties in the cosmological parameters, redshift
distributions, shot noise, and the bias prescription used to relate the quasar
clustering to the underlying dark matter. These are the strongest constraints
obtained to date on PNG using a single population of large-scale structure
tracers, and are already at the level of pre-{\it Planck} constraints from the
cosmic microwave background. A conservative forecast for a {\it Large Synoptic
Survey Telescope}-like survey incorporating mode projection yields
σ(fNL)∼5 -- competitive with the {\it Planck} result --
highlighting the power of upcoming large scale structure surveys to probe the
initial conditions of the universe.Comment: 5 pages, 5 figures. v3: version accepted in PRL. v4: systematics
templates publicly available on www.earlyuniverse.org/code, no change to
pape