Open paulmartini opened 1 year ago
@paulmartini I'd like to make some progress on this ticket but I could use your help with the templates.
First, you've suggested three possible template sets, V01, T06, and S07, but if we had to select one, which one would you select? I can fit subsets of objects with all three (in separate calls to FastSpecFit), but I suggest we select a fiducial set (or derive our own from DESI observations themselves??).
Could you (or a student) put the templates into a simple FITS format that I can use (and put them somewhere I can access on NERSC)? For example, here are the current templates-- https://data.desi.lbl.gov/public/external/templates/fastspecfit/README.txt https://data.desi.lbl.gov/public/external/templates/fastspecfit/1.0.0/
Second, can you provide more input regarding the power-law fitting you would like me to perform? I did play around with this previously and used a simple F_nu \propto \nu^(-alpha)
or F_lambda \propto \lambda^(alpha-2)
model, with alpha
ranging between 0.5 and 2.0. Is this what you had in mind, too? For reference, this is the paper I've been looking at, but let me know if you have other favorites--
https://www.aanda.org/articles/aa/pdf/2017/08/aa30378-16.pdf
Finally, do you want me to try to fit the host-galaxy continuum, too, or can we safely ignore it?
I'm looking forward to hearing your thoughts.
Hi @moustakas --
Great that you plan to work on this soon!
I've made a fits file with four templates and placed them in a single fits file here:
/global/cfs/cdirs/desi/users/martini/qsospecfit/fetemplates.fits
The first (EXTNAME = DEFAULT) is the variant of the Vestergaard & Wilkes (2001) template we've used the most. The three other templates are another Vestergaard & Wilkes template and the two others mentioned above. The EXTNAMEs should be clear, and there is more information in a COMMENT in each extension. Let me know if you'd like these arranged differently.
I suggest just using the one I called 'DEFAULT.' Zhefu Yu has done experiments with the others and not found big differences. The big disadvantage of the Vestergaard & Wilkes template is that there is not empirical data at the location of MgII, and ignoring an FeII contribution at those wavelengths probably leads to an overestimate of the MgII line flux and FWHM. Yet for scaling relations this may not matter too much, given the calibration process.
The power-law fit you mention is a good match to observations of the ultraviolet continuum, namely alpha
= 0.5 to 2.0 where F_nu \propto \nu^(-\alpha)
.
For quasars we can ignore the host-galaxy continuum. I think it is okay to ignore it.
@arjundey has recommended adding CII] 2326 to the line-list.
Whittle+85 (https://ui.adsabs.harvard.edu/abs/1985MNRAS.213...33W/abstract) recommend a few non-parametric line-profile measurements that may be useful.
Some additional AGN/QSO modeling refs (courtesy @bd-j!)--
Lyu+22 added semi-empirical AGN optical-FIR templates, and different SF dust emission templates - this is built off a very old prospector version, and I think uses a custom FSPS: https://ui.adsabs.harvard.edu/abs/2022ApJ...941..191L/abstract
This adds a power law continuum with attenuation independent from the stars - I hope it will be added to prospector soon. https://ui.adsabs.harvard.edu/abs/2024arXiv240302304W/abstract https://ui.adsabs.harvard.edu/abs/2024arXiv240302304W/abstract
This uses Lyu+22 and also a version with just a scaled Polletta 2006 QSO template added on: https://ui.adsabs.harvard.edu/abs/2024ApJ...968....4P/abstract https://ui.adsabs.harvard.edu/abs/2024ApJ...963..128B/abstract
This adds a scaled empirical AGN emission line template from Richardson 2014 - this is actually in prospector as AGNSpecModel, but was really just used to explore the potential impact of AGN emission lines on photo-as https://ui.adsabs.harvard.edu/abs/2022arXiv220802794N/abstract
There's been some other work about including emission lines more flexibly, including produced by AGN-like ionizing continuua, in https://ui.adsabs.harvard.edu/abs/2024arXiv240504598L/abstract https://ui.adsabs.harvard.edu/abs/2024ApJ...969L...5L/abstract
@paulmartini, @yzfkieran817 @wei-jian
Below is the VW01 template (originally provided by @paulmartini), resampled to constant-velocity pixels and then smoothed by a Gaussian kernel ranging in sigma between 150 and 10,000 km/s.
A couple questions:
Hi @moustakas -- Here are my answers (and @yzfkieran817's):
This purpose of this ticket is to organize discussion on the desired contents of a QSO VAC, which would be constructed from a QSO catalog.
The desired contents are emission line and continuum measurements that could be useful for the LyA, GQC, and GQP WGs. For example, we should be able to improve quasar redshifts with information about BH mass and luminosity relative to Eddington. This in turn requires measurements of continuum luminosities, parametric and non-parametric line shape information, removal of FeII emission, and power-law fits to the continuum.
The parent QSO catalog would be the official QSO catalog produced by the LSS WG. Therefore this does not include low-luminosity AGN in galaxies whose continua include significant stellar emission. The rationale is they have quite secure redshifts from the stellar features, and they are usually sufficiently low redshift and luminosity that issues like Fe template subtraction do not apply.
Here are some quantities to start a discussion of the contents:
Some references on Fe templates are:
For reference, here is a figure with a fit to the Fe emission around Mg II from a recent paper led by Zhefu Yu:
Caption: The four colored lines are a PL continuum and the three Fe templates referenced above. Also okay to ignore the single-epoch spectrum.