Faint white dwarf flux standards: data and models

Abhijit Saha; Annalisa Calamida; Armin Rest; Gautham Narayan; Ivan Hubeny; Jay B. Holberg; Karl D. Gordon; Ralph C. Bohlin; Susana Deustua; Thomas Matheson

REVIEW 2 major objections 5 minor 15 references

Faint white dwarf flux standards: data and models

T0 review · 2 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper establishes 19 faint white dwarfs as new CALSPEC flux standards, with model spectra from 900 Å to 30 μm fitted to HST STIS and WFC3 data.

desk verdict Solid, transparent release of 19 faint WD flux standards with new STIS UV data, but the 1% accuracy claim only holds where data exist; the mid-IR is model extrapolation at 2–3%. read the letter →

arxiv 2411.09049 v1 pith:HH5CE7O3 submitted 2024-11-13 astro-ph.IM

Ralph C. Bohlin , Susana Deustua , Gautham Narayan , Abhijit Saha , Annalisa Calamida , Karl D. Gordon , Jay B. Holberg , Ivan Hubeny

show 2 more authors

Thomas Matheson Armin Rest

This is my paper · ORCID

classification astro-ph.IM

keywords whitedwarffluxstandardsCALSPECSTISWFC3modelatmospheresJWSTcalibrationspectralenergydistribution

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Fainter standard stars are needed because the current bright standards saturate many JWST detector modes. The paper adds 19 white dwarfs with V magnitudes between 16.5 and 18.7, scattered across the sky, and fits pure-hydrogen model atmospheres to new STIS ultraviolet spectra and six-band WFC3 photometry. The resulting spectral energy distributions run continuously from 900 Å to the JWST limit of 30 μm and are placed on the current 2020 CALSPEC flux scale. The authors judge these SEDs adequate for an initial delivery to CALSPEC, with a goal of 1% accuracy in the observed range.

What carries the argument

The load-bearing device is the Hubeny tlusty grid of pure-hydrogen non-local-thermodynamic-equilibrium white-dwarf model atmospheres: 132 models with effective temperatures 20,000–95,000 K and surface gravities log g from 7.0 to 9.5. The fitting procedure minimizes a reduced chi-square formed from binned STIS spectra and synthetic WFC3 photometry, with log g held at the values from the prior Axelrod analysis and interstellar extinction treated with average extinction curves, usually R(V)=3.1. The best-fit model provides the entire 900 Å–30 μm spectrum, including wavelengths that no instrument observed, so the machinery must carry the argument past the 1.6 μm edge of the data.

What would settle it

Observe all 19 stars with JWST MIRI at roughly 5, 10, and 20 μm and compare those fluxes to the delivered model SEDs; if the mean offset exceeds about 3%, the model extrapolation that carries the paper's central claim is not adequate.

Watch

Extended reading notes

Core claim

The paper's claim is that 19 faint white dwarfs can serve as practical flux standards for JWST and other large-aperture telescopes. For each star, the authors fit Hubeny tlusty pure-hydrogen NLTE model atmospheres to new STIS spectrophotometry from 1150 to 3000 Å and to six-band WFC3 photometry from 0.28 to 1.6 μm, varying only effective temperature and selective extinction while holding surface gravity fixed. The best-fit model then supplies the complete spectral energy distribution from 900 Å to 30 μm. The paper argues that these predicted SEDs are already adequate for an initial delivery to CALSPEC on the 2020 HST flux scale, with agreement near 1% where data exist and an estimated 2–3% uncertainty in the extrapolated mid-infrared.

Load-bearing premise

The load-bearing premise is that the model atmospheres correctly predict the infrared brightness of each star from the temperature and reddening fit to data that stop at 1.6 μm, even though the 2–30 μm range is never directly measured.

Editorial extensions

If this is right

If the SEDs are correct, JWST can observe these fainter standards in its normal science modes, avoiding the small subarray modes used for bright standards.
The all-sky placement of the 19 white dwarfs gives Euclid, Roman, and Rubin an accessible extension of the CALSPEC flux scale.
Agreement between JWST calibrations based on these white dwarfs and those based on A- and G-type standards would validate the model extrapolation into the mid-infrared.
The SEDs can be posted to CALSPEC immediately as an initial faint-standard delivery, with revisions to follow as data and models improve.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Editorial extension: the seven stars whose 1350–1600 Å bins were successfully fit are the safer ultraviolet anchors; the other twelve carry unexplained FUV residuals that may reflect nonstandard interstellar reddening rather than model error.
Editorial extension: if the G140L mismatch is interstellar in origin, high-resolution ultraviolet spectra of the twelve discrepant stars could map the anomalous extinction and potentially restore them as full-wavelength standards.
Editorial extension: because the photometry of the 19 stars is internally consistent at the few-millimagnitude level, the same SEDs could serve as a faint transfer network tying ground-based optical surveys to the CALSPEC scale.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Tim, here’s the short version: this is a solid, transparent data release that gives you 19 faint white dwarf SEDs with new STIS UV spectra, on the current CALSPEC scale, with all-sky coverage. Worth having. But the 1% accuracy in the abstract applies only where the data constrain the fit—roughly 0.25–1.6 μm for most stars. The 2–30 μm part is model extrapolation, and the paper’s own numbers cap it at 2–3%. So treat the mid-IR as provisional.

What’s genuinely new: the 19 HST STIS G140L/G230L spectra between 1150 and 3000 Å, and the fitted model SEDs tied to the 2020 Bohlin et al. CALSPEC scale. That’s a real step beyond Axelrod et al. (2023), which used only WFC3 photometry on the 2014 scale. The fitting is the standard CALSPEC chi-square routine; nothing exotic, but it’s reproducible and the WFC3 residuals are small (Table 5 offsets are all within 0.006 mag). The paper also gives clear extinction handling and flags where it deviates from R(V)=3.1. Good.

Where it gets soft: the mid-IR. The longest observed point is F160W at 1.6 μm; everything beyond is pure tlusty model prediction. Section 3.4 is candid: the two independent grids (tlusty vs tmap) disagree by up to 2% at 30 μm, and preliminary JWST checks are good to only ~3% in the 2–30 μm band. So the abstract’s “goal of 1% accuracy” is a goal, not a delivered error bar for the full range. Also, for 12 of the 19 stars the two shortest G140L bins are dropped because no average extinction curve fits; residuals are many sigma. That doesn’t invalidate the fits, but it means those Teff/E(B-V) values carry extra systematic uncertainty that propagates straight into the Rayleigh-Jeans tail. The paper says so, but a reader skimming the abstract could easily over-trust the 1% figure.

Bottom line: for calibration folks this is a useful, honest addition. I’d use the SEDs for relative calibration shortward of 1.6 μm and treat the IR with caution until independent validation lands. It deserves a serious referee and, with minor wording changes to the accuracy claim, should be publishable. Bring it to the group if you’re doing JWST or Rubin calibration work; otherwise it’s a maybe.

Referee Report

2 major / 5 minor

Summary. The paper presents new HST/STIS UV spectrophotometry (1150–3000 Å) for 19 faint white dwarfs, combines it with six-band WFC3 photometry, and fits pure-hydrogen tlusty NLTE model SEDs to derive Teff and E(B-V). The resulting model SEDs are intended as an initial delivery to CALSPEC, providing continuous coverage from 900 Å to 30 μm with the stated goal of 1% accuracy. The fitting procedure is the standard CALSPEC χ2 approach, with WFC3 residuals at the few-milli-mag level and an explicit discussion of model-grid and JWST consistency uncertainties. The paper is transparent about the fact that the mid-IR portion is an extrapolation, but it frames the 1% goal as applying to the full wavelength range even though the direct constraints end at 1.6 μm and the UV constraints are omitted for a majority of the stars.

Significance. If the delivered SEDs are reliable, this work fills a real need: faint, all-sky flux standards for JWST and other large-aperture facilities, usable in standard science modes without subarray complications. The paper's strengths are the clear description of the data reduction, the reproducible and standard fitting methodology, the small WFC3 residuals, and the explicit comparison with an independent model grid and with preliminary JWST calibrations. The main limitation is that the 2–30 μm SEDs are uncalibrated model extrapolations, with the paper's own uncertainty estimates giving a 2% grid-agreement floor and about 3% JWST consistency. The manuscript is honest about many of these caveats, but the abstract and discussion overstate the degree to which 1% accuracy is demonstrated over the full 900 Å–30 μm range. The contribution is nevertheless potentially valuable for the calibration community, provided the accuracy claims are revised and the impact of the omitted FUV bins on the extrapolated fluxes is quantified.

major comments (2)

[Abstract; §3.4; Table 4] The abstract states a goal of 1% accuracy for SEDs covering 900 Å to 30 μm, but the paper's own uncertainty analysis in §3.4 gives a 2% floor from the tlusty/tmap grid comparison in the mid-IR and only ~3% consistency with preliminary JWST calibrations over 2–30 μm. Direct observational constraints end at F160W (1.6 μm), and for 12 of the 19 stars the shortest STIS bins are omitted (Table 3, Table 4). Thus the 1% goal is established, at best, only over the 0.27–1.6 μm region. Please revise the abstract and Section 4 to state separately the demonstrated accuracy in the directly constrained region (~1%) and the estimated accuracy in the extrapolated 2–30 μm region (~2–3%), rather than implying 1% over the full range.
[§3.2; Table 4] The FUV residuals for the 12 stars with omitted G140L bins reach -7.8% (WDFS1302+10) and +5.4% (WDFS2351+37), and two stars cannot be fit with any average extinction curve. Because Teff and E(B-V) are degenerate with UV extinction, the fitted parameters and the extrapolated Rayleigh-Jeans tail may be biased by the omitted FUV data. The paper lists possible causes but does not quantify the impact. Please add a sensitivity test that perturbs Teff and E(B-V) within ranges consistent with the observed residuals and reports the resulting changes in predicted 2–30 μm fluxes; this would provide a concrete systematic uncertainty budget for the delivered SEDs.

minor comments (5)

[Section 4] The text says the SEDs provide coverage 'at a resolution R=5,000,' while the STIS observations are at R~500. Clarify that R=5,000 refers to the model SED grid, not the observations.
[Table 4] The final column is labeled 'Resid (%)' but the text in §3.2 describes the residuals as 'many sigma'; please state both the percentage and the assumed sigma (e.g., the ~1% repeatability) so the reader can judge the significance.
[Table 3] The notes (a), (b), and (c) are informative, but they are easy to misread because the listed star names are not repeated in the table body; consider adding explicit star-name columns or a separate table of excluded bins per star.
[References] Gordon et al. (2024) is cited as 'submitted'; please update the reference if it has been accepted or published.
[ORCID list] The ORCID entry for Ivan Hubeny reads 'Hubenby'; this typo should be corrected.

Circularity Check

0 steps flagged · score 1.0 of 10

No self-definitional circularity: the SEDs are fits to external STIS/WFC3 data with model extrapolation to 30 um; the 2-30 um accuracy limitation is empirical, not a forced equivalence.

full rationale

The paper's derivation chain is a standard model-fitting calibration: observed STIS spectrophotometry and WFC3 photometry are external constraints, and the Hubeny tlusty model grid is fit by varying Teff and E(B-V) to minimize chi-squared (Section 3.1). The resulting SEDs are the fitted models, not a re-statement of the fitted parameters under another name. The long-wavelength SED from 1.6 to 30 um is a model extrapolation, explicitly acknowledged in Section 3.4 as having a minimum 2% uncertainty from grid-to-grid disagreement and 'actual errors in the IR can be larger' if systematics affect both grids. This is an empirical accuracy limitation, not a circular reduction. The conversion of Axelrod et al. (2023) photometry from the 2014 to the 2020 CALSPEC flux scale uses the three primary WD model SEDs, but those primaries are independent of the 19 target stars, so the scale transfer does not define the target SEDs by construction. The paper also compares against an independent model grid (Rauch tmap) and against external JWST flux calibrations, providing non-circular checks. The only mild concern is that the calibration infrastructure (CALSPEC procedures, model grid, flux scale) originates from the same group's prior work, but this is normal continuity in calibration science and is not load-bearing in the sense that the target SEDs are forced by the inputs. No specific equation or fitted parameter is renamed as a prediction, and no uniqueness claim is imported from self-citations. Therefore no significant circularity is found; the score reflects minor self-citation and inherited calibration scale rather than any self-definitional step.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities. Its parameters are the standard stellar and extinction parameters fitted to data, plus empirically chosen fitting bins. The model grid, flux scale, and extinction curves are taken from prior literature, mostly from the same collaboration.

free parameters (4)

Effective temperature Teff per star = 19810 to 57830 K (Table 4)
Adjusted in the chi-square fit to match STIS and WFC3 data; the main free parameter of each SED.
Selective extinction E(B-V) per star = A(V) values 0.003 to 0.323 (Table 4)
Adjusted in the chi-square fit along with Teff; critical for UV SED shape.
Total-to-selective extinction R(V) per star = Default 3.1; exceptions 3.4, 3.5, 4.2 (Table 4)
R(V) is fixed at 3.1 for most stars but changed by hand for four stars to improve residuals.
STIS fitting wavelength bins = Table 3 ranges, with per-star omissions
Wavelength bins are chosen manually to avoid interstellar lines and quasi-molecular features, and some bins are omitted for specific stars.

assumptions (5)

domain assumption Hubeny tlusty pure-hydrogen NLTE model atmospheres accurately predict DA white dwarf SEDs from 900 Å to 30 μm.
The predicted SEDs are the model fits; this is the central modeling assumption, invoked in Section 3 and relied on for all IR extrapolation.
domain assumption Gordon et al. (2023) average extinction curves with the adopted R(V) describe interstellar extinction along each sightline.
Extinction corrections are applied using these curves in Section 3.2; the paper notes two stars where no average curve works.
domain assumption Axelrod et al. (2023) log g values for these white dwarfs are correct and can be held fixed.
Section 3.1 fixes log g from Axelrod rather than fitting it, so errors in those values propagate into the SEDs.
domain assumption The CALSPEC 2020 absolute flux scale defined by Bohlin et al. (2020) is correct.
All photometry and spectra are tied to this scale in Section 2.2; if the scale has systematics, the new SEDs inherit them.
domain assumption WFC3 and STIS calibration repeatabilities stated in Section 3.1 correctly quantify the measurement uncertainties.
The chi-square weights and quoted uncertainties depend on these repeatability figures.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Faint white dwarf flux standards: data and models." pith.science (2026). https://pith.science/paper/HH5CE7O3

@misc{pith2026241109049,
  author       = {Pith},
  title        = {Pith review of: Faint white dwarf flux standards: data and models},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HH5CE7O3}},
  note         = {Machine review of arXiv:2411.09049}
}

abstract

Fainter standard stars are essential for the calibration of larger telescopes. This work adds to the CALSPEC (calibration spectra) database 19 faint white dwarfs (WDs) with all-sky coverage and V magnitudes between 16.5 and 18.7. Included for these stars is new UV (ultraviolet) HST (Hubble Space Telescope) STIS (Space Telescope Imaging Spectrometer) spectrophotometry between 1150 and 3000~\AA\ with a resolution of $\sim$500. Pure hydrogen WD models are fit to these UV spectra and to six-band HST/WFC3 (Wide Field Camera 3) photometry at 0.28 to 1.6~\micron\ to construct predicted model SEDs (spectral energy distributions) covering wavelengths from 900~\AA\ to the JWST (James Webb Space Telescope) limit of 30~\micron\ using well-established CALSPEC procedures for producing flux standards with the goal of 1\% accuracy.

Figures

Figures reproduced from arXiv: 2411.09049 by the authors.

**Figure 2.** Model fit and residuals for WDFS2351+37 as in [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 7 canonical work pages

[1]

2023, ApJ, 951, 78

Axelrod, T., Saha, A., Matheson, T., et al. 2023, ApJ, 951, 78

2023
[2]

C., & Deustua, S

Bohlin, R. C., & Deustua, S. E. 2019, AJ, 157, 229

work page 2019
[3]

C., Deustua, S

Bohlin, R. C., Deustua, S. E., & de Rosa, G. 2019, The Astronomical Journal, 158, 211

work page 2019
[4]

C., Gordon, K

Bohlin, R. C., Gordon, K. D., & Tremblay, P.-E. 2014, PASP, 126, 711

2014
[5]

C., Hubeny, I., & Rauch, T

Bohlin, R. C., Hubeny, I., & Rauch, T. 2020, AJ, 160, 21 4 http://www.stsci.edu/hst/instrumentation/reference-data-for-calibration-and-tools/astronomical-catalogs/calspec

work page 2020
[6]

C., M´ esz´ aros, S., Fleming, S

Bohlin, R. C., M´ esz´ aros, S., Fleming, S. W., et al. 2017, AJ, 153, 234

2017
[7]

2019, ApJ, 872, 199

Calamida, A., Matheson, T., Saha, A., et al. 2019, ApJ, 872, 199

2019
[8]

W., et al

Calamida, A., Matheson, T., Olszewski, E. W., et al. 2022, ApJ, 940, 19

work page 2022

Show all 15 references

[9]

D., Clayton, G

Gordon, K. D., Clayton, G. C., Decleir, M., et al. 2023, ApJ, 950, 86

2023
[10]

D., Bohlin, R., Sloan, G

Gordon, K. D., Bohlin, R., Sloan, G. C., et al. 2022, AJ, 163, 267

2022
[11]

1995, ApJ, 439, 875 13

Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 13

1995
[12]

S., & Cardelli, J

Mathis, J. S., & Cardelli, J. A. 1992, ApJ, 398, 610

1992
[13]

2019, ApJS, 241, 20

Narayan, G., Matheson, T., Saha, A., et al. 2019, ApJS, 241, 20

2019
[14]

A., Clayton, G

Valencic, L. A., Clayton, G. C., & Gordon, K. D. 2004, ApJ, 616, 912

2004
[15]

N., Bohlin, R

Witt, A. N., Bohlin, R. C., & Stecher, T. P. 1984, ApJ, 279, 698

1984

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [1]

2023, ApJ, 951, 78

Axelrod, T., Saha, A., Matheson, T., et al. 2023, ApJ, 951, 78

2023

[2] [2]

C., & Deustua, S

Bohlin, R. C., & Deustua, S. E. 2019, AJ, 157, 229

work page 2019

[3] [3]

C., Deustua, S

Bohlin, R. C., Deustua, S. E., & de Rosa, G. 2019, The Astronomical Journal, 158, 211

work page 2019

[4] [4]

C., Gordon, K

Bohlin, R. C., Gordon, K. D., & Tremblay, P.-E. 2014, PASP, 126, 711

2014

[5] [5]

C., Hubeny, I., & Rauch, T

Bohlin, R. C., Hubeny, I., & Rauch, T. 2020, AJ, 160, 21 4 http://www.stsci.edu/hst/instrumentation/reference-data-for-calibration-and-tools/astronomical-catalogs/calspec

work page 2020

[6] [6]

C., M´ esz´ aros, S., Fleming, S

Bohlin, R. C., M´ esz´ aros, S., Fleming, S. W., et al. 2017, AJ, 153, 234

2017

[7] [7]

2019, ApJ, 872, 199

Calamida, A., Matheson, T., Saha, A., et al. 2019, ApJ, 872, 199

2019

[8] [8]

W., et al

Calamida, A., Matheson, T., Olszewski, E. W., et al. 2022, ApJ, 940, 19

work page 2022