--- title: "Delayed sources and imperfect mixing" description: "Read a fixed 2026 dwarf-galaxy study through its figures: Ia delay, r-process source degeneracies, yield sensitivity and failures of a mixed one-zone model." canonical_url: "https://sandboxphysics.com/stories/galactic-chemical-evolution/delayed-sources-and-imperfect-mixing/" language: "en" page_type: "story" alternate_languages: zh-CN: "https://sandboxphysics.com/zh/stories/galactic-chemical-evolution/delayed-sources-and-imperfect-mixing/" last_modified: "2026-11-05" --- [Galactic chemical evolution](https://sandboxphysics.com/stories/galactic-chemical-evolution/) Galactic chemical evolution · L4 # Delayed sources and imperfect mixing Sandbox Physics · 2026-11-05 Similar magnesium-to-iron tracks can correspond to different formation times. Do barium, strontium, yttrium and europium then identify the source of heavy elements and its waiting time? The fixed 2026 preprint by Henderson and collaborators makes this question concrete. Its figures show both what the data constrain under a specified model and how those conclusions depend on formation history, yields and mixing. ## Before you begin Delay-time distributions, birth weighting and conditional measurement densities. This guide uses the complete 24-page arXiv:2610.15946v1 preprint. It does not independently rerun the authors' multi-element sampler. ## Identify the sample and the model first The paper uses literature Keck/DEIMOS samples for Sculptor and Draco. They are not the Hill FLAMES 99-star table in the previous tutorial. Its multi-element model links star formation to gas evolution and adopts a particular stellar mass distribution and yield prescriptions. Our calculation keeps only fixed-history, effective-magnesium-and-iron evolution; the additional elements and fitted gas-linked history belong to the paper’s model. The r-process is rapid neutron capture: nuclei capture neutrons faster than the relevant decays. Neutron-star mergers can provide such conditions; some prompt stellar explosions may also contribute. The paper represents a prompt channel using magnetorotational-supernova yields and an occurrence fraction. The occurrence fraction is conditional on that proposed yield model; identifying the astrophysical sources remains the research question. Reading: Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1 ## Figure 2: similar chemistry, different clocks Figure 2 varies the minimum Type Ia delay among 0.05, 0.1, 0.2 and 0.5 gigayears, refitting other parameters in each case. It is not simply a time translation of one fixed formation history. The magnesium-to-iron tracks on the right can look similar while the formation and cumulative-formation histories on the left differ. Check the lookback-time axis before deciding which direction means later. Table 2 gives corresponding minimum neutron-star-merger delays of 12, 15, 38 and 64 million years. These four refitted results show how the inferred merger delay depends on the Ia-delay and formation-history assumptions. The adopted Ia power-law slope is minus 1.1 and the merger slope minus 1.5. These distributions include event-rate amplitudes and are not interchangeable with our unit-normalized yield kernel. The listed minimum negative log likelihoods, 1521, 1559, 1656 and 1704, summarize those conditional fits. They are not Bayes factors obtained by integrating all parameter and model priors. Reading: Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1 ## Figures 5 and 6: is a prompt channel required? Compare merger-only production, a minimum merger delay fixed at 50 million years plus a prompt source, and prompt-only production. The actual figure legends use black dotted, black solid and blue dashed lines respectively. Their grayscale caption wording is not fully consistent with those legends, so identify the curves by model name and line pattern. The merger-only model returns conditional minimum delays around 14 million years for Sculptor and 16 million years for Draco. That short interval must include the formation of both neutron stars, their inspiral and merger, and mixing of the ejecta into star-forming gas. Sections 4.3–4.5 explain why the authors see a physical tension in that time budget and investigate an additional prompt channel. The mixed model fixes the merger minimum at 50 million years, described by the authors as a realistic but arbitrary test value. It is not a measured universal lower bound. Under the stated likelihood, the mixed model improves on merger-only production for Sculptor but not for Draco. A preference in one galaxy cannot be promoted into a universal source identification. These channels can produce partly similar abundance trends. The interpretation also depends on yields, occurrence rates, detection limits and selection. For an abundance upper limit, the likelihood uses the cumulative probability of a measurement below that threshold. This gives a distinct contribution for a star without a precise abundance detection. Reading: Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1 ## Event rates are not element-production rates Figure 8 uses the merger-only solution in the first column of Table 3. The blue solid curve is the s-process contribution from asymptotic-giant-branch stars, a late stage of lower- and intermediate-mass stars where neutron captures are relatively slow. The black dashed curve is the merger r-process contribution. The upper panel gives event rates and the lower panel gives fractions of barium production. Yield per event connects those quantities. Its time axis comes from the fitted history, rather than an independent set of measured stellar ages. Likewise, the 100 curves in Figure 7 are parameter realizations, not 100 independent data sets or a coverage test. The yield-rate degeneracy has a striking scale. The discussion notes that europium yields around one ten-millionth of a solar mass per event can require thousands of mergers, whereas yields around one hundred-thousandth can reduce the required count to tens. A chemical total does not separately measure both an event count and its yield. Figure 9 varies particular strontium and yttrium yield scales with other parameters held fixed. It exposes sensitivity to the adopted yields and stellar mass distribution. The suggested rescaling remains conditional on those choices. Reading: Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1 ## A calculable limit of instantaneous mixing A one-zone model assigns a single gas composition to each time. To see what instantaneous mixing omits, consider two equal-mass gas reservoirs. Initially one has an element fraction 0.002 and the other zero. There is no production or loss, only equal opposing exchange at rate $\kappa$ per unit time. Adding their rate equations conserves the sum. Subtracting them makes the concentration contrast decay at twice the one-way exchange rate. With $\kappa$ one per gigayear, half a gigayear leaves a finite contrast about 0.0007358. This is an original conservative two-reservoir example, not the paper's europium transport solution or a random-scatter channel in the one-zone model. $$ \frac{d(Z_1-Z_2)}{dt}=-2\kappa(Z_1-Z_2),\qquad Z_1-Z_2=0.002e^{-2\kappa t} $$ The equal masses make the exchange terms equal and opposite. The product $\kappa$ times elapsed time is dimensionless. Reading: Independent one-zone chemical-evolution Python materials ## Finish with Appendix A, including the failed galaxy The Fornax comparison in Appendix A is an essential counterexample. The one-zone fit concentrates formation into roughly two gigayears, unlike the more extended, episodic history supported by photometric work. Figures 10 and 11 and the metallicity distribution reveal the pressure on that model. Returning a merger delay around 11 million years does not repair a mismatched formation history. Three questions remain useful after closing the paper. Does a heavy-element delay survive a changed Ia clock? Can event rates and event yields be determined separately? And if stars form in different places, is instantaneous mixing sufficient? The fixed magnesium-and-iron model calculates gas flows and element production with effective yields. It contains no europium, r-process production or spatial abundance scatter. Reading: Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1 ## Exercises and explained answers 1. After one gigayear of the stated two-reservoir exchange, what happens to the total element amount and the contrast? Read the explanation The total remains unchanged because exchange terms cancel. The contrast is 0.002 times the natural exponential of minus two, about 0.0002707. Exchange redistributes the element without creating or destroying it. 2. Does the changing fitted merger delay in Table 2 mean that Type Ia explosions physically trigger neutron-star mergers? Read the explanation No. The delay parameters and the refitted formation history are linked within a chemical fit. A parameter correlation does not establish a causal trigger between astrophysical events. 3. A mixed source gives a higher likelihood in one galaxy. Has a unique r-process origin been found? Read the explanation No. Other galaxies, event yields, occurrence assumptions, upper limits, selection and mixing still matter. Conditional model preference is a research constraint, not event-by-event source identification. ## Try the prediction ### [Compare histories in one gas reservoir](https://sandboxphysics.com/labs/galactic-chemical-evolution/) 1. Compare rapid and extended formation at fixed yields. Find nearby abundance points that correspond to different times. 2. Read the gas and birth-population compositions at one time. Explain why the current one-zone export has no spatial abundance-scatter field. ## Standalone scientific materials [Download the scientific ZIP](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/chemical-evolution-science-v1.zip) · [Read setup and calculation scope](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/README.md) The offline calculation needs Python, NumPy and SciPy. It solves four histories and compares saved synthetic measurements. Optional fixed Hill extraction runs locally with pdftotext; the ZIP has no original observations, paper or browser Engine. Individual scripts, sizes and SHA-256 - [README.md](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/README.md) (5754 B) `5a638f0c540b5d489ab149d2af7b8713a376a85a66f6a1a241a57052fee5ea1a` - [LICENSE.txt](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/LICENSE.txt) (1271 B) `c8d56106a2c4041b440e94509790c6403bbad8908db7efdbf09bbfe4c893ce82` - [requirements.txt](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/requirements.txt) (27 B) `39b93fa1b8397265839a510ca7601608248a7f3524067b61667e977bf4351475` - [source-identities.json](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/source-identities.json) (1618 B) `f7d7dc4c91ee350ab7f809a2c4c54c7f2e6c650a87e3b4b7476cdaa1c59ed127` - [fixed-inputs.json](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/fixed-inputs.json) (533 B) `98aaf54949a88e6ee5e4697b5ff25050bfff5790faa17d6d6ee5b7b02f9f1ac2` - [model-family.json](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/model-family.json) (537 B) `58dac14f9a01ed44a43c4d45dc88e12aa76c9a26c5fd977940099811242751a4` - [own-observations.json](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/own-observations.json) (10029 B) `3c2f318702cf26a370d4d298b8e4ed37da26e19c5052d51ef9a0daee06dba2f0` - [compare.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/compare.py) (3237 B) `785b41dd356c57e90165197b70ec11adf6271daafe85e096e793f237bcc2c860` - [evolution.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/evolution.py) (4572 B) `4deabefc6e7087e882bf44fe43b582762ba461d08c634d7e5237615760fcf108` - [extract\_hill.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/extract_hill.py) (2549 B) `59aa4cd821574c65110549c41e513fe7161210789c319eeb9342b26d659598d3` - [fetch\_hill.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/fetch_hill.py) (1705 B) `d31ccb38d10b8cd8a90194ba6a6d5b3e97af81b41c3dbef903c48ecc9b413696` - [io\_helpers.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/io_helpers.py) (1894 B) `2a07f3ab2f0792f19d5575e5fbb38a6d1370df6e01c20c1cde0a3e2fa072a959` - [make\_observations.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/make_observations.py) (1327 B) `5498e13a03606dc0adf97488f0f1d1bef92f04ed275aa16f7fe111a0d5683e79` - [make\_records.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/make_records.py) (1461 B) `3f25e9ee127836c6e2ade3f922de38941ca941bf26ebd9c85bf21538655b611a` - [make\_reference.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/make_reference.py) (4034 B) `cb0b779f168fe6e8ec93614da0eae1de7f82c37d61a632f05b5e987fe4e14f69` - [observation\_quadrature.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/observation_quadrature.py) (1662 B) `5b1ad8917ef50bff8a63248c883e5b56f72875eaa49470b41d98db44ca21186a` - [plot.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/plot.py) (3976 B) `55f9242160e7bd86321117136b5fc12a735c39cb1d5ae2ccd57481f76fb2ec6d` - [validate.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/validate.py) (7153 B) `1e089465eae28b0bcbf8bf9e696c3ccd7e07a645bcf7173cb536c72e878e5c00` - [verify\_package.py](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/verify_package.py) (926 B) `16d0c1b6185e7107b33fd33639456aa4806f3984269ca671f9b1b70bf487cd3d` - [SHA256SUMS.json](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/SHA256SUMS.json) (3112 B) `28a1797407ef24a767fde830473a49ca5d8ddb675f13ce03135ba1a2fe4b567a` - [chemical-evolution-science-v1.zip](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/chemical-evolution-science-v1.zip) (26888 B) `63dc82bb31bcd276c9abc76191e3201804d13d068811c5ff27bd244b17dcc06d` [Previous: Recompute a fixed Sculptor abundance sample](https://sandboxphysics.com/stories/galactic-chemical-evolution/recompute-a-fixed-sculptor-sample/) ## Sources and calculation scope 1. [Henderson, Kirby, de los Reyes & Surman, Constraints on the r-process Delay Time from Dwarf Spheroidal Galaxies, arXiv:2610.15946v1](https://arxiv.org/abs/2610.15946v1) Full fixed 24-page preprint, both appendices, Figures 1–12 and Tables 1–4 read; source figures and tables visually inspected. Posted 11 October 2026, internal draft 14 October. Actual PDF SHA-256 is 25e460e527705c40f64feaf0bc06e0a840f8967b67fed681c3acec5ad71ebcdc; server ETag differs and is not treated as a content signature. This is a primary-paper guide, not a rerun of its multi-element sampling or an identification of r-process sites. Original PDF and figures are not redistributed. 2. [Independent one-zone chemical-evolution Python materials](https://sandboxphysics.com/stories/galactic-chemical-evolution/materials/chemical-evolution-science-v1.zip) Original continuous convolution/DOP853 producers, saved fourteen-field inventories and birth populations, finite refinement, normalized conditional mixtures and mathematical limits. Public inputs and figures are original synthetic data. Default is offline; optional fixed Hill source is fetched and analyzed on the reader’s computer. No App/Engine, original paper, observed rows or measured-derived plot is included.