Fine-structure constant with BOSS

Constraint on the time variation of the fine-structure constant
with the SDSS-III/BOSS DR12 quasar sample

Franco D. Albareti, Johan Comparat, Carlos M. Gutiérrez, Francisco Prada, Isabelle Pâris, David Schlegel, Martín López-Corredoira, Donald P. Schneider, Arturo Manchado, D. A. García-Hernández, Patrick Petitjean and Jian Ge
Instituto de Física Teórica (UAM/CSIC), Universidad Autónoma de Madrid, Cantoblanco, E-28049 Madrid, Spain
Instituto de Astrofísica de Canarias (IAC), La Laguna, E-38205 Tenerife, Spain
Departamento de Astrofísica, Universidad de La Laguna, La Laguna, E-38206 Tenerife, Spain
Campus of International Excellence UAM+CSIC, Cantoblanco, E-28049 Madrid, Spain
Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía, E-18080 Granada, Spain
INAF, Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34131 Trieste, Italy
Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA, 94720, USA
Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA
Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA
CSIC, Spain
Institut d’Astrophysique de Paris, CNRS-UPMC, UMR7095, 98bis bd Arago, 75014 Paris, France
Department of Astronomy, University of Florida, Gainesville, FL 32611-2055, USA
‘la Caixa’-Severo Ochoa Scholar.E-mail:
Accepted for publication in MNRAS.

From the Sloan Digital Sky Survey (SDSS) Data Release 12, which covers the full Baryonic Oscillation Spectroscopic Survey (BOSS) footprint, we investigate the possible variation of the fine-structure constant over cosmological time-scales. We analyse the largest quasar sample considered so far in the literature, which contains 13 175 spectra (10 363 from SDSS-III/BOSS DR12 + 2812 from SDSS-II DR7) with redshift . We apply the emission-line method on the [O iii] doublet () and obtain for the relative variation of the fine-structure constant. We also investigate the possible sources of systematics: misidentification of the lines, sky OH lines, H and broad line contamination, Gaussian and Voigt fitting profiles, optimal wavelength range for the Gaussian fits, chosen polynomial order for the continuum spectrum, signal-to-noise ratio and good quality of the fits. The uncertainty of the measurement is dominated by the sky subtraction. The results presented in this work, being systematics limited, have sufficient statistics to constrain robustly the variation of the fine-structure constant in redshift bins () over the last 7.9 Gyr. In addition, we study the [Ne iii] doublet () present in quasar spectra and discuss the systematic effects on using these emission lines to constrain the fine-structure constant variation. Better constraints on using the emission-line method would be possible with high-resolution spectroscopy and large galaxy/qso surveys.

line: profiles – quasars: emission lines – cosmology: observations – surveys – large-scale structure of Universe.
pagerange: Constraint on the time variation of the fine-structure constant with the SDSS-III/BOSS DR12 quasar sampleReferencespubyear: 2002

1 Introduction

Reference Quasar spectra SDSS release Time ago (Gyr)
Bahcall EDR (EDR) 0.16 0.80 7.0
Gutierrez DR6 (DR6) 0.00 0.80 7.0
Rahmani DR7 (DR7) 0.02 0.74 6.7
This work (2015) DR12 (DR12) 0.04 1.00 7.9

For a CDM cosmology with , and from Planck+WMAP-9+BAO (Planck).

Note: Since we have a larger sample than Gutierrez, we expect a factor of improvement in the error just from purely statistical reasons. In Figs 9 and 10, it is shown that the error is dominated by the sky subtraction algorithm, which suggests that the performed analysis have reached the maximum precision with the available data.

(a) *


Table 1: Summary of the results obtained by recent works based on the [O iii] emission line method for the possible variations of the fine-structure constant.

Since Dirac’s philosophical argument (Dirac) against the fixed value of fundamental constants of Nature, several experiments have been performed to constrain possible variation on dimensionless constants of physical theories. Fundamental constants of physics could be thought of as parameters which enter in our description of Nature but they cannot be predicted with our current theories and should be measured. Dirac’s idea is based on the unlikely fact that the most fundamental constants of the Universe have a certain fixed value (at a given energy) with no apparent relation with the real world. It is more likely that their present values are the result of a dynamical process, which had yielded the fundamental constants as they are measured today. Therefore, they should be considered as characterizing the state of the Universe (Uzan2003). There are many current theoretical frameworks which allow for such variation of the fundamental constants, for instance, string theory (Maeda), modified gravity and theories with extradimensions (e.g. Clifton). Moreover, the experimental bounds on their variation have become a stringent test for those theoretical models (e.g. Thompson; Leal). The most studied fundamental constants are the fine-structure constant , the Newton gravitational constant and the electron-to-proton mass ratio (Uzan2003; Uzan2011; Garcia-Berro).

The fine-structure constant governs the electromagnetic coupling between photons and charged particles . The current constraint on its relative variation , over geological time-scales, is up to (2 Gyr ago); obtained from the Oklo phenomenon (e.g. Petrov). It has also been reported up to ( Gyr ago) from meteorites (Olive); which also excludes possible variations on the scales of the Solar system. On the other hand, there are also constraints, , based on the cosmic microwave background (CMB; Landau; Planck) at and from big bang nucleosynthesis, the latter being model-dependent. By measuring fine-structure multiplets at different redshift in the absorption or emission spectra of galaxies and quasars, located at different directions in the sky, one can measure an estimate of the variation of with time or space over cosmological scales.

The first measurements on the variation of from astronomical observations reached an accuracy of (Savedoff; BahcallSalpeter; BahcallSchmidt; BahcallSargent). Since then, the methodology and understanding of systematics has dramatically improved. Current measurements of absorption multiplets along the line of sight of three quasars around redshift 1.5, observed with spectral resolving power at UVES/ESO-VLT, reached the level (Evans). Using emission lines, an accuracy of was achieved analysing quasar spectra at (Gutierrez; Rahmani), taken with the Sloan Digital Sky Survey (SDSS) spectrograph.

The measurements on absorption features on a quasar spectrum are currently limited by the precision in the absolute wavelength calibration of the spectra, i.e., m s using spectra with (Molaro; Evans; ESO). Furthermore, the so-called many-multiplet (MM) method used in Evans, although more precise, remains controversial as several assumptions are made, the most important one being ionization and chemical homogeneity. These assumptions may induce systematic biases on the value of .

In this article, we use the method based on the [O iii] emission lines, first proposed by BahcallSalpeter, which is less affected by systematics. In particular, there is no need for assuming ionization and chemical homogeneity, since the studied lines have the same profile (the transitions originate at the same upper energy level). Furthermore, the emission-line method suffers of much less spectral distortion, since the measurements of are done on a spectral window  Å  as compared to  Å when the MM method is used. With a large ensemble of quasars and/or using high-resolution spectroscopy, the uncertainty can be reduced significantly, and will compete with the absorption method when using high-resolution spectroscopy.

Figure 1: Left-hand panel: sky distribution of the full SDSS-III/BOSS DR12Q quasars (297 301) in J2000 equatorial coordinates. Right-hand panel: number of quasars with [O iii] emission lines in our fiducial sample (10 363 quasars) in bins. S/N (10 363 quasars), black solid line; S/N (4015 quasars), blue dashed line; and S/N (1498 quasars), red dotted line.

The beginning of the SDSS survey opened a new era of precision, allowing us to use big samples of quasars; thus, reducing the statistical uncertainty of the measurement of (see Table 1). Here, we extend these works by using the SDSS-III/BOSS Data Relase 12 (SDSS-DR12; DR12), which covers the full Baryonic Oscillation Spectroscopic Survey (BOSS) survey footprint with an area coverage of 10 000 deg. In contrast to these previous investigations, we use spectra obtained with the current BOSS spectrograph (Smee) instead of the previous SDSS-I/II instrument, making our BOSS sample totally independent from previous works. Moreover, the spectral range of the BOSS spectrograph allows an extension of the redshift interval for the [O iii] doublet from to . The number of quasar spectra is increased by a factor of 5 with respect to SDSS-DR7. All these spectra have been visually inspected and classified as quasars by the BOSS collaboration, and their products are provided in the SDSS-III/BOSS Data Release 12 Quasar catalogue (DR12Q; see Paris12). For the final constraint on , we combine in this work the BOSS sample with the previously studied SDSS-II DR7 quasar sample.

There are several emission doublets, in addition to [O iii] (), that can be used to measure as noted by Bahcall, and first used by Grupe. Gutierrez analysed different doublets and found that the [Ne iii] () and [Si ii] () doublets appear in quasar spectra with sufficient frequency to have a meaningful sample. Results for [Si ii] are consistent with no variation of the fine-structure constant, although the uncertainty is an order of magnitude bigger than for [O iii] , and this doublet can only be used at low redshift for optical spectra. However, they obtained a positive variation of the fine-structure constant, , when the [Ne iii] lines are used. No explanation was found for this positive variation. In this work, we also analyse the [Ne iii] lines to check whether the same effect is present in our BOSS quasar sample.

There are investigations which use Si iv absorption lines () to obtain a precision of (Chand). This method also avoids the assumption of ionization and chemical homogeneity. However, since the separation between both lines is only , the wavelength precision needed in the laboratory for the separation between both lines is five times higher than using [O iii] lines. Nevertheless, these constraints apply to the redshift interval , which does not overlap with our range, thus they are complementary to the ones reported in this research.

Finally, in the light of the upcoming large galaxy surveys, like eBOSS and DESI, that will provide millions of high-redshift galaxy spectra, we also discuss using galaxies instead of quasars to set constraints on the fine-structure constant.

The paper is organized as follows. First, in Section 2, we describe the data set used for our analysis. Next, in Section 3, the methodology is presented, the emission-line method is explained, and the code and simulations to analyse the spectra are described. In Section 4, we study several samples to check for systematics. Then, our results are presented in Section 5. Finally, we provide in Section 6 a summary of the main conclusions achieved with this research project.

2 Sample description

All the spectra used in this investigation were downloaded from the SDSS Database. This survey (York), which began taking observations in 1998, consists of a massive collection of optical images and spectra from astronomical objects including stars, galaxies and quasars. For this purpose, there is a dedicated 2.5-m wide-angle optical telescope at Apache Point Observatory in New Mexico (USA; for more details, see Gunn). The third phase of this project (SDSS-III; Eisenstein) includes BOSS (Dawson) among its four main surveys. The data analysed in this research were provided by BOSS, and it is used for measuring for the first time. The SDSS-III/BOSS pipeline (Bolton) classifies the objects as quasars with a minimization procedure to fit the observed spectrum to multiple galaxy and quasar spectrum templates for all allowed redshifts. Then, a visually-inspected quasar catalogue is built from these objects. Our fiducial sample is obtained from the DR12Q catalogue version (Paris12).

The wavelength coverage of the SDSS-III/BOSS spectrograph is 3600-10 400 Å and that of the SDSS-II spectrograph is 3800-9200 Å. The BOSS sample is homogeneous since all the spectra have been obtained with the same instrument, and it is independent from previous investigations. The wider coverage of the new spectra allows consideration of higher redshifts (up to for [O iii] doublet) than in the previous SDSS-II analysis based on the same method (see Table 1). The BOSS spectrograph has two channels (blue and red) whose wavelength coverage is - and -, respectively. The resolving power ranges from 1560 at to 2270 at (blue channel) and from 1850 at to 2650 at (red channel). For our sample, the [O iii] lines fall in the red channel for of the quasars. The number of pixels of each spectrum is about for the BOSS spectra and for the SDSS-I/II spectra. The pixel spacing is uniform in log-wavelengths (). More complete information about the SDSS-I/II and BOSS spectrographs can be found in Smee.

2.1 Data selection

Figure 2: Left-hand panel: data points (1416) for which , namely the absolute value of the difference between the measured line separation at redshift in rest frame and the local one, is bigger than 1 Å plotted as a function of redshift (and the wavelength observed for [O iii] 4960). We compare with a typical sky spectrum: the [O iii] positions for these spectra correlate with the sky emission lines. Hence, these high values of are due to bad sky subtractions and/or low S/N. These spectra are removed. Right-hand panel: a spectrum removed from the sample because of the sky emission-line criteria. For this quasar, we get  Å. The weak [O iii] line is affected by the subtraction of the 7995 Å OH sky emission line, indicated by the vertical red dashed line.

The SDSS-III/BOSS DR12Q catalogue contains 297 301 objects. Fig. 1 (left-hand panel) shows the quasar distribution in the sky. We summarize below the main selection criteria in order to define our fiducial sample from this catalogue.

  1. Redshift . This limitation is imposed by the wavelength range of the BOSS optical spectrograph and the position of the [O iii] lines. This criterion decreases the sample down to 45 802 quasars.

  2. S/N. We impose a mild constraint on the signal-to-noise ratio (S/N) of the stronger [O iii] line (5008 Å) in order to preserve a large number of spectra. Constraints on the expected width and amplitudes of the lines help in avoiding misidentifications of the [O iii] doublet (see Section 4). This selection reduces the sample from 45 802 to 13 023 objects.

  3. Non-converging fits. Since we analyse spectra with low S/N, there are some cases where the Gaussian fit to the lines does not converge. 1 244 spectra are discarded, leaving us with 11 779 spectra.

  4. Sky emission lines. Strong atmospheric lines, for instance the O i 5578 Å line, are poorly or not completely removed by the SDSS sky subtraction algorithm. This may lead to a wrong identification of the [O iii] lines and to include low S/N spectra (Gutierrez). Both effects will produce outliers. We use the SDSS sky mask for Ly forest studies which contains 872 lines (see Delubac, for more details) to remove spectra whose [O iii] lines lie within a particular distance from the strongest sky lines. Even though we vary the distance [O iii] – sky lines, use different set of sky lines (according to their intensity), or evaluate other conditions (S/N, fit errors, etc.) to remove affected spectra; we usually eliminate good spectra for each bad spectra eliminated. Thus, these tests decrease significantly the number of quasars while not being very effective: typically of the outliers are not removed. Thus, we decided to eliminate all spectra for which the separation between both lines differ by more than 1 Å from the local value (see the last paragraph in Section 3.3). Fig. 2 (left-hand panel) shows that the distribution of these outliers is correlated with a typical sky spectrum. From a visual inspection, we observed that these spectra have low S/N, and they are in fact contaminated by sky emission line subtraction (see right-hand panel of Fig. 2). This effect causes us to discard 1 416 spectra ( of the previous 11 779 quasars). Finally, we have 10 363 quasar spectra (our ‘fiducial sample’).

The presence of broad H  emission line (4861 Å) near the weak [O iii] line 4960 Å could produce a blueshift in the determination of the [O iii] line position. This could mimic a positive variation on the fine-structure constant. Therefore, a constraint on the strength and/or width of the H  emission line has been imposed on previous investigations (Bahcall; Gutierrez; Rahmani). However, we do not restrict any characteristic of the H  line in our fiducial sample. We obtain a weighted mean for using as weights the uncertainty in computed with the standard errors for the position of the lines derived from the Gaussian fits. The contamination of H  is automatically taken into account. For instance, a broad H  line near the [O iii] 4960 line means a bad Gaussian fit. Thus, we obtain larger errors in the position of the line centroids and, consequently, in . In Section 4, we analyse several samples where the S/N is constrained to check that the H  contamination has little weight on the final constraint value.

An electronic table is published along with the paper which contains all the information of each spectrum from our fiducial sample of 10 363 quasars (see Appendix A).

The distribution of the selected quasars in redshift according to their selected S/N is plotted in Fig. 1 (right-hand panel). Fig. 3 (left-hand panel) displays a composite image built with all the spectra from our fiducial sample sorted by redshift. The right-hand panel shows the [O iii] doublet in rest frame.

Figure 3: Composite image with our fiducial sample of 10 363 BOSS quasar spectra sorted by redshift. Left-hand panel: the whole range of wavelengths is shown. From right to left, the strongest emission lines are H  6565 Å; [O iii; H  4861 Å; H  4341 Å; [Ne iii; [O ii] 3730 Å; [Ne v] 3426 Å; Mg ii 2796 Å and C iii] 1906 Å. The narrow straight line at 5579 Å is the strong [O i] atmospheric line. Right-hand panel: wavelength interval centred at the [O iii] doublet in rest frame.

3 Methodology

3.1 Measurement method

To first order, the difference between the energy levels of an atom is proportional to . Transitions between energy levels of the same atom at a given ionization level, with the same principal quantum number and different total angular momentum , have an energy difference proportional to . These groups of transitions are called fine-structure multiplets. Savedoff first realized that the fine structure of these energy levels could be used to break the degeneracy between the redshift effect and a possible variation of .

The value of the fine-structure constant can be measured through the separation between absorption or emission multiplets in the spectra of distant quasars (Uzan2003) as


where are the wavelengths of the transitions and subscript 0 and stand for their value at redshift zero (theoretical/laboratory values) and at redshift , respectively. For illustrative purposes, expression (1) can be approximated by


where is the local separation between both wavelengths, and is the difference between the measured line separation at redshift in rest frame and the local one. Thus, in principle, the larger the difference between the pair of lines, the better the precision for measuring .

Concerning emission lines, the most suitable pair of lines is the [O iii] doublet, which is often present in quasar spectra with relatively high-S/N. The vacuum values for the [O iii] doublet wavelengths are


which are published in the NIST Atomic Spectra Database.111 These transitions are forbidden (they correspond to magnetic dipole and electric quadrupole transitions), and they are not observed in the laboratory. The wavelength experimental values are obtained indirectly by first computing the energy levels from observed wavelengths using a theta-pinch discharge (Pettersson). The wavelength separation has directly been measured in the infrared from H ii regions using a balloon-borne telescope and Michelson interferometer (Moorwood). Both measurements of the wavelength separation, indirectly with the theta-pinch discharge and directly with the Michelson interferometer, are in good agreement, being the Michelson interferometer more accurate with an error Å.

From equation (2), a determination of with a precision of allows for an uncertainty of in when using the [O iii] doublet. The precision from the NIST atomic data allows for a determination of up to , which is a bit less than the uncertainty in our result. One could perform a blind analysis in order to search for a possible variation on , where the absolute wavelength values are not required, if one had a large enough sample distributed in redshift. However, the precision on the absolute wavelengths limits the usefulness of high-resolution spectroscopy until better measurements of the [O iii] lines (or just their separation) are available.

3.2 Implementation

Figure 4: Seventh-order polynomial fits (red) to the continuum spectrum with their residuals for three typical quasar spectra at different redshift. The gaps in the residuals are the masked regions corresponding to (from right to left) , the [O iii] doublet, , , and Mg ii (black dashed lines).

The code developed for the analysis of the quasar spectra follows the one described in Gutierrez, although there are some modifications and more information has been extracted from the analysis. We describe the main characteristics of our code below.

3.2.1 Wavelength sampling

We consider only the experimental data together with their errors as processed by the SDSS pipeline to obtain the constraint on the possible variation of . We do not resample the wavelength range by using an interpolation method. Since the pixel spacing is uniform in log-wavelengths, a given range of wavelengths in rest frame has the same number of pixels , i.e.


and is independent of the redshift of the object. All the wavelength intervals with the same width in rest frame will have the same number of experimental points.

3.2.2 Fit of the continuum spectrum

First, we fit a seventh-order polynomial to subtract the continuum spectrum while masking regions where strong and wide emission lines are present (, , , , Mg ii and the [O iii] doublet). Our method differs from Gutierrez in that they use a cubic local spline to fit the continuum masking strong emission lines. The chosen order of the polynomial provides enough degrees of freedom to reproduce different continuum features. In Section 3, we test how our measurement for is affected by changing the polynomial order. Hundreds of continuum spectra fits were checked by eye. The residuals from the fits are smaller than the errors on the flux densities. Fig. 4 shows three different spectra with their continuum fit and residuals.

3.2.3 Signal-to-noise ratio

We follow Gutierrez for the determination of S/N. Hence, we compute the standard deviation of the flux between and (where is the redshift of the quasar) where there are no strong emission or absorption lines. Then, we search for the maximum of the [O iii] 5008 line, and determine S/N as the ratio between the maximum of the line and the previously computed standard deviation. Although for a more reliable determination of the S/N, it is better to use a Gaussian fit to the line. This procedure avoids possible issues related when fitting data with very low S/N. This S/N is used in the criterion ii (Section 2) to build our fiducial sample.

3.2.4 Measurement of the emission-line wavelengths

To measure the wavelengths of the [O iii] doublet, our fitting code needs as input an accurate estimate of the redshift of the quasar, at least with an error . This allows a search for the emission lines in a window around the expected location of the [O iii] lines. The SDSS pipeline provides a determination of the redshift based on a fit to different templates; we refer to Bolton for more details. These redshift estimates have errors between and , which are sufficient for our purposes. Moreover, there is also a visual redshift estimation which can be found in the quasar catalogue DR12Q (Paris12). The difference between both redshift estimates (if any) is usually . We decided to adopt the visual redshifts.

Figure 5: [O iii] (left-hand panel) and [Ne iii] (right-hand panel) lines for SDSS-J121417.80+293143.4, at redshift . The measured for this quasar is , and , . The measured for [Ne iii] is not consistent with zero regardless of the profile; see Fig. 14 and last paragraph of Section 4 for discussion. Each panel shows the flux density for each pixel with their respective error bars (solid symbols), together with the Gaussian fit (dotted red curve) and the pseudo-Voigt profile (thick grey curve) to each of the lines. The fitting procedure (described in the text) only takes into account the experimental data (solid symbols) weighted by their error bars. Notice how the deviation of the line centroid position derived from our Gaussian fit (vertical solid line) with respect to the expected position of the line (vertical dashed line) according to the visual redshift provided by the DR12Q catalogue are well correlated for the same pair of [O iii] and [Ne iii] lines, and for the different set of lines. The green shaded vertical areas highlight the uncertainty for the expected position of the lines due to the quasar redshift error (). Also shown is a fourth-order spline interpolation to the spectrum after subtracting the continuum (thin solid line). The [Ne iii] lines are weaker by one order of magnitude than the [O iii] lines, which is usually the case for all the spectra showing both pair of lines. The weak line near the stronger [Ne iii] line is blended with He i and H .

The centroid positions of the [O iii] emission lines are determined by four different methods.

  1. Gaussian profile method.

    First, we search for the maximum flux value in an window around the expected position of the line (according to the redshift provided by the DR12Q catalogue). This procedure automatically erases any bias produced by the redshift value. Then, we make an initial Gaussian fit around the position of the maximum flux value using a fixed width of . From this first fit, we obtain a new position for the line centroid and a Gaussian width. These values are used as initial parameters for the final fit of the lines; namely, the wavelength range considered to perform the final fit is centred around the position of the line centroid, and it is four times the Gaussian width of the lines. This approach means that we consider pixels up to away from the centre of the line. Hence, some lines are fitted using pixels, while others with pixels depending on the line width. The fit takes into account the flux errors for each pixel, i.e., we use the column found in each spectrum as weights for the fit. Our final centroid measurement for each considered line corresponds to the centroid of the Gaussian fit done in the last step of the adopted procedure. We also derive an error for using the standard errors for the centre position of the Gaussians. This is our main method for measuring .

  2. Voigt profile method.

    Following the same procedure than when using a Gaussian profile, we make the fit with a Voigt profile instead of a Gaussian. More precisely, we use a pseudo-Voigt profile which is a linear combination of a Gaussian and a Lorentzian profile. Then, we have one more parameter, i.e. the amplitude of the Lorentzian function, while its width and its position are the same as those for the Gaussian profile.

    In Fig. 5, we depict the [O iii] and [Ne iii] lines for the same quasar spectrum to illustrate the Gaussian and Voigt fitting methods.

  3. Integration method.

    Here, the centroids of the lines are obtained by integrating around from the position of the fitted Gaussian profile. This technique provides indications of whether there is contamination. However, due to the mid-resolution of the spectra , this method is not very accurate.

  4. Modified Bahcall method.

    In Bahcall the authors used a different approach to compute the line positions. They performed a third-order spline interpolation to the stronger [O iii] 5008 line, then fitted this interpolation to the weaker 4960 line by adjusting the amplitude and separation of the profile. We have modified this method by using a Gaussian fit to the stronger line rather than a third-order spline.

Although we have described four different methods, the main results for presented in this work are based on the Gaussian fitting method, while the other three are used only for comparison (see Section 4).

Finally, our final result for and its error is obtained in the same way as in Chand, namely we compute a weighted mean and a weighted standard deviation, where the errors for of each spectrum are used as weights.

3.3 Simulated spectra

Figure 6: Left-hand panel: errors for obtained from the simulations (standard deviation of the measurements on 100 realizations of each real spectrum) and standard errors from the Gaussian fits for our fiducial sample. The solid line represents a one-to-one correspondence, while the dashed lines have slopes of 2 and 0.5. Only the simulation and fit errors smaller than are shown. Right-hand panel: errors estimated from the simulations as a function of redshift. Spectra with are shown as red crosses ( of the total). There is a clear division between two different set of spectra which correlates with the sky emission lines (see discussion in the main text).

Figure 7: Left-hand panel: measurements of using Gaussian and Voigt fitting profiles. Non-compatible measurements at 1 are shown as red crosses ( of the total). Right-hand panel: errors from the Gaussian and Voigt fitting. Non-compatible measurements at 1 are shown as red crosses.

In order to test the robustness and accuracy of our methodology, we generate realizations of quasar spectra using as noise a normal distribution centred at the flux value, and taking the error in each pixel as the standard deviation. From our fiducial sample (10 363 quasars), we simulate 100 realizations for each spectrum ( a million in total). This number of realizations provides reasonable statistics to derive an error from the standard deviation of the measurements on the realizations of each real spectrum, while the computation time remains reasonable ( d) using a standard-size computer. The estimated error derived from the simulations includes


where is the error derived from the Gaussian fits, which is our error estimate for each real spectrum; is the error from different continuum subtraction due to the Gaussian noise, and is the systematic error of our code. Then, we expect and their difference will be an indication of the continuum and systematic errors.

Fig. 6 (left-hand panel) shows the correlation between the error in from the Gaussian fits of each real spectrum and the standard deviation for of its 100 realizations. The standard deviations from the simulations are within a factor of of the standard errors from the fits for () of the cases when both quantities are (). This shows that our code and the continuum subtraction do not introduce noticeable systematic errors compared to the Gaussian fitting. However, there is a set of spectra ( of the total) for which the simulations provide much larger errors . Fig. 6 (right-hand panel) shows the errors from the simulations as a function of redshift for our fiducial sample. Red crosses stand for spectra whose Gaussian fit error (). The errors are distributed in two clouds of points. For the cloud with , the virtual realizations of each spectrum seem to differ significantly from the real spectrum. Since we use the error in each pixel to build the realizations, the relative error is large for these spectra, which is an indication of a low S/N ratio or large absolute errors in the pixels, for instance in wavelength regions with sky emission lines. In fact, the cloud with bigger errors mimics the sky spectrum. Note also the strong correlation between this cloud of points and the spectra with large Gaussian fitting errors (red crosses). The other set of points with are close to our error estimation on the measurement of based on the Gaussian fits.

As a further proof, we also simulate realizations of the 1416 dropped spectra because of sky emission lines (criterion iv, see Section 2). We found that more than 80 of the spectra have . This confirms that these spectra have very low S/N and/or large pixels error due to the proximity of the lines to strong sky emission lines.

3.4 Gaussian versus Voigt fitting profiles

The results obtained when using Voigt profiles instead of Gaussian ones are compared in Fig. 7. The Voigt and Gaussian measurements are -compatible for the of the cases ( at ). Regarding the errors, there is no clear improvement when using either of both methods. However, Voigt profiles have one more parameter and restrict the number of degrees of freedom. Due to the spectral mid-resolution and the fact that the [O iii] lines are very narrow, there are often only a few pixels to fit, which frequently lead to non-convergent fits. This reduces the quasar sample in quasars. Further discussion about both profiles can be found in Section 4.

4 Systematics

In this section, we examine the possible unnoticed systematic errors by analysing different quasar samples. Table 2 summarizes all the samples considered together with their mean redshifts and the measured value for .

We consider the following sources of systematic errors.

  1. Misidentification of the lines. The expected line widths and amplitudes are useful to avoid misidentification of the [O iii] emission lines. (a) Line widths: since both lines originate on the same upper energy level, their width must coincide. We check that this is the case by considering quasars whose [O iii] line widths are the same within a relative fraction. For more than half of our fiducial sample, the [O iii] line widths differ by less than 10 (see Table 2). (b) Amplitude ratio: atomic physics states that the amplitude ratio between the [O iii] 5008 and [O iii] 4960 lines is 2.98 (Storey) (as quoted in Section 5, we obtain ). Thus, we consider different samples where this ratio differs by less than a certain amount from (see Table 2). All the samples considered in this test yield results for compatible with zero. Fig. 8 displays the Gaussian widths and fluxes of both [O iii] emission lines for our fiducial sample.

  2. Windows for the Gaussian fits. We use a wavelength range of 2 around each [O iii] line in order to obtain the final Gaussian fit to the line profiles. We study how our results depend on this choice. By considering a larger wavelength interval, the results are more affected by the H  contamination and possible asymmetries on the line wings. The differences in the number of spectra for these samples [which are obtained by applying the selection criteria (i)(iv) discussed in Section 2.1] arise because of the criteria concerning the non-converging fits and the sky emission lines described in Section 2.

  3. H  contamination. We analyse samples where the ratio between S/N and S/N is constrained. Despite the fact that the value for decreases as we place more stringent constraints on H , it is always consistent with no variation in within the errors. This analysis demonstrates that the strength and/or width of the H  line do not affect substantially the result for when a weighted mean is adopted.

  4. Continuum subtraction. We use a seventh-order polynomial to subtract the continuum spectrum. We examine if the polynomial order has important effects on our measurements. Our values for and their errors are only slightly affected by the chosen polynomial order.

  5. Goodness of Gaussian fits. We quantify the quality of the Gaussian fits by the coefficient. All the considered samples show values for consistent with no variation in .

  6. Broad lines. We also study samples where the width of both lines is less than a certain value (in km s). These samples are consistent with no variation of . Samples built from narrow lines km s may be more affected by misidentification of [O iii] lines as sky lines.

    No. of quasar spectra Redshift
    No. of quasar spectra Redshift
    Fit width No. of quasar spectra Redshift
    S/N No. of quasar spectra Redshift
    Pol. order (cont.) No. of quasar spectra Redshift
    (both fits) No. of quasar spectra Redshift
    [O iii (km s) No. of quasar spectra Redshift
    Method No. of quasar spectra Redshift
    Gaussian (weighted)
    Modified Bahcall
    Gauss versus Voigt No. of quasar spectra Redshift
    Gaussian profiles
    Voigt profiles
    Mixed profiles
    Table 2: Results for considering several samples with different constraints. The number of quasar spectra, the mean and standard deviation of the redshift and the value for are shown.

    Figure 8: Left-hand panel: Gaussian widths (in km s) for both [O iii] lines. Both lines originate on the same upper energy level, then their widths must coincide (red dashed line). Right-hand panel: fluxes for both [O iii] lines. The theoretical flux ratio is (red dashed line). The entire fiducial quasar sample is shown.
  7. Different methods for measuring the [O iii] line position. We compare the results obtained by the methods to measure the position of the [O iii] lines described in Section 3.2.4. Since not all the methods provide an error for the measurement, we cannot calculate a weighted mean, and it is necessary to select a more restricted sample. Then, we consider a sample where the difference between the widths of the lines is less than , the amplitude ratio is constrained to differ from the theoretical value 2.98 (Storey) by less than 0.5, and the S/N is smaller than half the S/N.

  8. Gaussian versus Voigt profiles. We compare the results for 8485 quasars from our fiducial sample after dropping 1878 spectra with non-converging Voigt fits (this reduction increases the statistical error). We also compute a ‘mixed’ value for where for each spectrum we use the value for the variation of the fine-structure constant with smaller error, either or .

We have also analysed the standard deviation and errors of the results for as a function of redshift (Fig. 9). Even though we have imposed a constraint on our initial sample based on the sky emission lines, the standard deviation and errors still correlate with the sky. In particular, for the correlation with the moving standard deviation, this means that the precision in our measurement of along the whole redshift interval is limited by the sky subtraction algorithm.

Figure 9: Errors for derived from the Gaussian fits (grey points) for our fiducial sample, moving mean of the these errors (blue line) using overlapping bins (100 spectra per bin, ), moving standard deviation of measurements using the same bins (red line) and a typical sky spectrum.

5 Results

5.1 [O iii] lines

We used a total of 10 363 quasar spectra, drawn from the SDSS-III/BOSS DR12Q catalogue, after applying the selection criteria (i)(iv) (see Section 2), to measure the possible variation of the fine-structure constant. The following measurement is obtained:

This value is consistent with the previous results reported in different investigations based on the same method: Bahcall, Gutierrez, and Rahmani. The redshift dependence of the measurements is shown in Fig. 10 (left-hand panel), where several bins have been made taking into account the redshift intervals affected by the sky (shaded zones). In the right-hand panel, we show the results obtained from the simulations described in Section 3, using the same redshifts intervals for the bins. The main differences between the real results and the simulations are in the regions where there are strong sky lines (shaded regions), while being in agreement in the remaining zones. Detailed information about each bin for the real data can be found in Table 3.

Figure 10: Left-hand panel: versus redshift (real data). Details about each bin are listed in Table 3. Right-hand panel: versus redshift (simulations). A typical sky spectrum and shadowed regions where the sky contamination is expected to be large, are shown as reference.

Our results are little affected by the specific constraints imposed in our sample as discussed in Section 4. For instance, we vary the width for the Gaussian fits, the contamination of H , the polynomial order used to fit the continuum spectrum, the quality of the Gaussian fits and test different methods to measure . The most important effect found is that by considering broader widths for the Gaussian fits, the results are more affected by the contamination from H  and possible asymmetries in the line wings. We have also checked for possible misidentifications of the [O iii] emission lines using their expected widths and amplitude ratio.

Table 4 contains the results for when the lower bound on the S/N is increased. All the results remain consistent with no variation of the fine-structure constant. In Fig. 11, the measured for our fiducial sample as a function of the S/N are plotted together with their errors.

Redshift interval No. of quasar spectra Redshift
Table 3: Detailed information about the bins in Fig. 10.
S/N No. of quasar spectra Redshift
Table 4: Results for several samples with different constraints on the S/N. For each sample, the number of quasar spectra, the mean redshift, together with its standard deviation and the value for are shown.

Figure 11: versus S/N (top panel) with the moving standard deviation (black lines) using overlapping bins (100 spectra per bin) and the error on (bottom panel) in linear-log scale for our fiducial sample. The deviation of from zero and its error steadily decreases as the S/N increases.

The distribution of BOSS quasars in the sky (see Fig. 1, left-hand panel) suggests to divide the sample into two, one for the North galactic cap and one for the South galactic cap. Table 5 describes the results for each galactic hemisphere, and no statistical meaningful difference is found. In order to look for a spatial variation, we also carried out a more precise analysis by fitting a dipole. First, we fixed the direction on the sky of the dipole and performed a linear fit () of the measurements of the variation of the fine-structure constant as a function of , where is the angle between the dipole axis and a quasar in the sky. Different fits were done for the dipole axis lying in a grid in RA and Dec . (). However, there is not statistical significance to determine the dipole axis with a meaningful error, i.e. smaller than the whole sky. There has been a claim on a significant deviation of from being a constant as a function of space (King), with a dipole amplitude in the direction RA h and Dec.. Fixing the dipole in that direction, we get for the dipole amplitude and for the monopole term, which are not precise enough to compare with that work.

Galactic hemisphere No. of quasar spectra Redshift
Table 5: Results for the North and South galactic hemispheres.

We are inclined to parametrize the possible time variation of with redshift . This is justified since any possible variation on must be dominated by the local geometry of space-time (at least if we consider the dynamics of the Universe as the main reason for such variation). Therefore, one is led to consider the possible variation of as a function of redshift () or the Ricci scalar (), where is the scale factor, the Hubble parameter and is the deceleration parameter. Since the Ricci scalar is not known for each quasar, it is straightforward to consider a possible variation with redshift. In contrast, for a time parametrized model of the variation of the analysis depends on the particular cosmology considered. Since there is no significant clear dependence, we use a linear model in redshift. Then, for


we obtain


which do not show any dependence of with redshift.

From this sample, we also obtain a value for the line ratio , where is the maximum flux density of the line, and is the Gaussian width. The value reported is a weighted mean where the S/N is used as weights. The quoted systematic error is computed from the analysis of samples with different polynomial orders for the continuum fit and different widths for the line fitting (see Table 6), since this quantity is more affected by these two parameters. The value we obtain is in agreement with the best current theoretical value, i.e., (Storey).

Polynomial order No. of quasar spectra Redshift
3 2.96
5 2.94
7 2.96
9 2.98
Fit width No. of quasar spectra Redshift
Table 6: Results for the line ratio when polynomials of different orders are used to subtract the continuum and different range for the Gaussian fits are used. For each sample, the number of quasar spectra, the mean redshift together with its standard deviation and the value for are shown.
Sample No. of quasar spectra Redshift
DR7 (SDSS cont.)

Table 7: Results for SDSS-II/DR7 and BOSS (SDSS-III/DR12) samples. For each sample, the number of quasar spectra, the mean redshift, together with its standard deviation and value for are shown.
Redshift interval No. of quasar spectra Redshift

Table 8: Values of using the combined sample BOSS+DR7 for the same redshift bins as in Table 3.

Figure 12: Difference between and measurements for the 41 quasars observed by SDSS-I/II that were re-observed by BOSS. Both values for are consistent.

Finally, we have also considered quasar spectra from the SDSS-II/DR7, which were observed using the previous spectrograph instead of the upgraded BOSS spectrograph (see Section 2). From the DR7 quasar catalogue (Schneider), which contains 105 783 quasars, we select a sample of 2853 quasars up to redshift using the same criteria described in Section 2. This number is similar to the quasar spectra considered by Rahmani. We re-analyse this sample using the methodology presented in this work, and we find . By combining this DR7 sample with our fiducial BOSS (DR12) quasar sample (after eliminating 41 spectra which were re-observed by BOSS), we obtain our final sample which contains a total of 13 175 quasars. The value obtained for this combined sample is reported as a final result of this investigation:

Table 7 shows the results for DR7, DR7 using the continuum fit provided by the SDSS pipeline,222The SDSS pipeline provides a continuum fit for the DR7 spectra. The good agreement between the value for obtained with the SDSS continuum fit and our continuum fit is a good test for our code. BOSS (DR12) and the combined BOSS+DR7. It can be seen that the mean redshift for the DR7 sample is lower than that for BOSS. Note that there is also a big difference on the mean S/N of both samples: S/N and S/N, which also explains why the statistical errors for do not reflect the expected reduction due to the increase in the number of quasars of our BOSS sample. Table 8 shows the results of using the combined sample in the same redshift bins considered for our fiducial sample. Fig. 12 shows the difference of the values obtained for for the 41 re-observed quasars. Both BOSS and DR7 measurements are in perfect agreement within the error bars. This test is a good check for the reliability of our code and the consistency of the SDSS spectra obtained with different spectrographs.

There are massive galaxy surveys to be carried out during the next decade. For instance, eBOSS and DESI will take spectra from millions of galaxies. Therefore, it is interesting to give an estimation of the accuracy that will be obtained when using galaxy spectra instead of quasars. For this, we have analysed the galaxy spectra collected by the DEEP2 survey (DEEP2) taken with resolving power . From this sample, we found 4056 galaxies with strong [OIII] lines. Naively, one would expect that the error on should be


where takes into account the different number of objects in each sample, stands for the different resolution of the spectra and is an extra factor due to the different characteristics of quasar and galaxy emission lines which affect the uncertainty of the line positions. This last factor is proportional to the line widths and inversely proportional to the line fluxes:


These numbers for the [O iii] 5008 line are approximately FWHM km s, FWHM km s, Flux and Flux in units of erg cm s, obtained from the DEEP2 sample and from our fiducial sample. Thus, the expected error is , where we have considered the error of our fiducial sample . Applying the same criteria described in Section 2 to the DEEP2 galaxy sample, we get . Thus, the upcoming future galaxy surveys will be quite competitive for constraining the variation of the fine-structure constant at low redshift . Fig. 13 shows the error on for the DEEP2 and BOSS samples.

Figure 13: Errors on as a function of the FWHM of [O iii] 5008 for our fiducial BOSS quasar sample (blue points) and for the DEEP2 galaxy sample (red crosses).

5.2 [Ne iii] lines

Figure 14: [Ne iii] and [O iii] measurements for . Empty symbols stand for spectra with S/N and solid squares for spectra with S/N. [Ne iii] measurements have a clear tendency to a positive variation of , which is due to a systematic effect affecting the [Ne iii] measurement. The same effect has already been noticed by Gutierrez, and it is explained in Section 5.2.

Figure 15: Left-hand panel: [Ne iii] 3869 and [Ne iii] 3968 together with our Gaussian fits (solid line) from a high-resolution () spectrum of the Planetary Nebula IC 418 obtained with the FIES spectrograph at the NOT telescope. Right-hand panel: IC 418 spectrum centred at [Ne iii] 3968 line. The two close lines are H  3971 Å and He i 3965 Å.

Figure 16: Left-hand panel: median-stacked quasar spectra with broad [Ne iii] lines (increasing line width from bottom to top). Both [Ne iii] lines are shown (). The weak [Ne iii] line is blended with the two lines H  3971 Å and He i 3965 Å. Right-hand panel: shift produced by the H  line in the [Ne iii] 3968 line as a function of the line intensity ratio of both lines as measured from the Planetary Nebula convolved spectrum.

We also measure from 462 quasar spectra with [Ne iii] emission lines the following constraint on the fine-structure constant:


to be compared with


obtained by Gutierrez. The analysis of the [Ne iii] lines reveals the same systematic effect previously observed, namely a clear tendency for a positive variation of . Fig. 14 compares the results obtained for for spectra where both [O iii] and [Ne iii] lines are present. To account for this effect, a shift on the theoretical or observed values of the wavelengths for the [Ne iii] lines is necessary. There are experimental (NeIIIlines) and indirect (Kramida) values for the wavelengths of the [Ne iii] lines which are in agreement with errors . We use the NIST values for the [Ne iii] lines


The results for the [O iii] doublet guarantee the good calibration of the SDSS spectra (and many more independent scientific results based on the SDSS spectra). Thus, we have measured the [Ne iii] lines using a high-resolution optical spectrum from the planetary nebula IC 418. The IC 418 optical spectrum (∼3600-7200 Å) was taken under service time at the Nordic Optical Telescope (NOT; Roque de los Muchachos, La Palma) in 2013 March with the FIES spectrograph. We used FIES in the low-resolution mode () with the 2.5 arsec fibre (centred at the central star of IC 418). Three exposures of 1200 s each were combined into a final IC 418 spectrum, reaching a S/N (in the stellar continuum) of 60 at 4000 Å and in excess of 150 at wavelengths longer than 5000 Å (see Anibal, for more observational details). To measure , we need to know the ratio


which is independent of the peculiar velocity of the planetary nebula. From our data, we obtain


compared to the one using NIST values for the wavelengths


The difference between the two values translates into a variation on . Thus, the measured wavelength separation for the [Ne iii] doublet does not account for the positive variation on observed using these lines. Fig. 15 (left-hand panel) shows the Gaussian fit to the [Ne iii] line profiles present in the IC 418 spectrum.

The IC 418 spectrum shows two different lines near the [Ne iii] 3968 Å line (see Fig. 15, right-hand panel). The stronger one is H  3971 Å, the other one is He i 3965 Å. Hence, we search for a possible blending of the [Ne iii] line 3968 with these two lines in our much lower spectral resolution quasar spectra. Fig. 16 (left-hand panel) shows stack quasar spectra with broad [Ne iii] emission lines. It can be seen that the weak [Ne iii] line is blended.

To quantify the displacement produced by the blending with H  line, we did a Gaussian convolution of the Planetary Nebula spectrum to lower the resolution down to . Since the line intensity ratio of [Ne iii] and H  may differ in the quasar narrow emission-line region and the Planetary Nebula, we show in Fig. 16 (right-hand panel) the shift produced by the H  line as a function of the ratio [Ne iii] /H . We get a shift 0.6 Å when H /[Ne iii] is 0.5. This explains the systematic found when using [Ne iii] lines to measure the variation of the fine-structure constant in previous studies (Gutierrez; Grupe).

6 Summary

The main conclusions of this work are as follows.

  1. From 45 802 objects at classified as quasars in the SDSS-III/BOSS DR12 quasar catalogue, we have extracted a sample of 10 363 quasars with [O iii] emission lines. Combining this fiducial sample with a sample of 2853 previously studied SDSS-II/DR7 quasars, we got a final sample of 13 175 after eliminating 41 re-observed quasars.

  2. With this combined sample, we have estimated a value for the possible variation of the fine-structure constant of , which represents the most accurate result obtained with this methodology.

  3. We have also studied how much our results change when analysing the fiducial sample according to different properties (width, amplitude, S/N and coefficient of the [O iii] lines), and when modifying some parameters of the analysis (polynomial order for the continuum subtraction, different methods to determine the line position, e.g. Gaussian/Voigt profiles). We conclude that our results are quite robust, and they are consistent with no variation of the fine-structure constant.

  4. From over one million simulated realizations of quasar spectra, we conclude that the precision of our emission-line method is dominated by the error from the Gaussian fits. Hence, the error from the continuum subtraction and any possible systematics from our code are small.

  5. The standard deviation of the results as a function of redshift correlates with the sky. This result suggests that our main source of uncertainty is determined by the sky subtraction algorithm.

  6. We have determined the ratio of the [O iii] transition lines to be , which is in good agreement with previous experimental and theoretical values.

  7. The same systematic effect previously noticed by Gutierrez has been found on the [Ne iii] lines measurement. Incorrect measurement for the separation of the [Ne iii] has been excluded as a possible explanation, and a blending of the H  and the [Ne iii] 3968 has been identified as the source of this effect.

  8. The measurement of using SDSS-III/BOSS spectra has reached the maximum precision unless better sky subtraction algorithms are developed. To obtain better constraints () using the emission-line method, high-resolution spectroscopy () is mandatory.

  9. We note that future large galaxies survey like eBOSS or DESI could provide quite stringent constraint for at low redshift, following our analysis of galaxy spectra taken from the DEEP2 survey.


FDA, JC, and FP acknowledge support from the Spanish MICINNs Consolider-Ingenio 2010 Programme under grant MultiDark CSD2009-00064, MINECO Centro de Excelencia Severo Ochoa Programme under grant SEV-2012-0249, and MINECO grants AYA-2012-31101 and AYA2014-60641-C2-1-P. FDA also acknowledges financial support from ‘la Caixa’-Severo Ochoa doctoral fellowship, UAM+CSIC Campus of International Excellence and Instituto de Astrofísica de Canarias for a summer stay where this work began. AM and DAGH acknowledge support provided by the Spanish Ministry of Economy and Competitiveness (MINECO) under grant AYA-2011-27754.

Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is

SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, University of Cambridge, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofísica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

This article is also partially based on service observations made with the Nordic Optical Telescope operated on the island of La Palma by the Nordic Optical Telescope Scientific Association in the Spanish Observatorio del Roque de Los Muchachos of the Instituto de Astrofísica de Canarias.

Appendix A

We publish along with this paper an electronic table with the combined SDSS-III/BOSS DR12 and SDSS-II/DR7 sample of 13 175 quasars used in this work. Table A1 describes the information and format of each column. The table is available in the following link 08/11/stv1406.DC1/

Column Name Format Description
1 SDSS_NAME STRING SDSS-DR12 designation (J2000)
2 RA DOUBLE Right Ascension in decimal degrees (J2000)
3 DEC DOUBLE Declination in decimal degrees (J2000)
5 PLATE INT32 Spectroscopic plate number
6 MJD INT32 Spectroscopic MJD (55 000 SDSS-III/BOSS spectra, 55 000 SDSS-II spectra)
7 FIBER INT32 Spectroscopic fiber number
8 Z_VI DOUBLE Redshift from visual inspection
9 Z_PIPE DOUBLE Redshift from BOSS pipeline
10 ERR_ZPIPE DOUBLE Error on BOSS pipeline redshift
11 ALPHA FLOAT from the Gaussian fits
12 ERR_ALPHA FLOAT Standard error for from the Gaussian fits
13 SN_O1 FLOAT S/N for the [O iii] 4960 line
14 SN_O2 FLOAT S/N for the [O iii] 5008 line
15 O1_FIT FLOAT Line centroid for the [O iii] 4960 line
16 O2_FIT FLOAT Line centroid for the [O iii] 5008 line
17 ERR_O1 FLOAT Error on the line centroid for the [O iii] 4960 line
18 ERR_O2 FLOAT Error on the line centroid for the [O iii] 5008 line
19 O1_AMPLITUDE FLOAT Gaussian amplitude at the centre for the [O iii] 4960 line
20 O2_AMPLITUDE FLOAT Gaussian amplitude at the centre for the [O iii] 5008 line
21 O1_WIDTH FLOAT Gaussian width for the [O iii] 4960 line
22 O2_WIDTH FLOAT Gaussian width for the [O iii] 5008 line
23 FILE_NAME STRING File name to download from the SDSS server
Table A1: Description of the electronic table with the combined sample (13 175 quasars) published along with the paper.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description