Step-by-Step ETC Guide for MIRI MRS Spectroscopy of a Late M Star

A walk-through of the JWST Exposure Time Calculator (ETC) for the MIRI MRS Spectroscopy of a Late M Star example science program is provided, demonstrating how to select exposure parameters for this observing program.

Dated material

This example was created pre-launch, and the ETC has been updated since its creation. You may see differences in the details of the results from the ETC, the information provided, or the appearance of the ETC GUI from what is shown herein.

Please refer to JWST Example Science Programs for more information.

On this page

See also: MIRI MRS Spectroscopy of a Late M Star, JWST ETC Exposure Time Calculator OverviewProposal Planning Video Tutorials

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

For the MIRI MRS Spectroscopy of a Late M Star example program, we focus on selecting exposure parameters for MIRI MRS. We start by defining a scene with a point source that is a late M star. We show how to run ETC calculations to achieve the desired SNR for MIRI MRS. 

The optimal exposure specifications (e.g., the number of groups and integrations) are the inputs needed for the Astronomer's Proposal Tool (APT) MIRI MRS APT Template, which is used to specify the observing program.

The ETC workbook associated with this example science program is called #28: MIRI MRS Spectroscopy of a Late M Star and can be selected from the Example Science Program Workbooks drop-down tab in the ETC workbooks page. The nomenclature and reported SNR values are based on ETC 1.5. There may be subtle differences if using a different version of ETC.



Define scenes and sources in the ETC

See also: JWST ETC Defining a New Source, JWST ETC Source Spectral Energy Distribution

Words in bold are GUI menus/
panels or data software packages; 
bold italics are buttons in GUI
tools or package parameters.

We first set up a scene with a single late M star. We define the following source in ETC:

Star: a point source star with a spectrum described by a PHOENIX Stellar Model for a M5V star of effective temperature 3,500 K and log(g) = 5.0. The star has zero redshift and zero extinction, and it is normalized to be 100 mJy at 10.0 μm wavelength.

Note that the position of the source in the scene can be viewed in the lower left Scene Sketch pane. By selecting the checkbox in the Plot column in the Select a Source pane, the spectral energy distribution (SED) of the selected source can be plotted. (Note: it may be helpful to limit the wavelength axis to the range relevant to the MIRI MRS mode, 5–28 µm).  This source is named Late M Star in the Select a Source pane and assigned to Late M Star Scene in the Select a Scene pane.

Define sources for MIRI simultaneous imaging

We assume the same type of star used as the MRS target for simultaneous imaging as late M stars are common stars. The only difference is that, under the Renorm tab in the Source Editor pane, at the Normalize in Bandpass parameter, we set Johnson K-band magnitude to 11.5 (vegamag). This source is named Imager Star in the Select a Source pane and assigned to Simultaneous Imaging in the Select a Scene pane.



Select MIRI medium resolution spectroscopy (MRS) Calculation

See also: JWST ETC Creating a New CalculationJWST ETC Backgrounds

MIRI MRS observations are taken in each of 3 grating settings: SHORT (A), MEDIUM (B), and LONG (C). For the observation at each MRS grating setting, all 4 channels are observed simultaneously. By obtaining the observations at each of the 3 grating settings, you can obtain the complete 5–28 μm spectrum.

Our goal is to detect the late M star at a SNR > 25 for all wavelengths shortward of 26.0 μm, so we run ETC Calculations for MIRI MRS spectroscopy for all 3 grating settings of all 4 MRS channels to determine the exposure parameters we need to achieve this SNR.

Since the JWST background is position dependent, fully specifying background parameters is important for the most accurate SNR Calculation. After selecting the Calculation in the Calculations pane, the adjacent Calculation editor pane becomes active: Under the Background tab, select Medium for Background configuration, which corresponds to the 50th percentile of the sky background. Though we don't provide RA and Dec for a specific star, we provide generic RA and Dec coordinates of 04:30:00 and +25:00:00, respectively.



Select instrument parameters

See also: JWST ETC IFU Strategies

Here, we set up the ETC workbook for MRS Calculations by editing parameters in the Calculations editor pane (adjacent to the Calculations pane):

  • Under the Instrument Setup tab:
    Calculation #1 shows our initial Calculation to assess the SNR for Band 3C. Here, Channel is set to Channel 3 (MRS_Long), and Wavelength Range is set to Long (C) for Band 3C.  For calculations for other bands, one would choose the appropriate combination of Channel and Wavelength Range.

  • Under the Detector Setup tab:
    • Subarray is set to FULL (only full frame readout is supported for the MRS).
    • Number of Exposures per specification is set to 2 because the IFU Nod In Scene option chosen under the Strategy tab assumes a 2-point dither. By so doing, we simulate a 4-point dither in ETC.
    • We chose the Readout Pattern as FAST.

  • Under the Strategy tab:
    • We selected the IFU Nod in Scene option. This allows for background subtraction using the other dither in the 2-point dither (assumed by ETC for IFU Nod in Scene). The source is detected in both dithers, and the IFU Nod in Scene option assumes the source detected in both dithers is co-added. This makes it easier to simulate the 4-point dither pattern (i.e., by setting Exposures per specification to 2, as noted above).
    • Aperture location was set to Centered on source with the only choice in the drop-down menu being "Late M Star".
    • Aperture radius was set to
      • 0.5 (arcsec) for channel 1 (Calculation # 8, 10, and 12) and channel 2 (Calculations # 7, 9, and 11);
      • 0.75 (arcsec) for channel 3 (Calculations # 1, 2, and 4) and channel 4 (Calculations # 3, 5, and 6). The radius is larger for the longer wavelength channels because the PSF grows larger at longer wavelengths.
    • Nod Position in Scene was set to the values listed below. 
      As with the aperture radius, we increase the x and y coordinates of the nod position in scene with wavelength, in order to account for the larger PSF at longer wavelengths.
      • 1 for X and (arcsec) for channel 1 (Calculations # 8, 10, and 12);
      • 1.5 for X and Y (arcsec) for channel 2 (Calculations # 7, 9, and 11);
      • 1.75 for X and Y (arcsec) for channel 3 (Calculations # 1, 2, and 4);
      • 2.25 for X and Y (arcsec) for channel 4 (Calculations # 3, 5, and 6).



Adjust exposure parameters to obtain desired signal-to-noise ratio

See also: JWST ETC Batch Expansions

Here we assess the SNR by editing parameters in the Calculations editor pane (adjacent to the Calculations pane):

  • Under the Detector Setup tab:
    • We chose the Readout Pattern as FAST because the same SNR could be achieved with slightly less total exposure time than for the SLOW readout pattern. We tested this by fixing the number of exposures and groups, keeping the same Subarray (FULL), and finding the number of groups needed for each of SLOW and FAST readout pattern to achieve the same SNR at the longest wavelength of the band.

    • Because the stellar photosphere spectrum monotonically decreases with wavelength, in general, we begin our ETC experiment by determining exposure parameters for each grating setting of channels 3 and 4, since we want SNR > 25 for all 3 grating settings of channels 1, 2, and 3 and the "Short" and "Medium" grating settings of channel 4, and SNR > 25 for all wavelengths <26 μm in channel 4 with the "Long" grating setting. Attaining this minimal SNR for channels 3 and 4 first, we should obtain even greater SNR for channels 1 and 2 because of the shape of the stellar photosphere spectrum. The number of integrations and groups for channels 1 and 2 of the MRS do not have to be the same as those for MRS channels 3 and 4. The only requirement on this is that the total exposure time should be the same for channels 1/2 versus channels 3/4. This is useful in this example science program because the star is so much brighter at shorter wavelengths in the MRS than longer wavelengths. The ETC shows that, when the same number of groups and integrations are used for all 4 MRS channels, pixels begin to saturate in MRS channel 1 when the number of groups is set so that SNR = 25 at 26.0 μm wavelength in channel 4 is achieved. By using more integrations and less groups for channels 1/2 than for channels 3/4 of MRS, saturation at shorter wavelengths can be avoided.

    • We wanted to avoid saturation warnings in the ETC, though we note that for MIRI, as long as saturation does not happen within 5 groups in an integration, the data is still usable even if the groups past the 5th group are saturated (for more, please see the MIRI Cross-Mode Recommended Strategies article).  So we took the following approach:
      • First, we set the Integrations per exposure to 1 to see if the desired SNR could be achieved by increasing the Groups per integration but without returning any saturation warnings in ETC.
      • For a given grating setting, if we could set the Groups per integration for all 4 MRS channels to achieve the desired S/N for Integrations per exposure equal to 1, without the ETC returning a saturation warning (though the end of the ramp being saturated does not mean the data is unusable; please see above), then we used these Groups per integration and Integrations per exposure values. If not, then, if the saturation happened in the shorter wavelength channels, we increased the number of Integrations per exposure for channels 1/2 by 1 and tried again (leaving Integrations per exposure and Groups per integration for channels 3/4 alone).
      • This iterative process found that, by setting the number of Integrations per exposure to 1, the number of Groups per integration needed to be 27 for the SHORT (A) grating setting (Calculations # 3, 4, 11, and 12) and 59 for the MEDIUM (B) grating setting (Calculations # 2, 5, 9, and 10) to achieve the desired S/N. For the LONG (C) grating setting (Calculation # 1, 6, 7, and 8), by setting the Integrations per exposure to 2 and the Groups per integration to 49 for channels 1/2 and Integrations per exposure to 1 and Groups per integration to 98 for channels 3/4, we were able to avoid saturation warnings and achieve the desired S/N > 25 for wavelengths < 26 microns.

Table 1 shows the exposure parameters used in the ETC for the Late M Star program found from the above-described experimentation in ETC.


Table 1. Exposure parameters used in the ETC for the Late M Star example program 

BandCalculation ID #

Exposures per

specification

Integrations

per exposure

Groups per

integration

1A122127
1B102159
1C82249
2A112127
2B92159
2C72249
3A42127
3B22159
3C12198
4A32127
4B52159
4C62198

Select MIRI simultaneous imaging Calculation

Along with the MIRI MRS observations, we include simultaneous imaging with the MIRI imager.  We wish to determine the brightest star that could be observed with the longest exposure time setting of the imager without saturation. We use the F1130W filter, as this allows observing the brightest stars without saturation, since most stars will not have significant excess emission above that from the stellar photosphere.

Similar to the MIRI MRS Calculations, we choose Medium for Background Configuration.

Select instrument parameters

Here we assess the SNR in Calculation #14, by editing parameters in the Calculations editor pane:

  • Under the Instrument Setup tab, set the Filter parameter to F1130W

  • Under the Detector Setup tab:
    • Subarray is set to FULL;
    • Readout pattern is set to FAST because it needs to be the same readout pattern as the MRS observations;
    • Exposures per specification is set to 1 because we are only interested in determining the magnitude of the field star that saturates for the same exposure time as the MRS observation
    • Integrations per exposure is set to 1 (matching the MRS channel 3/4 value for its Long" grating setting)
    • Groups per integration is set to 98 (again, matching the MRS channel 3/4 value for its "Long" grating setting)

  • Under the Strategy tab:
    • Imaging Aperture Photometry is the only option at the top, so it is chosen;
    • Under Aperture location, the Specify offsets in scene is selected, with X and Y set to 0;
    • Aperture radius is set to 0.3 (arcsec), which is the default value for MIRI imaging in ETC.
    • For Sky annulus, the inner radius and outer radius are set to 0.45 and 0.7 (arcsec), respectively.

The magnitude under the Renorm tab of the Source Editor is then varied until the Calculation returns a saturation warning. This happens for Johnson K magnitude less than 11.5 magnitudes. So K = 11.5 is the saturation limit for the simultaneous imaging obtained for the "Long" grating setting. Somewhat brighter stars could therefore be observed in simultaneous imaging without saturation for the "Medium" and "Short" grating settings, for which the Groups per integration is less.


Target acquisition Calculation

See also: MIRI MRS Target Acquisition, JWST ETC MIRI Target Acquisition

A target aquisition (TA) is required in order to position the target precisely on the detector.  TA for the MRS improves pointing precision to 90 mas.

Calculation #13, where we selected Target Acquisition under the Strategy tab in the Calculations editor pane, shows our initial calculation to determine which parameters to specify for TA for the target star.

  • Under the Backgrounds tab, we entered the coordinates for the generic star (04:30:00 +25:00:00) and selected Medium for Background configuration;

  • Under the Instrument Setup tab, we updated the selection of Acq Mode to TA for LRS-Slit, MRS and LRS-Slitless. For Filter, we chose FND (Neutral Density);

  • Under the Detector Setup tab, the only option for Subarray for MRS TA is FULL.  We choose Readout pattern to be FAST, as the FASTGRPAVG pattern results in longer exposure times. The only option for Integrations and Exposures are, for each, 1. We set the number of Groups to the relatively low value of 6;

  • Under the Strategy tab, the only field under this tab for target acquisition is Aperture centered on source, which for our case is set to the Late M Star.

From this calculation, we see the SNR is ~ 175.



Run ETC Calculations

Running the Calculations with these parameters gives SNR of 25 at the longest wavelengths of Bands 4A and 4B and at 26.0 μm in Band 4C. At wavelengths shortward of these wavelengths in Bands 4A, 4B, and 4C, the SNR is greater than 25, and at all wavelengths in Bands 1A, 1B, 1C, 2A, 2B, 2C, 3A, 3B, and 3C, the SNR is much greater than 25. These SNR values are provided graphically in the plot displayed under the 2D SNR tab in the Images pane at the lower left of the screen when the Calculations tab of the workbook is selected.

With the exposure parameters now determined for this program, we can populate the observation template in APT.  See the Step-by-Step APT Guide for MIRI MRS Spectroscopy of a Late M Star to complete the proposal preparation for this example science program.




Latest updates
Originally published