pith. machine review for the scientific record. sign in

arxiv: 2605.12721 · v1 · submitted 2026-05-12 · 🌌 astro-ph.GA

Recognition: 1 theorem link

· Lean Theorem

The Panchromatic Hubble Andromeda Southern Treasury (PHAST). II. The Spatially Resolved Recent Star Formation History in M31

Tobin M. Wainer , Benjamin F. Williams , Zhuo Chen , Margaret Lazzarini , Julianne J. Dalcanton , Eric F. Bell , Kameron Goold , Andrew Dolphin
show 14 more authors
Meredith J. Durbin Stefany L. Fabian Dub\'on Karoline M. Gilbert Puragra Guhathakurta Francois Hammer L. Clifton Johnson Eric W. Koch Kristen. B. W. McQuinn Ekta Patel Vaishnav V. Rao Julia Roman-Duval Adam Smercina Debby Tran Daniel R. Weisz
Authors on Pith no claims yet

Pith reviewed 2026-05-14 19:52 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords M31Andromeda galaxystar formation historycolor-magnitude diagramsHubble Space Telescopespatially resolvedgalactic ringsrecent star formation
0
0 comments X

The pith

M31's star formation rate has declined globally with a pronounced drop in the last 40 million years.

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

The paper uses Hubble Space Telescope imaging of the southern disk of the Andromeda galaxy to measure how star formation has changed over the past 500 million years. Researchers divided the area into more than 6500 small patches and fitted color-magnitude diagrams in each to reconstruct the history of star birth. The resulting maps show a clear overall slowdown, with the steepest decline occurring in the most recent 40 million years and appearing strongest near the companion galaxy M32. When these southern measurements are combined with earlier northern data, they cover two thirds of M31's star-forming disk and give average rates of 0.445 solar masses per year over the last 100 million years and 0.285 over the last 20 million years. The pattern is interpreted as the final phase of a longer-term reduction in activity that has been underway for several billion years and is driven largely by fading star formation inside the disk's ring structures.

Core claim

The authors establish that the recent star formation rate in M31 has undergone a global decline, with a pronounced drop in the last ~40 Myr that is most evident near M32. Fitting color-magnitude diagrams across over 6500 individual 0.01 kpc² regions yields mean SFRs of 0.445 ± 0.006 M⊙ yr⁻¹ over the last 100 Myr and 0.285 ± 0.014 M⊙ yr⁻¹ over the last 20 Myr in the joint PHAST+PHAT footprint, implying total disk SFRs of ~0.67 and ~0.43 M⊙ yr⁻¹. The decline is interpreted as the late stage of a multi-Gyr wind-down from a more active past state and is driven mainly by decreasing activity within the ring features. CMD-based SFHs produce a synthetic GALEX FUV image matching observations, while F

What carries the argument

Spatially resolved recent star formation histories obtained by fitting color-magnitude diagrams in over 6500 individual 0.01 kpc² regions using stellar population synthesis models.

If this is right

  • The global decline is driven mainly by decreasing activity within the ring features of the disk.
  • The observed trend represents the late stage of a multi-Gyr wind-down from a previously more active state.
  • CMD-derived SFHs accurately describe recent star formation, as shown by the matching synthetic FUV image.
  • FUV+24 μm prescriptions underestimate the CMD-based 100-Myr average SFR by a factor of ~2.1.
  • Star-formation tracers that implicitly average over ~100 Myr are unreliable when the recent SFR is evolving.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the decline continues, M31's total star formation rate may fall further in the coming tens of millions of years.
  • The localized drop near M32 suggests that interactions with the companion could be accelerating the slowdown in that sector of the disk.
  • The mismatch between CMD and FUV tracers indicates that time-resolved methods are needed for galaxies whose star formation rate is changing on short timescales.
  • Applying the same small-region CMD-fitting technique to other nearby spirals could identify whether similar multi-Gyr wind-down phases are common.

Load-bearing premise

The CMD-fitting procedure using standard stellar population synthesis models recovers the true recent star formation history without large systematic biases from crowding, variable extinction, or inaccuracies in the adopted isochrones and initial mass function for M31's stellar populations.

What would settle it

An independent measurement of the star formation rate in the last 20-40 million years, for example from counting very young star clusters or from H-alpha emission in the same regions, that yields values close to the 100-Myr average would falsify the claimed sharp recent drop.

Figures

Figures reproduced from arXiv: 2605.12721 by Adam Smercina, Andrew Dolphin, Benjamin F. Williams, Daniel R. Weisz, Debby Tran, Ekta Patel, Eric F. Bell, Eric W. Koch, Francois Hammer, Julianne J. Dalcanton, Julia Roman-Duval, Kameron Goold, Karoline M. Gilbert, Kristen. B. W. McQuinn, L. Clifton Johnson, Margaret Lazzarini, Meredith J. Durbin, Puragra Guhathakurta, Stefany L. Fabian Dub\'on, Tobin M. Wainer, Vaishnav V. Rao, Zhuo Chen.

Figure 1
Figure 1. Figure 1: Left: Stellar density in the PHAST survey, as defined by stars per square arcsecond with 21.5 < F814W < 23. Right: 50% completeness limits as determined by ASTs as a function of local stellar density, estimated from bootstrapped sampling of across the PHAST footprint. Completeness declines at higher stellar densities due to crowding. The gray shaded region in the right panel, and illustrated on the color b… view at source ↗
Figure 2
Figure 2. Figure 2: Spatial distribution of recovered extinction parameters across the PHAST footprint, as determined by the CMD-fitting procedure in each analysis region. Left panel shows the foreground extinction AV , while the middle panel shows the differential extinction dAV , and the right panel shows the total extinction AV +dAV . Regions are color-coded by extinction, with darker values indicating higher total extinct… view at source ↗
Figure 3
Figure 3. Figure 3: Padova stellar isochrones spanning log(Age/yr) = 6.6–10.2 plotted in CMD space (F475W−F814W vs. F814W), color-coded by age. Black stars mark the main-sequence turnoff point for each isochrone, defined as the transition from core hydro￾gen burning to the subgiant phase. The gray shaded region shows the excluded region described in Section 3.1. The secondary axis to the right shows the absolute F814 magnitud… view at source ↗
Figure 4
Figure 4. Figure 4: Overview of three representative PHAST regions spanning a range of stellar densities and environments. Upper Left: A color mosaic of the southern disk of M31, and the three selected regions (6068, 835, and 2284) highlighted (not perfectly to scale). For each region, we show: (upper left) the observed CMD used in SFH fitting, with excluded areas shaded in gray; (upper right) a color cutout of the region fro… view at source ↗
Figure 5
Figure 5. Figure 5: PHAST SFH maps over the indicated timescales in terms of lookback time. The top rows are the most recent epochs, and the subsequent rows look at older lookback times. Each pixels represents one 100 pc2 region, colored by the star formation rate surface density ΣSFR in units of M⊙ yr−1 kpc−2 . In each panel, we denote the location of M32 with a blue star. In the bottom right panel, we also include the avera… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison between the observed and synthetic FUV emission in the PHAST footprint. Left: observed GALEX FUV image, reprojected to the PHAST spatial grid. Middle: synthetic FUV map constructed by forward-modeling the CMD-derived SFHs with FSPS following the methodology of Lewis et al. (2017), including a two-component extinction model based on the Cardelli et al. (1989) law. Right: residual map showing the … view at source ↗
Figure 7
Figure 7. Figure 7: Full M31, PHAT and PHAST combined SFH maps over the indicated timescales. Five time bins are shown moving to older lookback times downward in the plot. The left side of the images is the PHAST region, shown in [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The star formation history of M31 from the PHAT+PHAST as a function of lookback time, where the most recent time is to the right. The black line is the mean SFR per given time bin, while the red lines show the uncertainty for each bin. In each panel, the gray dashed line shows the average SFR over the plotted time baseline. We show both the SFR, and the ΣSFR on the right axis. X-axis ranges were chosen to … view at source ↗
Figure 9
Figure 9. Figure 9: the spatial distribution of recent star formation in M31 as a fractional offset from the disk-wide mean SFR surface density in each age interval. In each panel we plot (ΣSFR(t)i − ⟨ΣSFR(t)⟩)/⟨ΣSFR(t)⟩, where ΣSFR(t)i is the SFR surface density in an individual spatial cell and ⟨ΣSFR(t)⟩ is the mean over all cells with ⟨ΣSFR⟩0−500 Myr > 6×10−4 M⊙ yr−1 kpc−2 . This mask restricts the comparison to regions co… view at source ↗
Figure 10
Figure 10. Figure 10: Shown is the radial dependence of the recent star formation history in M31. Bottom: map of the mean SFR surface density, ΣSFR, averaged over the last 100 Myr across the combined PHAT+PHAST footprint. The blue star marks the center of M32, and the red cross marks the center of M31. Two matched comparison regions are overlaid on the northern and southern disk, each subdivided into four bins of equal size. T… view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of the SFHs for the North and South regions of the M31 disk. Left panel shows the SFR, while the right shows the SFR surface density ΣSFR as a function of lookback time, where the most recent time is to the right. Blue shows the SFH from the PHAT survey in the North from L15, while the purple are the results presented here from the PHAST survey in the South. Both surveys are declining over the … view at source ↗
Figure 12
Figure 12. Figure 12: Subregions used to investigate local variations in the recent SFH. Left: SFR surface density map for the 0–25 Myr time bin from [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Comparison between SFR surface densities inferred from FUV+24 µm emission and from the CMD-derived SFHs. Top left: SFR surface density map, ΣSFRFUV+24, computed using the Leroy et al. (2008) calibration applied to the GALEX FUV and Spitzer/MIPS 24 µm data (Martin et al. 2005; Gordon et al. 2006). Top right: mean SFR surface density over the last 100 Myr from the CMD-based SFHs, ΣSFRCMD,100, evaluated on t… view at source ↗
Figure 14
Figure 14. Figure 14: Comparison between the observed FUV emission and synthetic FUV maps constructed under different assumptions for the recent SFH. Top left: observed GALEX FUV flux, background-subtracted and measured on the PHAST spatial grid. Top middle: synthetic FUV map generated by setting the SFR in the last 100 Myr in each region with a constant SFR equal to SFRFUV+24µm. Top right: synthetic FUV map generated by by se… view at source ↗
Figure 15
Figure 15. Figure 15: The SFH for PHAST regions with a density less than 1.5 stars per arcsecond2 , where our completeness limits allow us to probed older ages. Here we extend back to 1.2 Gyr, where we can reliably detect the MS turnoff above our completeness limit. We see a long decline in the SFR over this period of time. First, we consider whether a simple recalibration of the FUV+24 µm tracer could reconcile it with the CM… view at source ↗
read the original abstract

We use Hubble Space Telescope optical imaging from the Panchromatic Hubble Andromeda Southern Treasury (PHAST) to measure the spatially resolved recent star formation history (SFH) across the southern disk of M31. We fit color-magnitude diagrams (CMDs) of over 6500 individual 0.01 kpc$^2$ regions to measure SFHs over the last $\sim$500 Myr. The resulting maps show coherent structure that traces the ringed morphology of the disk. We find a clear global decline in the recent SFR, with a pronounced drop in the last $\sim$40 Myr that is most evident in the region closest to M32. Combining PHAST and PHAT measurements, we now cover two thirds of M31's star-forming disk with homogeneous SFHs, yielding the highest-resolution spatially resolved SFHs of M31. Inside the joint footprint, we measure mean SFRs of $0.445 \pm0.006$ M$_\odot$ yr$^{-1}$ over the last 100 Myr and $0.285 \pm 0.014$ M$_\odot$ yr$^{-1}$ over the last 20 Myr, implying total disk SFRs of $\sim$0.67 and $\sim$0.43 M$_\odot$ yr$^{-1}$, respectively. The observed decline is interpreted as the late stage of a multi-Gyr wind-down from a previously more active state. Because recent star formation in M31 is concentrated primarily in the rings, the global decline is driven mainly by decreasing activity within those features. We also compare the CMD-based SFR surface densities to those inferred from FUV+24 $\mu$m prescriptions and find that the FUV-based calibration underestimates the CMD-based 100 Myr average by a factor of $\sim$2.1. However, the PHAST SFHs produce a synthetic GALEX FUV image that agrees well with observations, indicating that the CMD-derived SFHs provide an accurate description of recent star formation. The mismatch with the FUV+24 $\mu$m estimates underscores that tracers implicitly averaged over $\sim$100 Myr are not reliable when the recent SFR is evolving.

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

0 major / 2 minor

Summary. The paper measures the spatially resolved recent star formation history (SFH) in the southern disk of M31 by fitting color-magnitude diagrams (CMDs) from PHAST HST optical imaging in over 6500 individual 0.01 kpc² regions over the last ~500 Myr. It reports a global decline in SFR with a pronounced drop in the last ~40 Myr (most evident near M32), mean SFRs of 0.445 ± 0.006 M⊙ yr⁻¹ (last 100 Myr) and 0.285 ± 0.014 M⊙ yr⁻¹ (last 20 Myr) within the joint PHAST+PHAT footprint covering two-thirds of the star-forming disk, and validates the SFHs via a synthetic GALEX FUV image that matches observations while noting that FUV+24 μm prescriptions underestimate the CMD-based 100 Myr average by a factor of ~2.1.

Significance. If the results hold, this delivers the highest-resolution spatially resolved SFHs yet for two-thirds of M31's star-forming disk, backed by >6500 independent CMD fits with reported uncertainties and direct external validation through the synthetic FUV image agreement with GALEX data. The homogeneous PHAST+PHAT coverage and the explicit quantification of the FUV tracer mismatch in an evolving SFR system represent clear strengths for studies of nearby galaxy evolution.

minor comments (2)
  1. [Abstract] Abstract: the total-disk SFR extrapolations (~0.67 and ~0.43 M⊙ yr⁻¹) are stated without specifying the exact area scaling or footprint fraction used, which would aid reproducibility.
  2. [Methods] The time-binning scheme for the SFH recovery (last ~500 Myr) is referenced but the precise bin widths and their justification for resolving the ~40 Myr drop are not detailed enough to evaluate temporal resolution.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, accurate summary of our methods and results, and recommendation to accept. The referee's comments highlight the strengths of the PHAST SFH maps and the FUV validation, which align with our own interpretation of the data.

Circularity Check

0 steps flagged

No significant circularity in the SFH derivation chain

full rationale

The paper's central results are obtained by directly fitting observed CMDs in >6500 spatial regions to standard stellar population synthesis models, then computing time-binned averages of the resulting SFHs to report mean SFRs (0.445 ± 0.006 M⊙ yr⁻¹ over 100 Myr and 0.285 ± 0.014 M⊙ yr⁻¹ over 20 Myr). These quantities are outputs of the fit rather than inputs redefined by the same equations. Validation is supplied externally via synthetic GALEX FUV images generated from the derived SFHs that match observed photometry, plus consistency checks against independent PHAT data over the joint footprint. No self-citation chain, ansatz, or uniqueness theorem is invoked to force the reported decline or mean rates; the derivation remains self-contained against the observed CMDs and external tracers.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central measurements rest on standard stellar population synthesis models and CMD-fitting techniques drawn from prior literature rather than derived in this work.

free parameters (2)
  • SFH time binning
    The choice of age bins and time resolution for the SFH recovery is set by the authors to balance signal and resolution.
  • Region-specific extinction and distance
    Extinction and distance parameters are fitted or adopted per region using standard M31 values.
axioms (1)
  • domain assumption Stellar evolution isochrones and population synthesis models accurately predict the locations of stars of given age and metallicity in the CMD.
    This assumption underpins the CMD fitting used to extract SFHs from the observed photometry.

pith-pipeline@v0.9.0 · 5831 in / 1487 out tokens · 68938 ms · 2026-05-14T19:52:21.088527+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages · 2 internal anchors

  1. [1]

    Star Formation Rate Indicators

    Alexander, D. M., Hickox, R. C., Aird, J., et al. 2025, NewAR, 101, 101733 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, Astronomy and Astrophysics, 558, A33, aDS Bibcode: 2013A&A...558A..33A. https://ui.adsabs.harvard.edu/abs/2013A&A...558A..33A 30 WAINER ET AL. Astropy Collaboration, Price-Whelan, A. M., Sip˝ocz, B. M., et al. ...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1208.2997 2025

  2. [2]

    2002, Monthly Notices of the Royal Astronomical Society, 332,

    http://adsabs.harvard.edu/abs/2000PASP..112.1383D —. 2002, Monthly Notices of the Royal Astronomical Society, 332,

    work page 2002

  3. [3]

    2013, The Astrophysical Journal, 775, 76, aDS Bibcode: 2013ApJ...775...76D

    http://adsabs.harvard.edu/abs/2002MNRAS.332...91D —. 2013, The Astrophysical Journal, 775, 76, aDS Bibcode: 2013ApJ...775...76D. https://ui.adsabs.harvard.edu/abs/2013ApJ...775...76D Dolphin, A. E., Saha, A., Skillman, E. D., et al. 2003, The Astronomical Journal, 126, 187, aDS Bibcode: 2003AJ....126..187D. https://ui.adsabs.harvard.edu/abs/2003AJ....126....

    work page doi:10.48550/arxiv.2403.19798 2013

  4. [4]

    M., Dalcanton, J

    http://adsabs.harvard.edu/abs/2010ApJ...724.1030G Gogarten, S. M., Dalcanton, J. J., Williams, B. F., et al. 2010, The Astrophysical Journal, 712, 858, aDS Bibcode: 2010ApJ...712..858G. https://ui.adsabs.harvard.edu/abs/2010ApJ...712..858G Gommers, R., Virtanen, P., Haberland, M., et al. 2024, scipy/scipy: SciPy 1.14.0, vv1.14.0, Zenodo, doi:10.5281/zenod...

    work page doi:10.5281/zenodo.12522488 2010

  5. [5]

    2009, The Astrophysical Journal, 703, 614, aDS Bibcode: 2009ApJ...703..614K

    http://adsabs.harvard.edu/abs/2017ApJ...839...78J Kang, Y ., Bianchi, L., & Rey, S.-C. 2009, The Astrophysical Journal, 703, 614, aDS Bibcode: 2009ApJ...703..614K. https://ui.adsabs.harvard.edu/abs/2009ApJ...703..614K Kennicutt, R. C., & Evans, N. J. 2012, Annual Review of Astronomy and Astrophysics, 50,

    work page 2009

  6. [6]

    http://adsabs.harvard.edu/abs/2012ARA%26A..50..531K Kennicutt, Jr., R. C. 1989, The Astrophysical Journal, 344, 685, aDS Bibcode: 1989ApJ...344..685K. https://ui.adsabs.harvard.edu/abs/1989ApJ...344..685K —. 1998, Annual Review of Astronomy and Astrophysics, 36, 189, aDS Bibcode: 1998ARA&A..36..189K. https://ui.adsabs.harvard.edu/abs/1998ARA&A..36..189K K...

    work page doi:10.48550/arxiv.2602.13114 1989

  7. [7]

    The Origin and Evolution of Multiple Star Systems

    https: //ui.adsabs.harvard.edu/abs/2008A%26A...482..883M/abstract Martin, D. C., Fanson, J., Schiminovich, D., et al. 2005, The Astrophysical Journal, 619, L1, aDS Bibcode: 2005ApJ...619L...1M. https://ui.adsabs.harvard.edu/abs/2005ApJ...619L...1M Martin, D. C., Wyder, T. K., Schiminovich, D., et al. 2007, The Astrophysical Journal Supplement Series, 173,...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2303.08774 2005

  8. [8]

    https: //ui.adsabs.harvard.edu/abs/2005NewAR..49..453S/abstract Smith, M. W. L., Eales, S. A., Gomez, H. L., et al. 2012, The Astrophysical Journal, 756, 40, aDS Bibcode: 2012ApJ...756...40S. https://ui.adsabs.harvard.edu/abs/2012ApJ...756...40S Somerville, R. S., & Davé, R. 2015, Annual Review of Astronomy and Astrophysics, 53, 51, aDS Bibcode: 2015ARA&A...

    work page 2012

  9. [9]

    M., Sandstrom, K., Leroy, A., & Smith, J.-D

    http://adsabs.harvard.edu/abs/2015ApJ...806..198W Whitcomb, C. M., Sandstrom, K., Leroy, A., & Smith, J.-D. T. 2023, The Astrophysical Journal, 948, 88, aDS Bibcode: 2023ApJ...948...88W. https://ui.adsabs.harvard.edu/abs/2023ApJ...948...88W Williams, B. F. 2002, Monthly Notices of the Royal Astronomical Society, 331, 293, aDS Bibcode: 2002MNRAS.331..293W....

    work page 2023

  10. [10]

    F., Dalcanton, J

    http://stacks.iop.org/1538-3881/126/i=3/a=1312 Williams, B. F., Dalcanton, J. J., Stilp, A., et al. 2013, The Astrophysical Journal, 765, 120, aDS Bibcode: 2013ApJ...765..120W. https://ui.adsabs.harvard.edu/abs/2013ApJ...765..120W Williams, B. F., Garcia, M. R., Kong, A. K. H., et al. 2004, The Astrophysical Journal, 609, 735, aDS Bibcode: 2004ApJ...609.....

    work page 2013

  11. [11]

    Shown in Figure15, are the SFH results for all PHAST regions with density lower than 1.5 stars per arcsecond 2, extending back to 1.2 Gyr

    From Figure 3, we can see that in these less dense regions, we can therefore reliably probe to somewhat older ages than what is considered in the main body of the manuscript. Shown in Figure15, are the SFH results for all PHAST regions with density lower than 1.5 stars per arcsecond 2, extending back to 1.2 Gyr. These results are in striking agreement wit...

    work page 2017

  12. [12]

    Here we extend back to 1.2 Gyr, where we can reliably detect the MS turnoff above our completeness limit

    Figure 15.The SFH for PHAST regions with a density less than 1.5 stars per arcsecond 2, where our completeness limits allow us to probed older ages. Here we extend back to 1.2 Gyr, where we can reliably detect the MS turnoff above our completeness limit. We see a long decline in the SFR over this period of time. First, we consider whether a simple recalib...

    work page 2008

  13. [13]

    We find best-fit values ofa FUV = 0.36 anda 24 = 3.1×10 −3. This refit substantially reduces the global normalization offset and improves the slope of the region-to-region residual, where the best-fit slope shifts to 0.42, closer to the 1.0 expected for the one-to-one relation, but still dramatically off base. Additionally, the improvement is mostly in th...

    work page 2008