A user-friendly package and workflow for generating effective homogeneous rheologies for the study of the long-term orbital evolution of multilayered planetary bodies
Pith reviewed 2026-05-15 11:29 UTC · model grok-4.3
The pith
A package automates reduction of layered planetary interiors to equivalent homogeneous generalized Voigt rheologies for tidal evolution studies.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The package identifies the secular relaxation poles and residues of a layered Maxwell body's tidal Love number response and inverts them to obtain the parameters of an equivalent homogeneous generalized Voigt body, delivering elastic, gravitational, viscous, and Voigt-element contributions in a form ready for numerical tidal-evolution codes.
What carries the argument
Numerical identification of secular relaxation poles and residues from the propagator-matrix Love number, followed by inversion into generalized Voigt compliance.
If this is right
- The returned parameters can be inserted directly into existing time-domain orbital and spin-evolution codes without repeated layered boundary-value solves.
- Dominant-mode selection yields reduced models that retain tidal response accuracy over a user-specified frequency interval while using fewer relaxation elements.
- The method supplies a reproducible bridge between detailed stratified interior models and the homogeneous rheologies required by long-term planetary evolution simulations.
- Application to a five-layer lunar model produces both full and reduced representations suitable for downstream orbital studies.
Where Pith is reading between the lines
- The same pole-residue workflow could be tested on other bodies such as icy satellites to generate rheologies for their tidal-heating calculations.
- If the equivalence proves robust, existing N-body or spin-evolution integrators could incorporate more realistic interior layering through simple parameter substitution.
- The approach might be extended to compressible layers or non-Maxwell rheologies by adapting the propagator matrix and inversion steps accordingly.
Load-bearing premise
The numerical extraction of poles and residues from the layered Love-number response accurately reproduces the original viscoelastic behavior over the frequency range of orbital evolution without significant truncation or artifacts.
What would settle it
Comparison of the frequency-dependent Love numbers or dissipation rates produced by the equivalent homogeneous model against those of the original layered model at several frequencies matching typical orbital periods would falsify the equivalence if the values differ beyond a chosen tolerance.
read the original abstract
We present a user-friendly, open-source Wolfram Language package that automates the construction of an effective homogeneous generalized Voigt rheology for a spherically symmetric, incompressible layered body with Maxwell solid layers. It provides a practical bridge between layered interior models and time-domain simulations of tidal evolution. The package combines three components: (i) a forward computation of the degree-2 tidal Love number based on the propagator-matrix formulation for incompressible stratified viscoelastic bodies; (ii) numerical identification of the secular relaxation poles and residues of the layered model; and (iii) inversion of the resulting response into the compliance of an equivalent homogeneous generalized Voigt body. The implementation is based on the equivalence established for multilayer Maxwell bodies and includes an optional dominant-mode selection procedure for obtaining reduced rheological models over a prescribed frequency range. The package returns the parameters of the equivalent homogeneous model, including elastic, gravitational, viscous, and Voigt-element contributions, in a format suitable for downstream numerical applications. As a case study, we apply the package to a five-layer lunar interior model and obtain its equivalent generalized Voigt representation, together with a reduced model that preserves the tidal response over the frequency interval relevant for orbital evolution while using fewer relaxation elements. This package makes the reduction from stratified viscoelastic interiors to effective homogeneous rheologies reproducible and accessible. It allows physical tidal dissipation models to be used in long-term orbital and spin-evolution studies without having to repeatedly solve the full layered boundary-value problem.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a Wolfram Language package that automates construction of an effective homogeneous generalized Voigt rheology for spherically symmetric incompressible layered Maxwell viscoelastic bodies. It combines propagator-matrix computation of degree-2 Love numbers, numerical extraction of secular relaxation poles and residues, and inversion to equivalent Voigt compliance parameters, with an optional dominant-mode reduction procedure. The package is demonstrated on a five-layer lunar interior model, returning elastic, gravitational, viscous, and Voigt-element parameters in a form suitable for time-domain orbital evolution codes.
Significance. If the numerical pole identification and inversion steps faithfully reproduce the layered response, the package supplies a reproducible bridge between stratified interior models and long-term tidal simulations, allowing physical dissipation to be incorporated without repeated solution of the full boundary-value problem.
major comments (1)
- [Lunar case study] Lunar case study: the manuscript reports only the resulting Voigt parameters and a reduced-mode variant but supplies no tabulated or plotted comparison of k2(ω) (real and imaginary parts) between the original layered propagator-matrix model and the reconstructed homogeneous model across the 10^{-12}–10^{-4} rad s^{-1} interval. Without such a metric (e.g., RMS error or maximum deviation), truncation of weak poles or conditioning issues in the residue fit could produce systematic deviations in dissipation, undermining the asserted practical utility for orbital evolution studies.
minor comments (1)
- [Methods] The statement that the implementation 'is based on the equivalence established for multilayer Maxwell bodies' should be accompanied by an explicit citation to the prior reference establishing this equivalence, placed in the methods section describing the inversion step.
Simulated Author's Rebuttal
We thank the referee for their constructive comment, which identifies a valuable opportunity to strengthen the validation of our package. We address the point below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Lunar case study] Lunar case study: the manuscript reports only the resulting Voigt parameters and a reduced-mode variant but supplies no tabulated or plotted comparison of k2(ω) (real and imaginary parts) between the original layered propagator-matrix model and the reconstructed homogeneous model across the 10^{-12}–10^{-4} rad s^{-1} interval. Without such a metric (e.g., RMS error or maximum deviation), truncation of weak poles or conditioning issues in the residue fit could produce systematic deviations in dissipation, undermining the asserted practical utility for orbital evolution studies.
Authors: We agree that a direct quantitative comparison of k2(ω) is necessary to confirm the accuracy of the reduction and to rule out systematic errors from pole truncation or numerical conditioning. In the revised manuscript we will add a new figure showing the real and imaginary parts of k2(ω) computed from the original five-layer propagator-matrix model and from the reconstructed homogeneous generalized Voigt model over the full frequency interval 10^{-12} to 10^{-4} rad s^{-1}. The figure will be accompanied by a table reporting the root-mean-square error and maximum absolute deviation in both real and imaginary parts. These additions will be placed in the lunar case-study section and will not change the underlying methodology or results, but will provide the explicit metric requested. revision: yes
Circularity Check
No significant circularity; workflow implements established equivalence
full rationale
The paper's chain consists of (i) propagator-matrix computation of degree-2 Love numbers for the layered incompressible Maxwell model, (ii) numerical extraction of secular poles/residues, and (iii) inversion to generalized Voigt compliance parameters. These operations are presented as a direct implementation of a prior-established mathematical equivalence for multilayer Maxwell bodies, not derived or fitted within the present manuscript. No step equates a fitted quantity to a prediction by construction, invokes a self-citation as the sole justification for uniqueness, or renames an input as output. The central deliverable is a reproducible software workflow whose numerical fidelity can be checked externally against the original layered response; the derivation therefore remains self-contained and non-circular.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Spherically symmetric incompressible stratified viscoelastic body composed of Maxwell solid layers
- domain assumption Equivalence of the tidal response of multilayer Maxwell bodies to that of a homogeneous generalized Voigt body
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
numerical identification of the secular relaxation poles and residues of the layered model; and (iii) inversion of the resulting response into the compliance of an equivalent homogeneous generalized Voigt body
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
generalized Voigt rheology ... C(s)=1/α+γ + Cv(s)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.