Step-by-Step ETC Guide for NIRCam and MIRI Coronagraphy of HR8799 b

A walk-through of the JWST ETC for the "NIRCam and MIRI Coronagraphy of HR8799 b" Example Science Program is provided, demonstrating how to select exposure parameters for this observing program.

On this page

Introduction

Main articles: NIRCam and MIRI Coronagraphy of HR8799 bJWST Exposure Time Calculator Overview
See also: JWST High-Contrast ImagingHCI RoadmapMIRI Coronagraphic ImagingNIRCam Coronagraphic Imaging,

The JWST Exposure Time Calculator (ETC) performs signal-to-noise (SNR) calculations for the JWST observing modes. Sources of interest are defined by the user and assigned to scenes which are used by the ETC to run calculations for the requested observing mode. 

This JWST Example Science Program, "NIRCam and MIRI Coronagraphy of HR8799 b", is based on the on observations in the Guaranteed Time Observer (GTO) Program #1194: "Characterization of the HR 8799 planetary system and planet search" (PI: Charles A. Beichman) with great simplification. The goals of this GTO program are two-fold: to search for previously unknown planets using NIRCam coronagraphy and the physical characterization of the known planets, HR8789bcde, using NIRCam and MIRI multi-filter photometry. This example science program presents a simplified version of the GTO program, with the aim of showcasing the workflow of building a coronagraphic imaging observation of a point source using both NIRCam and MIRI.

Following the instructions in Stage 5 of the HCI Roadmap, we focus on selecting the exposure parameters to detect the Beta Pictoris debris disk at the desired signal-to-noise ratio (SNR) for the MIRI Coronagraphic Imaging and NIRCam Coronagraphic Imaging observing modes. The optimal exposure specifications (e.g., number of groups and integrations) are the input needed for the Astronomer's Proposal Tool (APT) observation template, which is used to specify an observing program and submit proposals.

The ETC workbook associated with this Example Science Program is called "#36: NIRCam and MIRI Coronagraphy of HR 8799 b" and can be selected from the Get a Copy of an Example science program dropdown on the ETC Workbooks page to get the Read-only version. The nomenclature and reported SNR values in this article are based on ETC v.1.5. There may subtle differences if using a different version of ETC.


Define Sources and Scenes 

Main article: JWST ETC Scenes and Sources Page Overview
See also: JWST ETC Defining a New Source, JWST ETC Source Spectral Energy Distribution

In order to perform our intended calculations, we must first build a library of relevant sources and scenes for our ETC calculations. In the ETC, scenes are idealized representations of spatial (two angular coordinates) and spectral brightness distributions, before being observed by a telescope. For MIRI and NIRCam Coronagraphic Imaging modes, scenes are composed of relatively small areas (8.19" and 6.36", respectively) centered on the coronagraphic masks. When defined, a scene should contain the source targets of an observation, as well as all other nearby sources that could contribute to the observed target and background fluxes within that area.

The source library for this program consists of five sources: the central star, HR8799, its planet b, orbiting at a distance of ~1.73'' and three reference stars to be used for PSF subtraction in the ETC calculations. The reference stars differ by brightness and spectral type, and such differences will be exploited to illustrate separate use cases. 

To create a new source we go to the Scenes and Sources page and click the NEW button within the Select a Source pane. Once created, the source properties can be modified using various tabs within the Source Editor pane. Creating five new sources, we change the source properties as described below (note that the changes made using the source editor under each tab are applied to the selected source only when the SAVE button is activated):

Table 1. Input Source Parameters

ID #ID tabContinuum tabRenorm tabShape tab
  • we update the Source identity information to "HR8799"
  • from the library of "Pheonix Stellar Modelswe Select the "F0V 7250 4.0" template stellar spectrum1
  • we normalize the spectral energy distribution for the source by providing the source's integrated magnitude of "5.24" "Vegamag" in the "Johnson" "K" bandpass
  • we keep the default option for the Shape of source as Point
2
  • we update the Source identity information to "planet b"
  • we assign the source a Blackbody Spectrum with an affective temperature (a Tb) of "1100" K2
  • we normalize the spectral energy distribution for the source by providing the source's integrated magnitude of "14.25" "Vegamag" in the "Johnson" "K" bandpass
  • we keep the default option for the Shape of source as Point
3
  • we update the Source identity information to "Reference (HD218261)"
  • from the library of "Pheonix Stellar Modelswe Select the "F8V 6250 4.5" template stellar spectrum3
  • we normalize the spectral energy distribution for the source by providing the source's integrated magnitude of "5.14" "Vegamag" in the "Johnson" "K" bandpass
  • we keep the default option for the Shape of source as Point
4
  • we update the Source identity information to "Reference (HD218381)"
  • from the library of "Pheonix Stellar Modelswe Select the "K0V 5250 4.5" template stellar spectrum4
  • we normalize the spectral energy distribution for the source by providing the source's integrated magnitude of "4.41" "Vegamag" in the "Johnson" "K" bandpass
  • we keep the default option for the Shape of source as Point
5
  • we update the Source identity information to "Reference (HD220657)"
  • from the library of "Pheonix Stellar Modelswe Select the "G0III 5750 3.0" template stellar spectrum5
  • we normalize the spectral energy distribution for the source by providing the source's integrated magnitude of "3.04" "Vegamag" in the "Johnson" "K" bandpass
  • we keep the default option for the Shape of source as Point

Table notes:

1 Values from SIMBAD

2 This is a simplifying assumption for illustrative purposes. It is understood that a more realistic spectral shape, e.g. from an existing observed spectrum and/or model spectrum should be used if available

3 Values from SIMBAD

4 Values from SIMBAD

5 Values from SIMBAD, the spectral type used here is the closest available in the ETC library of PHOENIX spectra to the SIMBAD value of F8III. Such difference is of little importance for imaging at the MIRI wavelengths


With our source properties now defined, we will create the astronomical scenes in which to place them. To perform our ETC calculations we require 3 Science scenes for our science calculations, and 3 Reference scenes on which to perform our Target Acquisition calculations on the three reference sources. To create a new scene, we need to click on the NEW button in the Select a Scene pane. To add a source to a scene we select the desired source in the Select a Source pane, select the scene to which the source will be associated from the Select a Scene, and click the ADD SOURCE button. We define our scenes and associate them with the appropriate sources, as follows: 

Table 2. Scene configurations, as seen in the Select a Scene pane:

IDName6Sources
1

Science - NIRCam 

1,2,3
2Science - NIRCam bright1,2,4
3Science - MIRI1,2,5
4Target Acq - HD2182613
5Target Acq - HD2183814
6Target Acq - HD2206575

Table Notes:

6To generate the scene Name we go into the Source Editor pane and update the Source Identity Information 


With our sources now associated with the relevant scenes, we can use the Offsets tab in the Source Editor pane to place the sources at their desired locations. The source's spatial offsets are defined with respect to the center of the scene. Selecting each scene and source in turn (refer to the Sources and Scene Tables section for a more in-depth explanation) we provide the following source offsets:

  • In each Science Scene (i.e. scene ID: 1, Name: "Science - NIRCam"; ID: 2, Name: "Science - NIRCam bright"; and ID: 3, Name: "Science - MIRI", respectively):
    • we place the target star (i.e. source ID: 1, Name: "HR8799") at an X Offset and Y Offset of "0" arcsec
    • we place the planet (i.e. source ID1, Name: "planet b") at an X Offset of "-1.6" arcsec and Y Offset of "0.65" arcsec
    • we place the reference star source (i.e. source ID: 3, Name: "Reference (HD218261)"ID4, Name: "Reference (HD218381)"; and ID5, Name: "Reference (HD220657)", respectively) at an X Offset and Y Offset of "10" arcsec7.
  • In each Reference Scene (i.e. scene ID4Name: "Target Acq - HD218261"; ID5Name: "Target Acq - HD218381"; and ID5Name: "Target Acq - HD220657", respectively)
    • we position the reference PSF target (i.e. scene ID: 3, Name: "Reference (HD218261)", ID4, Name"Reference (HD218381)" and ID5, Name: "Reference (HD220657)", respectively) at the center of the scene, defining an X Offset and Y Offset of "0" arcsec.

7In order to perform ETC calculations using either one of the reference stars, they have to be added to the ETC scene and be placed at a sufficiently large distance, as to not contaminate the target and thus not interfere with the flux of the target star or its planet. This is just a trick that allows to perform simulations of the PSF subtraction process. The "10'' arcsec Offset from the target in each of the scenes is a convenience value and does not reflect the true angular distance between target and reference star).

Once each scene has been defined an idealized idealized 2D representation of the scene can be viewed in the Scene Sketch pane.


Define Calculations

Main article: JWST ETC Creating a New Calculation
See also: JWST ETC Calculations Page Overview

The Calculations page in the ETC is where we specify the instrument mode, background, instrument and detector configuration, observing setup and extraction strategy for a given calculation. The aim of these calculations is to make a comparative analysis of locations in the scene and determine the exposure parameters required to obtain the desired signal-to-noise ratio (SNR). 

With JWST, a coronagraphic observation consists of several steps. First, the telescope must slew towards the desired target. After the slew, a target acquisition (TA) observation is performed to place the desired target behind the selected coronagraphic mask. A science observation follows. This observation sequence (slew→TA→science observation) is repeated twice: once for the actual science target and once for the PSF reference star. The latter is used to build a coronagraphic PSF model that is subtracted from the coronagraphic target observation.

The ETC allows the user to simulate all 4 observations: 2 targets x 2 modes (TA or coronagraphic exposure). The ETC also allows the user to simulate the PSF subtraction process for the primary target coronagraphic observation. In this article we will illustrate how.

We note that the Example Science Program #36 ETC workbook associated to this page only exemplifies 3 of the 4 observations; the calculation to derive the SNR in a given coronagraphic mode for the PSF reference star is not shown. The reason behind this is as follows: the PSF reference star that should be used for any coronagraphic observation should be at least as bright as the main target—this way, by using an exposure time that is equal to the exposure time for the main target, the SNR of the PSF reference image will be at least as high as that of the main target's observation. Thus the reference PSF subtraction step will not contribute substantially to the total noise budget of the measurement.

NIRCam Coronagraphic Imaging Calculations

Main article: NIRCam Coronagraphic Imaging 
See also: HCI RoadmapJWST ETC BackgroundsJWST ETC Coronagraphy Strategy

For the NIRCam coronagraphic imaging observations, we subdivide the calculations in two groups: one using "Scene 1: Science - NIRCam" (with PSF reference star HD218261); the other using "Scene 3: Science - NIRCam bright" (with the brighter reference star HD218381). The underlying rationale for the setup of all the NIRCam coronagraphic observations in this article is to reach an integrated SNR of ~100.

Select Instrument Parameters

For all of the NIRCam coronagraphic imaging observations, the following applies: 

  • under the Detector Setup tab we set the Groups per Integration equal to "10", the Readout pattern to RAPID and the Subarray to SUB320 (nominal for LW coronagraphy with NIRCam). This readout pattern allows all the frames to be saved. Data volume is not an issue when using the small 320x320 subarray. Moreover, a 10 arcsec field of view is large enough to have HR8799b fit with a great margin. We also keep the Exposures per specification equal to "1" to simulate having no dithers for the science target observation (note that dithers are recommended for only PSF reference star coronagraphic observations). With all the parameters set in the Detector Setup tab defined, we use the Integrations per Exposure menu to achieve the desired SNR. With this readout we are using only ~8% of the saturation dynamic of the detector. For challenging (fainter) sources, it is recommended to use more but in this case, using this RAPID setup guarantees a SNR of ~100 on HR 8799 b in 1 to 10 minutes depending on the NIRCam filter as well as giving the flexibility to switch to a brighter star (calculation 9)  while keeping the same readout and saving all the frames.


  • Under the Strategy tab/ Observation sub-tab, we set the Scene rotation equal to "0" and the PSF subtraction source equal to either "2: HD218261" or "7: HD218381"—depending on which scene is being used. From the PSF Subtraction menu we select the Optimal (PSF autoscaling) option, allowing for the reference star to brightness to be re-scaled to that of the science target for PSF subtraction. Variations to this strategy are discussed below, under Advanced Strategies.


  • Under the Strategy tab/ Extraction sub-tab, we select "3: planet b" as the SNR Source; set the Aperture radius equal to "0.08" arcseconds; Sky annulusInner radius equal to "0.1" arcseconds and Sky annulusOuter Radius equal to "0.2" arcseconds.

Adjust Exposure Time to Achieve Desired SNR

Scene "1: Science - NIRCam"

This scene contains HR 8799, HR 8799b and reference star HD218261. There are 8 calculations defined for this scene corresponding to the 8 medium band filters supported for NIRCam coronagraphic observations in the LW channel using the MASKLWB occulter. The number of Integrations per Exposure were adjusted to reach SNR of ~100, as summarized in Table 3.

Table 3. NIRCam Coronagraphy calculation exposure parameters for scene ID: 1, Name: "Science - NIRCam"

Calc ID #Instrument Setup tabDetector Setup tab
CoronagraphFilterGroups per integration8Integrations per exposure9
1MASKLWBF250M108
2MASKLWBF300M107
3MASKLWBF335M105
4MASKLWBF360M104
5MASKLWBF410M105
6MASKLWBF430M1012
7MASKLWBF460M1060
8MASKLWBF480M1060

Table Notes:

8Number of groups: set to the maximum (10) as getting as many frames as possible is an advantage for post-processing in high contrast imaging. Also, in the event of partial saturation, one can recover the slope from the first few unsaturated groups. Finally, one can reject bad frames associated with an eventual cosmic ray.

9Integrations per Exposure chosen to reach SNR of ~100.

Scene "2: Science - NIRCam bright"

This science scene contains HR 8799, HR 8799b and reference star HD218381. This PSF reference star is brighter (4.41 Johnson K, Vegamag) than that used in Scene ID: 1, (5.14 Johnson K, Vegamag). We only define one calculation for this Scene, calculation ID #9, using the MASKLWB Coronagraph and F335M Filter. In the Detector Setup tab we set the number of Groups per integration to "10" and Integrations per exposure to "5". 

The Instrument and Detector Setup of calculation #9 is identical to that of calculation #3; however the selected PSF subtraction source in the Strategy tab is different. The SNR for this calculation is ~108, compared to ~106 for calculation #3. The small difference in SNR is due to the lower noise introduced in the subtraction step when using a brighter reference star.


MIRI Coronagraphic Imaging Calculations

See also:  MIRI Coronagraphic Imaging,  HCI RoadmapJWST ETC Backgrounds

For MIRI coronagraphic imaging we only use scene ID: 3, Name: "Science - MIRI" (whereby the PSF reference star target is ID: 6, Name: "Reference (HD220657)"). Similarly to the NIRCam case, the underlying rationale for the setup of all the MIRI coronagraphic observations in this example science program is to reach a SNR ~ 100.

Select Instrument Parameters

The following set-up is defined for all of the MIRI coronagraphic imaging calculations:

  • Under the Detector Setup tab, because of the brightness of our targets we keep the default Readout Pattern of FASTFAST1 mode provides short MULTIACCUM exposures to maximize dynamic range and minimize noise in a background-dominated regime. We set the number of Integrations per exposure equal to "5" and keep the Exposures per specification equal to "1"—this choice corresponds to having no dithers for the science target observation (note that dithers are instead recommended for PSF reference star coronagraphic observations). With all the other parameters set in the Detector Setup tab, we use the Groups per Integration menu to achieve the desired SNR.  We are following the MIRI Cross-Mode Recommended Strategies for bright stars. Since planet b is really easily detected and we are using less than 10% of the dynamic range of the detector (before saturation), it would be also ok to increase the number of integrations and keep the number of groups/integrations fixed as we did for NIRCam. Since MIRI allows to use a number of groups/integrations larger than 10, we prefer to raise this parameter and largely meet the flux calibrations requirements (not guaranteed with less than 5 groups/integrations). 
  • Under the Strategy tab/ Observation sub-tab, we set the Scene rotation equal to "0" and the PSF subtraction source equal to "5: Reference (HD220657)". From the PSF subtraction menu we select "Optimal (PSF autoscaling)", which allows for re-scaling of the reference star brightness to that of the science target. Variations to this strategy are discussed below, under Advanced Strategies.
  • Under the Strategy tab/ Extraction sub-tab, we set the SNR source to be "3: planet b"; the Aperture radius equal to "0.3" arcseconds, and Sky annulusInner radius equal to "0.45" arcseconds and Sky annulus Outer Radius equal to "0.7" arcseconds.

Adjust Exposure Time to Achieve Desired SNR

There are 3 MIRI coronagraphic imaging calculation in the ETC workbook which utilize the three different 4QPM masks. The Groups per Integration were increased to achieve a SNR of ~100, as summarized in Table 4.

Table 4. MIRI Coronagraphy calculation exposure parameters for scene ID: 3, Name: "Science - MIRI"

Calc ID #Instrument Setup:
Mask/Filter
Detector Setup:
Groups per integration1
104QPM/F1550C102
114QPM/F1140C27
124QPM/F1065C26

1Groups per Integration chosen to reach SNR of ~100.


Target Acquisition Calculations

See Also: NIRCam Coronagraphic Imaging Target AcquisitionJWST ETC NIRCam Target Acquisition, MIRI Coronagraphic Imaging Target Acquisition, JWST ETC NIRCam Target Acquisition

Two types of target acquisition calculations are used in the Example Science Program #36: one for the main science target, HR 8799 and its planet b, the other for the PSF reference stars.

Science Target TA

The TA calculations for the science targets in the Example Science Program #36 ETC workbook are configured to reach SNR of 100. With this goal in mind we adjust the Instrument and Detector Setup to reach approximately the target SNR.

NIRCam Target Acquisition on Science Target

For the NIRCam Target Acquisition calculation (Calculation ID13) we configure the scene to be "1: Science - NIRCam". This scene has 3 objects in it (Sources ID: 1,2,3, i.e. HR8799, HR8799b and HD218216, respectively). In the Strategy tab we need to specify that Aperture centered on source should be set to "1: HR 8799".

Note that setting the Aperture centered on source to the reference star HD218216 would result in an error: the PSF reference star in this scene has been placed at ~10 arcseconds from the main target; however the calculation domain for TA SNR in the ETC is smaller than this area and the calculation would fail—this is why we require Reference Scenes for TA on the PSF reference stars.

 Given the K=5.24 mag brightness of the source, in the Instrument Setup tab we set the Filter menu to "F335M + ND square (bright)" (see this page). In the Detector Setup tab we set the Groups per integration to "65", and the Readout Pattern to BRIGHT2. Note that these are the only 2 parameters that can be adjusted for this mode, and that only certain fixed values of NGROUPS are allowed. The resulting SNR for this calculation is ~112 while the minimum recommended value to guarantee a good acquisition and centering is 30.

MIRI Target Acquisition on Science Target

For the MIRI Target Acquisition calculation (Calculation IDs: 14, only 4QPM/1550 setup):

  • Under the Scene tab, we configure the Scene for calculation to be "3: Science - MIRI".
  • In the Strategy tab we need to specify that "Aperture centered on source" should be set to "1: HR 8799". 
  • Under the Instrument Setup tab, we set the filter to "FND (Neutral density)", required because of the source brightness (see this page). Given that we will be using the 4QPMs optical elements for our science observations, the Acq Mode menu needs to be set in turn to each one of the "4QPM" values, "TA for 4QPM/1065", "TA for 4QPM/1140", and "TA for 4QPM/1550".
  • In the Detector Setup tab, we set the Readout Mode to "FAST" and the number of Groups per Integration to "10". With MIRI, TA can be performed using two different readout modes: FAST and FASTGRPAVG. In general, users should consider using the FAST readout mode; however, for TA with the Lyot coronagraph concerning fainter stars (where longer than the shortest integration times are needed), FASTGRPAVG can be used. Given the brightness of our target source and use of the 4QPMs, the FAST readout mode is the appropriate choice.

    The resulting SNR is ~100 for all the 4QPM TA calculations.

Reference Targets TA

Similarly to the science target case, the TA calculations for the PSF reference stars in the Example Science Program #36 ETC workbook are configured to reach SNR of ~100.

  • Calculation ID 15: NIRCam Target Acquisition on reference star HD218261.
    For this calculation we configure the Scene to be "4: Reference - HD218261". This scene has 1 object in it (Source ID3, Name: "Reference (HD218216)"). Given the very similar brightness of the targets, the setup for this calculation is identical to calculation ID #16 (NIRCam TA on the science target). The resulting SNR is ~119.

  • Calculation ID 16: NIRCam Target Acquisition on reference star HD218381.
    For this calculation we configure the Scene to be "5: Reference - HD218381". This scene has 1 object in it (Source ID4, Name"Reference (HD218381)"). This brighter target is used to demonstrate possible trade-offs in the science strategy. The brighter target allows much faster Reference PSF observations (see the APT page for this Example Science Program). The gain in execution time for TA is minimal given that TA observations are already quite fast. In order to reach SNR ~100 for this brighter target, the calculation setup is identical to that of Calculation ID #17, except for the Readout pattern in the Detector Setup tab, which is set to RAPID since this star is brighter than HR 8799 The resulting SNR is ~102. 


  • Calculation ID 17: MIRI Target Acquisition on reference star HD220657.
    For this calculation we configure the Scene to be "6: Reference - HD220657". This scene has 1 object in it (Source ID2, Name:
    "Reference (HD220657)"). This source is brighter than HR8799 and thus we set number of Groups per integration to "4" in the Detector Setup panel. 4 is the minimum recommended number of groups for MIRI TA—a lower number would not allow a robust fit of the accumulated charge vs. time slope. Even with the smallest number of reads, the resulting SNRs are ~201 - 202 when selecting the Acq Mode to be TA for 4QPM/1065, TA for 4QPM/1140, and TA for 4QPM/1550, under the Instrument Setup tab.



Advanced strategies

The above example has been simplified for the purpose of the training. Since the objective of the program is to study the atmosphere of the exoplanet HR 8799 b, one should in principle use a synthetic model spectrum which has been optimized with existing observations and which account for the presence of an atmosphere rather than a perfect black body. For GTO program 1194, the team has calculated the corresponding integrated flux in each of the NIRCam and MIRI filters using the most up-to-date models. Nevertheless, Figure 1 below shows that using a synthetic model spectrum with clouds (generated by Travis Barman as in Rajan et al. 2015), the SNR previously adjusted to ~100 for a 1100K black body spectrum normalized at K=14.25 (Vega mag) is comprised between ~50 an ~110 using a ~1100K cloudy synthetic model spectrum also at 1100K but normalized at K=14.05 to preserve the same energy and match the mid-IR tail.

Because HR8799 b can be easily detected and characterized  with JWST, the approximation is fair.

For a more challenging object, fainter and/or closer to the coronagraph's IWA (typically inside 1 arcsecond for these coronagraphs), it is recommended to upload a realistic model spectrum and be careful about its normalization, using previous photometric points, often taken from the ground with different filters or filter systems.

Figure 1. HR 8799 b's SNR (ETC Plots) across wavelength for both NIRCam and MIRI. 


Left: the input spectrum is a 1100 K black body and the SNR is roughly constant at ~100. Right: the input spectrum is now a modeled synthetic spectrum taken from Rajan et al. 2015 (Teff = 1100K, log(g) = 3.00, radius = 0.90 RJup) with a similar normalization. With the same readout parameters, the SNR now varies from ~50 to ~110 for NIRCam and is still close to ~100 for MIRI (because our normalization makes the long wavelengths match well as shown on the top-right plot). 


With the exposure parameters now determined for this program, we can populate the observation template in APT. See the Step-by-Step APT Guide to complete the proposal preparation for this example science program. 


Links

SIMBAD entry for HR8799

SIMBAD entry for HD218261

SIMBAD entry for HD218381

SIMBAD entry for HD220657


References

A. Rajan, T. S. Barman, R. Soummer et al. arXiv:1508.02395
Characterizing the Atmospheres of the HR8799 Planets with HST/WFC3