Lie symmetry classification and invariant solutions of time-fractional telegraph systems with variable coefficients

This paper presents a complete Lie symmetry classification of time-fractional telegraph systems with variable coefficients, identifying three distinct symmetry classes and deriving exact invariant solutions in terms of Mittag-Leffler, generalized Wright, and Fox $H$-functions to model transport phenomena with memory and nonlocal effects.

The Gist

Imagine you are trying to predict how a signal travels through a very strange, bumpy road. In the real world, signals (like heat moving through a material or electricity moving through a chip) don't just zip along instantly. They have "memory." If the road was bumpy yesterday, the signal might still be wobbling today because of that past experience. They also don't just move in a straight line; they spread out like a wave and diffuse like a drop of ink in water.

Mathematicians use a special tool called a Telegraph Equation to describe this kind of movement. But when the material is complex (like a semiconductor with uneven properties) and the "memory" effect is strong, the standard math isn't enough. That's where this paper comes in.

Here is a simple breakdown of what the authors did, using some everyday analogies:

1. The Problem: A Road with Changing Rules

The authors are studying a specific type of equation (a "time-fractional telegraph system") that models these signals.

• The "Road" (Coefficients): Imagine the road isn't flat. Some parts are slippery, some are sticky, and the rules change depending on where you are (spatially varying coefficients).
• The "Memory" (Fractional Derivative): Unlike a normal car that only cares about the road right under its tires, this "signal car" remembers the road it drove on for the last hour. The math uses something called the Riemann–Liouville fractional derivative to keep track of this history.

2. The Tool: The "Symmetry" Detective

To solve these messy equations, the authors used a method called Lie Symmetry Analysis.

• The Analogy: Imagine you have a complex, tangled knot of string. You want to untangle it to see the pattern. You look for "symmetries"—ways you can rotate, stretch, or shift the knot without changing its fundamental shape.
• What they did: They acted like detectives, looking for these hidden symmetries in their equations. They asked: "If I change the time or the position in a specific way, does the equation still look the same?"
• The Discovery: They found that the answer depends entirely on the relationship between two things: the transport coefficient (how fast the signal moves, like the road's smoothness) and the potential function (external forces pushing the signal).

3. The Three "Families" of Solutions

The authors discovered that depending on how the road and the forces relate to each other, the equations fall into three distinct families (or classes).

• Family 1: The most general case. The road and forces are related in a specific, complex way.
• Family 2: A slightly simpler relationship where the forces are tied to the road's shape in a specific formula.
• Family 3: The most special case where the forces are perfectly balanced with the road's shape.

For each family, they built an "Optimal System."

• The Analogy: Think of this as a master keyring. Instead of trying every single key to open a door, they found the smallest, most efficient set of keys (symmetries) that can open every possible door in that family.

4. The Result: Cracking the Code

Once they found the right keys (symmetries), they could simplify the complex equations.

• The Reduction: They turned a difficult problem with two variables (time and space) into a simpler problem with just one variable (a "fractional ordinary differential equation").
• The Solution: They solved these simpler problems and wrote down the exact answers. These answers aren't simple numbers; they are expressed using special mathematical "super-functions" named after famous mathematicians:
  • Mittag-Leffler functions: The "fractional cousins" of the standard exponential functions we use in basic physics.
  • Generalized Wright functions and Fox H-functions: Even more complex tools needed to describe the "memory" and "non-local" behavior of the system.

Why Does This Matter?

The paper claims these solutions are benchmarks.

• The Analogy: Imagine engineers are building a new computer simulation to design better car brakes or faster microchips. They need a "gold standard" answer to check if their computer is working correctly.
• Because the authors found exact, closed-form solutions (the "gold standard"), engineers can run their complex computer models and compare them to these exact answers. If the computer model matches the paper's solution, the engineers know their model is accurate.

Summary

In short, this paper is a mathematical map. It tells us exactly how to navigate a specific type of complex, memory-filled signal transport problem. By finding the hidden symmetries, the authors turned a messy, unsolvable-looking puzzle into a set of clear, exact formulas. These formulas act as a "truth check" for scientists and engineers trying to model real-world systems like heat flow in special materials or electricity in uneven semiconductors.

Note: The paper focuses strictly on the mathematical classification and the derivation of these exact formulas. It does not claim to have solved a specific industrial problem yet, nor does it discuss clinical uses; it provides the mathematical tools (the exact solutions) that others can use to validate their own models.

Technical Summary

1. Problem Statement

The paper addresses the mathematical modeling of transport processes in industrial and physical systems that exhibit memory effects and nonlocal behaviors. While classical telegraph equations model signal propagation with finite speed, they fail to capture the hereditary properties found in viscoelastic media, anomalous heat conduction, and charge transport in disordered semiconductors.

The authors investigate a specific class of time-fractional telegraph systems with spatially varying coefficients, defined by the system:
$$\begin{cases} \frac{\partial^\alpha u}{\partial t^\alpha} = v_x \\ \frac{\partial^\alpha v}{\partial t^\alpha} = f(x)u_x + g(x)u \end{cases}$$
where:

• $\alpha > 0$ is the fractional order (using the Riemann–Liouville derivative).
• $f(x) > 0$ represents spatially dependent transport properties (e.g., diffusivity).
• $g(x)$ represents external potential or forcing effects in inhomogeneous media.

The Core Challenge: Previous studies on fractional telegraph equations were limited to specific coefficient forms or restricted conditions. A complete Lie symmetry classification for this system with general differentiable functions $f(x)$ and $g(x)$ was an open problem. The goal was to identify all functional forms of $f(x)$ and $g(x)$ that admit symmetry extensions and to derive exact invariant solutions.

2. Methodology

The authors employed Lie symmetry analysis extended to fractional partial differential equations (FPDEs). The methodology involved the following steps:

1. Infinitesimal Generator Construction: They defined the general infinitesimal generator $X = \xi \partial_x + \tau \partial_t + \mu \partial_u + \phi \partial_v$ and its prolongation, including the specific formulas for fractional extended infinitesimals $\mu^{(\alpha)}$ and $\phi^{(\alpha)}$ involving the Riemann–Liouville derivative.
2. Determining Equations: By applying the invariance criterion $\tilde{X}(\Delta)|_{\Delta=0} = 0$, they derived an overdetermined system of determining equations.
3. Linearity Proof: A crucial lemma proved that the infinitesimals $\mu$ and $\phi$ must be linear with respect to the dependent variables $u$ and $v$. This simplified the system by eliminating higher-order supplementary terms ($\mu_1, \phi_1$).
4. Classification: By analyzing the constraints on $f(x)$ and $g(x)$ derived from the determining equations, the authors classified the system into distinct symmetry cases.
5. Optimal Systems: For each symmetry class, they constructed optimal systems of one-dimensional Lie subalgebras to ensure a non-redundant classification of group-invariant solutions.
6. Symmetry Reduction: Using the characteristic equations of the optimal systems, the original FPDEs were reduced to fractional ordinary differential equations (FODEs).
7. Exact Solutions: The reduced FODEs were solved analytically, expressing solutions in terms of special functions: Mittag-Leffler functions, generalized Wright functions, and Fox H-functions.

3. Key Contributions

• Complete Lie Symmetry Classification: The paper provides the first complete classification of Lie symmetries for time-fractional telegraph systems with general spatially varying coefficients.
• Identification of Symmetry Classes: The study reveals that the symmetry structure depends fundamentally on the relationship between $f(x)$ and $g(x)$. The authors identified three distinct symmetry classes (plus a generic case) based on specific functional forms of $g(x)$ relative to $f(x)$:
  1. Generic Case: Only trivial scaling and linear superposition symmetries exist.
  2. Case 1: $g(x) = \frac{\lambda_2 \sqrt{f(x)}}{\omega_{\lambda_1}(x)} + \frac{f'(x)}{2}$, admitting an additional scaling symmetry.
  3. Case 2: $g(x) = \lambda_2 \sqrt{f(x)} + \frac{f'(x)}{2}$, admitting a translation-like symmetry.
  4. Case 3: $g(x) = \frac{f'(x)}{2}$, admitting the largest symmetry algebra (including a rotation-like symmetry in the $u-v$ plane).
• Reduction to FODEs: The authors successfully reduced the complex FPDE system to solvable FODEs for each symmetry class.
• Closed-Form Solutions: They derived explicit invariant solutions expressed via advanced special functions (Fox H-functions, Wright functions), which are essential for describing fractional transport phenomena.

4. Key Results

The study yielded specific invariant solutions for the identified cases:

• Case 1 & 2 (Reduced to FODEs with variable coefficients):
  • Solutions were expressed using Fox H-functions ($H_{p,q}^{m,n}$) for $0 < \alpha < 1$ and generalized Wright functions ($\Psi$) for $\alpha \geq 1$.
  • The similarity variables involved combinations of $t$ and integrals of $1/\sqrt{f(x)}$.
• Case 3 ($g(x) = f'(x)/2$):
  • This case allowed for the most significant reduction. The system could be decoupled into separable equations.
  • Solutions were expressed as sums of products of spatial functions and temporal functions involving Mittag-Leffler functions ($E_{\alpha, \beta}$) and Wright functions.
  • Notably, if $f(x) = x^{2m}$, the system reduces to a separable form for any $m > 0$, a generalization of previous results which were limited to $f(x)=x^2$.

5. Significance and Implications

• Theoretical Advancement: This work bridges the gap between classical telegraph equation symmetry analysis and fractional calculus, extending the theory to systems with arbitrary spatial heterogeneity.
• Benchmarking: The derived exact analytical solutions serve as critical benchmarks for validating numerical methods used in industrial transport modeling. Numerical schemes for fractional PDEs often struggle with accuracy; these exact solutions provide a ground truth for verification.
• Physical Insight: The solutions, characterized by Mittag-Leffler and Fox H-functions, mathematically capture the "memory" and "nonlocal" nature of transport in complex media (e.g., viscoelasticity, anomalous diffusion).
• Future Directions: The authors suggest that this framework can be extended to nonlinear fractional telegraph equations, multidimensional systems, and other fractional operators (e.g., Caputo derivative), further strengthening the mathematical foundation for modeling complex industrial systems.

In summary, this paper provides a rigorous mathematical framework for analyzing and solving time-fractional telegraph systems with variable coefficients, offering both a complete symmetry classification and a rich set of exact solutions that are vital for understanding and modeling memory-dependent transport phenomena.