pith. sign in

arxiv: 2512.12434 · v2 · submitted 2025-12-13 · 🌌 astro-ph.IM

SIMLA: The Spitzer Infrared Spectrograph Mapping Legacy Archive

Grant P. Donnelly , Cory M. Whitcomb , Lindsey Hands , Sara E. Duval , J.-D. T. Smith , Karin Sandstrom , David Carroll , McKenna Dowd
show 4 more authors
Brandon S. Hensley Leslie K. Hunt Edward Walsh Julie Watson
This is my paper

Pith reviewed 2026-05-16 22:30 UTC · model grok-4.3

classification 🌌 astro-ph.IM
keywords Spitzer IRSmid-infrared spectroscopyspectral cubesmapping observationsdata archiveinfrared astronomydata reduction
0
0 comments X

The pith

The Spitzer/IRS Mapping Legacy Archive supplies complete mid-infrared spectral cubes for hundreds of objects after novel pipeline processing.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper presents the Spitzer/IRS Mapping Legacy Archive (SIMLA), a collection of spectral cubes derived from low-resolution mapping-mode observations by the Spitzer Infrared Spectrograph spanning 5.2 to 38 microns at resolutions of roughly 60 to 130. The authors apply a custom pipeline to subtract foregrounds, backgrounds, and detector artifacts, using the CUBISM code to assemble the cubes and automatically flag bad pixels. Validation comes from showing that photometry synthesized from the cubes typically matches independent WISE measurements to within a few percent. The archive covers galaxies, molecular clouds, supernova remnants, HII regions, and other targets, with products slated for public release at the NASA/IPAC Infrared Science Archive.

Core claim

We present the Spitzer/IRS Mapping Legacy Archive (SIMLA); a complete set of mid-infrared spectral cubes built from low-resolution mapping-mode fixed-target observations from Spitzer/IRS (5.2-38 micron, R~60-130). Each cube has been carefully treated to remove astronomical foregrounds and backgrounds as well as detector effects using a novel pipeline. Cube assembly was facilitated by the CUBISM code, which included automatic detection and removal of bad pixels.

What carries the argument

The novel pipeline that removes astronomical foregrounds, backgrounds, and detector effects while assembling cubes with the CUBISM code for automatic bad-pixel detection and removal.

Load-bearing premise

The novel pipeline correctly removes astronomical foregrounds, backgrounds, and detector effects without introducing artifacts that would affect scientific use.

What would settle it

Independent verification showing that synthetic photometry extracted from the SIMLA cubes deviates systematically from WISE photometry by more than a few percent, or the presence of visible processing artifacts in the released spectral cubes.

Figures

Figures reproduced from arXiv: 2512.12434 by Brandon S. Hensley, Cory M. Whitcomb, David Carroll, Edward Walsh, Grant P. Donnelly, J.-D. T. Smith, Julie Watson, Karin Sandstrom, Leslie K. Hunt, Lindsey Hands, McKenna Dowd, Sara E. Duval.

Figure 1
Figure 1. Figure 1: The distribution of SIMLA maps on-sky. Red and blue squares indicate SL and LL cubes, respectively. The size of these markers do not represent the real fields of view, but are logarithmically scaled to the covered fields of view with an arbitrary scale factor applied to all. In this IRS sample, there is a total area on sky of 2.94 deg2 . The total unique area covered by SL is 0.9 deg2 and by LL, 2.34 deg2 … view at source ↗
Figure 2
Figure 2. Figure 2: Three-color images made using mosaics of SIMLA spectral cubes of the Eagle Nebula (top), the Antennae galaxies (bottom left), and the Cassiopeia A supernova remnant (bottom right). For each image, the pixel values of each color correspond to the slice of a SIMLA cube mosaic at the listed wavelength [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Example spectra from SIMLA spectral cubes of the Galactic Hii region G43 (blue), the star-forming galaxy M101 (orange), the luminous infrared galaxy NGC 7469 (green), and the binary Herbig Ae stars PDS144 (red). 0 20 40 60 80 100 120 Pixel x 0 20 40 60 80 100 120 Pixel y Inter order SL1 SL2 SL3 Red PU Blue PU 0 20 40 60 80 100 120 Pixel x 0 20 40 60 80 100 120 Pixel y Inter order LL1 LL2 LL3 S p e c t r al… view at source ↗
Figure 4
Figure 4. Figure 4: Basic calibrated data (BCD) example images for Spitzer/IRS with major sections labeled. The left and right panels show BCDs for the short-low (SL) and long-low (LL) modules, respectively. For both, the red and blue sections highlight the wavesamp regions for the sub-slits within either module. The green sections are the “bonus orders” (SL3, LL3), on which light from the 2nd slit (SL2, LL2) is sampled withi… view at source ↗
Figure 5
Figure 5. Figure 5: Schematic of the SIMLA pipeline to produce a cube. Blue boxes indicate SIMLA intermediate products derived from more basic inputs shown as yellow rounded boxes. Green hexagons are operations. Note that most of the SIMLA pipeline is oriented towards producing backgrounds. but nearby in time and in sky location. Thus, both the pixel behavior and the zodiacal emission captured by the background would be simil… view at source ↗
Figure 6
Figure 6. Figure 6: Examples of the three main components of SIMLA backgrounds are shown for one SL AOR: the 2D zodiacal emission model adjusted for the CVZ spectrum is in the top left (Section 3.1.1), the interpolated baseline frame scaled for the zodiacal intensity of this AOR is in the top right (Section 3.1.2), and the average of the stack of dark shards is in the bottom left (Section 3.2). In the bottom right panel, the … view at source ↗
Figure 7
Figure 7. Figure 7: Observations at the continuous viewing zone (CVZ) track the time-dependent offsets in the zodiacal emission model. Here, we show these offsets as the 28.75 µm intensities (center of LL1) of observations at the CVZ for one RAMPTIME as colored hexagons. The orange curve shows the modeled zodiacal intensity at 28.75 µm of a given time at the CVZ minus the modeled mean intensity at the CVZ. This pattern oscill… view at source ↗
Figure 8
Figure 8. Figure 8: Example spectrum for a typical off-source ob￾servation with a decomposition of the CVZ-adjusted zodi￾acal emission model. The black spectrum is extracted off of the baseline frame-subtracted (see Section 3.1.2) median of all BCDs within a dark AOR that only contains zodia￾cal emission. The dashed red line is the unchanged zodiacal emission model for the AOR (Zk(λ)), and the dotted red line is the scaled CV… view at source ↗
Figure 9
Figure 9. Figure 9: Top: SL binned baseline frames for a certain RAMPTIME (60.95 s), with each panel in the row corresponding to a different zodiacal emission intensity. The baseline frames capture the pixel excess artifact/dark current. As the intensity increases, there are more negative (darker) pixels as the zodiacal emission model tends to overestimate. Bottom: spectra extracted from each of the above binned baseline fram… view at source ↗
Figure 11
Figure 11. Figure 11: Shard selection diagram for a SL AOR on NGC 1512. Colored rectangles show the sky positions of shards on a WISE image, with the colors indicating whether they passed either the WISE or BCD spectrum cut, both, or neither. Only green shards qualify for use as part of a background. The full color key is given in the top left, and the cut values in the lower right. In this example, there are no blue-colored c… view at source ↗
Figure 10
Figure 10. Figure 10: Top: WISE image showing an example of the LL slit positions during an observation of a star. The red rectangles show the full sky positions for the LL1 (upper) and LL2 (lower) slits. The blue rectangles illustrate the positions of the shards on the sky. Note that the shards do not span the full length because of the edge trimming. Bottom: the corresponding BCD for the observation shown above. The thick bl… view at source ↗
Figure 12
Figure 12. Figure 12: Results from the experiment described in Sec￾tion 3.2.4 to find the best values for ∆zmax and dtarg. The pixel values in both panels show the median σbgsub/σ0 across all test images (see text) for the combination of cuts. The top panel shows a “coarse” parameter search, and the bot￾tom panel shows a “fine” search localized around the best region in the top panel, represented by the red dashed box. In both… view at source ↗
Figure 13
Figure 13. Figure 13: A two-dimensional histogram of the spectra extracted from dark regions in SIMLA cubes (see Section 4.1). The color bar on the right shows the density of spectral points. The vertical breaks at 7.7, 14.5, and 22 µm align with the edges of spectral segments (SL2 and SL1, etc). By far the highest densities of SIMLA-derived surface brightnesses for dark regions are near 0 MJy sr−1 (black dashed line), as desi… view at source ↗
Figure 14
Figure 14. Figure 14: Synthetic WISE W3 photometry derived from spectra extracted from 5′′ radius apertures from SIMLA cubes, compared with matched extractions of WISE W3 pho￾tometry. All apertures in this figure have full wavelength coverage. The black line indicates the identity line, and the green line indicates the median (1.19) of points with observed brightness greater than 0.1 MJy sr−1 . spectra were extracted from SIML… view at source ↗
Figure 15
Figure 15. Figure 15: The difference between spectra with no correction for the inter-order (IO) signal (left) and spectra with a correction (right). The correction is only applied to SL cubes. Each spectrum is the median of extractions from SIMLA cubes binned by surface brightness; the gray regions in the left panel show the bin ranges, and the bins are the same in both panels. The spectra have not been stitched together. Col… view at source ↗
Figure 17
Figure 17. Figure 17: Summary of noise levels in the dark spectra extracted from SIMLA cubes (see Section 4.1), separated by RAMPTIME and suborder. The noise is quantified as the standard deviation of emission-free spectra. The box heights show the inter-quartile ranges, the orange lines represent the median values, and the whiskers indicate 1.5 times the in￾ter-quartile range. by CUBISM. This can result in striping artifacts … view at source ↗
read the original abstract

We present the Spitzer/IRS Mapping Legacy Archive (SIMLA); a complete set of mid-infrared spectral cubes built from low-resolution mapping-mode fixed-target observations from Spitzer/IRS (5.2-38 micron, R~60-130). Contained in this dataset are spectral maps for several hundred spatially-resolved and unresolved objects, including galaxies, molecular clouds, supernova remnants, HII regions, and more. Each cube has been carefully treated to remove astronomical foregrounds and backgrounds as well as detector effects using a novel pipeline. Cube assembly was facilitated by the CUBISM code, which included automatic detection and removal of bad pixels. We describe the SIMLA pipeline for reducing and validating the cubes, and we show that synthetic photometry derived from SIMLA spectra and corresponding WISE photometry typically agree within a few percent. SIMLA products and documentation related to their use will soon be available at the NASA/IPAC Infrared Science Archive (DOI:10.26131/IRSA655).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript presents the Spitzer/IRS Mapping Legacy Archive (SIMLA), a complete set of mid-infrared spectral cubes (5.2-38 μm, R~60-130) constructed from low-resolution mapping-mode fixed-target observations of several hundred objects including galaxies, molecular clouds, supernova remnants, and HII regions. A novel pipeline removes astronomical foregrounds, backgrounds, and detector effects, with cube assembly performed via the CUBISM code that includes automatic bad-pixel detection and removal. Validation consists of showing that synthetic photometry extracted from the cubes agrees with independent WISE photometry to within a few percent. The products are intended for public release via the NASA/IPAC Infrared Science Archive.

Significance. If the pipeline demonstrably produces spectral cubes free of processing artifacts at the level needed for scientific analysis, SIMLA would constitute a substantial legacy resource for mid-IR spectroscopy. It would enable uniform, ready-to-use access to Spitzer/IRS mapping data across a wide range of astrophysical targets, reducing the barrier for studies that require spatially resolved spectra without requiring each user to re-reduce the raw data.

major comments (3)
  1. [Abstract / Validation] The only quantitative validation provided is the few-percent agreement between synthetic photometry from the cubes and WISE broadband fluxes. Because WISE supplies only four integrated measurements, this test cannot detect narrow spectral residuals, continuum tilts, or localized bad-pixel artifacts that would remain visible in the R~60-130 cubes; a more sensitive validation (e.g., direct comparison to higher-resolution spectra or residual maps) is required to support the claim that the novel pipeline removes foregrounds, backgrounds, and detector effects without introducing artifacts.
  2. [Pipeline Description] No detailed error budget, uncertainty propagation, or quantitative residual assessment is presented for the processed cubes. Without these, it is impossible to evaluate whether the cubes meet the precision needed for typical scientific applications such as line-flux measurements or continuum fitting on faint sources.
  3. [Data Selection] The criteria used to select or exclude individual mapping observations from the final archive are not specified. This omission affects the claimed completeness of the dataset and the reproducibility of the sample for future users.
minor comments (2)
  1. [Abstract] Clarify whether 'spectral cubes' and 'spectral maps' are used interchangeably throughout the text.
  2. [Abstract] Confirm the final DOI (10.26131/IRSA655) and provide a direct link or accession number once the archive is deposited.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their detailed and constructive review of our manuscript on the SIMLA archive. We address each major comment below and have revised the manuscript to strengthen the validation, documentation, and reproducibility aspects.

read point-by-point responses

  1. Referee: [Abstract / Validation] The only quantitative validation provided is the few-percent agreement between synthetic photometry from the cubes and WISE broadband fluxes. Because WISE supplies only four integrated measurements, this test cannot detect narrow spectral residuals, continuum tilts, or localized bad-pixel artifacts that would remain visible in the R~60-130 cubes; a more sensitive validation (e.g., direct comparison to higher-resolution spectra or residual maps) is required to support the claim that the novel pipeline removes foregrounds, backgrounds, and detector effects without introducing artifacts.

    Authors: We agree that broadband WISE photometry alone provides limited sensitivity to narrow spectral features or localized artifacts. In the revised manuscript we have added direct comparisons of extracted spectra from SIMLA cubes against archival high-resolution IRS staring-mode observations for a representative subset of targets (galaxies, HII regions, and supernova remnants). These comparisons confirm that line ratios, continuum slopes, and overall spectral shapes are preserved to within the expected calibration uncertainties. We have also included example residual maps after background subtraction to illustrate the level of foreground removal achieved. revision: yes

  2. Referee: [Pipeline Description] No detailed error budget, uncertainty propagation, or quantitative residual assessment is presented for the processed cubes. Without these, it is impossible to evaluate whether the cubes meet the precision needed for typical scientific applications such as line-flux measurements or continuum fitting on faint sources.

    Authors: We acknowledge the value of a quantitative error budget. The revised manuscript now contains a dedicated subsection describing the uncertainty contributions from each pipeline stage (background modeling, bad-pixel flagging, and cube assembly). We propagate these uncertainties into per-voxel error maps that are delivered alongside the science cubes. Typical 1-sigma uncertainties for continuum and line measurements are quantified for representative source brightness levels. revision: yes

  3. Referee: [Data Selection] The criteria used to select or exclude individual mapping observations from the final archive are not specified. This omission affects the claimed completeness of the dataset and the reproducibility of the sample for future users.

    Authors: We thank the referee for highlighting this omission. The revised Section 2 now explicitly states the selection criteria: all publicly available low-resolution mapping-mode observations of fixed targets as of the archive cutoff date, excluding only those with insufficient integration time, severe telemetry gaps, or known instrument anomalies. We report the total number of observations examined and the fraction retained, along with a table summarizing excluded cases and the rationale for each exclusion. revision: yes

Circularity Check

0 steps flagged

No circularity: validation uses independent external WISE photometry

full rationale

The paper is a data-release description of spectral cubes produced by a novel reduction pipeline applied to Spitzer/IRS mapping observations. The only quantitative validation step is a comparison of synthetic broadband photometry extracted from the cubes against independent WISE photometry, which lies outside the pipeline inputs and is not used to fit any parameter. No equations, self-definitions, or self-citation chains are presented that reduce any claimed result to a fitted quantity defined by the same data. The central product (the archive cubes) is therefore not forced by construction from its own validation metric.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the novel pipeline and CUBISM code produce scientifically usable cubes; no free parameters or invented entities are quantified in the abstract, and the only explicit axiom is the reliability of automatic bad-pixel removal.

axioms (1)
  • domain assumption CUBISM code accurately assembles cubes and removes bad pixels automatically from Spitzer/IRS mapping data.
    Invoked in the description of cube assembly and bad-pixel handling.

pith-pipeline@v0.9.0 · 5508 in / 1232 out tokens · 55065 ms · 2026-05-16T22:30:36.217820+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Nowhere Left to Hide: Uncovering All of the Massive Young Embedded Star Clusters in the Antennae with JWST

    astro-ph.GA 2026-05 unverdicted novelty 7.0

    JWST observations identify all massive young embedded star clusters in the Antennae, revealing they are extremely young, heavily obscured, and account for ~60% of the ionizing luminosity.

  2. JWST Observations of Starbursts: Dust Processing in the M82 Superwind

    astro-ph.GA 2026-04 conditional novelty 7.0

    PAH abundance remains constant at ~1% throughout the M82 superwind to 5 kpc, indicating shielding within cool cloud surfaces rather than destruction by the hot outflow.

  3. JWST Observations of Starbursts: Dust Processing in the M82 Superwind

    astro-ph.GA 2026-04 conditional novelty 7.0

    JWST observations find constant ~1% PAH abundance in the M82 superwind to 5 kpc, consistent with shielding in surface layers of cool clouds and possible replenishment.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages · cited by 2 Pith papers · 1 internal anchor

  1. [1]

    Armus, L., Charmandaris, V., & Soifer, B. T. 2020, Nature Astronomy, 4, 467, doi: 10.1038/s41550-020-1106-3

    work page doi:10.1038/s41550-020-1106-3 2020

  2. [2]

    D., Allamandola, L

    Boersma, C., Bregman, J. D., Allamandola, L. J., Temi, P., & Maragkoudakis, A. 2024, ApJ, 975, 177, doi: 10.3847/1538-4357/ad7d08

    work page doi:10.3847/1538-4357/ad7d08 2024

  3. [3]

    Genzel, R., & Cesarsky, C. J. 2000, ARA&A, 38, 761, doi: 10.1146/annurev.astro.38.1.761

    work page doi:10.1146/annurev.astro.38.1.761 2000

  4. [4]

    , author =

    Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18, doi: 10.1086/423134 IRS Instrument Team, & Science User Support Team. 2021, IRS Instrument Handbook, NASA IPAC DataSet, IRSA487, doi: 10.26131/IRSA487

    work page doi:10.1086/423134 2004

  5. [5]

    L., Franz, B

    Kelsall, T., Weiland, J. L., Franz, B. A., et al. 1998, ApJ, 508, 44, doi: 10.1086/306380

    work page doi:10.1086/306380 1998

  6. [6]

    F., Steinz, J

    Kessler, M. F., Steinz, J. A., Anderegg, M. E., et al. 1996, A&A, 315, L27

    work page 1996

  7. [7]

    E., Glaccum, W

    Krick, J. E., Glaccum, W. J., Carey, S. J., et al. 2012, ApJ, 754, 53, doi: 10.1088/0004-637X/754/1/53

    work page doi:10.1088/0004-637x/754/1/53 2012

  8. [8]

    Law, D. R., E. Morrison, J., Argyriou, I., et al. 2023, AJ, 166, 45, doi: 10.3847/1538-3881/acdddc

    work page doi:10.3847/1538-3881/acdddc 2023

  9. [9]

    2020, Nature Astronomy, 4, 339, doi: 10.1038/s41550-020-1051-1

    Li, A. 2020, Nature Astronomy, 4, 339, doi: 10.1038/s41550-020-1051-1

    work page doi:10.1038/s41550-020-1051-1 2020

  10. [10]

    2011, title Preliminary Results from NEOWISE: An Enhancement to the Wide-field Infrared Survey Explorer for Solar System Science , , 731, 53, 10.1088/0004-637X/731/1/53

    Mainzer, A., Bauer, J., Grav, T., et al. 2011, ApJ, 731, 53, doi: 10.1088/0004-637X/731/1/53

    work page doi:10.1088/0004-637x/731/1/53 2011

  11. [11]

    T., Morris, P., Boulanger, F., & Okumura, K

    Reach, W. T., Morris, P., Boulanger, F., & Okumura, K. 2003, Icarus, 164, 384, doi: 10.1016/S0019-1035(03)00133-7

    work page doi:10.1016/s0019-1035(03)00133-7 2003

  12. [12]

    M., Bolatto, A

    Sandstrom, K. M., Bolatto, A. D., Bot, C., et al. 2012, ApJ, 744, 20, doi: 10.1088/0004-637X/744/1/20

    work page doi:10.1088/0004-637x/744/1/20 2012

  13. [13]

    Smith, J. D. T., Armus, L., Dale, D. A., et al. 2007, PASP, 119, 1133, doi: 10.1086/522634

    work page doi:10.1086/522634 2007

  14. [14]

    T., Helou, G., & Werner, M

    Soifer, B. T., Helou, G., & Werner, M. 2008, ARA&A, 46, 201, doi: 10.1146/annurev.astro.46.060407.145144

    work page doi:10.1146/annurev.astro.46.060407.145144 2008

  15. [15]

    Starkey, C. A. 2016, PhD thesis, University of Toledo, Ohio van Dishoeck, E. F. 2004, ARA&A, 42, 119, doi: 10.1146/annurev.astro.42.053102.134010

    work page doi:10.1146/annurev.astro.42.053102.134010 2016

  16. [16]

    L., et al., 2010, @doi [ ] 10.1088/0004-6256/140/6/1868 , http://adsabs.harvard.edu/abs/2010AJ....140.1868W 140, 1868

    Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868

    work page internal anchor Pith review doi:10.1088/0004-6256/140/6/1868 2010