ALMA2014 LBC SVDATA: Difference between revisions
| Line 59: | Line 59: | ||
*'''Observing narrow ''thermal'' line emission at high resolution is challenging''': Except for Mira, which has the advantage of having compact line emission, as well as emission from highly excited or even masing lines, the imaging of line emission at the very high resolution afforded by the longest baselines is challenging. For example, at 89 GHz with 0.25 km/s channels, 30 antennas, 3.5 hours on-source, and an angular resolution of 70 mas, the expected spectral rms noise per channel is 2.6 mJy/beam but only 83 K in surface brightness sensitivity. As a result, it was necessary to taper the spectral line images from HL Tau and SDP.81 to achieve good detections. | *'''Observing narrow ''thermal'' line emission at high resolution is challenging''': Except for Mira, which has the advantage of having compact line emission, as well as emission from highly excited or even masing lines, the imaging of line emission at the very high resolution afforded by the longest baselines is challenging. For example, at 89 GHz with 0.25 km/s channels, 30 antennas, 3.5 hours on-source, and an angular resolution of 70 mas, the expected spectral rms noise per channel is 2.6 mJy/beam but only 83 K in surface brightness sensitivity. As a result, it was necessary to taper the spectral line images from HL Tau and SDP.81 to achieve good detections. | ||
*'''Combining bands can produce amazing results''': For sources that have significant spatial changes in continuum spectral index, and wide observed fractional bandwidth (>15%) for which the signal-to-noise is high throughout the frequency range, it is possible to use a technique in {{clean}} to simultaneously account for the changing spectral index (<math>S_{\nu}\propto \nu^{\alpha}</math>) to produce a high fidelity image, and a map of the spectral index along with its estimated error (see Rau, U., & Cornwell, T.J. 2011, A&A, 532, AA71). However, the caveat that the signal to noise must be very good, especially at the frequency extrema is a strong one. The results of this technique on HL Tau Band 6+7 are dramatic, and it also provides some improvement for Mira at Band 3. However, our test show that for the weaker emission at Band 3 for HL Tau and all bands for SDP.81 at native resolution, this technique does not work well. Bottom line: use with caution, but it can produce amazing results. Since this method automatically masks the resulting spectral index map where it estimates <math>\alpha/\alpha_{error}< 4</math> a good rule of thumb is that if the spectral index map is mostly masked in the image signal regions, there is inadequate S/N for this technique has not been successful. | *'''Extended emission at high angular resolution can best be imaged with multiscale''': Standard clean techniques that use a delta function formulation for the clean beam, tend to make extended emission look pixelated, due to the so-called "clean-instability". For the LBC SV data we had good success using the ''multiscale'' parameter in {{clean}} with typically ''multiscale=[0,5,15]'' where these scales are defined as multipliers of the cellsize. | ||
*'''Combining bands can produce amazing results''': For sources that have significant spatial changes in continuum spectral index, and wide observed fractional bandwidth (>15%) for which the signal-to-noise is high throughout the frequency range, it is possible to use a technique in {{clean}} to simultaneously account for the changing spectral index (<math>S_{\nu}\propto \nu^{\alpha}</math>) to produce a high fidelity image, and a map of the spectral index along with its estimated error (see Rau, U., & Cornwell, T.J. 2011, A&A, 532, AA71). However, the caveat that the signal to noise must be very good, especially at the frequency extrema is a strong one. The results of this technique on HL Tau Band 6+7 are dramatic, and it also provides some improvement for Mira at Band 3. However, our test show that for the weaker emission at Band 3 for HL Tau and all bands for SDP.81 at native resolution, this technique does not work well. Bottom line: use with caution, but it can produce amazing results. Since this method automatically masks the resulting spectral index (<math>\alpha</math>) map where it estimates <math>\alpha/\alpha_{error}< 4</math> a good rule of thumb is that if the spectral index map is mostly masked in the image signal regions, there is inadequate S/N for this technique has not been successful. | |||
Revision as of 18:47, 15 February 2015
This page is currently under construction.
To navigate the CASAguides pages, visit [http://casaguides.nrao.edu/ casaguides.nrao.edu ]
Overview
This casaguide describes the imaging used for the ALMA 2014 Long Baseline Campaign (LBC) Science Verification (SV) data. The purpose of this campaign was to commission baseline lengths up to 15 km in preparation for this capability to be offered for Cycle 3. As part of the overall commissioning plan, five science targets were chosen to span the range of properties that would be needed to offer the capability and to demonstrate to the community both the power and caveats associated with long baseline observing. The SV portion of the campaign was carried out from October 12 to November 14, 2014. As described below the data products for four of the five targets can be dowloaded from the regional ALMA Science Portals, while the imaging scripts and salient details are provided on this page (the final target 3C138, will be released separately with other SV polarization data in the near future).
Obtaining the Data
To download the data, go to the ALMA Science Portal (SP) Science Verification data page for the region closest to your location:
Choose a LBC SV science target: 13. Juno, 14. Mira, 15. HL Tau, or 16. SDP.81
There you will find gzipped tar files of the form (these files can be unpacked with the tar -xvzf command):
- Target_Band_Freq_UnCalibratedData: Raw, uncalibrated data in ALMA Science Data Model (asdm) format
- Target_Band_Freq_CalibrationScripts: Calibration scripts that go along with the UnCalibratedData. These scripts can be executed once placed in the same directory as the UnCalibratedData
- Target_Band_Freq_CalibratedData: Calibrated and combined uv-data. In the case of continuum datasets, an appropropriate amount of channel averaging has been applied to reduce the data size. Line data sets are provided at their full resolution. If self-calibration was possible, you will also find the final phase-only and amplitude tables in a subdirectory called "calibration". This gives the choice of either performing your own self-calibration or applying the supplied tables. If you are going to apply these tables, be sure to also run any flagdata commands in the imaging scripts. If the self-calibration tables are not present, we deemed self-calibration unfeasible. Note that the imaging scripts are provided on this page, not at the science portal.
- Target_Band_Freq_ReferenceImages: Final reference images with self-calibration applied where possible. For the few cases where primary beam correction is warranted these files are also provided. Continuum images with .pbcor in the name have already been corrected. The line files with .flux in the name are the primary beam images themselves and can be applied to the corresponding line data if desired using the CASA task impbcor.
- README: Summarizes the data available for that science target.
The data have been packaged in this manner to allow a wide range of interaction with the data, ranging from complete re-reduction, to simply having a look at the final reference images. Please only download the data products you expect to use -- we expect the initial load on the servers
could be high and the UnCalibratedData, and in the case of targets/bands with spectral line data the CalibratedData, files are quite large.
CASA Version
The data were calibrated and imaged using CASA 4.2.2. At this time, we have no reason to believe that they would not work in CASA 4.3 with a few caveats: (1) The calibration scripts check that the version they were made in matches the version being used. If you want to try them in CASA 4.3 you will need to change this section of each calibration script. (2) Be aware that data imported in CASA 4.3 cannot then be imaged in CASA 4.2.2 due to the necessary addition of a new keyword in the MS format. Please note, however, that only CASA 4.2.2 has been tested.
CASA can be downloaded from: http://casa.nrao.edu/casa_obtaining.shtml The CASA 4.2.2 version is called ALMA pipeline CASA 4.2.2
LBC SV Imaging
For each LBC SV target there is a brief description of the data and imaging, along with links to the imaging scripts. For spectral line targets, a table of the lines observed is also given. The "detection" columns are merely descriptive to give a sense of what the data show, they are not meant to be taken as scientific assessments. The continuum was detected for all four of the science targets at each observed band. For all of the targets the number of antennas employed was of order 30.
The most important take home messages:
- UV-Coverage is critically important: For each target/band a scheduling block was designed to run for 60-70 minutes, with about 30 minutes on source per execution. The number of executions carried out depended on the complexity of the imaging required, from short snapshots of the bright compact emission from Juno to to the many executions required for the complex morphologies of HL Tau and SDP.81.
- Observing narrow thermal line emission at high resolution is challenging: Except for Mira, which has the advantage of having compact line emission, as well as emission from highly excited or even masing lines, the imaging of line emission at the very high resolution afforded by the longest baselines is challenging. For example, at 89 GHz with 0.25 km/s channels, 30 antennas, 3.5 hours on-source, and an angular resolution of 70 mas, the expected spectral rms noise per channel is 2.6 mJy/beam but only 83 K in surface brightness sensitivity. As a result, it was necessary to taper the spectral line images from HL Tau and SDP.81 to achieve good detections.
- Extended emission at high angular resolution can best be imaged with multiscale: Standard clean techniques that use a delta function formulation for the clean beam, tend to make extended emission look pixelated, due to the so-called "clean-instability". For the LBC SV data we had good success using the multiscale parameter in clean with typically multiscale=[0,5,15] where these scales are defined as multipliers of the cellsize.
- Combining bands can produce amazing results: For sources that have significant spatial changes in continuum spectral index, and wide observed fractional bandwidth (>15%) for which the signal-to-noise is high throughout the frequency range, it is possible to use a technique in clean to simultaneously account for the changing spectral index ([math]\displaystyle{ S_{\nu}\propto \nu^{\alpha} }[/math]) to produce a high fidelity image, and a map of the spectral index along with its estimated error (see Rau, U., & Cornwell, T.J. 2011, A&A, 532, AA71). However, the caveat that the signal to noise must be very good, especially at the frequency extrema is a strong one. The results of this technique on HL Tau Band 6+7 are dramatic, and it also provides some improvement for Mira at Band 3. However, our test show that for the weaker emission at Band 3 for HL Tau and all bands for SDP.81 at native resolution, this technique does not work well. Bottom line: use with caution, but it can produce amazing results. Since this method automatically masks the resulting spectral index ([math]\displaystyle{ \alpha }[/math]) map where it estimates [math]\displaystyle{ \alpha/\alpha_{error}\lt 4 }[/math] a good rule of thumb is that if the spectral index map is mostly masked in the image signal regions, there is inadequate S/N for this technique has not been successful.
Juno - Asteroid
Asteroid (3) Juno was observed in continuum in Band 6 (233 GHz) to demonstrate the ephemeris capability. The observations were taken over 60% of Juno's rotation period of 7.2 hours. The data are split into approximately 18min intervals (10min on source) spread across 10 individual datasets and imaged and self-calibrated separately in order to prevent smearing due to the rapid rotation. The angular resolution achieved ranges from 32 x 24 mas to 42 x 36 mas, differing primarily due to the level of phase stability as the observations span the transition from night through dawn and into daytime, but changing uv-coverage also plays a role.
Juno Band 6 (5 executions)
Mira - AGB Star
Mira (omicron Cet) was observed in continuum and spectral lines in Bands 3 and 6. Angular resolutions achieved for the continuum data are 70 x 60 mas and 34 x 24 mas in Bands 3 (94 GHz) and 6 (231 GHz), respectively. In order to achieve the best possible dynamic range, the continuum and SiO v=1 line data were self-calibrated independently. The maser emission from the SiO v=1 line was used to derive self-calibration for the line data. In order to obtain high dynamic range on the band 3 continuum, which has a wide fractional bandwidth, it was also necessary to take into account the spectral index with nterms=2.
Mira Band 3 (3 executions)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| SiO v=0 | 2-1 | 2.1 | 86.8470 | weak maser? + thermal |
| SiO v=1 | 2-1 | 1771 | 86.2434 | strong maser + thermal |
| SiO v=2 | 2-1 | 3523 | 85.6405 | thermal + absorption |
| 29SiO v=0 | 2-1 | 2.1 | 85.7590 | weak maser? + thermal |
Mira Band 6 (2 executions)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| H30alpha | - | - | 231.90093 | none |
| H2O v2=1 | 5(5,0)-6(4,3) | 3450.7 | 232.68670 | thermal |
| SiO v=0 | 5-4 | 20.84 | 217.10498 | weak maser? + thermal |
| SiO v=1 | 5-4 | 1789.8 | 215.59595 | strong maser + thermal |
| SiO v=2 | 5-4 | 3541.7 | 214.08854 | thermal + absorption |
| 29SiO v=0 | 5-4 | 20.58 | 214.38576 | weak maser? + thermal |
HL Tau - Protoplanetary Disk
HL Tau was observed in continuum and spectral lines in Band 3, and (only) continuum in Bands 6 & 7. Angular resolutions achieved for the continuum data are 85 x 61 mas, 35 x 22 mas, and 30 x 19 mas in Bands 3 (101.9 GHz), 6 (233.0 GHz) and 7 (343.5 GHz), respectively. It was necessary to uv-taper the spectral line data so the resulting angular resolution is coarser (ranging from 0.25arcsec to 1.1 arcsec). In addition to the more standard image products we also provide an example of combining the Band 6 and Band 7 continuum data using nterms=2 to derive the spectral index and create a high fidelity image.
HL Tau Band 3 Lines (7 executions for each of two tunings: 89 and 115 GHz)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| HCN | 1-0 | 0.0 | 88.63185 | absorption |
| HCO+ | 1-0 | 0.0 | 89.18853 | outflow, disk, absorption |
| CN | 1-0 | 0.0 | 113.49097 | absorption |
| 12CO | 1-0 | 0.0 | 115.2712 | outflow, confused |
Note only the frequency of the strongest hyperfine (HF) transition is listed for HCN (3 HFs) and CN (4 HFs).
HL Tau Band 3 Continuum (combined continuum from both line tunings)
HL Tau Band 6 (continuum only, 9 executions)
HL Tau Band 7 (continuum only, 10 executions)
HL Tau Band 6+7 (nterms=2) (continuum only)
SDP.81 - Lensed Galaxy
The lensed galaxy SDP.81 was observed in continuum and spectral lines in Bands 4, 6 & 7. Angular resolutions achieved for the continuum data are 60 x 54 mas, 39 x 30 mas and 31 x 23 mas in Bands 4 (151 GHz), 6 (236 GHz), and 7 (290 GHz), respectively. It was necessary to uv-taper the spectral line data so the resulting angular resolution is coarser (up to 169 mas). This source lies at a red-shift of ~3, accounting for the large offset in rest frequency compared to the observing band. In addition to the more standard image products we also provide an example of combining the Band 6 and Band 7 continuum data using nterms=2 to derive the spectral index and create a high fidelity image. However it is important to note that this procedure did not work well at full resolution, as evidenced by the almost complete masking of the resulting spectral map -- this is an indication that the spectral index and hence resulting flux-densities are not well constrained. However, using a bit of taper (larger beam) improved the signal enough for the technique to work.
SDP.81 Band 4 (12 executions)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| CO | 5-4 | 55.3 | 576.26793 | detection |
SDP.81 Band 6 (9 executions)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| CO | 8-7 | 154.9 | 921.79970 | detection |
| H2O v=0 | 2(0,2)-1(1,1) | 53.4 | 987.92676 | weak detection |
SDP.81 Band 7 (11 executions)
| Line | Transition | El (K) | Rest Freq | detection |
|---|---|---|---|---|
| CO | 10-9 | 248.9 | 1151.98544 | detection |
SDP.81 Band 6+7, (tapered nterms=2) (continuum only)
