iSpec is an integrated spectroscopic software framework suitable for the creation of spectral libraries such as the Benchmark Stars library and the determination of atmospheric parameters (i.e. effective temperature, surface gravity, metallicity) and individual chemical abundances.
The framework can be used in automatic massive analysis, but it also includes a user-friendly visual interface that can easily interoperate with other astronomical applications such as TOPCAT, VOSpec and splat which provide access to the Virtual Observatory.
iSpec: An integrated software framework for the analysis of stellar spectra]iSpec: An integrated software framework for the analysis of stellar spectra
S. Blanco-Cuaresma et al.] S. Blanco-Cuaresma,
P. Jofré and
Univ. Bordeaux, LAB, UMR 5804, F-33270, Floirac, France.
CNRS, LAB, UMR 5804, F-33270, Floirac, France.
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, U.K.
Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden.
Received — ; accepted —
Keywords : spectral libraries – spectral analysis – spectroscopic surveys
On-going high resolution spectroscopic surveys such as the Gaia ESO survey (GES, ) and the future HERMES/GALAH () provide an enormous quantity of high quality spectra. This represents a unique opportunity to unravel the history of our Galaxy by studying the chemical signatures of large samples of stars.
This huge quantity of spectroscopic data challenges us to develop automatic processes to perform the required analysis. Numerous automatic methods have been developed over the past years (e.g ; ; ; ; ; ; , to name a few) to treat the spectra and derive atmospheric parameters from large datasets, where each of them uses different reference models and minimization strategies. A consequence of this diversity is the inhomogeneity of the published measurements. This can be seen when we compare atmospheric parameters from different sources in the bibliographical compilation PASTEL ().
To improve the understanding of the systematics between the various characterization methods, we have assembled a library of the 34 Gaia FGK Benchmark Stars, covering different regions of the HR diagram and spanning a wide range in metallicity (see in this proceedings). Their temperature and surface gravity were determined from fundamental relations and [Fe/H] from a dedicated spectroscopic analysis ().
In Sect. [, we describe the iSpec software framework111We plan to freely distribute iSpec under the Open Source terms of the GNU Affero General Public License. that was used to prepare the spectra of the library (Blanco-Cuaresma et al, submitted). In Sect. [, we show an application of iSpec to chemical tagging by using its capabilities for the determination of atmospheric parameters and abundances (Blanco-Cuaresma et al, in preparation).
The main functionalities for spectra treatment and library creation that are integrated into iSpec are the following:
Cosmic rays removal: Spectra may contain residuals from cosmics rays which are automatically detected by using a median filter to smooth out single flux deviations.
Continuum normalization: The continuum points of a spectrum are found by applying a median and maximum filter with different window sizes. Afterwards, a polynomial or group of splines can be use to model the continuum and finally normalize the spectrum by dividing all the fluxes by the model.
Resolution degradation: The spectral resolution can be degraded by convolving the fluxes with a Gaussian of a given FWHM (km/s).
Additionally, iSpec is capable of determining atmospheric parameters (i.e effective temperature, surface gravity, metallicity, micro/macroturbulence, rotation) and individual chemical abundances. Two different techniques can be applied:
Synthetic spectra fitting technique: iSpec compares an observed spectrum with synthetic ones generated on the fly, in a similar way as Spectroscopy Made Easy (SME) tool does (). A least-squares algorithm (Levenberg-Marquardt implementation) minimizes the difference between the synthetic and the observed spectra in order to converge towards the closest synthetic spectrum. In each iteration, the algorithm will vary each one of the free parameters at a time and prognosticate in what direction it should move to find the solution. Specific regions of the spectrum can be selected to minimize the computation time and focus onto the most relevant parts to better identify stars (i.e. wings of H-/MgI triplet, Fe I/II lines).
Equivalent widths method: Gaussian models are fitted for a given list of Fe I/II lines, from their integrated area we can derive their respective equivalent width and thus their abundances. The algorithm for determining the atmospheric parameters is a least-squares technique (as explained above) based on the assumption of excitation and ionization equilibrium, similar to GALA () and FAMA ().
“SPECTRUM: a Stellar Spectral Synthesis Program” () is used to generate synthetic spectra and derive abundances from equivalent widths. iSpec offers an wide range of options such as different atomic line lists, model atmospheres (i.e. ATLAS or MARCS ) and solar abundances.
It is worth noting that the framework can be used in automatic massive analysis through python scripts, but it also includes a user-friendly visual interface that can easily interoperate with other astronomical applications such as TOPCAT222http://www.star.bris.ac.uk/~mbt/topcat/, VOSpec333http://www.sciops.esa.int/index.php?project=ESAVO&page=vospec and splat444http://star-www.dur.ac.uk/~pdraper/splat/splat.html facilitating a indirect way to access the Virtual Observatory555http://www.ivoa.net/
The derivation of atmospheric parameters has been tested using the Gaia FGK Benchmark Stars library and associated parameters (see Sect. [). The results are very consistent as shown in Fig. 1, 2 and 3.
Stars born from the same giant molecular cloud are expected to share the same homogeneous chemical patterns as shown by . These findings open the possibility of using the technique of chemical tagging to identify common formation sites in the disk as proposed by . This requires to apply automatic and massive spectral analysis to handle huge quantities of spectra.
We have tested and validated the performance of iSpec for this purpose by automatically analysing three different open clusters observed by the NARVAL spectropolarimeter mounted on the 2m Telescope Bernard Lyot () (located at Pic du Midi, France). After deriving their atmospheric parameters and the individual abundances for a list of given elements, we have shown that stars belonging to the cluster share a common chemical pattern, with a very low dispersion for certain elements (below 0.05 dex). Additionally, the synthetic spectra fitting technique and equivalent width method are able to recover similar chemical patterns as shown in Fig. 4, 5 and 6. A detailed discussion of the abundances of individual elements will be presented in a forthcoming paper.
iSpec has been proven to be a suitable framework for the creation of spectral libraries such as the Gaia FGK Benchmark Stars (Blanco-Cuaresma et al., submitted), together with the derivation of atmospheric parameters and individual chemical abundances (Blanco-Cuaresma et al., in preparation). The code will be released on http://www.blancocuaresma.com
This work was partially supported by the Gaia Research for European Astronomy Training (GREAT-ITN) Marie Curie network, funded through the European Union 7th Framework Programme [FP7/2007-2013] under grant agreement n. 264895. All the software used in the data analysis has been provided by the Open Source community.
- Allende Prieto, C., 2007, AJ, 134, 1843-1848
- Aurière, M., 2003, EAS Publications Series, 9, 105
- Bertaux, J. L., Lallement, R., Ferron, S., Boone, C., Bodichon, R., 2013, arXiv eprint
- Block, D. L., Freeman, K. C., Puerari, I., 2010, Galaxies and their Masks, 319
- De Silva, G. M., Freeman, K. C., Asplund, M., Bland-Hawthorn, J., Bessell, M. S., Collet, R., 2007, AJ, 133, 1161-1175
- Freeman, K., Bland-Hawthorn, J., 2002, ARAA, 40, 487-537
- Gilmore, G., Randich, S., Asplund, M., Binney, J., Bonifacio, P., Drew, J., Feltzing, S., Ferguson, A., Jeffries, R., Micela, G., Negueruela, I., Prusti, T., Rix, H.-W., Vallenari, A., Alfaro, E., Allende-Prieto, C., Babusiaux, C., Bensby, T., Blomme, R., Bragaglia, A., Flaccomio, E., François, P., Irwin, M., Koposov, S., Korn, A., Lanzafame, A., Pancino, E., Paunzen, E., Recio-Blanco, A., Sacco, G., Smiljanic, R., Van Eck, S., Walton, N., 2012, The messenger, 147, 25-31
- Gray, R. O., Corbally, C. J., 1994, AJ, 107, 742-746
- Gustafsson, B., Edvardsson, B., Eriksson, K., Jørgensen, U. G., Nordlund, Å., Plez, B., 2008, AAP, 486, 951-970
- Jofré, P., Panter, B., Hansen, C. J., Weiss, A., 2010, A&A, 517
- Jofré, P., Heiter, U., Soubiran, C., Blanco-Cuaresma, S., Pancino, E., Bergemann, M., Cantat-Gaudin, T., Gonzalez Hernandez, J. I., Hill, V., Lardo, C., de Laverny, P., Lind, K., Magrini, L., Masseron, T., Montes, D., Mucciarelli, A., Nordlander, T., Recio-Blanco, A., Sobeck, J., Sordo, R., Sousa, S. G., Tabernero, H., Vallenari, A., Van Eck, S., Worley, C. C., 2013, arXiv eprint
- Jofré, P., Heiter, U., Blanco-Cuaresma, S., Soubiran, C., 2013, IWSSL Proceedings
- Katz, D., Soubiran, C., Cayrel, R., Adda, M. Cautain, R., 1998, AAP, 338, 151-160
- Koleva, M., Prugniel, P., Bouchard, A., Wu, Y., 2009, A&A, 501, 1269-1279
- Kurucz, R. L., 2005, Memorie della Societa Astronomica Italiana Supplementi, 8, 14
- Magrini, L., Randich, S., Friel, E., Spina, L., Jacobson, H., Cantat-Gaudin, T., Donati, P., Baglioni, R., Maiorca, E., Bragaglia, A., Sordo, R., Vallenari, A., 2013, AAP, 558, 38
- Mucciarelli, A., Pancino, E., Lovisi, L., Ferraro, F. R., Lapenna, E., 2010, APJ, 766, 78
- Recio-Blanco, A., Bijaoui, A., de Laverny, P., 2006, MNRAS, 370, 141-150
- Soubiran, C., Le Campion, J.-F., Cayrel de Strobel, G., Caillo, A., 2010, A&A, 515, 111
- Valenti, J. A, Piskunov, N., 1996, AAPS, 118, 595-603