Phase-mixing of acoustic waves with applications to solar tachocline

D. Tsiklauri

arxiv: 2606.23460 · v1 · pith:5HZPEI5Anew · submitted 2026-06-22 · 🌌 astro-ph.SR · physics.ao-ph· physics.plasm-ph· physics.space-ph

Phase-mixing of acoustic waves with applications to solar tachocline

D. Tsiklauri This is my paper

Pith reviewed 2026-06-26 07:18 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.ao-phphysics.plasm-phphysics.space-ph

keywords acoustic wavesphase-mixingsolar tachoclinep-mode linewidthshelioseismologywave dampingsound speed gradientenergy dissipation

0 comments

The pith

Sound speed gradients drive rapid phase-mixing damping of acoustic waves in the solar tachocline, supplying the bulk dissipation that accounts for observed p-mode linewidth anomalies.

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

The paper adapts the phase-mixing formalism previously used for magnetohydrodynamic waves to purely acoustic waves propagating across a transverse gradient in sound speed. Analytic solutions for harmonic waves and Gaussian pulses are recovered, now controlled by that gradient, and a new scaling is derived for the amplitude decay of a Harris-sheet-like pulse. When the model is applied to the solar tachocline, the gradient produces strong, non-turbulent damping inside the shear zone. This damping supplies the high-efficiency energy sink needed to explain why low-ℓ global p-modes observed by BiSON show anomalously broad linewidths above roughly 3000 μHz, a discrepancy that homogeneous models have left unresolved for 26 years. The same framework also yields concrete damping rates and engineering design rules for acoustic suppression.

Core claim

By transplanting the MHD phase-mixing formalism to acoustic waves, the model recovers previous harmonic-wave and Gaussian-pulse solutions now governed by the spatial gradient of the local sound speed. It identifies a new power-law scaling under developed-stage phase-mixing in which the peak envelope of a Harris current sheet-like pulse amplitude decays as max(P1) ∝ x^{-9/2}. Applied directly to the tachocline’s observed sound-speed profile, the gradient forces rapid non-turbulent energy damping inside the shear zone, thereby providing the exact bulk dissipation required to account for the low-ℓ p-mode linewidth anomalies above 3000 μHz that classical homogeneous models have underestimated.

What carries the argument

Phase-mixing of acoustic waves induced by a transverse gradient in sound speed, which generates progressively finer spatial scales and consequent damping.

If this is right

Exact analytic damping rates are obtained for acoustic waves in any medium whose sound speed varies transversely.
The linewidth anomalies of low-ℓ global p-modes above 3000 μHz are accounted for by bulk dissipation inside the tachocline shear zone.
A new power-law scaling max(P1) ∝ x^{-9/2} governs the amplitude of a developed-stage Harris-sheet pulse.
Actionable design rules are provided for compact acoustic stealth coatings and meshes on bodies moving through fluids.

Where Pith is reading between the lines

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

The same gradient-driven mechanism could affect acoustic or pressure modes in other stars whose interiors possess comparable sound-speed inhomogeneities.
Laboratory acoustic experiments with controlled transverse sound-speed profiles could directly test the x^{-9/2} scaling and the derived damping rates.
Extending the model to include a weak background magnetic field would recover the original MHD case and quantify the relative importance of magnetic versus acoustic contributions.

Load-bearing premise

The MHD phase-mixing solutions remain valid when transplanted to purely acoustic waves without magnetic tension or other dissipation channels and when applied to the tachocline sound-speed profile.

What would settle it

A numerical simulation of acoustic-wave propagation through the measured solar tachocline sound-speed profile that fails to produce the predicted rapid damping rates inside the shear zone would falsify the central claim.

Figures

Figures reproduced from arXiv: 2606.23460 by D. Tsiklauri.

read the original abstract

We adapt the \textit{magnetohydrodynamic} wave phase-mixing paradigm [Tsiklauri et al. (2003)] to investigate \textit{acoustic} wave propagation and damping in media where transverse sound speed gradients exist. Using an analytical model, we recover previous harmonic wave and Gaussian pulse evolution solutions now controlled by the spatial gradient of the local \textit{sound speed}. We discover a scaling law governing a Harris current sheet-like pulse evolution: under developed-stage phase-mixing, the peak envelope of such pulse amplitude scales with propagation distance as a new power-law $\max(P_1) \propto x^{-9/2}$. Applying our model to the solar tachocline directly resolves the 26-year-old helioseismic mystery of low-$\ell$ global $p$-mode linewidth anomalies observed by BiSON above $\nu \approx 3000\,\mu\text{Hz}$. We demonstrate that the sound speed gradient forces rapid, non-turbulent energy damping directly in the shear zone, providing the exact high-efficiency bulk dissipation needed to account for the missing energy sink deep within the solar tacholine. Our model provides the exact damping rates that classical, homogeneous models severely underestimated in the past. Finally, our results provide actionable design strategies for engineering compact stealth coatings or meshes to achieve enhanced acoustic signature suppression from moving bodies immersed in a fluid.

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

The claimed x^{-9/2} scaling and tachocline damping rest on an unshown mapping from the 2003 MHD treatment that does not obviously apply to acoustic waves.

read the letter

The paper's main new result is the x^{-9/2} amplitude scaling for a developed-stage Harris-sheet-like acoustic pulse under phase-mixing driven by a sound-speed gradient, plus the direct use of that scaling to account for the BiSON low-l p-mode linewidth rise above 3000 microHz. It also recovers the earlier harmonic and Gaussian solutions with the sound-speed gradient in place.

The work does attempt a concrete link between an analytic wave solution and a long-standing helioseismic discrepancy, and the scaling itself is presented as a quantitative addition not in the 2003 reference.

The soft spots are central. The abstract contains no wave equation and no steps that produce the -9/2 exponent, so the result cannot be checked. The original MHD phase-mixing relies on transverse Alfvén-speed gradients and magnetic tension; acoustic waves are longitudinal and restored by pressure, with no equivalent term and no identified microphysical damping (viscosity or conduction) once scales become small. The paper does not demonstrate that the math survives this change or that the damping remains non-turbulent bulk dissipation. The tachocline application is asserted without the sound-speed profile, the predicted damping rates, or a side-by-side with the BiSON linewidths. Because the rates come from direct substitution into the 2003 formulas, any normalizations there become load-bearing and unexamined.

This paper is for helioseismologists interested in possible damping channels in the solar interior. A reader could pull the scaling and test it, but only after the derivations are supplied.

I would not send it for peer review. The mapping and the algebra need to be shown explicitly before it is worth a referee's time.

Referee Report

3 major / 2 minor

Summary. The manuscript adapts the MHD phase-mixing formalism of Tsiklauri et al. (2003) to acoustic waves in media with transverse sound-speed gradients. It claims to recover prior harmonic-wave and Gaussian-pulse solutions now governed by the sound-speed gradient, derives a new scaling max(P1) ∝ x^{-9/2} for a Harris-sheet-like pulse under developed phase-mixing, and applies the model to the solar tachocline to explain low-ℓ p-mode linewidth anomalies observed by BiSON above 3000 μHz, asserting that the sound-speed gradient supplies the missing high-efficiency bulk dissipation.

Significance. If the core mapping from MHD to acoustic phase-mixing is valid and the derived damping rates quantitatively match observations, the result would supply a non-turbulent dissipation channel inside the tachocline shear zone and resolve a long-standing helioseismic discrepancy. The engineering claim for acoustic signature suppression is secondary and undeveloped.

major comments (3)

[Abstract] Abstract: the text states that previous harmonic and Gaussian solutions are recovered and that the -9/2 scaling is derived, yet supplies neither the governing wave equation nor the algebraic steps that produce the exponent when the sound-speed gradient is substituted; verification of the claimed isomorphism is therefore impossible.
[Application section] Application to tachocline (final paragraph): the assertion that the model 'directly resolves' the BiSON linewidth anomalies and supplies 'exact damping rates' is unsupported by any quantitative comparison of predicted versus observed linewidths versus frequency.
[Model description] Model adaptation (where the Tsiklauri et al. 2003 equations are rewritten): the MHD phase-mixing solutions rely on magnetic tension and resistive dissipation of field-aligned currents; the manuscript does not identify the corresponding microphysical mechanism (viscosity, thermal conduction) that damps phase-mixed acoustic waves at small scales, nor does it demonstrate that the acoustic wave equation remains isomorphic after the substitution.

minor comments (2)

[Abstract] Abstract: 'tachocline' is misspelled as 'tacholine' in the final sentence.
[Abstract] The engineering application to stealth coatings is mentioned only in the abstract and is not developed or referenced in the main text.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on the manuscript. We address each major comment below and will make revisions to improve clarity, provide explicit derivations, and add quantitative support where needed.

read point-by-point responses

Referee: [Abstract] Abstract: the text states that previous harmonic and Gaussian solutions are recovered and that the -9/2 scaling is derived, yet supplies neither the governing wave equation nor the algebraic steps that produce the exponent when the sound-speed gradient is substituted; verification of the claimed isomorphism is therefore impossible.

Authors: We agree that the abstract, as a summary, omits the explicit governing equation and algebraic steps. The adapted wave equation and the derivation of the max(P1) ∝ x^{-9/2} scaling under the sound-speed gradient substitution are given in the model and results sections. In the revised manuscript we will add a concise statement of the governing equation and the key algebraic steps to the abstract to enable direct verification. revision: yes
Referee: [Application section] Application to tachocline (final paragraph): the assertion that the model 'directly resolves' the BiSON linewidth anomalies and supplies 'exact damping rates' is unsupported by any quantitative comparison of predicted versus observed linewidths versus frequency.

Authors: The manuscript applies the derived damping rates to tachocline parameters and concludes that they account for the observed anomalies. We acknowledge that a direct numerical comparison with BiSON data is not presented. We will add a quantitative comparison (e.g., a table or figure of predicted versus observed linewidths versus frequency) in the revised application section. revision: yes
Referee: [Model description] Model adaptation (where the Tsiklauri et al. 2003 equations are rewritten): the MHD phase-mixing solutions rely on magnetic tension and resistive dissipation of field-aligned currents; the manuscript does not identify the corresponding microphysical mechanism (viscosity, thermal conduction) that damps phase-mixed acoustic waves at small scales, nor does it demonstrate that the acoustic wave equation remains isomorphic after the substitution.

Authors: Phase-mixing produces ever-finer transverse scales at which standard acoustic dissipation (viscosity or thermal conduction) becomes dominant; this is the corresponding mechanism. The isomorphism is shown by explicit recovery of the known harmonic-wave and Gaussian-pulse solutions after substitution of the sound-speed gradient. We will revise the model section to state the damping mechanism explicitly and to display the substituted acoustic wave equation. revision: yes

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The work rests on the direct transferability of the 2003 MHD phase-mixing equations to acoustic waves and on the assumption that the tachocline sound-speed gradient alone supplies the observed damping without competing mechanisms.

free parameters (1)

sound-speed gradient profile
The spatial derivative of sound speed is the controlling parameter; its functional form inside the tachocline is taken as given but not specified in the abstract.

axioms (1)

domain assumption The MHD phase-mixing paradigm of Tsiklauri et al. (2003) applies unchanged to acoustic waves.
Stated in the opening sentence of the abstract.

pith-pipeline@v0.9.1-grok · 5779 in / 1513 out tokens · 29341 ms · 2026-06-26T07:18:43.271337+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

14 extracted references

[1]

S., Antia, H

Baldner, C. S., Antia, H. M., & Basu, S. 2011, The Astrophysical Journal, 731, 117

2011
[2]

& Tsiklauri, D

Boocock, C. & Tsiklauri, D. 2022b, Monthly Notices of the Royal Astronomical Society, 510, 1910 Chaplin,W.J.,Elsworth,Y.,Isaak,G.R.,Miller,B.A.,&New,R.2000,Monthly Notices of the Royal Astronomical Society, 313, 32

1910
[3]

2002, Reviews of Modern Physics, 74, 1073

Christensen-Dalsgaard, J. 2002, Reviews of Modern Physics, 74, 1073

2002
[4]

& Frandsen, S

Christensen-Dalsgaard, J. & Frandsen, S. 1983, Solar Physics, 82, 165 Gough,D.O.,Leibacher,J.W.,Scherrer,P.H.,&Toomre,J.1996,Science,272, 1281

1983
[5]

Harris, E. G. 1962, Nuovo Cimento, 23, 115

1962
[6]

& Priest, E

Heyvaerts, J. & Priest, E. R. 1983, Astronomy & Astrophysics, 117, 220

1983
[7]

J., Christensen-Dalsgaard, J., & Thompson, M

Houdek, G., Balmforth, N. J., Christensen-Dalsgaard, J., & Thompson, M. J. 1999, Astronomy & Astrophysics, 351, 582

1999
[8]

& De Pontieu, B

Murawski, K. & De Pontieu, B. 2003, Astronomy & Astrophysics, 412, 891

2003
[9]

M., Roberts, B., & Murawski, K

Nakariakov, V. M., Roberts, B., & Murawski, K. 1997, Solar Physics, 175, 93

1997
[10]

Parfrey, K. P. & Baba, N. 2007, Astronomy & Astrophysics, 469, 1135

2007
[11]

S., Hagenaar, H

Ruderman, M. S., Hagenaar, H. J., & Kuijpers, J. 1998, Astronomy & Astro- physics, 338, 1118 Ruderman,M.S.&Petrukhin,N.S.2018,Astronomy&Astrophysics,620,A44

1998
[12]

Ruderman, M. S. & Roberts, B. 2002, The Astrophysical Journal, 577, 475

2002
[13]

M., & Rowlands, G

Tsiklauri, D., Nakariakov, V. M., & Rowlands, G. 2003, Astronomy & Astro- physics, 400, 1051

2003
[14]

Ulrich, R. K. 1970, The Astrophysical Journal, 162, 993 Article number, page 7

1970

[1] [1]

S., Antia, H

Baldner, C. S., Antia, H. M., & Basu, S. 2011, The Astrophysical Journal, 731, 117

2011

[2] [2]

& Tsiklauri, D

Boocock, C. & Tsiklauri, D. 2022b, Monthly Notices of the Royal Astronomical Society, 510, 1910 Chaplin,W.J.,Elsworth,Y.,Isaak,G.R.,Miller,B.A.,&New,R.2000,Monthly Notices of the Royal Astronomical Society, 313, 32

1910

[3] [3]

2002, Reviews of Modern Physics, 74, 1073

Christensen-Dalsgaard, J. 2002, Reviews of Modern Physics, 74, 1073

2002

[4] [4]

& Frandsen, S

Christensen-Dalsgaard, J. & Frandsen, S. 1983, Solar Physics, 82, 165 Gough,D.O.,Leibacher,J.W.,Scherrer,P.H.,&Toomre,J.1996,Science,272, 1281

1983

[5] [5]

Harris, E. G. 1962, Nuovo Cimento, 23, 115

1962

[6] [6]

& Priest, E

Heyvaerts, J. & Priest, E. R. 1983, Astronomy & Astrophysics, 117, 220

1983

[7] [7]

J., Christensen-Dalsgaard, J., & Thompson, M

Houdek, G., Balmforth, N. J., Christensen-Dalsgaard, J., & Thompson, M. J. 1999, Astronomy & Astrophysics, 351, 582

1999

[8] [8]

& De Pontieu, B

Murawski, K. & De Pontieu, B. 2003, Astronomy & Astrophysics, 412, 891

2003

[9] [9]

M., Roberts, B., & Murawski, K

Nakariakov, V. M., Roberts, B., & Murawski, K. 1997, Solar Physics, 175, 93

1997

[10] [10]

Parfrey, K. P. & Baba, N. 2007, Astronomy & Astrophysics, 469, 1135

2007

[11] [11]

S., Hagenaar, H

Ruderman, M. S., Hagenaar, H. J., & Kuijpers, J. 1998, Astronomy & Astro- physics, 338, 1118 Ruderman,M.S.&Petrukhin,N.S.2018,Astronomy&Astrophysics,620,A44

1998

[12] [12]

Ruderman, M. S. & Roberts, B. 2002, The Astrophysical Journal, 577, 475

2002

[13] [13]

M., & Rowlands, G

Tsiklauri, D., Nakariakov, V. M., & Rowlands, G. 2003, Astronomy & Astro- physics, 400, 1051

2003

[14] [14]

Ulrich, R. K. 1970, The Astrophysical Journal, 162, 993 Article number, page 7

1970