Multi-component modelling of galaxies is a valuable tool in the effort to quantitatively understand galaxy evolution, yet the use of the technique is plagued by issues of convergence, model selection and parameter degeneracies. These issues limit its application over large samples to the simplest models, with complex models being applied only to very small samples. We attempt to resolve this dilemma of "quantity or quality" by developing a novel framework, built inside the Zooniverse citizen science platform, to enable the crowdsourcing of model creation for Sloan Digitial Sky Survey galaxies. We have applied the method, including a final algorithmic optimisation step, on a test sample of 198 galaxies, and examine the robustness of this new method. We also compare it to automated fitting pipelines, demonstrating that it is possible to consistently recover accurate models that either show good agreement with, or improve on, prior work. We conclude that citizen science is a promising technique for modelling images of complex galaxies, and release our catalogue of models.

Dust attenuation in star-forming spiral galaxies affects stars and gas in different ways due to local variations in dust geometry. We present spatially resolved measurements of dust attenuation for a sample of 232 such star-forming spiral galaxies, derived from spectra acquired by the SDSS-IV MaNGA survey. The dust attenuation affecting the stellar populations of these galaxies (obtained using full spectrum stellar population fitting methods) is compared with the dust attenuation in the gas (derived from the Balmer decrement). Both of these attenuation measures increase for local regions of galaxies with higher star formation rates; the dust attenuation affecting the stellar populations increases more so than the dust attenuation in the gas, causing the ratio of the dust attenuation affecting the stellar populations to the dust attenuation in the gas to decrease for local regions of galaxies with higher star formation rate densities. No systematic difference is discernible in any of these dust attenuation quantities between the spiral arm and inter-arm regions of the galaxies. While both the dust attenuation in the gas and the dust attenuation affecting the stellar populations decrease with galactocentric radius, the ratio of the two quantities does not vary with radius. This ratio does, however, decrease systematically as the stellar mass of the galaxy increases. Analysis of the radial profiles of the two dust attenuation measures suggests that there is a disproportionately high concentration of birth clouds (incorporating gas, young stars and clumpy dust) nearer to the centres of star-forming spiral galaxies.

Bars in galaxies are thought to stimulate both inflow of material and radial mixing along them. Observational evidence for this mixing has been inconclusive so far however, limiting the evaluation of the impact of bars on galaxy evolution. We now use results from the MaNGA integral field spectroscopic survey to characterise radial stellar age and metallicity gradients along the bar and outside the bar in 128 strongly barred galaxies. We find that age and metallicity gradients are flatter in the barred regions of almost all barred galaxies when compared to corresponding disk regions at the same radii. Our results re-emphasize the key fact that by azimuthally averaging integral field spectroscopic data one loses important information from non-axisymmetric galaxy components such as bars and spiral arms. We interpret our results as observational evidence that bars are radially mixing material in galaxies of all stellar masses, and for all bar morphologies and evolutionary stages.

The exact nature of the arms of spiral galaxies is still an open question. It has been widely assumed that spiral arms in galaxies with two distinct symmetrical arms are the products of density waves that propagate around the disk, with the spiral arms being visibly enhanced by the star formation that is triggered as the passing wave compresses gas in the galaxy disk. Such a persistent wave would propagate with an approximately constant angular speed, its pattern speed Op. The quasi-stationary density wave theory can be tested by measuring this quantity and showing that it does not vary with radius in the galaxy. Unfortunately, this measurement is difficult because Op is only indirectly connected to observables such as the stellar rotation speed. Here, we use the detailed information on stellar populations of the grand-design spiral galaxy UGC 3825, extracted from spectral mapping, to measure the offset between young stars of a known age and the spiral arm in which they formed, allowing the first direct measure of Op at a range of radii. The offset in this galaxy is found to be as expected for a pattern speed that varies little with radius, indicating consistency with a quasi-stationary density wave, and lending credence to this new method.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 63 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.

We present quantified visual morphologies of approximately 48,000 galaxies observed in three Hubble Space Telescope legacy fields by the Cosmic And Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and classified by participants in the Galaxy Zoo project. 90% of galaxies have z < 3 and are observed in rest-frame optical wavelengths by CANDELS. Each galaxy received an average of 40 independent classifications, which we combine into detailed morphological information on galaxy features such as clumpiness, bar instabilities, spiral structure, and merger and tidal signatures. We apply a consensus-based classifier weighting method that preserves classifier independence while effectively down-weighting significantly outlying classifications. After analysing the effect of varying image depth on reported classifications, we also provide depth-corrected classifications which both preserve the information in the deepest observations and also enable the use of classifications at comparable depths across the full survey. Comparing the Galaxy Zoo classifications to previous classifications of the same galaxies shows very good agreement; for some applications the high number of independent classifications provided by Galaxy Zoo provides an advantage in selecting galaxies with a particular morphological profile, while in others the combination of Galaxy Zoo with other classifications is a more promising approach than using any one method alone. We combine the Galaxy Zoo classifications of "smooth" galaxies with parametric morphologies to select a sample of featureless disks at 1 < z < 3, which may represent a dynamically warmer progenitor population to the settled disk galaxies seen at later epochs.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We identify 885,503 type 1 quasar candidates to i<22 using the combination of optical and mid-IR photometry. Optical photometry is taken from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III/BOSS), while mid-IR photometry comes from a combination of data from the Wide-Field Infrared Survey Explorer (WISE) "ALLWISE" data release and several large-area Spitzer Space Telescope fields. Selection is based on a Bayesian kernel density algorithm with a training sample of 157,701 spectroscopically-confirmed type-1 quasars with both optical and mid-IR data. Of the quasar candidates, 733,713 lack spectroscopic confirmation (and 305,623 are objects that we have not previously classified as photometric quasar candidates). These candidates include 7874 objects targeted as high probability potential quasars with 3.5<z<5 (of which 6779 are new photometric candidates). Our algorithm is more complete to z>3.5 than the traditional mid-IR selection "wedges" and to 2.2<z<3.5 quasars than the SDSS-III/BOSS project. Number counts and luminosity function analysis suggests that the resulting catalog is relatively complete to known quasars and is identifying new high-z quasars at z>3. This catalog paves the way for luminosity-dependent clustering investigations of large numbers of faint, high-redshift quasars and for further machine learning quasar selection using Spitzer and WISE data combined with other large-area optical imaging surveys.

We explore the extinction/reddening of ~35,000 uniformly selected quasars with 0<z<5.3 in order to better understand their intrinsic optical/ultraviolet spectral energy distributions. Using rest-frame optical-UV photometry taken from the Sloan Digital Sky Survey's (SDSS) 7th data release, cross-matched to WISE in the mid-infrared, 2MASS and UKIDSS in the near-infrared, and GALEX in the UV, we isolate outliers in the color distribution and find them well described by an SMC-like reddening law. A hierarchical Bayesian model with a Markov Chain Monte Carlo sampling method was used to find distributions of powerlaw indices and E(B-V) consistent with both the broad absorption line (BAL) and non-BAL samples. We find that, of the ugriz color-selected type 1 quasars in SDSS, 2.5% (13%) of the non-BAL (BAL) sample are consistent with E(B-V)>0.1 and 0.1% (1.3%) with E(B-V)>0.2. Simulations show both populations of quasars are intrinsically bluer than the mean composite, with a mean spectral index (${\alpha}_{\lambda}$) of -1.79 (-1.83). The emission and absorption-line properties of both samples reveal that quasars with intrinsically red continua have narrower Balmer lines and stronger ionizing spectral lines, the latter indicating a harder continuum in the extreme-UV and the former pointing to differences in black hole mass and/or orientation.

We explore the mid-infrared (mid-IR) through ultraviolet (UV) spectral energy distributions (SEDs) of 119,652 luminous broad-lined quasars with 0.064<z<5.46 using mid-IR data from Spitzer and WISE, near-infrared data from Two Micron All Sky Survey and UKIDSS, optical data from Sloan Digital Sky Survey, and UV data from Galaxy Evolution Explorer. The mean SED requires a bolometric correction (relative to 2500A) of BC=2.75+-0.40 using the integrated light from 1um-2keV, and we further explore the range of bolometric corrections exhibited by individual objects. In addition, we investigate the dependence of the mean SED on various parameters, particularly the UV luminosity for quasars with 0.5<z<3 and the properties of the UV emission lines for quasars with z>1.6; the latter is a possible indicator of the strength of the accretion disk wind, which is expected to be SED dependent. Luminosity-dependent mean SEDs show that, relative to the high-luminosity SED, low-luminosity SEDs exhibit a harder (bluer) far-UV spectral slope, a redder optical continuum, and less hot dust. Mean SEDs constructed instead as a function of UV emission line properties reveal changes that are consistent with known Principal Component Analysis (PCA) trends. A potentially important contribution to the bolometric correction is the unseen extream-UV (EUV) continuum. Our work suggests that lower-luminosity quasars and/or quasars with disk-dominated broad emission lines may require an extra continuum component in the EUV that is not present (or much weaker) in high-luminosity quasars with strong accretion disk winds. As such, we consider four possible models and explore the resulting bolometric corrections. Understanding these various SED-dependent effects will be important for accurate determination of quasar accretion rates.

In the restframe UV, two of the parameters that best characterize the range of emission-line properties in quasar broad emission-line regions are the equivalent width and the blueshift of the CIV line relative to the quasar rest frame. We explore the connection between these emission-line properties and the UV through X-ray spectral energy distribution (SED) for radio-quiet (RQ) quasars. Our sample consists of a heterogeneous compilation of 406 quasars from the Sloan Digital Sky Survey and Palomar-Green survey that have well-measured CIV emission-line and X-ray properties (including 164 objects with measured Gamma). We find that RQ quasars with both strong CIV emission and small CIV blueshifts can be classified as "hard-spectrum" sources that are (relatively) strong in the X-ray as compared to the UV. On the other hand, RQ quasars with both weak CIV emission and large CIV blueshifts are instead "soft-spectrum" sources that are (relatively) weak in the X-ray as compared to the UV. This work helps to further bridge optical/soft X-ray "Eigenvector 1" relationships to the UV and hard X-ray. Based on these findings, we argue that future work should consider systematic errors in bolometric corrections (and thus accretion rates) that are derived from a single mean SED. Detailed analysis of the CIV emission line may allow for SED-dependent corrections to these quantities.

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

The properties of Galactic molecular clouds tabulated by Solomon etal (1987) (SRBY) are re-examined using the Boston University-FCRAO Galactic Ring Survey of 13CO J=1-0 emission. These new data provide a lower opacity tracer of molecular clouds and improved angular and spectral resolution than previous surveys of molecular line emission along the Galactic Plane. We calculate GMC masses within the SRBY cloud boundaries assuming LTE conditions throughout the cloud and a constant H2 to 13CO abundance, while accounting for the variation of the 12C/13C with Galacto-centric radius. The LTE derived masses are typically five times smaller than the SRBY virial masses. The corresponding median mass surface density of molecular hydrogen for this sample is 42 Msun/pc^2, which is significantly lower than the value derived by SRBY (median 206 Msun/pc^2) that has been widely adopted by most models of cloud evolution and star formation. This discrepancy arises from both the extrapolation by SRBY of velocity dispersion, size, and CO luminosity to the 1K antenna temperature isophote that likely overestimates the GMC masses and our assumption of constant 13CO abundance over the projected area of each cloud. Owing to the uncertainty of molecular abundances in the envelopes of clouds, the mass surface density of giant molecular clouds could be larger than the values derived from our 13CO measurements. From velocity dispersions derived from the 13CO data, we find that the coefficient of the cloud structure functions, vo=sigma_v/R^{1/2}, is not constant, as required to satisfy Larson's scaling relationships, but rather systematically varies with the surface density of the cloud as Sigma^{0.5} as expected for clouds in self-gravitational equlibrium.