A spatially resolved study of ionized regions in galaxies at different ...

A spatially resolved study of ionized regions in galaxies at different ... A spatially resolved study of ionized regions in galaxies at different ...

cab.inta.csic.es
from cab.inta.csic.es More from this publisher
02.01.2015 Views

58 3 • IFS of a GEHR in NGC 6946 the spectrophotometric standard star frames. All the corrections performed to the images, mentioned in precedent sections, do not correct for the differences in the shape of the emission lines along the cross-dispersion axis. This effect is also present in sky-subtracted spectra from long-slit spectroscopy, showing characteristic residuals in the location of the sky emission lines. This effect is more severe in fiber-based spectrographs, since the effective dispersion and shape of the emission lines varies fiber-to-fiber. As mentioned in Section 3.2, PPak includes additional fibers that probe the sky far enough from the science FOV to avoid contamination by the astronomical object. One simple method to perform the sky subtraction is to extract all the 36 sky fibers, create a median sky spectra and substract it from the science exposures. Nevertheless, this does not take into account the distortion in the wavelength solution along the cross-dispersion. To increase the accuracy of the sky substraction, a second method can be applied. The routine create sky ppak determines the sky spectra corresponding to any science fiber by an interpolation of the spectra obtained through the sky lines. As the set of Hii regions in our images presents an inhomogeneous distribution through the field of view of the IFU, we followed another method to derive a representation of the sky. The routine create sky clip creates a sky spectra by obtaining the median between a certain number of adjacent spectra, clipping those ones with a flux over a certain threshold of the standard deviation. The results from the last two methods are very close from each other in our case, with differences less than 3%. The use of one or other method depends on the nature of the science observations, so it is advisable to check both solutions. Both results create a better representation of the sky than the simple method, with less deep residuals, but still present. These residuals are specially strong in the red part of the spectrum (8000-11000 Å). As seen in Section 3.1, the Hii complex is located at the extreme of the NE arm of the galaxy. As the sky fibers surround the object, some of them might have some diffuse nebular emission from the end of the arm contaminating the sky spectra. We compared the flux of the Hα line in representative spectra with and without sky subtraction. The difference in the flux was negligible 8 . This reliable, free of nebular emission, sky spectra was subtracted from the science frames. 3.3.8 Mosaics As explained before, adjacent spectra in a RSS file may not correspond to spatially near positions, so it is important to reorder the spectra to their original position in the sky. A position table is required, relating the spectra to their locations in the sky. Then, it is possible to create a regular gridded data cube by interpolating the data spatially, and reconstruct the original image of the target at any wavelength. This step is needed if DAR correction is 8 The Hα contribution of the sky spectra is about 10% for the lowest value in the FOV.

3.3. Data Reduction 59 compulsory. If it is not the case, this step can be avoided and work with the RSS file and its corresponding position table. In the case of PMAS in the PPak mode it is necessary to interpolate at each wavelength the intensity at each spatial location in order to create a data cube in an artificial squared grid pattern. This can be done by using E3D or its complementary routines. Nevertheless, when performing dithering exposures to fill the holes between adjacent fibers to get a filling factor of 1, the procedure is slightly different and interpolation can be avoided. First, all the dithering exposures must be reassembled in a single RSS file. In the case of covering all the FOV, the exposures have to have certain overlap and the offsets between adjacent exposures had to be selected in a way that the overlap is complete between spectra. We used the routine Mosaic rss overlap of R3D to build a final RSS for each night which covered the entire FOV. As mentioned before, it is possible to construct the datacube without interpolating the data. This is accomplished by means of the regularization scheme, in which all the spectra result from the combination of measured values. A grid of a certain pixel scale, smaller than that of the original fibers 9 is constructed over the FOV, and for each pixel the spectrum that enters in it is averaged, obtaining a final datacube. In the second night, when getting the exposures in the 7000-10100 Å range (GROT = -72), the de-rotator and guiding stopped at the second pointing of the dithering exposures 10 , loosing the pointing for the second and the third dithering exposures. This prevents us to perform a mosaic covering all the FOV. We tried to recover the positions of the last dithering, but no good results where achieved in the last mosaic, so we ended with only the first dithering for the red part of the spectrum. Nevertheless, this is not important if only ratios are involved and absolute flux from emission lines are not needed. In section 3.4.8, we will see that this problem will not affect our results when dealing with integrated properties. The mosaic in the blue range was built without any problems. The final RSS blue mosaic contained 993 spectra (3 × 331), while the red RSS frame 331, since it is only the pointing for the first dithering exposure. 3.3.9 Correction for atmospheric absorption In some cases, the absorption of the atmosphere is so strong that it can depress some emission lines, altering fluxes and ratios needed to obtain the physical properties of the target. This is specially important in the red range of the spectrum, namely from 8000 Å. We have followed the procedure described by Diaz et al. (1985) to remove the atmospheric water-vapour absorption bands in order to determine the strengths of the sulphur lines 9 A pixel size of 1/3 of the original spaxel size is recommended. This value should also be used when interpolating to create datacubes when there are no dithering exposures or building maps for emission lines. In this case, a nearest neighbour interpolation is also recommended. 10 It has to be taken into account that the dithering method in PPak is experimental, so the telescope does not record in the header any information of the offsets. These values have to be kept in the log of observations in order to build the mosaic.

3.3. D<strong>at</strong>a Reduction 59<br />

compulsory. If it is not the case, this step can be avoided and work with the RSS file and its<br />

correspond<strong>in</strong>g position table.<br />

In the case <strong>of</strong> PMAS <strong>in</strong> the PPak mode it is necessary to <strong>in</strong>terpol<strong>at</strong>e <strong>at</strong> each wavelength<br />

the <strong>in</strong>tensity <strong>at</strong> each sp<strong>at</strong>ial loc<strong>at</strong>ion <strong>in</strong> order to cre<strong>at</strong>e a d<strong>at</strong>a cube <strong>in</strong> an artificial squared<br />

grid p<strong>at</strong>tern. This can be done by us<strong>in</strong>g E3D or its complementary rout<strong>in</strong>es. Nevertheless,<br />

when perform<strong>in</strong>g dither<strong>in</strong>g exposures to fill the holes between adjacent fibers to get a fill<strong>in</strong>g<br />

factor <strong>of</strong> 1, the procedure is slightly <strong>different</strong> and <strong>in</strong>terpol<strong>at</strong>ion can be avoided. First, all the<br />

dither<strong>in</strong>g exposures must be reassembled <strong>in</strong> a s<strong>in</strong>gle RSS file. In the case <strong>of</strong> cover<strong>in</strong>g all the<br />

FOV, the exposures have to have certa<strong>in</strong> overlap and the <strong>of</strong>fsets between adjacent exposures<br />

had to be selected <strong>in</strong> a way th<strong>at</strong> the overlap is complete between spectra. We used the rout<strong>in</strong>e<br />

Mosaic rss overlap <strong>of</strong> R3D to build a f<strong>in</strong>al RSS for each night which covered the entire<br />

FOV. As mentioned before, it is possible to construct the d<strong>at</strong>acube without <strong>in</strong>terpol<strong>at</strong><strong>in</strong>g the<br />

d<strong>at</strong>a. This is accomplished by means <strong>of</strong> the regulariz<strong>at</strong>ion scheme, <strong>in</strong> which all the spectra<br />

result from the comb<strong>in</strong><strong>at</strong>ion <strong>of</strong> measured values. A grid <strong>of</strong> a certa<strong>in</strong> pixel scale, smaller than<br />

th<strong>at</strong> <strong>of</strong> the orig<strong>in</strong>al fibers 9 is constructed over the FOV, and for each pixel the spectrum th<strong>at</strong><br />

enters <strong>in</strong> it is averaged, obta<strong>in</strong><strong>in</strong>g a f<strong>in</strong>al d<strong>at</strong>acube.<br />

In the second night, when gett<strong>in</strong>g the exposures <strong>in</strong> the 7000-10100 Å range (GROT =<br />

-72), the de-rot<strong>at</strong>or and guid<strong>in</strong>g stopped <strong>at</strong> the second po<strong>in</strong>t<strong>in</strong>g <strong>of</strong> the dither<strong>in</strong>g exposures 10 ,<br />

loos<strong>in</strong>g the po<strong>in</strong>t<strong>in</strong>g for the second and the third dither<strong>in</strong>g exposures. This prevents us<br />

to perform a mosaic cover<strong>in</strong>g all the FOV. We tried to recover the positions <strong>of</strong> the last<br />

dither<strong>in</strong>g, but no good results where achieved <strong>in</strong> the last mosaic, so we ended with only the<br />

first dither<strong>in</strong>g for the red part <strong>of</strong> the spectrum. Nevertheless, this is not important if only<br />

r<strong>at</strong>ios are <strong>in</strong>volved and absolute flux from emission l<strong>in</strong>es are not needed. In section 3.4.8, we<br />

will see th<strong>at</strong> this problem will not affect our results when deal<strong>in</strong>g with <strong>in</strong>tegr<strong>at</strong>ed properties.<br />

The mosaic <strong>in</strong> the blue range was built without any problems. The f<strong>in</strong>al RSS blue mosaic<br />

conta<strong>in</strong>ed 993 spectra (3 × 331), while the red RSS frame 331, s<strong>in</strong>ce it is only the po<strong>in</strong>t<strong>in</strong>g<br />

for the first dither<strong>in</strong>g exposure.<br />

3.3.9 Correction for <strong>at</strong>mospheric absorption<br />

In some cases, the absorption <strong>of</strong> the <strong>at</strong>mosphere is so strong th<strong>at</strong> it can depress some<br />

emission l<strong>in</strong>es, alter<strong>in</strong>g fluxes and r<strong>at</strong>ios needed to obta<strong>in</strong> the physical properties <strong>of</strong> the<br />

target. This is specially important <strong>in</strong> the red range <strong>of</strong> the spectrum, namely from 8000 Å.<br />

We have followed the procedure described by Diaz et al. (1985) to remove the <strong>at</strong>mospheric<br />

w<strong>at</strong>er-vapour absorption bands <strong>in</strong> order to determ<strong>in</strong>e the strengths <strong>of</strong> the sulphur l<strong>in</strong>es<br />

9 A pixel size <strong>of</strong> 1/3 <strong>of</strong> the orig<strong>in</strong>al spaxel size is recommended. This value should also be used when<br />

<strong>in</strong>terpol<strong>at</strong><strong>in</strong>g to cre<strong>at</strong>e d<strong>at</strong>acubes when there are no dither<strong>in</strong>g exposures or build<strong>in</strong>g maps for emission l<strong>in</strong>es.<br />

In this case, a nearest neighbour <strong>in</strong>terpol<strong>at</strong>ion is also recommended.<br />

10 It has to be taken <strong>in</strong>to account th<strong>at</strong> the dither<strong>in</strong>g method <strong>in</strong> PPak is experimental, so the telescope does<br />

not record <strong>in</strong> the header any <strong>in</strong>form<strong>at</strong>ion <strong>of</strong> the <strong>of</strong>fsets. These values have to be kept <strong>in</strong> the log <strong>of</strong> observ<strong>at</strong>ions<br />

<strong>in</strong> order to build the mosaic.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!