🔢 mathematics

Explicit form of relaxation tensor for isotropic extended Burgers model and its spectral inversion

This paper presents a systematic analytical method to derive the explicit relaxation tensor for the inhomogeneous isotropic quasi-static extended Burgers model (q-EBM) and utilizes it to compute and invert the eigenvalue clusters of Earth's free oscillations.

Original authors: Youjun Deng, Ching-Lung Lin, Gen Nakamura

Published 2026-02-12

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Youjun Deng, Ching-Lung Lin, Gen Nakamura

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to understand how the Earth "rings" like a giant bell after a massive earthquake.

If the Earth were a solid, perfect steel ball, it would ring with very clear, sharp notes. But the Earth isn't steel; it’s "squishy" and "gooey" in places (like the mantle). This "squishiness" is what scientists call anelasticity. Because the Earth is a bit like a mix of a solid and a thick liquid, it doesn't just ring; it absorbs some of the sound, making the notes fade and change in complex ways.

This paper is a mathematical "instruction manual" for describing that complex, squishy ringing. Here is the breakdown:

1. The Model: The "Springs and Honey" Machine

The researchers use something called the Extended Burgers Model (EBM).

Think of the Earth’s interior as a massive, invisible machine made of millions of tiny springs and jars of honey (dashpots) connected together.

The Springs: Represent the "elastic" part. When you push them, they snap back instantly. This is the part that creates the "ringing" sound.
The Honey: Represents the "viscous" part. When you push it, it resists, but it also flows slowly, absorbing energy and making the movement "sluggish."

By adding more springs and more jars of honey in a specific sequence, the scientists can create a model that perfectly mimics how the Earth’s mantle behaves—capturing both the quick snaps and the slow, gooey flows.

2. The Problem: The Math was "Too Messy"

Before this paper, scientists knew this "springs and honey" model was useful, but the math required to calculate the Earth's specific "notes" (called eigenvalues) was a nightmare. It was like trying to calculate the exact melody of a piano played underwater; the "water" (the gooiness) makes the math incredibly tangled and difficult to solve directly.

3. The Breakthrough: The "Volume vs. Shape" Trick

The authors found a clever way to untangle the mess. They realized that when the Earth vibrates, it does two things at once:

The Squeeze (Volumetric): It gets slightly denser/thicker (like squeezing a sponge).
The Stretch (Deviatoric): It changes shape without changing its volume (like squeezing a ball of dough between your fingers).

In the past, these two movements were mathematically "tangled" together. The authors developed a new method to separate the squeeze from the stretch. By treating them as two different mathematical tracks, they were able to find an "explicit form"—essentially a clean, usable formula—for how the Earth's "squishiness" affects its ringing.

4. The "C-ev" Clusters: Finding the Hidden Pattern

Because the Earth is so complex, the "notes" it plays aren't single, lonely tones. Instead, they come in clusters (which the authors call C-ev’s).

Imagine instead of a single note on a flute, you hear a "chord" of several notes played at once. The authors proved that these clusters aren't random; they follow a very specific mathematical pattern. They showed that if you know the "notes" in these clusters, you can actually work backward to figure out exactly how much "spring" and how much "honey" is inside the Earth.

Summary: Why does this matter?

If we want to understand the Earth's interior—to better predict how seismic waves move or how the mantle flows—we need to know exactly how "squishy" it is.

This paper provides the mathematical lens that allows geophysicists to look at the "ringing" of the Earth and say: "Aha! Based on those specific clusters of notes, I can tell you exactly how much honey and how many springs are hidden deep beneath our feet."

This paper, titled "Explicit form of relaxation tensor for isotropic extended Burgers model and its spectral inversion," presents a mathematical framework for modeling the anelasticity of the Earth using the Extended Burgers Model (EBM). It addresses a significant gap in geophysical modeling by providing an analytical method to compute the eigenvalues of the Earth's free oscillations (FOE) when modeled as a viscoelastic medium.

1. The Problem

The Earth's mantle exhibits complex viscoelastic behavior that standard models often fail to capture fully. The Extended Burgers Model (EBM) is an integro-differential system capable of reconciling the mechanical behavior of upper mantle materials with the Earth as a whole. However, prior to this work, there was no systematic, analytical way to compute the eigenvalues of the FOE when using the EBM. Specifically, the researchers sought to:

Derive an explicit form of the relaxation tensor for inhomogeneous isotropic EBM.
Identify and analyze Clusters of Eigenvalues (C-ev’s)—specific groups of eigenvalues that characterize the system's oscillations.
Solve the inverse problem: determining the physical parameters of the EBM (viscosities and elasticities) from the observed C-ev’s.

2. Methodology

The authors employ several advanced mathematical techniques to bridge the gap between the integro-differential EBM and solvable eigenvalue problems:

Volumetric/Deviatoric Decomposition: To handle the complexity of the three-dimensional tensors, the authors decompose the strain and stress into volumetric (pressure-related) and deviatoric (shape-changing) parts. This allows the system to be broken down into two decoupled ordinary differential systems, making the derivation of the relaxation tensor tractable.
Quasi-Static EBM (q-EBM) Approach: The authors transform the Boltzmann-type viscoelastic system (BVS) into a quasi-static form. By integrating the system by parts, they separate the operator into an instantaneous term (elastic response) and a memory term (viscous response).
Simultaneous Eigenvalue Problem: To find the C-ev’s, the authors seek special radial eigenfunctions that are simultaneously eigenfunctions of both the volumetric and deviatoric operators while satisfying traction-free boundary conditions on the Earth's surface. This is achieved using spherical Bessel functions.
Augmented System & Characteristic Polynomials: The integro-differential equation is converted into a first-order system of partial differential equations (an "augmented system"). The eigenvalues are then found as the roots of a high-order characteristic polynomial.

3. Key Contributions

Explicit Relaxation Tensor: The paper provides the first explicit mathematical expression for the relaxation tensor of a 3D inhomogeneous isotropic EBM. This is expressed through a sum of exponential terms involving the eigenvalues of specific matrix blocks ( $L_0$ and $L_1$ ).
Analytical Characterization of C-ev’s: The authors prove that the eigenvalues of the FOE under the EBM model are not arbitrary but are roots of a specific polynomial. They demonstrate that these eigenvalues "cluster" around the eigenvalues of the instantaneous elastic part.
Spectral Inversion Formula: The paper establishes a rigorous method to perform spectral inversion. It proves that by knowing at least two clusters of eigenvalues, one can uniquely recover the relaxation parameters of the model.

4. Results

Mathematical Structure: The authors show that the C-ev’s consist of $2n+2$ real roots that interlace with the relaxation rates ( $\beta_j$ ), plus two additional complex roots.
Parameter Recovery: They derive specific formulas (e.g., equations 5.5 and 5.6 in the text) that allow for the calculation of the Lamé moduli ( $\lambda_0, \mu_0$ ) and the individual spring/dashpot constants from the spectral data.
Eigenfunction Identification: They successfully identified that the eigenfunctions of the C-ev’s share the same spatial structure as the eigenfunctions of the instantaneous elastic part, simplifying the search for solutions.

5. Significance

This research is highly significant for both geophysics and material science.

Geophysics: It provides a more accurate mathematical tool for interpreting the free oscillations of the Earth, allowing scientists to better understand the viscoelastic properties of the mantle.
Material Science: The methodology for spectral inversion and the use of the EBM can be applied to study the glassy states and anelasticity of various complex materials.
Theoretical Advancement: By providing an explicit form for the relaxation tensor, the authors have moved the EBM from a theoretical concept to a practical, computable tool for studying complex viscoelastic systems.