A new representation formula for the logarithmic corotational derivative -- a case study in application of commutator based functional calculus

This paper utilizes a newly developed commutator-based functional calculus to derive a new representation formula for the logarithmic spin tensor in continuum mechanics while also addressing related problems concerning matrix logarithms and stress-strain monotonicity, thereby demonstrating the calculus's broad utility in tensor and matrix analysis.

The Gist

The Big Picture: Why Do We Need This?

Imagine you are watching a video of a piece of rubber being stretched, twisted, and squashed. In physics and engineering (specifically Continuum Mechanics), we need to describe how that rubber changes over time.

The problem is that the rubber is moving. If you just measure how the shape changes from one second to the next, your measurement will be messed up by the fact that the rubber is also spinning or sliding around. It's like trying to measure how much a spinning top is wobbling while you are running alongside it.

To fix this, scientists use something called an "Objective Time Derivative." Think of this as a special camera filter that ignores the spinning and sliding, focusing only on the actual stretching and squashing.

For decades, scientists have used a specific filter called the Logarithmic Corotational Derivative. It is the "gold standard" because it perfectly translates the complex math of stretching into simple, usable formulas. However, there was a catch: to use this filter, you had to solve a very difficult puzzle involving the "eigenvalues" (the hidden numbers inside the shape) of the material. It was like trying to tune a radio by taking the radio apart and soldering the wires yourself every time you wanted to listen to music. It worked, but it was messy and hard to do on a computer.

The Breakthrough: A New "Universal Remote"

The authors of this paper (Bathory, Bulíček, Málek, and Průša) have invented a new way to calculate this filter.

Instead of taking the radio apart (solving for eigenvalues), they built a Universal Remote Control. They call this tool the "Commutator Based Functional Calculus."

Here is the analogy:

• The Old Way: To change the volume, you had to manually count the gears inside the machine.
• The New Way: You just press a button labeled "Volume," and the machine does the math for you using a special set of rules.

The Secret Ingredient: The "Commutator"

The magic behind their new remote is a concept called the Commutator.

In the world of matrices (grids of numbers), order matters. If you multiply Matrix A by Matrix B, you might get a different result than if you multiply B by A.

• $A \times B \neq B \times A$

This difference is called the Commutator ($AB - BA$). It's like the "friction" or "clash" that happens when two things don't get along in the same order.

The authors realized that this "clash" is actually the key to unlocking the math. Instead of trying to find the hidden numbers (eigenvalues) inside the material, they use the "clash" between the material's shape and its movement to calculate the answer directly.

The Main Result: A Simpler Formula

The paper derives a new, much simpler formula for the Logarithmic Spin (the part of the filter that handles the spinning).

The Old Formula:
Looked like a messy recipe requiring you to list every single ingredient (eigenvalue) and mix them in a specific, complicated way.
$$\Omega_{log} = W + \text{[Huge, messy sum involving eigenvalues]}$$

The New Formula:
Looks like a clean, elegant instruction:
$$\Omega_{log} = W - \sigma(\text{Commutator}) \times D$$

In plain English: "The spin is the basic rotation ($W$) minus a correction factor ($\sigma$) that depends on how much the shape and the movement 'clash' (the Commutator)."

This new formula is powerful because:

1. It's cleaner: No need to find hidden numbers.
2. It's faster: Computers can calculate "clashes" (matrix multiplications) much faster than they can find eigenvalues.
3. It's more robust: It works even when the math gets weird or the material is stretched to its limits.

Other Cool Tricks in the Box

The authors didn't just fix the spin problem; they showed that this "Universal Remote" (the Commutator Calculus) can solve other stubborn problems in physics:

1. The "Double-Check" Problem: They proved a rule about how the "Logarithm of a Shape" changes. Previously, people weren't sure if a certain equation was always true. Using their new remote, they proved it is true only if the shape and the movement are "friendly" (they commute). If they fight (don't commute), the rule breaks. This settled a long-standing debate in the field.
2. Stress and Strain: They used the tool to prove that certain materials behave predictably (monotonically) when you stretch them. This is crucial for designing safe bridges and tires.

The "Best of Both Worlds" Conclusion

The paper ends with a philosophical note. Usually, mathematicians have to choose between two methods:

• Power Series: Easy to write down, but only works if the numbers are small (like a radio that only works when the volume is low).
• Spectral Decomposition: Works for all numbers, but is messy and hard to write down (like taking the radio apart).

The authors say: "If we restrict ourselves to symmetric shapes (which most physical materials are), we can have the best of both worlds."

They define a new way to use their "Universal Remote" that works for any size of shape, without needing the messy "taking apart" method, while keeping the clean, simple algebra of the "Power Series" method.

Summary

This paper is like upgrading from a manual transmission car to an automatic one. The physics of the engine (continuum mechanics) hasn't changed, but the way we drive it (the math) is now smoother, faster, and much easier to handle. The authors replaced a complex, manual calculation with a sleek, algebraic tool that uses the "clash" between movement and shape to do the heavy lifting.

Technical Summary

1. Problem Statement

In continuum mechanics, particularly in the theories of elasto-plastic solids, hypo-elastic solids, and viscoelastic fluids, the concept of an objective time derivative is fundamental. While various objective derivatives exist (e.g., Oldroyd's upper-convected, Jaumann–Zaremba), the logarithmic corotational derivative is unique because it satisfies the identity $\overset{\circ}{\log} H = D$, where $H$ is the Hencky strain tensor and $D$ is the symmetric part of the velocity gradient.

The core difficulty addressed in this paper is the representation of the logarithmic spin tensor ($\Omega_{\log}$), which defines this specific derivative.

• Existing Limitations: The standard representation (Xiao et al., 1997) relies on the spectral decomposition of the left Cauchy–Green tensor $B$. This involves calculating eigenvalues and eigenvectors, which is computationally cumbersome and algebraically difficult for symbolic manipulations, especially when dealing with complex constitutive relations.
• Goal: To derive a new, closed-form representation of the logarithmic spin tensor that avoids explicit spectral decomposition and eigenvalue calculations, utilizing a newly developed commutator-based functional calculus.

2. Methodology

The authors employ a commutator-based functional calculus, treating the commutator operator $\text{ad}_A[X] = AX - XA$ as a central algebraic tool.

• Commutator Operator: The paper defines the action of functions on the commutator operator via formal power series. For a function $f$, the operator $f(\text{ad}_A)$ is defined such that it acts on matrices $X$ in the space of linear operators.
• Formal Power Series: The authors initially work with formal power series expansions for matrix functions (e.g., exponential, logarithm, hyperbolic cotangent) to derive identities. They acknowledge that while these series have convergence radii, the resulting algebraic structures can be rigorously extended to symmetric matrices via spectral decomposition (as detailed in their follow-up work, Bathory, 2025).
• Key Lemmas: The methodology relies on a series of lemmas establishing:
  • Derivatives of matrix exponentials and logarithms in terms of $\text{ad}_A$.
  • Identities relating the derivative of the matrix logarithm to commutators (e.g., $\frac{\partial \ln A}{\partial A}[\text{ad}_A[X]] = \text{ad}_{\ln A}[X]$).
  • Properties of isotropic tensorial functions regarding commutativity with the derivative operator.

3. Key Contributions and Results

A. New Representation Formula for Logarithmic Spin

The primary result is Theorem 1, which provides a new formula for the logarithmic spin tensor $\Omega_{\log}$:
$$\Omega_{\log} = W - \sigma(\text{ad}_H)[D]$$
Where:

• $W$ is the skew-symmetric part of the velocity gradient.
• $H = \frac{1}{2} \ln B$ is the Hencky strain.
• $D$ is the symmetric part of the velocity gradient.
• $\sigma(x) = \coth(x) - \frac{1}{x}$ is a specific scalar function.
• $\sigma(\text{ad}_H)$ is the operator defined by the power series of $\sigma$ applied to the commutator operator $\text{ad}_H$.

Significance: This formula expresses the spin tensor directly in terms of the velocity gradient and the Hencky strain using the commutator operator, eliminating the need for explicit eigenvalue decomposition of $B$.

B. Identities for Matrix Logarithm Derivatives

The authors derive and prove several new identities regarding the derivative of the matrix logarithm, refining existing literature (e.g., Neff et al., 2024):

• Commutativity Condition: They prove that $\frac{\partial \ln A}{\partial A}[AX + XA] = 2X$ if and only if $X$ commutes with $A$ (for symmetric positive definite $A$).
• Generalized Derivative Formulas: They provide closed-form expressions for derivatives of the form $\frac{\partial \ln A}{\partial A}[A^p X A^{-s} \pm A^{-s} X A^p]$ using hyperbolic functions of the commutator operator (e.g., $\sinh$, $\cosh$).

C. Monotonicity of Stress-Strain Relations

The paper addresses the corotational stability postulate. It proves the equivalence between two characterizations of monotonicity for isotropic tensorial functions:

1. A condition involving the derivative of a function composed with the logarithm of a tensor ($\ln A$).
2. A condition involving the derivative of the function with respect to the tensor itself ($G$).

Using the commutator calculus, the authors show that these two conditions are equivalent via a bijective transformation involving the operator $\sqrt{r_{p+s}}(\text{ad}_G)$. This simplifies the analysis of stability in constitutive models.

4. Significance and Implications

• Algebraic Efficiency: The commutator-based calculus allows for seamless symbolic manipulation of tensor-valued functions and their derivatives without resorting to the "messy" spectral decomposition (eigenvalues/eigenvectors) required by traditional methods.
• Generality: The framework demonstrates that the commutator operator is a universal tool in tensor/matrix analysis, applicable to a wide range of problems beyond just the logarithmic spin, including matrix logarithms and stress-strain monotonicity.
• Bridging Theory and Application: While the initial proofs rely on formal power series (which have convergence limits), the authors argue that for symmetric tensors (common in mechanics), these operators can be rigorously defined via spectral decomposition (Definition 1), preserving the algebraic simplicity of the commutator approach while extending validity beyond the radius of convergence.
• Foundation for Future Work: This paper serves as a "case study" for a broader functional calculus developed by the authors (Bathory, 2025), suggesting a new paradigm for handling objective rates and constitutive modeling in continuum mechanics.

In summary, the paper successfully replaces a complex, eigenvalue-dependent formula for the logarithmic spin with a compact, operator-based formula, demonstrating the power of commutator calculus in simplifying and generalizing problems in continuum mechanics.