Shifting, Interpolating, Integrating SED Data in Iris
IRIS 2.0 has added several science capabilities built upon the powerful SED building and modeling already implemented in the previous versions of the tool, making IRIS even more useful for astronomers. In this thread, we will discuss how IRIS can interpolate and smooth SEDs, shift SEDs to any cosmological redshift, and calculate the flux through a given interval or a photometric filter.
Last Update: 05 Aug 2013 - Updated to illustrate shifting, interpolation and integration features only. Previous version moved to "Iris 2.0 Features Tour".
- Interpolation and Smoothing
- Red- and Blue-shifting
- Integration of SEDs
For this thread, we use the SED of Arp 220, which is the closest ultra-luminous infrared galaxy (redshift z=0.018). The photometric data is imported from the NED SED service via the Load NED SED desktop icon. The more useful visualization of the SED of this source can be obtained in the Viewer window switching from flux density to flux in the "Units" field (here we use erg/s/cm2 against microns). In this visualization, we can clearly see that the infrared radiation dominates the SED.
Now that we have our SED loaded, we can explore the Science applications in IRIS. In the main IRIS window, click on the "Shift, Interpolate, Integrate" icon to open the Science window. The left column shows the list of SEDs available, while on the right the first tab, named "Redshift and Interpolation", allows users to move the SED of a source in redshift (upper panel), or produce new SEDs by interpolating the SED selected in the list on the left, using different techniques (lower panel). The second tab, called "Calculate Flux", allows users to evaluate integrated photometric quantities.
After checking that Arp 220 SED is selected in the left, we can choose the method for the interpolation, the interval limits, and the units of the final interpolated SEDs. In this case, the default values of the parameters ("Linear Spline" as the interpolation method, "Angstrom" as Units, and 1000 as the number of bins in the final interpolated SED) are good. In order to produce a smoothed version of the interpolated SED that would minimize the effect of any small scale scatter, we can also tick the "Smooth" box and choose the size of the box (expressed as a number of data points). Then, clicking on "Create New SED" will create the smoothed and interpolated version of the Arp 220 SED. The new interpolated SED will appear in the "open SED" frame in the "Science" "SED Builder" windows, and will automatically be displayed in the Iris Visualizer. Just for the sake of comparison, we overplot the original data over the interpolated SED using the "co-plot" function located under "Display" in the Iris Visualizer tool bar.
Iris is capable of shifting SEDs to and from any cosmological redshift under the "Redshift and Interpolation" tab in the "Science" window. Simply input the initial redshift of the SED in the field labeled "Initial redshift", input the new redshift into the "Move to redshift" field, and click "Create New SED". A new SED will appear in the list of SEDs on the "Open SEDs" box of the "Science" and SED Builder windows, containing the redshifted photometric points. The newly created SED will automatically appear on the Iris Visualizer.
Let's say that now we would like to understand how the SED of Arp 220 would look if it was located at a much higher redshift than the actual redshift (z=0.018), for example at z=0.5. We will enter 0.018 in the "Initial redshift" field and 0.5 in the "Move to redshift" field. The resulting change in the shape of the SED can be inspected qualitatively by coplotting the original SED and the redshifted SED at z=0.5 using the co-plot function (in Iris Visualizer tool bar, under "Display". See more about co-plotting in Visualizing SED Data in Iris).
We might also be interested in how the SED of this source would look like if it was in the vicinity of the Milky Way, say at redshift z=0. We can reconstruct the "local" version of the Arp 220 SED just setting to 0 the "Move to redshift" field and clicking again on "Create new SED". For Arp 220, there isn't much of a difference between the z=0 SED and the original SED. However, this method is useful for shifting high-redshift sources to rest frame.
Note: Any metadata associated with the original SED will be lost. This means if you shift a multi-segment SED, the result will be a one-segment SED; Only information from the Data tab in the Metadata Browser will remain.
Note: setting either one of the initial or final redshift values to a negative number will trigger an "Invalid redshift values" error message, since IRIS only takes into account the effect of the cosmological redshift.
Another useful science capability that has been added to IRIS 2.0 is the ability to evaluate integrated fluxes from the SEDs in spectral intervals that can be either defined by the users or selected from a repository of existing photometric systems. To calculate a meaningful flux, there must be enough photometric points in the spectral range we wish to integrate over. Therefore, it is a safe-bet to use an interpolated SED when calculating the flux, provided that the number of bins is large enough to densely sample the SED. We can calculate the flux directly from the original SED, but by using the interpolated SED, we don't have to worry about having enough data points in the spectral range we are integrating over.
We showcase the highlights of this capability using again the interpolated and smoothed SED of Arp 220 that we have produced previously.
In the "Science" window, select the "Calculate Flux" tab. Under the section "Add Passband", Iris provides two methods of calculating the flux: a simple user-defined passband or the use of pre-loaded photometry filters. This section describes the "Passband" option.
Iris can integrate the flux of an SED in a user-defined passband. The interval of integration can be over wavelength (in Angstroms), frequency (in Hz), or in energy (keV). The user selects the "Passband" option, enters the lower and upper limits in the empty fields adjacent to the "Passband" option, and selects the desired units in the drop-down menu on the right. By clicking on "Add", the passband, effective wavelength (in Angstroms) and the integrated flux (in Janskys) will automatically be added to the list in the "Results" table in the bottom of the window. Note that these values can be converted into any other units supported by Iris (i.e. the units listed in the Iris Visualizer "Units" browser).
Let's calculate some fluxes for our smoothed and interpolated Arp220 SED. Since Arp220 is a infrared luminous source, let's calculate the observed infrared flux on the interval [1, 1000] microns (or [1e4,1e7] Angstroms). Select "Passband" and enter the lower (1e4) and upper (1e7) limits in the blank fields, like in the figure below. In the drop-down menu on the right, select "Angstroms" as the unit of measure. Now, by clicking on "Add", the passband we have defined will be automatically added to the "Results" list (in this case, it will be the first), together with the calculated flux value and the effective wavelength measured in Angstrom. We can repeat the operation with any different spectral interval in any units, for example [1E9, 1E10] Hz, and the resulting integrated flux will appear in the "Results" panel after clicking on "Add".
The Photometry Filter Selector allows the user to easily calculate the flux of an object as it would be measured by a real photometric system used at a given observatory. It contains over 2100 photometric transmission curves provided by the Spanish Virtual Observatory, which span the full electromagnetic spectrum. You can open the Photometry Filter Selector from the "Calculate Flux" tab by selecting "Photometry filter" and clicking "Choose". You can also open this window from the Iris desktop toolbar, under Tools -->SED Builder --> Photometry Filter Browser.
In the browser, the catalog of facilities/observatories is on the left, and filter information appears on the right (see figure below). The facilities and corresponding filters are arranged in alphabetical order. Clicking on any of the arrows on the left of each name will toggle down the available filters for that observatory. When you select a filter, the fields on the right side of the window will fill up with basic information about the filter. There is a search "By String" bar in the top-right corner of the window, into which a facility or instrument name is typed to fetch a list of available filters. The user is allowed to select multiple filters at a time by holding down the Shift or Ctrl ("Command" on Mac) keys and clicking on the filter name. When all filters are selected, clicking the "Done" button at the bottom of the window will automatically close the Photometry Filter Selector, and calculate and display the integrated fluxes in the "Results" section in the "Science" window. Other previously calculated fluxes will remain in the list. Below is an example of the photometry browser, after selecting the WISE filters (which we will do in the example below).
For example, let's say that we want to evaluate what integrated fluxes from Arp 220 we would expect to observe in the four WISE bands (3.4, 4.6, 12 and 22 microns) from the WISE telescope. We open the Photometry Filter Selector by selecting the "Photometry filter" option in the "Calculate flux" tab of the "Science" window and clicking on "Choose". Now, we browse the list for the WISE filters, then click on the arrow next to the WISE folder to reveal the available filters. We select all four WISE filters by holding down the Shift key and clicking the top and bottom filters. We see the basic filter information appear in the fields on the right. Click "Done" to calculate the integrated fluxes. Now we see the flux from the WISE filters (see figure below), as well as the fluxes from our user-defined passbands. All the flux values can be saved to a ASCII file by clicking on "Save" (see next section), and can be extracted as a separate SED by clicking "Create SED" on the bottom of the window.
You can save the calculated fluxes into an ASCII file for external use by clicking the "Save" button on the bottom of the "Calculate Flux" tab. You have the option to save any or all of the calculated fluxes in ASCII format by highlighting the ones you wish to save; clicking "Save" without highlighting any of the Results will save all fluxes. A simple interface allow users to pick the location, name of the file, and the units for the effective wavelengths of the spectral intervals/photometric systems and fluxes, respectively, using the X and Y drop-down menus. The following is the file output for our Arp220 data, with the effective wavelength in microns and the flux in erg/s/cm2.
# Iris Flux Integration output # This file was generated by Iris, the VAO SED building and analysis tool # # # File created on Wed Jul 17 16:15:13 EDT 2013# # # Effective Wavelength in micron # Flux in erg/s/cm2 # # EffectiveWavelength Flux 3.315655859375 4.674049594152648E-11 4.56449921875 4.0619123056081297E-11 10.78684453125 2.227278813788558E-10 21.9149640625 9.096264184236841E-10 500.5 0.0013274507662364892 164885.875 1.3965920014327272E-5
|02 Jun 2013||First issue|
|05 Aug 2013||Updated to illustrate shifting, interpolation and integration features only. Previous version moved to "Iris 2.0 Features Tour".|