TWHydraBand7: Difference between revisions

From CASA Guides
Jump to navigationJump to search
Cbrogan (talk | contribs)
Swood (talk | contribs)
 
(21 intermediate revisions by 6 users not shown)
Line 5: Line 5:
[[Image:Roberge STIS.jpg|thumb|HST STIS image of TW Hya ([http://adsabs.harvard.edu/abs/2005ApJ...622.1171R Roberge et al. 2005]; Figure 7). The solid green line is the direction of maximum disk brightness at optical wavelengths.]]
[[Image:Roberge STIS.jpg|thumb|HST STIS image of TW Hya ([http://adsabs.harvard.edu/abs/2005ApJ...622.1171R Roberge et al. 2005]; Figure 7). The solid green line is the direction of maximum disk brightness at optical wavelengths.]]


TW Hya is a pre-main sequence classical T Tauri star at a distance of about 52+/-1 pc (Mamajek 2005,2010).  It is the most studied member of the TW Hydra association (TWA) of low mass stars. From a wide variety of previous observations from the infrared to submillimeter, TW Hya is known to have a hot inner disk extending to radii < 4 AU, which is optically thin in the IR, and a larger cold dust disk out to about 200 AU (see for example the introduction by Vacca & Sandell 2011, and references therein).  Recent optical interferometry finds that TW Hya also contains a hot optically thick disk on even smaller size scales of ~0.5 AU, and suggests that the optically thin disk could be due to gas clearing by a planet (Akeson et al. 2011).  TW Hya is apparently still accreting from its disk at a rate of about (4-20) x 10<sup>-10</sup> Msun/year and the most recent estimates of its spectral type, mass, and age are M2.5V, 0.4 Msun, and 3 Myr (Vacca & Sandell 2011).  
TW Hya is a pre-main sequence classical T Tauri star at a distance of about 59.5 pc (Gaia Collaboration et al. 2016).  It is the most studied member of the TW Hydra association (TWA) of low mass stars. From a wide variety of previous observations from the infrared to submillimeter, TW Hya is known to have a hot inner disk extending to radii < 4 AU, which is optically thin in the IR, and a larger cold dust disk out to about 200 AU (see for example the introduction by Vacca & Sandell 2011, and references therein).  Recent optical interferometry finds that TW Hya also contains a hot optically thick disk on even smaller size scales of ~0.5 AU, and suggests that the optically thin disk could be due to gas clearing by a planet (Akeson et al. 2011).  TW Hya is apparently still accreting from its disk at a rate of about (4-20) x 10<sup>-10</sup> Msun/year and the most recent estimates of its spectral type are K8 (Ducourant et al. 2014) or M0.5 (Herczeg & Hillenbrand 2014) or M2.5V (Vacca & Sandell 2011). The most recent estimates of its mass and age are 0.7 Msun and 10^6.9 Myr (Siess et al. 2000, Ducourant+14).  


Millimeter and submillimeter observations of the continuum and spectral lines are particularly useful for tracing in the outer cold disk.  Previous observations by the VLA at 7 mm (Wilner et al. 2000), ATCA at 3 mm (Wilner et al. 2003), and the SMA at 1.3, 0.87, and 0.45 mm (Qi et al. 2004, 2006, 2008 and Hughes et al. 2011) reveal Keplerian rotation in the disk and an inclination angle of about 7 degrees (i.e. almost face-on).  Detailed studies of the dust continuum properties from the SMA work suggest that there are centimeter sized particles within the cold proto-planetary disk.
Millimeter and submillimeter observations of the continuum and spectral lines are particularly useful for tracing in the outer cold disk.  Previous observations by the VLA at 7 mm (Wilner et al. 2000), ATCA at 3 mm (Wilner et al. 2003), and the SMA at 1.3, 0.87, and 0.45 mm (Qi et al. 2004, 2006, 2008 and Hughes et al. 2011) reveal Keplerian rotation in the disk and an inclination angle of about 7 degrees (i.e. almost face-on).  A study of the 870um dust continuum with ALMA long baselines traced millimeter-sized particles in the disk down to spatial scales of 1 AU. The high resolution map revealed bright and dark ring-shaped structures associated with the concentration of solids in the disk.


== ALMA Data Overview ==
== ALMA Data Overview ==
Line 14: Line 14:
[[Image:Qi_HCOp3_2.gif|thumb|SMA HCO+(3-2) emission from TW Hya (Qi et al. 2008; Figure 2).]]
[[Image:Qi_HCOp3_2.gif|thumb|SMA HCO+(3-2) emission from TW Hya (Qi et al. 2008; Figure 2).]]


ALMA Science Verification data at Band 7 (~345 GHz) was taken for TW Hya on April 22, 2011.  A scheduling block about 1.5 hours long was run three times in a row for a total of about 4.5 hours of observing time.  The names of the three ASDMs were:  uid://A002/X1d9d21/X3c1, uid://A002/X1d9d21/X5d8  and uid://A002/X1d9d21/X7df.  Nine antennas were available during these runs, but one has to be flagged.  All four available basebands were used, resulting in four spectral windows (spws) containing data.  Two basebands were placed in the Lower Sideband (LSB) and two basebands in the Upper Sideband (USB).  In the LSB the CO(3-2) line at a rest frequency of 345.79599 GHz is located in spw=2.  In the USB the HCO+(4-3) line at a rest frequency of 356.7342 GHz is located in spw=0.  The other two spectral windows do not contain strong spectral lines and are used for continuum.  Each spectral window is 0.5 GHz wide and the channel width is 122 kHz.  Because the ALMA correlator was configured to apply Hanning smoothing of the signal, the effective spectral resolution is about twice the channel width, which in this case is about 0.2 km/s. For the antenna configuration in use at the time, the angular resolution is expected to be about 1.5".  The median value of precipitable water vapor (PWV) for this period was 1.16 mm, as measured by the water vapor radiometers. This PWV corresponds to an opacity of 0.20 at the CO(3-2) line.  The mean wind speed was 6.2 m/s.
ALMA Science Verification data at Band 7 (~345 GHz) was taken for TW Hya on April 22, 2011.  A scheduling block about 1.5 hours long was run three times in a row for a total of about 4.5 hours of observing time.  The names of the three ASDMs were:  uid://A002/X1d9d21/X3c1, uid://A002/X1d9d21/X5d8  and uid://A002/X1d9d21/X7ef.  Nine antennas were available during these runs, but one has to be flagged.  All four available basebands were used, resulting in four spectral windows (spws) containing data.  Two basebands were placed in the Lower Sideband (LSB) and two basebands in the Upper Sideband (USB).  In the LSB the CO(3-2) line at a rest frequency of 345.79599 GHz is located in spw=2.  In the USB the HCO+(4-3) line at a rest frequency of 356.7342 GHz is located in spw=0.  The other two spectral windows do not contain strong spectral lines and are used for measuring the continuum.  Each spectral window is 0.5 GHz wide and the channel width is 122 kHz.  Because the ALMA correlator was configured to apply Hanning smoothing of the signal, the effective spectral resolution is about twice the channel width, which in this case is about 0.2 km/s. For the antenna configuration in use at the time, the angular resolution is expected to be about 1.5".  The median value of precipitable water vapor (PWV) for this period was 1.16 mm, as measured by the water vapor radiometers. This PWV corresponds to an opacity of 0.20 at the CO(3-2) line.  The mean wind speed was 6.2 m/s.


The ALMA CO(3-2) data presented here is similar to the Submillimeter Array data presented in Hughes et al. 2011 (ApJ, 727, 85), though the SMA data have ~3 times smaller channel width at 50 kHz.  
The ALMA CO(3-2) data presented here is similar to the Submillimeter Array data presented in Hughes et al. 2011 (ApJ, 727, 85), though the SMA data have ~3 times smaller channel width at 50 kHz.  
Line 20: Line 20:
HCO+(4-3) data has not previously been published, but SMA HCO+(3-2) data (at 267.55762 GHz) is presented in Qi et al. 2008 (ApJ, 681, 1396). These SMA data have comparable angular resolution but a wider 203 kHz channel width.
HCO+(4-3) data has not previously been published, but SMA HCO+(3-2) data (at 267.55762 GHz) is presented in Qi et al. 2008 (ApJ, 681, 1396). These SMA data have comparable angular resolution but a wider 203 kHz channel width.


[[Image:TWHya_HCOp4_3_moments.png|center|700px]] ''ALMA HCO+(4-3) moment maps from TW Hya,  with white continuum contours at 3 and 100 sigma. From left to right: integrated intensity, intensity weighted velocity field, intensity weighted velocity dispersion are shown.''
Cubes are available online for CO(2-1), CN(2-1), and CS(5-4) along with an analysis of the velocity dispersion from Teague et al. 2016.
 
'''Using the data for publication''': Please use the acknowledgement given at the bottom of the [https://almascience.nrao.edu/alma-data/science-verification Science Verification Data page].
 
[[Image:TWHya_HCOp4_3_moments.png|center|frame|700px]] ''ALMA HCO+(4-3) moment maps from TW Hya,  with white continuum contours at 3 and 100 sigma. From left to right: integrated intensity, intensity weighted velocity field, intensity weighted velocity dispersion are shown.''


==Obtaining the Data==
==Obtaining the Data==
Line 41: Line 45:
==TWHya Data Reduction Tutorial==
==TWHya Data Reduction Tutorial==


The tutorial (called a casaguide) for reducing these data '''using CASA version 3.4 (&ge; r19874)''' has been split into calibration and imaging pages:
The tutorial (called a casaguide) for reducing these data '''using CASA version 4.3.0''' has been split into calibration and imaging pages:


# '''[[TWHydraBand7_Calibration_3.4]]''' this page requires that you download the TWHYA_BAND7_UnCalibratedMSAndTablesForReduction directory
# '''[[TWHydraBand7_Calibration_4.3]]''' this page requires that you download the TWHYA_BAND7_UnCalibratedMSAndTablesForReduction directory
# '''[[TWHydraBand7_Imaging_3.4]]''' this page requires that you have either used [[TWHydraBand7_Calibration_3.4]] to obtain the TWHYA_BAND7_CalibratedData or that you have downloaded this directory.
# '''[[TWHydraBand7_Imaging_4.3]]''' this page requires that you have either used [[TWHydraBand7_Calibration_4.3]] to obtain the TWHYA_BAND7_CalibratedData or that you have downloaded this directory.


Alternatively you can just download the final images (TWHYA_BAND7_ReferenceImages directory) if you only want to see the final results.
Alternatively you can just download the final images (TWHYA_BAND7_ReferenceImages directory) if you only want to see the final results.


'''NOTE: CASA 3.4 or later is required to process the data using the guides above''' see http://casa.nrao.edu/casa_obtaining.shtml. The data products have not changed. The older CASA 3.3 casaguide versions are still available at: [[TWHydraBand7_Calibration_for_CASA_3.3]] and [[TWHydraBand7_Imaging_for_CASA_3.3]].
'''NOTE: CASA 4.3 is required to process the data using the guides above''' see http://casa.nrao.edu/casa_obtaining.shtml. The data products have not changed.  
* The older CASA 4.2 casaguide versions are still available at: [[TWHydraBand7_Calibration_4.2]] and [[TWHydraBand7_Imaging_4.2]].
* The older CASA 4.1 casaguide versions are still available at: [[TWHydraBand7_Calibration_4.1]] and [[TWHydraBand7_Imaging_4.1]].
* The older CASA 4.0 casaguide versions are still available at: [[TWHydraBand7_Calibration_4.0]] and [[TWHydraBand7_Imaging_4.0]].
* The older CASA 3.4 casaguide versions are still available at: [[TWHydraBand7_Calibration_3.4]] and [[TWHydraBand7_Imaging_3.4]].
* The older CASA 3.3 casaguide versions are still available at: [[TWHydraBand7_Calibration_for_CASA_3.3]] and [[TWHydraBand7_Imaging_for_CASA_3.3]].


NOTE: These guides are dynamic and will evolve as our understanding of how best to reduce ALMA data improves. Check back for updates periodically.
NOTE: These guides are dynamic and will evolve as our understanding of how best to reduce ALMA data improves. Check back for updates periodically.
Line 66: Line 75:
Tabs matter in python, make sure that commands that span more than one line and  
Tabs matter in python, make sure that commands that span more than one line and  
"for" loops keep their spacing. Sometimes (especially "for" loops) you may need to  
"for" loops keep their spacing. Sometimes (especially "for" loops) you may need to  
explicitly hit enter twice to get the command going.
explicitly hit enter twice to get the command going. Use the CASA command 'cpaste'
to enter block commands and preserve whitespace.
</source>
</source>



Latest revision as of 02:05, 10 January 2017


Science Target Overview

HST STIS image of TW Hya (Roberge et al. 2005; Figure 7). The solid green line is the direction of maximum disk brightness at optical wavelengths.

TW Hya is a pre-main sequence classical T Tauri star at a distance of about 59.5 pc (Gaia Collaboration et al. 2016). It is the most studied member of the TW Hydra association (TWA) of low mass stars. From a wide variety of previous observations from the infrared to submillimeter, TW Hya is known to have a hot inner disk extending to radii < 4 AU, which is optically thin in the IR, and a larger cold dust disk out to about 200 AU (see for example the introduction by Vacca & Sandell 2011, and references therein). Recent optical interferometry finds that TW Hya also contains a hot optically thick disk on even smaller size scales of ~0.5 AU, and suggests that the optically thin disk could be due to gas clearing by a planet (Akeson et al. 2011). TW Hya is apparently still accreting from its disk at a rate of about (4-20) x 10-10 Msun/year and the most recent estimates of its spectral type are K8 (Ducourant et al. 2014) or M0.5 (Herczeg & Hillenbrand 2014) or M2.5V (Vacca & Sandell 2011). The most recent estimates of its mass and age are 0.7 Msun and 10^6.9 Myr (Siess et al. 2000, Ducourant+14).

Millimeter and submillimeter observations of the continuum and spectral lines are particularly useful for tracing in the outer cold disk. Previous observations by the VLA at 7 mm (Wilner et al. 2000), ATCA at 3 mm (Wilner et al. 2003), and the SMA at 1.3, 0.87, and 0.45 mm (Qi et al. 2004, 2006, 2008 and Hughes et al. 2011) reveal Keplerian rotation in the disk and an inclination angle of about 7 degrees (i.e. almost face-on). A study of the 870um dust continuum with ALMA long baselines traced millimeter-sized particles in the disk down to spatial scales of 1 AU. The high resolution map revealed bright and dark ring-shaped structures associated with the concentration of solids in the disk.

ALMA Data Overview

SMA CO(3-2) emission from TW Hya (Hughes et al. 2011; Figure 2).
SMA HCO+(3-2) emission from TW Hya (Qi et al. 2008; Figure 2).

ALMA Science Verification data at Band 7 (~345 GHz) was taken for TW Hya on April 22, 2011. A scheduling block about 1.5 hours long was run three times in a row for a total of about 4.5 hours of observing time. The names of the three ASDMs were: uid://A002/X1d9d21/X3c1, uid://A002/X1d9d21/X5d8 and uid://A002/X1d9d21/X7ef. Nine antennas were available during these runs, but one has to be flagged. All four available basebands were used, resulting in four spectral windows (spws) containing data. Two basebands were placed in the Lower Sideband (LSB) and two basebands in the Upper Sideband (USB). In the LSB the CO(3-2) line at a rest frequency of 345.79599 GHz is located in spw=2. In the USB the HCO+(4-3) line at a rest frequency of 356.7342 GHz is located in spw=0. The other two spectral windows do not contain strong spectral lines and are used for measuring the continuum. Each spectral window is 0.5 GHz wide and the channel width is 122 kHz. Because the ALMA correlator was configured to apply Hanning smoothing of the signal, the effective spectral resolution is about twice the channel width, which in this case is about 0.2 km/s. For the antenna configuration in use at the time, the angular resolution is expected to be about 1.5". The median value of precipitable water vapor (PWV) for this period was 1.16 mm, as measured by the water vapor radiometers. This PWV corresponds to an opacity of 0.20 at the CO(3-2) line. The mean wind speed was 6.2 m/s.

The ALMA CO(3-2) data presented here is similar to the Submillimeter Array data presented in Hughes et al. 2011 (ApJ, 727, 85), though the SMA data have ~3 times smaller channel width at 50 kHz.

HCO+(4-3) data has not previously been published, but SMA HCO+(3-2) data (at 267.55762 GHz) is presented in Qi et al. 2008 (ApJ, 681, 1396). These SMA data have comparable angular resolution but a wider 203 kHz channel width.

Cubes are available online for CO(2-1), CN(2-1), and CS(5-4) along with an analysis of the velocity dispersion from Teague et al. 2016.

Using the data for publication: Please use the acknowledgement given at the bottom of the Science Verification Data page.

ALMA HCO+(4-3) moment maps from TW Hya, with white continuum contours at 3 and 100 sigma. From left to right: integrated intensity, intensity weighted velocity field, intensity weighted velocity dispersion are shown.

Obtaining the Data

To download the data, click on the region below that is closest to your location:

This will take you to a webpage with the following files:

  • TWHYA_BAND7_UnCalibratedMSAndTablesForReduction.tgz - This contains the uncalibrated data, already converted from raw data in ALMA Science Data Model (ASDM) format to CASA Measurement Sets (MS). We did this using the importasdm task in CASA. Along with the uncalibrated data, we also provide some tables that you will need for the calibration which cannot currently be generated inside of CASA (for Early Science, these tables will either be pre-applied or supplied with the data).
  • TWHYA_BAND7_CalibratedData.tgz - Contains only the calibrated uv-data for TWHya ready for imaging and self-calibration.
  • TWHYA_BAND7_ReferenceImages.tgz - The final spectral line and continuum images.

***Before*** you begin to download, read the details of TWHydraBand7#Data_Reduction_Tutorial below to see which files you want, these files are quite large and you may not want to download everything.

TWHya Data Reduction Tutorial

The tutorial (called a casaguide) for reducing these data using CASA version 4.3.0 has been split into calibration and imaging pages:

  1. TWHydraBand7_Calibration_4.3 this page requires that you download the TWHYA_BAND7_UnCalibratedMSAndTablesForReduction directory
  2. TWHydraBand7_Imaging_4.3 this page requires that you have either used TWHydraBand7_Calibration_4.3 to obtain the TWHYA_BAND7_CalibratedData or that you have downloaded this directory.

Alternatively you can just download the final images (TWHYA_BAND7_ReferenceImages directory) if you only want to see the final results.

NOTE: CASA 4.3 is required to process the data using the guides above see http://casa.nrao.edu/casa_obtaining.shtml. The data products have not changed.

NOTE: These guides are dynamic and will evolve as our understanding of how best to reduce ALMA data improves. Check back for updates periodically.

How to Use A casaguide

For tips on using CASA and ways CASA can be run, see EVLA_Spectral_Line_Calibration_IRC+10216#How_to_Use_This_casaguide page.

To learn how to extract executable Python scripts from the tutorial, see Extracting_scripts_from_these_tutorials.

In the guides

# In CASA
Regions of this color are CASA commands (or definitions) that need to be cut and 
pasted in sequence. Wait until one command is finished before pasting another. 
Tabs matter in python, make sure that commands that span more than one line and 
"for" loops keep their spacing. Sometimes (especially "for" loops) you may need to 
explicitly hit enter twice to get the command going. Use the CASA command 'cpaste' 
to enter block commands and preserve whitespace.
Information in this color shows excerpts from the CASA Logger output
This color shows you background information about the data or other types of reference material