arxiv: 2605.03810 · v1 · submitted 2026-05-05 · 🌌 astro-ph.SR · physics.flu-dyn

Recognition: unknown

Unifying Transport Models of Thermohaline Convection in Stars

Valentin A. Skoutnev

Authors on Pith no claims yet

Pith reviewed 2026-05-07 04:07 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.flu-dyn

keywords thermohaline convectionstellar mixingturbulent diffusionlinear stability analysisdouble-diffusive instabilitystellar interiorschemical transportunstable modes

0 comments

The pith

The full spectrum of unstable modes lets thermohaline diffusion interpolate between competing stellar mixing models.

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

The efficiency of thermohaline convection as a mixing process in stars remains unsettled because different theories predict turbulent diffusion coefficients that differ by orders of magnitude. This paper completes the linear stability theory for the instability and identifies two distinct classes of unstable modes: slow-growing ones at large length scales and fast-growing ones at small length scales. Reevaluating the diffusion coefficient D_μ from the entire spectrum shows that it interpolates between the previously proposed scalings depending on the local instability conditions. Resolving this would allow stellar evolution calculations to use a consistent mixing prescription rather than choosing between incompatible limits.

Core claim

By accounting for the complete set of linearly unstable modes in thermohaline convection, the associated turbulent diffusion coefficient can be shown to transition continuously between the different scalings advanced in earlier work as the controlling parameters vary.

What carries the argument

The effective turbulent diffusion coefficient D_μ constructed by integrating contributions across the spectrum of linear growth rates and wavelengths of both slow large-scale and fast small-scale modes.

If this is right

Stellar models gain a single formula for thermohaline mixing that automatically selects the appropriate scaling based on local conditions.
The dominant contribution to mixing shifts from slow to fast modes or vice versa as parameters change.
Predicted chemical profiles in stars become more reliable across different evolutionary stages.
Discrepancies between theory and observations of surface abundances can be addressed by using the interpolated value.

Where Pith is reading between the lines

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

High-resolution simulations could test whether the actual nonlinear transport follows the linear-mode prediction or requires additional saturation mechanisms.
The same mode-spectrum approach might apply to analogous double-diffusive instabilities in planetary oceans or atmospheres.
Observational tests could involve comparing models with and without the unified D_μ against asteroseismic or abundance data from specific stars.

Load-bearing premise

The actual turbulent transport rate is determined directly by the properties of the linearly unstable modes without needing to resolve nonlinear saturation or inter-mode coupling.

What would settle it

Numerical experiments that compute the mixing rate in the regime where both slow and fast modes are unstable and check whether the measured diffusion matches the interpolated expression or reverts to one of the limiting scalings.

Figures

Figures reproduced from arXiv: 2605.03810 by Valentin A. Skoutnev.

**Figure 1.** Figure 1: Thermohaline growth rate λ versus wavenumber k for representative parameters inside a red giant (NT = 0.1 s−1 , ν = 100 cm2 s −1 , P r = 10−6 , τ = 10−7 ). The growth rate is obtained by numerically solving Equation 1 for different density ratios R0 = N 2 T /N2 µ (colored curves). It attains a maximum value λ ≈ NT /R1/2 0 across a range of wavenumbers kT < k < kν when R0 < P r−1 . For larger density ratios… view at source ↗

**Figure 2.** Figure 2: Turbulent diffusion coefficient Dµ(k) = λ(k)/k2 for representative parameters inside a red giant (same as in view at source ↗

read the original abstract

Thermohaline convection is a standard chemical mixing process in stellar interiors, yet its mixing efficiency is not fully settled. Competing theories predict turbulent diffusion coefficients, $D_\mu$, that can differ by orders of magnitude, leading to uncertainties in stellar models and interpretations of observations. This paper explores a potential resolution to existing discrepancies. We first complete the linear stability theory and identify two types of unstable modes: slow growing modes at large length scales and fast growing modes at small length scales. We then reevaluate $D_\mu$ considering the full spectrum of unstable modes and find that it can self-consistently interpolate between previously proposed theoretical scalings across the instability parameter space. The question of thermohaline mixing efficiency in stars may be settled by future simulations that quantify the scale-dependent contributions of fast and slow modes to $D_\mu$ and determine how the modes dominating the transport change across parameter space.

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.

Desk Editor's Note private letter to a colleague

Completes the linear mode analysis with slow and fast branches but the D_μ interpolation still rests on an untested mapping from growth rates to transport.

read the letter

The main thing to know is that this paper finishes the linear stability analysis for thermohaline convection by identifying slow-growing modes at large scales and fast-growing modes at small scales, then argues that the full spectrum produces a D_μ that interpolates between earlier competing expressions. The linear part is the real advance here. The work systematically derives the dispersion relation and shows how both mode branches exist across the instability parameter space. Earlier theories each appear to have captured one limit, and this calculation makes the connection explicit. That clarification is useful and fills a gap in the linear theory. The citations to the prior models are appropriate and the approach stays within standard linear analysis. The soft spot is the step from the mode spectrum to D_μ. Linear growth rates and wavenumbers do not automatically yield a turbulent diffusivity; some closure is required for how the modes saturate and contribute to mixing. The abstract does not supply the explicit weighting or mapping rule, and the text itself states that future simulations are needed to determine which modes dominate the transport and how that changes with parameters. This makes the interpolation a plausible suggestion rather than a result that follows rigorously from the linear calculation alone. No nonlinear validation or saturation model is presented to support the claim that the linear scales set the effective transport across regimes. The paper is aimed at researchers who build stellar evolution models and need a consistent treatment of thermohaline mixing for chemical abundances. A reader working on the linear foundations will get value from the mode classification. Anyone hoping for a ready-to-use mixing prescription will still need the simulations the authors flag. I would send this to peer review. The linear completion is careful enough to deserve referee input, and the interpolation claim can be tested for the strength of its assumptions.

Referee Report

2 major / 2 minor

Summary. The manuscript completes the linear stability analysis of thermohaline convection in stellar interiors, identifying two classes of unstable modes: slow-growing modes at large scales and fast-growing modes at small scales. It then reevaluates the turbulent diffusion coefficient D_μ over the full spectrum of these modes and claims that the resulting expression self-consistently interpolates between previously proposed theoretical scalings for mixing efficiency across the relevant instability parameter space. The paper concludes that the question of thermohaline mixing efficiency may be settled by future simulations that determine the scale-dependent contributions of the modes.

Significance. If the interpolation is shown to follow rigorously from the linear mode spectrum, the result would unify disparate predictions for D_μ that currently differ by orders of magnitude, thereby reducing systematic uncertainties in stellar evolution calculations of chemical transport and surface abundance anomalies.

major comments (2)

[abstract and section on reevaluation of D_μ] The explicit mapping from the linear growth rates σ(k) and wavenumbers k of the unstable modes to the turbulent diffusivity D_μ is not derived or stated. The abstract and the reevaluation section assert that considering the full spectrum yields an interpolating D_μ, but no functional form (e.g., a weighted integral ∫ σ(k)/k² dk over the unstable band or a dominant-mode selection rule) is supplied. Without this step the interpolation claim cannot be verified and remains an extrapolation rather than a closed derivation. This is load-bearing for the central result.
[concluding paragraph] The manuscript acknowledges that nonlinear saturation, mode interactions, and secondary instabilities are not treated, yet proceeds to equate linear quantities directly to turbulent transport. The concluding paragraph states that future simulations are required to quantify which modes dominate D_μ, confirming that the interpolation rests on an untested assumption that the linear spectrum sets the effective diffusivity scaling across parameter space. This assumption is central to the unification claim and requires either an explicit closure or a demonstration that nonlinear effects do not alter the scaling.

minor comments (2)

[abstract] The abstract would be strengthened by including the explicit expression (or at least the functional dependence) for the interpolated D_μ rather than describing the procedure only in words.
[linear stability section] Notation for the instability parameter space and the definitions of the slow and fast mode branches should be introduced with a single equation reference early in the linear-stability section to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments correctly identify areas where the connection between the linear mode spectrum and the turbulent diffusivity requires greater explicitness, and where the limitations of the linear analysis should be emphasized. We address each major comment below and describe the revisions that will be incorporated in the next version of the manuscript.

read point-by-point responses

Referee: [abstract and section on reevaluation of D_μ] The explicit mapping from the linear growth rates σ(k) and wavenumbers k of the unstable modes to the turbulent diffusivity D_μ is not derived or stated. The abstract and the reevaluation section assert that considering the full spectrum yields an interpolating D_μ, but no functional form (e.g., a weighted integral ∫ σ(k)/k² dk over the unstable band or a dominant-mode selection rule) is supplied. Without this step the interpolation claim cannot be verified and remains an extrapolation rather than a closed derivation. This is load-bearing for the central result.

Authors: We agree that the explicit mapping from the linear growth rates and wavenumbers to D_μ was not stated with sufficient detail in the submitted manuscript. The reevaluation of D_μ in the paper follows the standard mixing-length approach used in prior thermohaline transport models, in which the contribution of each mode is estimated as proportional to σ(k)/k² and the effective diffusivity is obtained by integrating over the spectrum of unstable modes (both slow large-scale and fast small-scale). In the revised manuscript we will add a dedicated subsection that derives this functional form explicitly, shows the integral expression, and demonstrates analytically how the resulting D_μ interpolates between the previously discrepant scalings across the relevant range of the instability parameter. This addition will make the central claim directly verifiable from the linear analysis. revision: yes
Referee: [concluding paragraph] The manuscript acknowledges that nonlinear saturation, mode interactions, and secondary instabilities are not treated, yet proceeds to equate linear quantities directly to turbulent transport. The concluding paragraph states that future simulations are required to quantify which modes dominate D_μ, confirming that the interpolation rests on an untested assumption that the linear spectrum sets the effective diffusivity scaling across parameter space. This assumption is central to the unification claim and requires either an explicit closure or a demonstration that nonlinear effects do not alter the scaling.

Authors: We acknowledge that the present work is a linear stability analysis and does not treat nonlinear saturation or mode coupling. The manuscript unifies existing transport models, all of which themselves rest on linear-based estimates of D_μ; it does not claim to have derived a nonlinear closure. The concluding paragraph already states that future simulations are needed to determine the scale-dependent contributions of the fast and slow modes. In the revision we will strengthen this language to make explicit that the interpolation result holds within the linear framework and provides a consistent reference against which nonlinear simulations can be compared. We do not provide a nonlinear closure, as that lies outside the scope of the present study. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected in the derivation

full rationale

The paper completes the linear stability analysis to identify slow large-scale and fast small-scale unstable modes, then reevaluates D_μ over the full spectrum using standard linear growth-rate and wavenumber inputs. This produces an interpolation between prior scalings as a direct consequence of including both mode classes rather than by presupposing the target result in the inputs. No self-definitional equations, fitted parameters renamed as predictions, load-bearing self-citations, or ansatz smuggling are present; the mapping from linear quantities to diffusivity is an independent closure assumption common to the field, not a tautology. The derivation chain is therefore self-contained against external linear-analysis benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard fluid-dynamical assumptions for double-diffusive convection and on the premise that linear mode properties determine the turbulent diffusivity; no free parameters or new entities are introduced in the abstract.

axioms (2)

domain assumption The governing equations of incompressible, Boussinesq double-diffusive convection apply to thermohaline mixing in stellar interiors.
Invoked when completing the linear stability theory for the stellar context.
domain assumption The turbulent diffusion coefficient D_μ is determined by the spectrum of linearly unstable modes.
This is the key step that allows reevaluation of D_μ and the claimed interpolation.

pith-pipeline@v0.9.0 · 5449 in / 1480 out tokens · 61666 ms · 2026-05-07T04:07:32.682388+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

55 extracted references · 1 canonical work pages

[1]

C., & Mathis, S

Augustson, K. C., & Mathis, S. 2019, ApJ, 874, 83

2019
[2]

1969, JFM, 37, 289

Baines, P., & Gill, A. 1969, JFM, 37, 289

1969
[3]

J., Dempsey, A

Barker, A. J., Dempsey, A. M., & Lithwick, Y. 2014, ApJ, 791, 13

2014
[4]

B., & Bildsten, L

Bauer, E. B., & Bildsten, L. 2018, ApJL, 859, L19 —. 2019, ApJ, 872, 96

2018
[5]

2023, MNRAS, 518, 5465

Behmard, A., Sevilla, J., & Fuller, J. 2023, MNRAS, 518, 5465

2023
[6]

M., Garaud, P., & Stellmach, S

Brown, J. M., Garaud, P., & Stellmach, S. 2013, ApJ, 768, 34

2013
[7]

2010, A&A, 521, A9

Cantiello, M., & Langer, N. 2010, A&A, 521, A9

2010
[8]

2010, A&A, 522, A10

Charbonnel, C., & Lagarde, N. 2010, A&A, 522, A10

2010
[9]

2007, A&A, 467, L15

Charbonnel, C., & Zahn, J.-P. 2007, A&A, 467, L15

2007
[10]

2020, JFM, 903, A1

Cope, L., Garaud, P., & Caulfield, C. 2020, JFM, 903, A1

2020
[11]

G., Fraser, A

Cresswell, I. G., Fraser, A. E., Bauer, E. B., Anders, E. H., & Brown, B. P. 2025, ApJL, 986, L10

2025
[12]

K., Barker, A

Currie, L. K., Barker, A. J., Lithwick, Y., & Browning, M. K. 2020, MNRAS, 493, 5233

2020
[13]

Wachlin, F. C. 2013, A&A, 557, L12

2013
[14]

Woodward, P. R. 2024, MNRAS, 535, 1243

2024
[15]

Fraser, A. E. 2026, ApJL, 1001, L22

2026
[16]

2022, ApJ, 941, 164

Cantiello, M. 2022, ApJ, 941, 164

2022
[17]

E., Reifenstein, S

Fraser, A. E., Reifenstein, S. A., & Garaud, P. 2024, ApJ, 964, 184

2024
[18]

E., van Kan, A., Knobloch, E., Julien, K., & Liu, C

Fraser, A. E., van Kan, A., Knobloch, E., Julien, K., & Liu, C. 2025, JFM, 1020, R1

2025
[19]

2024, ApJL, 964, L15

Fuentes, J., Castro-Tapia, M., & Cumming, A. 2024, ApJL, 964, L15

2024
[20]

2011, ApJL, 728, L30 —

Garaud, P. 2011, ApJL, 728, L30 —. 2018, AnRFM, 50, 275

2011
[21]

2015, ApJ, 815, 42

Garaud, P., & Brummell, N. 2015, ApJ, 815, 42

2015
[22]

R., & Hurley, J

Glebbeek, E., Pols, O. R., & Hurley, J. R. 2008, A&A, 488, 1007

2008
[23]

2000, A&A, 354, 169

Gratton, R., Sneden, C., Carretta, E., & Bragaglia, A. 2000, A&A, 354, 169

2000
[24]

Z., & Garaud, P

Harrington, P. Z., & Garaud, P. 2019, ApJL, 870, L5

2019
[25]

1980, A&A, 91, 175

Kippenhahn, R., Ruschenplatt, G., & Thomas, H.-C. 1980, A&A, 91, 175

1980
[26]

1990, Stellar structure and evolution, Vol

Kippenhahn, R., Weigert, A., & Weiss, A. 1990, Stellar structure and evolution, Vol. 192 (Springer)

1990
[27]

2019, A&A, 621, A24

Lagarde, N., Reyl´ e, C., Robin, A., et al. 2019, A&A, 621, A24

2019
[28]

C., Siess, L., Church, R

Lattanzio, J. C., Siess, L., Church, R. P., et al. 2015, MNRAS, 446, 2673 Ligni` eres, F. 1999, A&A, 348, 933

2015
[29]

2021, A&A, 651, A84

Magrini, L., Lagarde, N., Charbonnel, C., et al. 2021, A&A, 651, A84

2021
[30]

Montgomery, M., & Dunlap, B. H. 2024, ApJ, 961, 197

2024
[31]

M., Frebel, A., Beers, T

Placco, V. M., Frebel, A., Beers, T. C., & Stancliffe, R. J. 2014, ApJ, 797, 21

2014
[32]

2013, Double-diffusive convection (Cambridge University Press)

Radko, T. 2013, Double-diffusive convection (Cambridge University Press)

2013
[33]

2021, ApJ, 923, 277

Renzo, M., & G¨ otberg, Y. 2021, ApJ, 923, 277

2021
[34]

2017, RSOS, 4

Salaris, M., & Cassisi, S. 2017, RSOS, 4

2017
[35]

2025, arXiv preprint arXiv:2509.18421

Schneider, F. 2025, arXiv preprint arXiv:2509.18421

work page arXiv 2025
[36]

2020, ApJL, 901, L18

Schwab, J. 2020, ApJL, 901, L18

2020
[37]

2018, ApJ, 862, 136

Sengupta, S., & Garaud, P. 2018, ApJ, 862, 136

2018
[38]

2022, MNRAS, 516, 3354

Sevilla, J., Behmard, A., & Fuller, J. 2022, MNRAS, 516, 3354

2022
[39]

A., et al

Shetrone, M., Tayar, J., Johnson, J. A., et al. 2019, ApJ, 872, 137

2019
[40]

Skoutnev, V. A. 2023, JFM, 956, A7

2023
[41]

A., & Veronis, G

Spiegel, E. A., & Veronis, G. 1960, ApJ, 131, 442

1960
[42]

Stancliffe, R. J. 2010, MNRAS, 403, 505 7

2010
[43]

J., Glebbeek, E., Izzard, R

Stancliffe, R. J., Glebbeek, E., Izzard, R. G., & Pols, O. R. 2007, A&A, 464, L57

2007
[44]

Stern, M. E. 1960, Tellus, 12, 172

1960
[45]

Stevenson, D. J. 1979, GAFD, 12, 139

1979
[46]

2022, ApJL, 935, L30 Th´ eado, S., & Vauclair, S

Tayar, J., & Joyce, M. 2022, ApJL, 935, L30 Th´ eado, S., & Vauclair, S. 2012, ApJ, 744, 123

2022
[47]

2011, ApJL, 728, L29

Traxler, A., Garaud, P., & Stellmach, S. 2011, ApJL, 728, L29

2011
[48]

Ulrich, R. K. 1972, ApJ, 172, 165

1972
[49]

2004, ApJ, 605, 874

Vauclair, S. 2004, ApJ, 605, 874

2004
[50]

C., Vauclair, G., Vauclair, S., & Althaus, L

Wachlin, F. C., Vauclair, G., Vauclair, S., & Althaus, L. G. 2022, A&A, 660, A30

2022
[51]

C., Vauclair, S., & Althaus, L

Wachlin, F. C., Vauclair, S., & Althaus, L. G. 2014, A&A, 570, A58

2014
[52]

2025, A&A, 693, A138

Wu, F., Song, H., Meynet, G., et al. 2025, A&A, 693, A138

2025
[53]

2019, JFM, 858, 228

Xie, J.-H., Julien, K., & Knobloch, E. 2019, JFM, 858, 228

2019
[54]

2017, Fluids, 2, 6

Xie, J.-H., Miquel, B., Julien, K., & Knobloch, E. 2017, Fluids, 2, 6

2017
[55]

2014, ApJ, 795, 118

Zemskova, V., Garaud, P., Deal, M., & Vauclair, S. 2014, ApJ, 795, 118

2014