i-Rheo-Tempo: A Model-Free, Quadrature-Free Reconstruction of the Shear Relaxation Modulus from Complex Viscosity

The paper introduces *i-Rheo-Tempo*, a robust, model-free, and quadrature-free method that accurately reconstructs the shear relaxation modulus from complex viscosity data by utilizing an exact second-derivative representation, thereby solving the frequency-to-time inverse problem without requiring parametric fitting or predefined spectra.

The Gist

Imagine you are trying to understand how a piece of chewing gum behaves. You can test it in two different ways:

1. The "Snap" Test (Time Domain): You pull the gum quickly and watch how it stretches and snaps back over time. This tells you exactly what the gum is doing right now.
2. The "Wiggle" Test (Frequency Domain): You wiggle the gum back and forth at different speeds (fast, slow, medium) and measure how stiff or sticky it feels at each speed. This tells you how the gum reacts to rhythms.

The Problem:
Scientists have known for decades that these two tests describe the same material. However, turning the "Wiggle" data into a "Snap" prediction is incredibly difficult. It's like trying to reconstruct a whole movie from a few scattered frames.

When scientists try to do this mathematically, the numbers get messy. The data is never perfect (it has noise), and we can't wiggle the gum infinitely fast or infinitely slow. Standard math methods often produce "ghosts"—fake wiggles and spikes in the prediction that don't actually exist in the real material. To fix this, scientists usually force the data into a pre-made "box" (a mathematical model), but that box might not fit the specific weirdness of your gum.

The Solution: i-Rheo-Tempo
The authors of this paper, Jorge Ramírez and Manlio Tassieri, have invented a new tool called i-Rheo-Tempo. Think of it as a "Magic Translator" that converts the "Wiggle" test directly into the "Snap" test without needing a pre-made box or doing messy, error-prone calculations.

Here is how it works, using simple analogies:

1. The "Slope" Trick (No More Guessing)

Imagine the "Wiggle" data is drawn as a winding mountain road.

• Old methods tried to calculate the total distance by adding up tiny steps (quadrature). If the map was a bit blurry, the steps added up to the wrong total.
• i-Rheo-Tempo doesn't care about the total distance. Instead, it looks at the steepness of the road. It asks: "Where does the road suddenly get steeper or flatter?"
• In math terms, it looks at the second derivative (the change in the slope). It turns the smooth, messy curve into a series of sharp "jumps" in steepness. Because it only cares about these jumps, it avoids the noise and errors that plague the old methods. It's like listening for the clunk of a gear shifting rather than trying to measure the smooth hum of the engine.

2. The "Bookends" Strategy (Handling the Edges)

Every experiment has a limit. You can't wiggle the gum at zero speed (eternity) or infinite speed (instant).

• The Problem: When you stop measuring, the math gets confused at the edges, creating fake "ghosts" in the final prediction.
• The Fix: i-Rheo-Tempo acts like a careful editor. Before it translates the data, it gently "anchors" the beginning and end of the story.
  • At the slow end, it uses physics to guess what the gum should be doing (like a rubber band settling down).
  • At the fast end, it smooths out the data just enough so the math doesn't trip over the edge.
  • Crucially, it only reports the final answer for the time range where it actually has data, so you never get a fake prediction for a time it hasn't measured.

3. The "Universal Translator"

The coolest part is that this tool isn't just for gum or rubber. The math behind it is so general that it works for any system where a "push" creates a "response" over time.

• It could translate the electrical "hum" of a battery into how it charges over time.
• It could translate the light "wiggle" of a crystal into how it vibrates.
• It works for everything from industrial tires to biological cells.

Why This Matters

Before this, if you wanted to know how a material would behave over time based on a frequency test, you had to guess the shape of the material's "skeleton" (the model) and hope you guessed right.

i-Rheo-Tempo removes the guesswork. It takes the raw experimental data and reconstructs the time-behavior directly, cleanly, and accurately. It's like taking a blurry, noisy photo and using a new algorithm to sharpen it into a crystal-clear picture without adding any fake details.

In a nutshell:
They found a way to listen to the "rhythm" of a material and instantly know exactly how it will "stretch and snap" in the future, without needing to force it into a pre-existing theory or getting lost in mathematical noise. It's a cleaner, smarter, and more honest way to understand how the world's materials behave.

Technical Summary

1. Problem Statement

In linear viscoelasticity, a central challenge is the reliable transformation of material functions between the frequency domain (dynamic moduli $G'(\omega)$, $G''(\omega)$) and the time domain (shear relaxation modulus $G(t)$).

• Limitations of Current Methods:
  • Direct Fourier Integration: Experimental data are finite in bandwidth, discretely sampled, and noisy. Direct numerical evaluation of Fourier integrals leads to oscillations, truncation artifacts, and high sensitivity to endpoint treatment.
  • Model-Dependent Approaches: The standard practice involves fitting data to a Generalized Maxwell model (a discrete sum of relaxation modes). While robust, this approach is inherently model-dependent; the reconstructed function reflects the imposed spectral structure as much as the experimental data.
  • Previous i-Rheo Framework: While the original i-Rheo method used a second-derivative representation to reduce artifacts, it was primarily applied to storage/loss moduli. Extending this to the complex viscosity domain to handle a wider range of experimental conditions required a new formulation.

2. Methodology: i-Rheo-Tempo

The authors introduce i-Rheo-Tempo, a method that reconstructs the shear relaxation modulus $G(t)$ directly from complex viscosity $\eta^*(\omega)$ without numerical quadrature or parametric fitting.

A. Theoretical Foundation

The method relies on an exact second-derivative representation of the inverse Fourier transform.

1. Complex Viscosity Definition: The relationship is reformulated using complex viscosity $\eta^*(\omega) = G^*(\omega)/(i\omega) = \eta'(\omega) - i\eta''(\omega)$.
2. Integration by Parts: By performing two successive integrations by parts on the inverse Fourier integral, the dependence shifts from the spectral functions themselves to their second derivatives (curvature) with respect to frequency.
   • For viscoelastic fluids ($G_e = 0$), the boundary terms vanish, yielding:
     $$G(t) = -\frac{2}{\pi t^2} \int_0^\infty \ddot{\eta}''(\omega) \sin(\omega t) d\omega$$
     (and a similar cosine form for $\eta'$).
3. Distributional Discretization: Since experimental data are discrete, the spectrum is approximated as piecewise linear.
   • In this representation, the first derivative is piecewise constant, and the second derivative is zero everywhere except at the nodes, where it manifests as Dirac delta functions representing "slope jumps" ($\Delta a_k$).
   • This converts the integral into a closed-form interval-slope sum:
     $$G(t) = \frac{2}{\pi t^2} \sum_{k} a_k [\text{Kernel}(\omega_{k+1}t) - \text{Kernel}(\omega_k t)]$$
   • Here, $a_k$ are the local slopes of the viscosity spectrum between frequency nodes. This formulation eliminates numerical quadrature entirely.

B. Practical Implementation & Boundary Conditioning

To handle the finite bandwidth of experimental data ($\omega_{min}$ to $\omega_{max}$), the method employs specific conditioning strategies:

• Low-Frequency Conditioning: The zero-frequency limit is treated explicitly. $\eta''(0)$ is set to 0 (for fluids), and $\eta'(0)$ is estimated via local fitting or a single-Maxwell fit of the terminal regime. This prevents unphysical boundary contributions that plague standard truncation methods.
• High-Frequency Completion: A modest local polynomial fit near $\omega_{max}$ is used to regularize the interpolation and suppress spurious curvature, though no physical meaning is assigned to this extension.
• Resampling: Data is interpolated in linear $\omega$ space and resampled on a logarithmic grid to ensure numerical stability.
• Validity Window: The reconstructed $G(t)$ is reported strictly within the reciprocal experimental window $[1/\omega_{max}, 1/\omega_{min}]$.
• Error Estimation: The method includes a diagnostic to estimate the "numerical fragility" of the long-time tail, which arises from the amplification of low-frequency uncertainties inherent in the second-derivative formulation.

3. Key Contributions

1. Quadrature-Free Inversion: The derivation of a closed-form, analytical expression for $G(t)$ based solely on interval slopes of the complex viscosity, removing the need for numerical integration.
2. Model Independence: The method does not assume a predefined relaxation spectrum (e.g., Maxwell modes) or parametric form, preserving the raw experimental information.
3. Robust Boundary Handling: A novel conditioning strategy that explicitly anchors the zero-frequency limit and regularizes the high-frequency end, avoiding the artifacts associated with artificial zero-slope extrapolation.
4. General Framework: While validated for rheology, the mathematical framework is applicable to any causal complex response function (e.g., dielectric permittivity, electrical impedance).

4. Results and Validation

The method was validated against a diverse set of synthetic and experimental datasets spanning nearly nine decades in frequency:

• Synthetic Benchmark (Burgers Model): Demonstrated exact reversibility. Transforming $G(t)$ back to the frequency domain recovered the original dynamic moduli without distortion, confirming internal consistency.
• Industrial Styrene-Butadiene Rubber (SBR): Showed quantitative agreement between the reconstructed $G(t)$ and direct step-strain measurements across the entire time window.
• Monodisperse Linear Polyisoprene & Polybutadiene Melts: Successfully reconstructed relaxation moduli from Time-Temperature Superposition (TTS) data. The method accurately captured the rubbery plateau and terminal regime, with reconstructed curves collapsing onto a common short-time curve for different molecular weights (as physically expected).
• Comb Polymers: Tested on hierarchical systems with branch retraction and backbone reptation. The method accurately resolved the multi-step stress relaxation spanning ten decades, capturing both fast and slow dynamics.
• Broadband Microrheology (DWS): Applied to data spanning nine decades, including high-frequency inertia-affected regimes. The method produced a smooth, physically consistent $G(t)$ without fitting, demonstrating robustness against noise and boundary artifacts.

Observation on Long-Time Behavior: In all cases, the reconstructed $G(t)$ exhibited apparent oscillations or noise at very long times (beyond the inverse of the lowest frequency). The authors clarify that these are numerical artifacts due to finite bandwidth and do not represent physical material properties. The method provides a threshold to identify where these fluctuations begin.

5. Significance

• Bridging the Gap: i-Rheo-Tempo provides a robust, model-free solution to the frequency-to-time inverse problem, a long-standing challenge in rheology.
• Experimental Fidelity: By avoiding parametric fitting, the method ensures that the reconstructed material function reflects the actual experimental observations rather than the assumptions of a chosen model.
• Broad Applicability: The framework extends the utility of rheological analysis to ultra-broadband datasets (like microrheology) and offers a general tool for analyzing linear response functions in other fields (dielectrics, electrochemistry, optics).
• Software Availability: The authors provide open-source MATLAB and Python implementations (via DOI 10.5525/gla.researchdata.2230), including a Graphical User Interface (GUI) for easy application to experimental data.

In conclusion, i-Rheo-Tempo establishes a new standard for converting dynamic frequency data into time-domain relaxation functions, offering a mathematically rigorous, artifact-reduced, and model-independent approach that is directly compatible with experimental constraints.