Atacama Large Millimeter/submillimeter Array (ALMA) ######################## Publications making use of these data must include the following statement in the acknowledgement: This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00006.E. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ." In addition, publications from NA authors must include the standard NRAO acknowledgement: "The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc." ################################################################################### This Readme file contains a description of the VY CMa Band 9 ALMA 2017 High Frequency Long Baseline Campaign data, scripts and data products. There is one ALMA execution block (ASDM) observed in a mixed continuum TDM (2 SPWs) and high sprectral resolution (FDM) (6 SPWs) spectral resolution mode. Spectral windows include continuum and H2O and SiO masers. The array included 49 antennas in a configuration with baseline lengths between ~113 and 13894 m. The resolution achieved is ~12 x 11 milli arcseconds. The RMS achieved in the delivered images are: ~2.0 mJy/beam (continuum) 0.1 Jy/beam / 0.44km/s bin on the 658-GHz water maser line image (in line-free regions) raising to ~0.8Jy/beam measured away from signal but in strong emission channels Additional information can be found in the comments section included in this README. The Calibrated Data and Reference Images provided here were produced using CASA version 5.6.1-8 ######################## COMMENTS: Data were obtained on 2017 November 3rd as a part of the high frequency long baseline EOC (Extention and Optimization of Capability) in Band 9. The array included 49 antennas in a configuration with baseline lengths between ~113 and 13894 m. The PWV was 0.56 mm. Note that data was taken for a new imaging capaility commissioning, so several experimental functions were tested: band-to-band (B2B) phase referecing with differentail gain calibration (DGC), 90deg-Walsh phase switch for the Band~9 sideband separation. The total observing time was 47 minutes, and the on-source time for VY CMa was 344 s. The LO1 frequencies for the target and the phase calibrator are 669 (Band 9) and 149 (Band 4) GHz, respectively. The bright QSO, J0522-3627, 28 deg away from the target, was used as the DGC source, as well as for bandpass and flux-scale calibration. The correlator was configured to have eight spectral windows in Band 9 with a bandwidth of 1.875 GHz each using the 90-deg Walsh phase switching to separate upper side band and lower side band. Spectral window 33 covers the 658 GHz H2O maser. Spectral window 39 covers the 676~GHz SiO maser line. Those two spectral windows have the frequency resolution of 15.625 kHz. The line-free continuum emission image using the averaged eight spectral windows with Briggs weighting (robust = 2) has the synthesized beam size of 12 × 11 mas and the image RMS noise is 1.5 mJy/beam. In the initial map, there is a bright and compact component (VY CMa) that cannot be resolved even with the longest baselines, so that the self-calibration can be applied to the data (such a compact component is ideal for performing self-calibration). In order to obtain the visibility data free from the atmospheric and instrumental phase errors, phase and amplitude self-calibrations were performed for the line-free continuum channels with the solution interval of 16 s to remove residual phase offsets between the spectral windows. The self-calibration solutions were then applied to all target data, and the bright maser peak, which has a higher S/N, was used for further self-calibration, now making and applying the solutions with an interval of 4 s to all data. Improvements were seen on the continuum peak, which increased from 42 to 135 mJy/beam, while the image rms noise was largely unchanged at 1.5 mJy/beam. This indicates that stochastic phase errors remained even after B2B phase referencing, which gave an image coherence loss of ~70%. The final self-calibrated image was synthesized with the eight spectral windows' line-free emission. In the course of the data reduction, there were steps in quite unsusual processes to deal with B2B phase referencing and troubles related to the above functions such as the 90-deg Walsh phase switching. One antenna (DA54) was found to lack the first DGC sequence in spectral windows 41 and 43 and required extra attention. ######################## VYCMa_B9_B2B_Scripts.tgz contains the calibration and imaging scripts: * uid___A002_Xc660ef_X2d26.ms.scriptForCalibration_noplots.py * VYCMa_scriptForSelfcalAndImaging.py VYCMa_B9_B2B_UncalibratedData.tgz contains the raw data: One ASDM datasets: * uid___A002_Xc660ef_X2d26 VYCMa_B9_B2B_CalibratedData.tgz contains the calibrated uv-data: * vycma.cont.ms * vycma.line.ms To re-create these calibrated datasets you would need to execute the calibration scripts (e.g uid___A002_Xc660ef_X2d26.ms.scriptForCalibration_noplots.py) provided in VYCMa_B9_B2B_Scripts.tgz VYCMa_Band9_ReferenceImages.tgz contains the line and continuum self-calibrated image products and primary beams: VYCMa_sci.spw29_31_33_35_37_39_41_43.cont.tt0.I.manual-selfcal-robust.image.fits VYCMa_sci.spw29_31_33_35_37_39_41_43.cont.tt0.I.manual-selfcal-robust.pb.fits VYCMa_sci.spw39.cube.I.manual-selfcal-29SiO-J16-15-676GHz-robust.image.fit VYCMa_sci.spw39.cube.I.manual-selfcal-29SiO-J16-15-676GHz-robust.pb.fits VYCMa_sci.spw33.cube.I.manual-selfcal-H20-658GHz-robust.image.fits VYCMa_sci.spw33.cube.I.manual-selfcal-H20-658GHz-robust.pb.fits The reference images can be re-created by executing the imaging script (VYCMa_scriptForSelfcalAndImaging.py). ################################################################################### HOW TO EXECUTE THE CALIBRATION SCRIPTS (1) Put an ALMA ASDM and the corresponding calibration script in the same directory (2) In CASA 5.6.1-8 e.g. execfile('uid___A002_Xc660ef_X2d26.ms.scriptForCalibration_noplots.py') HOW TO EXECUTE THE IMAGING SCRIPT WITH SELF-CALIBRATION (1) Put vycma.cont.ms and vycma.line.ms, and the imaging script in the same directory (2) In CASA 5.6.1-8 execfile('VYCMa_scriptForSelfcalAndImaging.py')