12539( 3) - 06-Feb-2012 18:03:33 - [ 1] HUBBLE SPACE TELESCOPE OBSERVING PROGRAM 12539 Version: 3 Check-in Time: 06-Feb-2012 18:03:33 STScI Edit Number: 1 Title A novel approach to find Lyman continuum leaking galaxies at z~0.3 with COS ------------------------------------------------------------------------------------ Type Cycle GO 19 ------------------------------------------------------------------------------------ Investigators Contact? PI: Prof. Nils Bergvall Uppsala Astronomical Observatory CoI: Ms. Elisabet Leitet Uppsala Astronomical Observatory Y CoI: Mr. Thomas Marquart Uppsala Astronomical Observatory N CoI: Dr. Matthew Hayes Observatoire Midi-Pyrenees N CoI: Prof. Goeran Oestlin Stockholm University N CoI: Dr. Erik Zackrisson Stockholm University N ------------------------------------------------------------------------------------ Abstract What drives the early phase of cosmic reionization is poorly known. Numerous low to intermediate-mass star forming galaxies are the most promising candidates, but the problem is to understand how the ionizing radiation can leak out into the intergalactic medium. Models indeed predict the escape of LyC radiation from dwarf starburst galaxies under certain conditions, but it has proven to be observationally very difficult to confirm this in the local galaxy sample. Except in very few reported cases the envelope of neutral gas seems to halt all Lyman continuum photons. To better understand this discrepancy, we here propose to observe with COS a small sample of local starburst galaxies selected from the SDSS for their very blue colours, young ages and weak H-alpha emission. This combination can, according to models, occur only when a large fraction of the ionizing radiation escapes from the galaxy. The redshifts of our four targets is z~0.3, which puts the LyC in a spectral region where the sensitivity of COS/G130M is high. Since it might be counterintuitive to select strong starbursts by weak H-alpha emission, it can be understood that all earlier searches for LyC emission at low redshifts have had a high rate of failure. Selection techniques utilizing star formation indicators like strong optical emission lines will bias the samples towards galaxies with a high neutral gas-fraction, and thereby against LyC escape. With the novel approach for target selection in this proposal, we are instead sure to optimize the success rate. ------------------------------------------------------------------------------------ Observations Description ------------------------ The scientific goal with this project is to search for escape of radiation in the Lyman continuum (LyC) of a sample of 4 galaxies with predicted high LyC escape fractions. For this purpose we will use the G130M grism of the COS FUV channel in TIME-TAG mode to directly measure the signal in the LyC. The redshifts of the targets span a range of z=0.243-0.353. At these redshifts the Ly-limit will fall where the throughput of COS is high, and is yet low enough that the data from SDSS allow us to properly model the galaxies. The central 12539( 3) - 06-Feb-2012 18:03:33 - [ 2] wavelength settings of G130M depend on the redshift of the galaxy, and have been chosen to fit the position of the Lyman break: Galaxy Lyman limit (A) Central wave (A) Visibility (min) f900 S/N J143004.07+022213.3 1233 1291 54 2.3E-17 11 J143832.04+390333.7 1133 1055 55 2.3E-17 7 J080250.09+190119.1 1193 1291 54 2.1E-17 8 J154459.67+115017.4 1134 1055 54 2.8E-17 9 In the table above are also shown the visibilities, estimated minimum flux densities in the 900 A continnum of the target galaxies and the corresponding S/N. Flux densities are given in units of erg/cm^2/s/A. The size of the galaxies, obtained from the SDSS g band data, were found to match the 2.5 arcsec COS aperture. NUV acquisition images from previous COS programes have however shown that the star-forming regions in UV are much more compact (GO 11522, 11727), thus we expect a resolution in the FUV continuum better than R~1000 (at 1150 A for G130M, 1055). This is well within the requirements for this project where we only need to resolve the Lyman break. The LyC flux is however assumed to escape from even smaller regions through channels in the ISM, and can likely be considered as a point source with better resolution for these wavelengths. Flux density predictions -------------------------------- Two of the galaxies have been observed by GALEX. The ratio between the observed GALEX fuv and SDSS g fluxes, after reddening corrections, was almost the same in the two cases. Therefore, in the two cases were no GALEX data existed, the fuv/g ratio was assumed to be the same. The stability of this ratio as long as the starburst is active is ensured also from the modelling. The flux density shortwards of the Lyman limit was derived from the observed f(lambda)_fuv (lambda ~ 1175 A in the z=0.3 frame) or the value corresponding to the g magnitude. The Starburst99 (Leitherer et al. 1999, ApjS 123, 3) prediction under 100% leakage conditions is f(lambda)_fuv/f(lambda)_LyC ~1.9 while our models (Zackrisson et al. 2001, A&A 375, 814) gives a lower value (~1.3). In our exposure time (ET) estimates we use the more conservative value (1.9). Dust extinction has been found to account for ~50% attenuation of LyC photons in starburst galaxies. Therefore we have used a conservative relative escape fraction of f_esc_rel=10%, even though our models predict much more. The predicted 900 A fluxes can be seen in the table above. S/N calculations -------------------------- The Lyman limit of our targets is redshifted to 1133-1233 A. In this spectral range the throughput of the FUV XDL detector changes dramatically with wavelength, and we have therefore limited our LyC S/N calculations to ~30 A baselines although the LyC flux is expected to be significant down to much shorter wavelengths. In the ET estimates we assume that the LyC is constant over the whole baseline, an assumption based on the coincidence of a rise towards shorter wavelengths for model spectra of young burst populations and 12539( 3) - 06-Feb-2012 18:03:33 - [ 3] that of a similar rise towards shorter wavelengths for the extinction curve. The spectra of the two single detections in Shapley et al. (2006, ApJ 651, 688) also seem to confirm this. For one of our targets, J143004.07+022213.3, the airglow Ly alpha emission line falls in the baseline, decreasing it further with 2 A. We chose to use the FP-POS strategy for the science exposures, using all 4 positions, 2 for each orbit. The total exposure times on target (5178-5304s) were derived from the alloted 2 orbits with visibility 54 (55) minutes each, minus total overhead (999-1025s). We used the online COS ETC to calculate the S/N in the LyC with settings point source, flat continuum normalized to the LyC flux given in the table, with zodiacal background, airglow and earth shine set to AVERAGE. The calculations were performed at the central wavelength of respectively baseline, and the S/N per resolution element was then multiplied by the square root of the total number of resolution elements in the baseline. From this we could derive a S/N of 7-11 (see table). These numbers have been calculated conservative in all steps, so any change in the S/N should be towards the better. All targets have a 3 sigma detection limit of f_esc, rel ~ 3%. Acquisition strategy -------------------------- The coordinates for the targets are obtained from the SDSS, with estimated errors < 0.2 arcsec (J.R. Pier et al. 2003, AJ 125, 1559; Sandor et el. Proceedings of the 8th European VLBI Network Symposium. September 26-29, 2006, Torun, Poland). The errors together with a comparison of the outcome of earlier COS programes (GO 12248), lead us to the conclusion that the ACQ/IMAGE method will give us the pointing precision we need. Integration times for the ACQ/IMAGE exposures were obtained from the COS Imaging Target Acquisition ETC, assuming the point source option and using the GALEX NUV fluxes when available or the corresponding model magnitudes. The GALEX NUV corresponds well to the peak sensitivity wavelength of the COS NUV detector at ~2300A. As an extra safety measure, to account for slight deviations from a pure point source structure, we added 1 magnitude to the NUV magnitudes before loading the ETC. A similar strategy has successfully been used by other teams when observing similar objects (e.g. GO 12248). We reach S/N=40 for exposure times 52-78 s. BUFFER-TIME estimates -------------------------- The BUFFER-TIME were in all FUV/G130M exposures found to be considerable longer than the exposure times, hence the BUFFER-TIME was set to the exposure time minus 100 s in all cases to minimize overheads (according to COS instrument handbook). Since the exposures were split in four (one for each FP-POS), we have given the value for the shortest exposure in each visit. Bright target notes -------------------------- The flux densities of the targets do not exceed the bright object limit. The conservative case of no extinction was applied to the flux above the Ly-limit of the galaxies, and the target Ly-alpha line does not fall in the wavelength window of the G130M setting for any of the targets. We found no violation of neither the global nor the local count rates in the FUV, and for NUV imaging 12539( 3) - 06-Feb-2012 18:03:33 - [ 4] acquisition the count rate was well below the limit at 360 cts/s in the 9x9 pixel boxes surrounding the targets. The BOT do return one health-and-safety issue for target 1 and its field. The warning is issued from the GSC II, while GALEX issues no warning. The coordinates of the object (assumed to be a **O5V** star by BOT), coincides with our target. Since our target itself is far from exceeding the bright object limit (GALEX FUV=20.1), we see no problem to observe target 1. Real Time Justification ----------------------- Calibration Justification ------------------------- Additional Comments ------------------- ------------------------------------------------------------------------------------ 12539( 3) - 06-Feb-2012 18:03:33 - [ 5] TARGET LIST Fixed Targets ------------------------------------------------------------------------------------------------------------------------------------ Tar| Target | Target | Target |Coord | Radial | Flux data No | Name | Description | Position |Eqnx | Vel. | ------------------------------------------------------------------------------------------------------------------------------------ 1 J143004.07+0 GALAXY, STARBURST RA=14H30M4.0613S +/- 0.2", J2000 V = 18.9 FUV=20.1 22213.3 DEC=+02D22'13.60" +/- 0.2" Reference Frame: ICRS Comments: This object was generated by the targetselector and retrieved from the NED database. The BOT gives a warning for a source whose coordinates coincides with the target. ------------------------------------------------------------------------------------------------------------------------------------ 2 J143832.04+3 GALAXY, STARBURST RA=14H38M32.0424S +/- 0.2", J2000 V = 19.0 90333.7 DEC=+39D03'33.77" +/- 0.2" Reference Frame: ICRS Comments: This object was generated by the targetselector and retrieved from the NED database. ------------------------------------------------------------------------------------------------------------------------------------ 3 J080250.09+1 GALAXY, STARBURST RA=08H02M50.0952S +/- 0.2", J2000 V = 19.0 90119.2 DEC=+19D01'19.20" +/- 0.2" Reference Frame: ICRS Comments: This object was generated by the targetselector and retrieved from the NED database. ------------------------------------------------------------------------------------------------------------------------------------ 4 J154459.66+1 GALAXY, STARBURST RA=15H44M59.6688S +/- 0.2", J2000 V = 18.5 FUV=19.9 15017.5 DEC=+11D50'17.56" +/- 0.2" Reference Frame: ICRS Comments: This object was generated by the targetselector and retrieved from the NED database. ------------------------------------------------------------------------------------------------------------------------------------ 12539( 3) - 06-Feb-2012 18:03:33 - [ 6] Visit: 01 Visit Requirements: On Hold Comments: Additional Comments: Exposures ------------------------------------------------------------------------------------------------------------------------------------ Exp | Target | Instr | Oper. | Aper |Spectral|Central| Optional |Num| Time | Special Num | Name | Config | Mode | or FOV |Element |Waveln.| Parameters |Exp| | Requirements ------------------------------------------------------------------------------------------------------------------------------------ 1 J143004.07+ COS/NUV ACQ/IMA PSA MIRRORA 1 57 S 022213.3 GE ------------------------------------------------------------------------------------------------------------------------------------ 2 J143004.07+ COS/FUV TIME-TA PSA G130M 1291 SEGMENT=BOTH, 1 1205 S 022213.3 G BUFFER-TIME=1085, EXTENDED=NO, see below FP-POS=ALL, FLASH=YES Comments: The BUFFER-TIME= exposure time - 100s differs between the splits. The given BUFFER-TIME is for the shortest exposure. ------------------------------------------------------------------------------------------------------------------------------------ Sub Exposures ------------------------------------------------------------------------------------------------------------------------------------ Target | Exp |Instr | Oper. | Aper |Spectral|Cent.|Primary |Secondary |Iteration |CR-SPLIT |Orbit |Duration Name | Num |Config | Mode | or FOV |Element |Wave.|Pattern Pos |Pattern Pos |Num |Num |Number | ------------------------------------------------------------------------------------------------------------------------------------ J143004.07+022213.3 1 COS/NUV ACQ/IMA PSA MIRRORA none none none none 1 N/A GE J143004.07+022213.3 2 COS/FUV TIME-TA PSA G130M 1291 none none none 1 1 1185.0 S G J143004.07+022213.3 2 COS/FUV TIME-TA PSA G130M 1291 none none none 2 1 1185.0 S G J143004.07+022213.3 2 COS/FUV TIME-TA PSA G130M 1291 none none none 3 2 1410.0 S G J143004.07+022213.3 2 COS/FUV TIME-TA PSA G130M 1291 none none none 4 2 1410.0 S G 12539( 3) - 06-Feb-2012 18:03:33 - [ 7] Visit: 02 Visit Requirements: ON HOLD On Hold Comments: Waiting the implementation of G130M 1055 setting Additional Comments: Exposures ------------------------------------------------------------------------------------------------------------------------------------ Exp | Target | Instr | Oper. | Aper |Spectral|Central| Optional |Num| Time | Special Num | Name | Config | Mode | or FOV |Element |Waveln.| Parameters |Exp| | Requirements ------------------------------------------------------------------------------------------------------------------------------------ 1 J143832.04+ COS/NUV ACQ/IMA PSA MIRRORA 1 67 S 390333.7 GE ------------------------------------------------------------------------------------------------------------------------------------ 2 J143832.04+ COS/FUV TIME-TA PSA G130M 1055 BUFFER-TIME=1108, EXTENDED=NO, 1 1210 S 390333.7 G FLASH=YES, FP-POS=ALL, see below SEGMENT=BOTH Comments: The BUFFER-TIME= exposure time - 100s differs between the splits. The given BUFFER-TIME is for the shortest exposure. ------------------------------------------------------------------------------------------------------------------------------------ Sub Exposures ------------------------------------------------------------------------------------------------------------------------------------ Target | Exp |Instr | Oper. | Aper |Spectral|Cent.|Primary |Secondary |Iteration |CR-SPLIT |Orbit |Duration Name | Num |Config | Mode | or FOV |Element |Wave.|Pattern Pos |Pattern Pos |Num |Num |Number | ------------------------------------------------------------------------------------------------------------------------------------ J143832.04+390333.7 1 COS/NUV ACQ/IMA PSA MIRRORA none none none none 1 N/A GE J143832.04+390333.7 2 COS/FUV TIME-TA PSA G130M 1055 none none none 1 1 1208.0 S G J143832.04+390333.7 2 COS/FUV TIME-TA PSA G130M 1055 none none none 2 1 1209.0 S G J143832.04+390333.7 2 COS/FUV TIME-TA PSA G130M 1055 none none none 3 2 1443.0 S G J143832.04+390333.7 2 COS/FUV TIME-TA PSA G130M 1055 none none none 4 2 1444.0 S G 12539( 3) - 06-Feb-2012 18:03:33 - [ 8] Visit: 03 Visit Requirements: On Hold Comments: Additional Comments: Exposures ------------------------------------------------------------------------------------------------------------------------------------ Exp | Target | Instr | Oper. | Aper |Spectral|Central| Optional |Num| Time | Special Num | Name | Config | Mode | or FOV |Element |Waveln.| Parameters |Exp| | Requirements ------------------------------------------------------------------------------------------------------------------------------------ 1 J080250.09+ COS/NUV ACQ/IMA PSA MIRRORA 1 78 S 190119.2 GE ------------------------------------------------------------------------------------------------------------------------------------ 2 J080250.09+ COS/FUV TIME-TA PSA G130M 1291 BUFFER-TIME=1071, EXTENDED=NO, 1 1200 S 190119.2 G FLASH=YES, FP-POS=ALL, see below SEGMENT=BOTH Comments: The BUFFER-TIME= exposure time - 100s differs between the splits. The given BUFFER-TIME is for the shortest exposure. ------------------------------------------------------------------------------------------------------------------------------------ Sub Exposures ------------------------------------------------------------------------------------------------------------------------------------ Target | Exp |Instr | Oper. | Aper |Spectral|Cent.|Primary |Secondary |Iteration |CR-SPLIT |Orbit |Duration Name | Num |Config | Mode | or FOV |Element |Wave.|Pattern Pos |Pattern Pos |Num |Num |Number | ------------------------------------------------------------------------------------------------------------------------------------ J080250.09+190119.2 1 COS/NUV ACQ/IMA PSA MIRRORA none none none none 1 N/A GE J080250.09+190119.2 2 COS/FUV TIME-TA PSA G130M 1291 none none none 1 1 1171.0 S G J080250.09+190119.2 2 COS/FUV TIME-TA PSA G130M 1291 none none none 2 1 1172.0 S G J080250.09+190119.2 2 COS/FUV TIME-TA PSA G130M 1291 none none none 3 2 1418.0 S G J080250.09+190119.2 2 COS/FUV TIME-TA PSA G130M 1291 none none none 4 2 1417.0 S G 12539( 3) - 06-Feb-2012 18:03:33 - [ 9] Visit: 04 Visit Requirements: ON HOLD On Hold Comments: Waiting the implementation of G130M 1055 setting Additional Comments: Exposures ------------------------------------------------------------------------------------------------------------------------------------ Exp | Target | Instr | Oper. | Aper |Spectral|Central| Optional |Num| Time | Special Num | Name | Config | Mode | or FOV |Element |Waveln.| Parameters |Exp| | Requirements ------------------------------------------------------------------------------------------------------------------------------------ 1 J154459.66+ COS/NUV ACQ/IMA PSA MIRRORA 1 52 S 115017.5 GE ------------------------------------------------------------------------------------------------------------------------------------ 2 J154459.66+ COS/FUV TIME-TA PSA G130M 1055 EXTENDED=NO, FLASH=YES, 1 1200 S 115017.5 G FP-POS=ALL, SEGMENT=BOTH, see below BUFFER-TIME=1095 Comments: The BUFFER-TIME= exposure time - 100s differs between the splits. The given BUFFER-TIME is for the shortest exposure. ------------------------------------------------------------------------------------------------------------------------------------ Sub Exposures ------------------------------------------------------------------------------------------------------------------------------------ Target | Exp |Instr | Oper. | Aper |Spectral|Cent.|Primary |Secondary |Iteration |CR-SPLIT |Orbit |Duration Name | Num |Config | Mode | or FOV |Element |Wave.|Pattern Pos |Pattern Pos |Num |Num |Number | ------------------------------------------------------------------------------------------------------------------------------------ J154459.66+115017.5 1 COS/NUV ACQ/IMA PSA MIRRORA none none none none 1 N/A GE J154459.66+115017.5 2 COS/FUV TIME-TA PSA G130M 1055 none none none 1 1 1196.0 S G J154459.66+115017.5 2 COS/FUV TIME-TA PSA G130M 1055 none none none 2 1 1195.0 S G J154459.66+115017.5 2 COS/FUV TIME-TA PSA G130M 1055 none none none 3 2 1406.0 S G J154459.66+115017.5 2 COS/FUV TIME-TA PSA G130M 1055 none none none 4 2 1425.0 S G 12539( 3) - 06-Feb-2012 18:03:33 - [ 10] Summary Form for Proposal 12539 Item Used in this proposal ------------------------------------------------------------------------------------------------------------------------------------ Apertures PSA ------------------------------------------------------------------------------------------------------------------------------------ Configurations COS/FUV, COS/NUV ------------------------------------------------------------------------------------------------------------------------------------ Opmodes ACQ/IMAGE, TIME-TAG ------------------------------------------------------------------------------------------------------------------------------------ Optional Parameters BUFFER-TIME=1071, BUFFER-TIME=1085, BUFFER-TIME=1095, BUFFER-TIME=1108, EXTENDED=NO, FLASH=YES, FP-POS=ALL, SEGMENT=BOTH ------------------------------------------------------------------------------------------------------------------------------------ Special Requirements ON HOLD ------------------------------------------------------------------------------------------------------------------------------------ Spectral Elements G130M, MIRRORA ------------------------------------------------------------------------------------------------------------------------------------ Target Names J080250.09+190119.2, J143004.07+022213.3, J143832.04+390333.7, J154459.66+115017.5 ------------------------------------------------------------------------------------------------------------------------------------