Next Article in Journal
Real-Time Identification of Strawberry Pests and Diseases Using an Improved YOLOv8 Algorithm
Next Article in Special Issue
A Uniqueness Theorem for Stability Problems of Functional Equations
Previous Article in Journal
Current Profiling Control for Torque Ripple Reduction in the Generating Mode of Operation of a Switched Reluctance Motor Drive
Previous Article in Special Issue
Some Properties on Normalized Tails of Maclaurin Power Series Expansion of Exponential Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order Correlation between Special Functions Inspired by Bone Fractal Operators

1
Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
2
Department of Civil Engineering, Chengdu Technological University, Chengdu 611730, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(10), 1279; https://doi.org/10.3390/sym16101279
Submission received: 29 August 2024 / Revised: 25 September 2024 / Accepted: 27 September 2024 / Published: 29 September 2024

Abstract

:
In recent years, our research on biomechanical and biophysical problems has involved a series of symmetry issues. We found that the fundamental laws of the aforementioned problems can all be characterized by fractal operators, and each type of operator possesses rich invariant properties. Based on the invariant properties of fractal operators, we discovered that the symmetry evolution laws of functional fractal trees in the physical fractal space can reveal the intrinsic correlations between special functions. This article explores the fractional-order correlation between special functions inspired by bone fractal operators. Specifically, the following contents are included: (1) showing the intrinsic expression in the convolutional kernel function of bone fractal operators and its correlation with special functions; (2) proving the following proposition: the convolutional kernel function of bone fractal operators is still related to the special functions under different input signals (external load, external stimulus); (3) using the bone fractal operators as the background and error function as the core, deriving the fractional-order correlation between different special functions.

1. Introduction

This article will use bone fractal operators as the backdrop, the error function as the link, and combine the symmetry evolution laws of functional fractal trees in the physical fractal space to reveal the fractional-order correlation between different special functions.
In recent years, our research on the self-similar structures of biomaterials and fractional-order mechanics has involved a series of symmetry and symmetry-breaking issues. For instance, when establishing the fractional-order constitutive model of tendon/ligament fibers, we found that the symmetry breaking of fractal cells introduces a time modulation factor in their relaxation response, further extending the relaxation time and successfully revealing the inevitable correlation between the fractional-order viscoelastic emergent behavior of biological fibers and their fractional-dimension self-similar structures [1]. When establishing the fractional-order impulse characteristics of neurons, we discovered that the symmetry breaking of step cells introduces an exponential modulation term in the dynamic operator, which also provides a functional explanation for why nerve fibers have spiny dendrites [2]. When studying the infinite-level physical fractal circuit model and fractional-order time-varying dynamic response of arterial blood flow, we found that the symmetry breaking of fractal cells is a necessary condition for the formation of the reflux time-varying response phenomena in real arterial blood flow [3].
In studying the aforementioned issues, we found that as the symmetry of the functional fractal tree evolves in the physical fractal space, the properties of the fractal operators become more enriched, as shown in Figure 1. Specifically, Figure 1a represents the abstract process of fractal cells in nerve fibers, forming a resistance-capacitance network with an infinitely expanding tree-like topology that is structurally perfectly symmetrical [2]. Figure 1b represents the abstract process of fractal cells in neuronal dendrites. Its assembly of components is consistent with that of nerve fibers, and it maintains a macroscopic binary bifurcation structure [2]. However, at this point, due to the introduction of an additional physical component, the functional fractal tree undergoes its first symmetry breaking. This phenomenon is caused by the unequal number of physical components at the two ends of the bifurcation. Physiologically, this multi-scale fractal tree model of dendrites can also be used to explain the electrophysiological changes associated with aging and neurodegenerative diseases [4,5]. If the topological structure of the functional fractal tree and the number of the physical components at the two ends of the bifurcation are changed simultaneously, a second symmetry breaking will occur, forming a new fractal step topology. Figure 1c represents the abstract process of fractal cells in compact bone fibers [6]; this soft/hard alternating multi-level self-assembly pattern is not unique to bone and is widely found in biological materials such as shell [7] and nacre [8]. In our study of the operator algebra method for viscoelastic theory, we found that fractal tree and fractal step operators share commonalities: both are quadratic radical operators, or in other words, both are non-rational operators. There are also differences: the fractal tree operator is extremely simple, while the fractal step operator is quite complex. The simplicity of the former comes from the symmetry of the fractal tree topology, while the complexity of the latter arises from the symmetry breaking of the fractal step topology. It is important to note that the methods for breaking symmetry are not unique. Figure 1b achieves symmetry breaking by changing the number of physical components at the two ends of the functional fractal tree. If the number of physical components is not changed, but only their functional properties are altered, symmetry breaking can also be achieved. Therefore, if the topological structure of the functional fractal tree and the functional properties of the physical components at the two ends of the bifurcation are changed simultaneously, the fractal step topology shown in Figure 1c can also be formed. We found that as the number and degree of symmetry breaking increase, the properties of the fractal operators become more enriched. For a perfectly symmetrical functional fractal tree (Figure 1a), it is difficult to characterize the various mechanical properties of complex organisms, such as the high modulus and toughness of bone [9,10,11,12,13], with just two types of physical components and the classical spring-dashpot model. The introduction of symmetry breaking increases the designability of fractal operators, enriches their basic functions, and broadens the applicability of the physical fractal space.
In studying the biomechanical and biophysical issues, we abstracted the “physical fractal space” and established the mechanics in the “physical fractal space [6]”, which is different from the traditional “geometric fractal model [14,15,16]”. Surprisingly, the mechanics in the physical fractal space exhibit common features: non-rational fractional-order fractal operators, mechanics governed by non-elementary functions, non-integer-order mechanics, and non-localized mechanics. These fractional-order fractal operators all possess rich invariance properties. Based on the invariant properties of fractal operators, we confirmed that the special functions from different disciplines and research backgrounds indeed have fractional-order correlations. With the help of fractal operators, the deep dependency relationships between special functions can be revealed. The fractal operators in the physical fractal space provide a new perspective for understanding non-elementary functions and fractional-order calculus [6].
In studying the convolutional kernel function of bone fractal operators, we confirmed that the error function is its core component [6]. We found that the fundamental laws in the physical fractal space can all be characterized by fractal operators and the convolution kernel functions of each type of the operators are related to the special functions. Despite the completely different objects of study, an astonishing consistency is displayed. How can we understand this consistency? Is there a correlation between these special functions? This study attempts to provide answers, including the following contents: (a) an analysis of the intrinsic expressions and characteristics in the convolutional kernel function of bone fractal operators; (b) a demonstration of the correlation between the convolutional kernel function of bone fractal operators and special functions under various input signals (external load, external stimulus); and (c) a revelation of the fractional-order correlation between different special functions based on bone fractal operators. The definition and operations of traditional fractional-order calculus often involve certain special functions, such as the Gamma function. In other words, the emergence of special functions has provided the theoretical foundation for the establishment and development of fractional-order calculus [17,18,19]. Now, inspired by the symmetry evolution laws of functional fractal trees in the physical fractal space, with the aid of bone fractal operators, we established a series of correlations between special functions, which represents a novel attempt to understand the commonalities across different disciplinary fields.
Section 2 discusses the intrinsic expression in convolutional kernel function of bone fractal operators. Section 3, Section 4 and Section 5 discuss the properties of the error function. Section 6, Section 7, Section 8, Section 9 and Section 10 discuss the correlation between error function and various special functions.

2. The Intrinsic Expression in Convolutional Kernel Function of Bone Fractal Operators

The previous article abstracted the physical fractal space from the compact bone fibers [6]. This section mainly discusses the intrinsic expression in the convolutional kernel function of bone fractal operators. For convenience, we first provide the algebraic expressions of bone fractal operators. When stress σ t is applied to bone fibers and strain ε t is generated, then we have
σ t = T p ε t ,
where T p is the fractal operator on the physical fractal space of the bone. The differential operator p is defined as follows:
p g t = d g d t .
Unlike the conventional approach, our research object is selected as a microfiber. Force and deformation are transmitted along the fiber chain, from which we can abstract a multi-level self-similar component tree. If the number of levels in the multi-level structure tends toward infinity, the self-similar component tree becomes a fractal component tree, as shown in Figure 1c. The “=” sign means that the intrinsic invariance in the physical fractal space is the equivalence of the fractal cell and fractal component [6].
Therefore, the bone fractal operator T needs to satisfy the following algebraic equation:
1 1 T 1 + 1 T 2 + 1 T + T 3 = T .
In mechanics, the right-hand side of Equation (3) can be regarded as the stiffness of the fractal component, and the left-hand side can be regarded as the stiffness of the fractal cell. It is further deduced that
T 1 + T 2 T 2 T 1 + T 2 T 3 T T 1 T 2 T 3 = 0 .
Equation (4) is the algebraic equation for the fractal operator T . This is a quadratic equation with a radical solution:
T = 1 2 T 3 ± T 3 2 + 4 T 1 T 2 T 3 T 1 + T 2 .
So far, the fractal operator T of compact bone is determined, as shown in Equation (5).
If the physical component operators are taken as
T 1 = G 1 , T 2 = G 2 , T 3 = G 3 μ 3 p ,
where both T 1 and T 2 are elastic elements and T 3 is a viscous element, then the main body of the fractal operator in Equation (5) becomes T ^ p = p ± p 2 + α p , where α is the physical constant, equivalent to the characteristic frequency [6].
Interestingly, T ^ p operators appeared in both hemodynamics and bone mechanics [6]. In hemodynamic research, Peng et al. [3] considered the operator T ^ p as a whole. Unlike Peng et al., in bone mechanics, Jian et al. separated the quadratic radical term p 2 + α p in the operator T ^ p and studied it separately [6]:
T p = p 2 + α p = p p + α .
The quadratic radical operator T p is decomposed into the product of two 1/2-order differential operators p and p + α . In Courant’s work [20], the operator p + α acting on the Heaviside unit step function η t generates a convolutional kernel function h t :
h t = p + α η ( t ) .
With the help of the displacement theorem [21], the expression for the kernel function h t can be written as follows:
h t = p + α η ( t ) = e α t p e α t η t = e α t p e α t = e α t π d d t 0 t e α τ t τ d τ
In Equation (7), e α t t is the characteristic function, 0 t e α τ t τ d τ is the characteristic integral of the convolutional kernel function of bone fractal operators, where e α τ is the exponential term brought by the displacement theorem. 1 t τ comes from the convolutional kernel function of the operator p .
If the input signal is f t and the fractal operator p + α acts on f t , then the output signal can be expressed as
p + α f ( t ) = e α t p e α t f t = e α t π d d t 0 t e α τ t τ f τ d τ .
Obviously, the operation of bone fractal operator p + α inevitably involves the integrals such as 0 t e α τ t τ d τ and 0 t e α τ t τ f τ d τ :
p e α t f t = 1 π d d t 0 t e α τ t τ f τ d τ ,
we call it the weighted integration of the characteristic function e α τ t τ . Obviously, the weighted integration is the core component of the 1/2-order fractal operator p + α . Moreover, the weighted integration is closely related to the error function. The previous article [6] has exported
0 t e α τ t τ η τ d τ = e α t 0 t e α ( τ t ) t τ d τ = e α t α t 0 e x 2 x α 2 x α d x = π e α t α erf α t ,
namely,
0 t e α τ t τ d τ = 0 t τ 0 e α τ t τ d τ = π e α t α erf α t .
In Equation (9), integration 0 t τ 0 e α τ t τ d τ can be considered as the 0 moment of the characteristic function e α τ t τ . Therefore, the kernel function can be written as
h t = p + α η ( t ) = e α t π d d t π e α t α erf α t = e α t α d d t e α t erf α t ,
and
p e α t = 1 π d d t 0 t e α τ t τ d τ = 1 α d d t e α t erf α t .
Combining the differential properties of error function,
p erf α t = α e α t π t ,
the final expression of the kernel function can be exported
h t = p + α η ( t ) = α erf α t + e α t π t ,
and
p e α t = α e α t erf α t + 1 π t = e α t h t = e α t p + α η ( t ) .
Equation (12) shows that the core component of characteristic integration 0 t e α τ t τ d τ is the error function erf α t . It can also be said that the core component of the 1/2-order fractal operator p + α is the error function erf α t . Therefore, the following analysis primarily focuses on the invariance properties of weighted integration and the error function.

3. The Moment of the Characteristic Function and Its Correlation with Error Function

This section considers the typical cases of the general input signal f τ . Common input signals f τ are typically elementary functions, such as the unit slope function τ , the unit acceleration function τ 2 , the sine function sin ( τ ) , the cosine function cos ( τ ) , the exponential function e β τ , and the logarithmic function ln τ , among others.
This section attempts to answer the following basic question: when the input signal f ( τ ) takes different types of elementary functions, is 0 t e α τ t τ f τ d τ still related to the error function? If the answer is yes, can this correlation reveal the correlation between the error function and different special functions? It should be emphasized that there are theoretically an infinite variety of elementary functions and their combinations. This study, however, will only discuss the following typical elementary functions.
When f ( τ ) is taken as a polynomial function, special cases can be considered first. The unit function f ( τ ) = τ 0 , corresponding to the 0 moment of the characteristic function, which has been discussed in the previous section. Here, when the input signal is the unit slope function f ( τ ) = τ ,
0 t τ e α τ t τ d τ = t α + 1 2 α 3 2 π e α t 1 + 2 α t erf α t .
In Equation (13), integration 0 t τ e α τ t τ d τ can be regarded as the first moment of the characteristic function e α τ t τ .
When the input signal takes the unit acceleration function f ( τ ) = τ 2 ,
0 t τ 2 e α τ t τ d τ = 1 4 α 5 2 2 α t ( 3 + 2 α t ) + e α t π ( 3 + 4 α t ( 1 + α t ) ) erf ( α t ) .
In Equation (14), integration 0 t τ 2 e α τ t τ d τ can be regarded as the second moment of the characteristic function e α τ t τ .
From Equations (13) and (14), the integration of the weighted characteristic function is related to the error function. To export more general expressions, let
0 t τ 0 e α τ t τ d τ = I 0 ,
take the derivative of the parameter variable α at both sides of Equation (15):
I 0 α = 0 t τ e α τ t τ d τ = I 1 .
Let
0 t τ n e α τ t τ d τ = I n ,
I n is the n-th moment of the characteristic function. It can be inferred that
I n = I n 1 α = n I 0 α n .
A very simple recursive relationship and general formula can be written using Equation (18):
I n + 1 = I n α , I n = n I 0 α n .
According to Equations (9) and (15), the 0 moment I 0 can also be written as
I 0 = π e α t α erf α t .
It can be inferred that the n-th moment I n is related to the error function. Therefore, we give the following proposition:
Proposition 1.
When the input signal f ( τ ) takes any n-degree polynomial function, 0 t e α τ t τ f ( τ ) d τ is related to the error function.

4. Weighted Integration and Its Correlation with Error Function When the Input Signal Is Exponential and Trigonometric Functions

Exponential functions are closely related to trigonometric functions. When the input signal f ( τ ) takes an exponential function e β τ , it is equivalent to taking the parameter variable α in the above equations as α + β :
0 t e β τ e α τ t τ d τ = π α + β e α + β t erf ( α + β t ) .
Although Equation (20) is trivial, it provides an opportunity to introduce sine and cosine functions. In Equation (20), let
β = ± i λ ,
we have
0 t e i λ τ e α τ t τ d τ = π α + i λ e α + i λ t erf ( α + i λ t ) ,
0 t e i λ τ e α τ t τ d τ = π α i λ e α i λ t erf ( α i λ t ) .
According to
cos ( λ τ ) = 1 2 e i λ τ + e i λ τ , sin ( λ τ ) = 1 2 i e i λ τ e i λ τ .
By combining Equations (21) and (22), we can derive
0 t cos λ τ e α τ t τ d τ = π 2 e α + i λ t erf α + i λ t α + i λ + e α i λ t erf α i λ t α i λ ,
0 t sin λ τ e α τ t τ d τ = π 2 i e α + i λ t erf α + i λ t α + i λ e α i λ t erf α i λ t α i λ .
Equations (23) and (24) have an integral on the real field at the left end and a complex number at the right end. But this is not contradictory. After operation on the right end, the complex numbers will disappear. It is important to note that the core of the right-hand term is the imaginary error function.
At this point, we can extract the following general proposition.
Proposition 2.
When the input signal f ( τ ) takes the exponential, sine, and cosine functions, 0 t e α τ t τ f ( τ ) d τ is related to the error function.

5. The Integral Property of Error Function

This section will confirm the following proposition.
Proposition 3.
When the input signal f ( τ ) takes different elementary functions, the integration operation results with the error function still possess correlation with the error function.
In the case of f ( τ ) taking the Heaviside unit function η t , using the partial integration method and combining it with the differential properties of the error function (Equation (11)), we can derive
0 t erf ( α τ ) η τ d τ = 0 t τ 0 erf ( α τ ) d τ = 0 t erf ( α τ ) d τ = 1 π α t e α t + 1 2 α + t erf ( α t ) .
In Equation (25), 0 t τ 0 erf ( α τ ) d τ is the 0 moment of error function. The 0 moment of the error function is still related to the error function.
The input signal takes the unit slope function—that is, f ( τ ) = τ . By using the partial integration method and combining with the differential properties of the error function (Equation (11)), we can derive
0 t τ erf ( α τ ) d τ = 1 8 α 2 2 α π ( 3 + 2 α t ) t e α t + ( 3 + 4 α 2 t 2 ) erf ( α t ) .
In Equation (26), 0 t τ erf ( α τ ) d τ is the first moment of error function. The first moment of the error function is still related to the error function.
The input signal takes the unit acceleration function—that is, f ( τ ) = τ 2 . By using the partial integration method and combining the differential properties of the error function (Equation (11)), we can derive
0 t τ 2 erf ( α τ ) d τ = 1 24 α 3 2 α π ( 15 + 10 α t + 4 α 2 t 2 ) t e α t + ( 15 + 8 α 3 t 3 ) erf ( α t ) .
In Equation (27), 0 t τ 2 erf ( α τ ) d τ is the second moment of the error function. The second moment of the error function is still related to the error function.
The 0, 1st, and 2nd moments of the error function are all related to the error function. We speculate that this phenomenon should be universal. There should be the following proposition: the n-th moment of the error function is still related to the error function.
Now we confirm the above proposition holds. We directly analyze the general situation. By using the partial integration method, we can obtain the n-th moment of the error function:
0 t τ n erf α τ d τ = 1 n + 1 0 t erf α τ d τ n + 1 = t n + 1 erf α t n + 1 1 n + 1 0 t τ n + 1 d erf α τ = t n + 1 erf α t n + 1 α π n + 1 0 t τ n + 1 / 2 e α τ d τ = t n + 1 erf α t n + 1 1 π α n + 1 n + 1 0 α t τ n + 1 / 2 e τ d τ .
The last term in Equation (28) involves the (n + 1)-order moment of the negative exponential function. The definition of the lower incomplete Gamma function is as follows [22]:
γ s , x = 0 x t s 1 e t d t ,
the (n + 1)-order moment of the negative exponential function can be written as
0 α t τ n + 1 / 2 e τ d τ = γ n + 3 / 2 , α t .
The general formula for the n-th moment of error function can be derived from Equation (29):
0 t τ n erf α τ d τ = t n + 1 erf α t n + 1 1 π α n + 1 ( n + 1 ) γ n + 3 / 2 , α t .
At this point, we confirm the following general proposition: the integration of the polynomial weighted error function is not only related to the error function but also to the Gamma function.
Let us take a different perspective and look at the general formula for the n-th moment. We focus on the recursive properties of the n-th moment. Equation (30) gives the moment of (n − 1)-order as follows:
0 t τ n 1 erf α τ d τ = t n erf α t n 1 π α n n γ n + 1 / 2 , α t ,
namely,
t n erf α t = n 0 t τ n 1 erf α τ d τ + 1 π α n γ n + 1 / 2 , α t .
Based on Equation (32), recalculating the n-th moment, we can derive the recursive formula:
0 t τ n erf α τ d τ = 0 t n 0 t 1 τ n 1 erf α τ d τ + γ n + 1 / 2 , α t 1 π α n d t 1 = 0 t n 0 t 1 τ n 1 erf α τ d τ d t 1 + 0 t γ n + 1 / 2 , α t 1 π α n d t 1 .
The recursive Formula (33) shows that the n-th moment of the error function can be expressed as an integral of the (n − 1)-order moment. Similarly, the moment of (n − 1)-order can be expressed as an integral of the (n − 2)-order moment. By repeating the above operations, we can conclude that
0 t τ n erf α τ d τ = n ! 0 t 0 t 1 0 t n 1 0 t n erf α τ d τ d t n d t 2 d t 1 + 0 t γ n + 1 / 2 , α t 1 π α n d t 1 + 0 t 0 t 1 0 t n 2 0 t n 1 γ 3 / 2 , α t n π α d t n d t n 1 d t 2 d t 1 = n ! p n 1 erf α t + k = 1 n 1 π α n + 1 k p k γ n + 3 / 2 k , α t .
By combining Equations (30) and (34), we can conclude that
0 t τ n erf α τ d τ = t n + 1 erf α t n + 1 γ n + 3 / 2 , α t π α n + 1 ( n + 1 ) = n ! p n 1 erf α t + k = 1 n p k γ n + 3 / 2 k , α t π α n + 1 k .
The last equation of Equation (35) is summarized as follows:
t n + 1 erf α t n + 1 ! γ n + 3 / 2 , α t π α n + 1 ( n + 1 ) ! = p n 1 erf α t + 1 n ! k = 1 n p k γ n + 3 / 2 k , α t π α n + 1 k ,
namely,
p n + 1 erf α t = t n + 1 n + 1 ! erf α t 1 π ( n + 1 ) ! k = 0 n p k γ n + 3 / 2 k , α t α n + 1 k .
Equation (36) indicates that the high-order integration of the error function is not only related to the error function but also to the high-order integration of the Gamma function. We can write the error function and Gamma function in Equation (36) on both sides of the equal sign:
t n + 1 n + 1 ! erf α t p n + 1 erf α t = 1 π ( n + 1 ) ! k = 0 n p k γ n + 3 / 2 k , α t α n + 1 k .
In Equation (37), the left end only has an error function, and the right end only has a Gamma function. Two unrelated functions can be given identity relationships through complex integration operations.
We can also understand Equation (37) from the perspective of integral equations. Consider the following integral equation:
t n + 1 n + 1 ! f t p n + 1 f t = 1 π ( n + 1 ) ! k = 0 n p k γ n + 3 / 2 k , α t α n + 1 k .
It is obvious that this integral equation has the following special solutions:
f t = erf α t .
That is to say, the error function is the solution of the nonlinear integral Equation (38).
It is important to note that the first term n ! p n 1 erf α t of the last equal sign in Equation (35) is closely related to the weighted feature integration formula 0 t τ n e α τ t τ d τ . According to the differential properties of the error function (Equation (11)), the weighted feature integration can be expressed as
0 t τ n e α τ t τ d τ = e α t 0 t τ n e α t τ t τ d τ = e α t t n e α t t = π α e α t t n p erf α t .
The convolution in the last term of Equation (40) can be calculated as follows (the proof can be found in Appendix A):
t n p erf α t = n ! p n + 1 p erf α t = n ! p n erf α t .
Therefore, the weighted feature integration can be transformed into
0 t τ n e α τ t τ d τ = n ! π α e α t p n erf α t .
Equation (42) indicates that when f ( τ ) taking τ n , 0 t e α τ t τ f ( τ ) d τ is still related to the n-th integration of the error function.
For ease of comparison, Equation (42) can be rewritten as
p n + 1 erf α t = α π e α t n + 1 ! 0 t τ n + 1 e α τ t τ d τ .
Comparing Equations (36) and (43),
α π e α t n + 1 ! 0 t τ n + 1 e α τ t τ d τ = t n + 1 n + 1 ! erf α t 1 π ( n + 1 ) ! k = 0 n p k γ n + 3 / 2 k , α t α n + 1 k .
Simplifying the above equation yields
0 t τ n + 1 e α τ t τ d τ = π α e α t t n + 1 erf α t 1 π k = 0 n p k γ n + 3 / 2 k , α t α n + 1 k .
Equation (44) indicates that when f ( τ ) taking τ n , 0 t e α τ t τ f ( τ ) d τ is not only related to error function but also to the n-th integration of Gamma function. Furthermore, according to Equation (44), we have
0 t τ n e α τ t τ d τ = π α e α t t n erf α t 1 π k = 0 n 1 p k γ n + 1 / 2 k , α t α n k ,
0 t τ n 1 e α τ t τ d τ = π α e α t t n 1 erf α t 1 π k = 0 n 2 p k γ n 1 / 2 k , α t α n 1 k .
By combining Equations (45) and (46) and eliminating the error function, we can obtain the recursive relationship:
0 t τ n e α τ t τ d τ t 0 t τ n 1 e α τ t τ d τ = π α e α t t π k = 0 n 2 p k γ n 1 / 2 k , α t α n 1 k 1 π k = 0 n 1 p k γ n + 1 / 2 k , α t α n k .
When the input signal f ( τ ) takes the sine sin ( τ ) and cosine cos ( τ ) functions, respectively, we can derive that
0 t sin ( τ ) erf ( α τ ) d τ = cos ( t ) erf ( α t ) + α 2 erf ( i + α ) t ( i + α ) + erf ( i + α ) t ( i + α ) ,
0 t cos ( τ ) erf ( α τ ) d τ = sin ( t ) erf ( α t ) + i α 2 erf ( i + α ) t ( i + α ) erf ( i + α ) t ( i + α ) .
When the input signal f ( τ ) takes the exponential function e β τ , we can derive that
0 t e β τ erf ( α τ ) d τ = 1 β e β t erf ( α t ) α erf ( i ( β α ) t ) β α .
We can introduce the 0 moment:
0 t τ 0 e β τ erf ( α τ ) d τ = 0 t e β τ erf ( α τ ) d τ = M 0 .
Taking the derivative of the parameter variable β at both ends yields, we can obtain the first moment:
M 0 β = 0 t τ e β τ erf ( α τ ) d τ = M 1 .
Taking the derivative of the parameter variable β at both ends yields again, we can obtain the second moment:
M 1 β = d 2 M 0 d β 2 = 0 t τ 2 e β τ erf ( α τ ) d τ = M 2 .
The n-th moment is defined as
M n = 0 t τ n e β τ erf ( α τ ) d τ .
Therefore, there are the recursive formulas and general formulas:
M n = M n 1 β , M n = n M 0 β n .
By exporting the 0-order moment M 0 , any n-th moment M n can be easily exported.
The above equations can also deduce the recursive formula for the n-th moment of the error function. Let β 0 :
N 0 = lim β 0 M 0 = lim β 0 0 t e β τ erf ( α τ ) d τ = 0 t erf ( α τ ) d τ ,
N n = lim β 0 M n = lim β 0 0 t τ n e β τ erf ( α τ ) d τ = 0 t τ n erf ( α τ ) d τ = lim β 0 M n 1 β = lim β 0 n M 0 β n .
Therefore, there exists an extremely simple recursive relationship:
N n + 1 = lim β 0 M n β .
For the first item, the following equation holds:
N 0 = 0 t erf ( α τ ) d τ = t erf ( α t ) + 1 2 α 2 π Γ 3 / 2 , α t 1 ,
where Γ s , x is the upper incomplete Gamma function Γ s , x = x t s 1 e t d t . Equation (59) can also be written as
N 0 = 0 t erf ( α τ ) d τ = t erf ( α t ) + 1 π α γ 3 / 2 , α t .
Γ s , x and γ s , x satisfy the following relationship [23,24]:
Γ s , x + γ s , x = Γ s
Equation (60) is a special case of Equation (30). Taking n = 0 , Equation (30) can degenerate into Equation (60).
When f ( τ ) takes any n-degree polynomial function, the weighted integration of the error function 0 t erf ( α τ ) f ( τ ) d τ is not only related to the error function but also to the Gamma function.

6. Algebraic Correlation between Error Function and Gamma Function

The authors of [6] involves the algebraic correlation between error function and Gamma function.
Equation (25) shows that the integration of the error function is related to the error function and exponential function. Equation (59) indicates that the integration of the error function is related to the error function and Gamma function. Simultaneous Equations (25) and (59) are as follows:
0 t erf ( α τ ) d τ = 1 π α t e α t + 1 2 α + t erf ( α t ) = t erf ( α t ) + 1 2 α 2 π Γ 3 2 , α t 1 .
By organizing the second equation, we can obtain that (the proof process can be found in Appendix B):
erf α t + 2 π Γ 3 2 , α t = 1 + 2 π α t e α t .
Equation (62) indicates that there is such a simple algebraic identity between the error function and Gamma function.

7. Fractional-Order Correlation between Error Function and Bessel Function

When studying the convolutional kernel function of bone fractal operators, we confirm that there is a fractional-order correlation between the error function and Bessel function. In [6], we confirm that
p 1 p + α η t = p erf α t α = e α t 2 J ¯ 0 α t 2 .
Equation (63) indicates that the 1/2-order fractional differential of the error function erf α t α , i.e., the operator p acting on erf α t α , generates an exponential weighted Bessel function e α t 2 J ¯ 0 α t 2 . It should be emphasized that the overbar notation of the Bessel functions is not conventional [25]. Equation (63) can also be understood in this way: the 0-order Bessel function can be represented by the fractional derivative of the error function.
Equation (63) reveals the fractional-order correlation between the error function and the 0-order corrected Bessel function.
By combining Equations (62) and (63), we can conclude that
p 1 p + α η t = 1 α p 1 + 2 π α t e α t 2 π Γ 3 2 , α t = e α t 2 J ¯ 0 α t 2 .
Equation (64) shows the fractional-order correlation between the Gamma function and Bessel function. According to Equation (63), we give the following proposition.
Proposition 4.
The correlation between the 0-order Bessel function and the error function can be established through a 1/2-order fractional differential operator p .

8. Fractional-Order Correlation between Error Function and Struve Function

When studying hemodynamics, Zhou et al. derived the kernel functions of fractal operators p 2 + α 2 p and p 2 α 2 p as follows [26]:
p 2 + α 2 p η t = 0 t α τ J 1 α τ d τ ,
p 2 α 2 p η t = 0 t α τ J ¯ 1 α τ d τ .
Zhou et al. also confirmed that the kernel function at the right end of Equations (65) and (66) can be integrated into an explicit form, namely [26]:
0 t 1 τ J 1 α τ d τ = 1 2 2 α t J 0 α t J 1 α t + π α t J 1 α t H 0 α t J 0 α t H 1 α t ,
0 t 1 τ J ¯ 1 α τ d τ = 1 2 2 α t J ¯ 0 α t J ¯ 1 α t + π α t J ¯ 0 α t H ¯ 1 α t J ¯ 1 α t H ¯ 0 α t ,
wherein H 0 and H 1 are 0-order and 1st-order Struve functions, respectively. H ¯ 0 and H ¯ 1 are 0th-order and 1st-order modified Struve functions, respectively [27]. The above equations show that there is a profound intrinsic correlation between the Bessel function and the Struve function, as well as between the modified Bessel function and the modified Struve function, and Equations (67) and (68) have similar forms, reflecting good symmetry.
In Section 7, we established the fractional-order correlation between the error function and Bessel function using the operator p . If we make the following transformation to Equation (66), we have
p 2 α 2 p η t = p α p + α η t p η t .
It is important to note the relationship p η t = δ t and substitute Equation (12) into Equation (69):
p 2 α 2 p η t = p + α α erf α t + e α t π t δ t = e α t p e α t α erf α t + e α t π t δ t = 1 π e α t d d t 0 t e α τ t τ α erf α τ + e α τ π τ d τ δ t .
According to Equations (66) and (70), we can conclude that
e α t p e α t α erf α t + e α t π t δ t = 1 2 2 α t J ¯ 0 α t J ¯ 1 α t + π α t J ¯ 0 α t H ¯ 1 α t J ¯ 1 α t H ¯ 0 α t = e α t π d d t 0 t e α τ t τ α erf α τ + e α τ π τ d τ δ t .
Equation (71) indicates that there is a fractional-order correlation between the error function and modified Struve function. Therefore, we give the following proposition.
Proposition 5.
The correlation between the modified Struve function and the error function can be established through a 1/2-order fractional differential operator p .
The above equation can be regarded as a 1/2-order fractional differential equation:
e α t p f t δ t = 1 2 2 α t J ¯ 0 α t J ¯ 1 α t + π α t J ¯ 0 α t H ¯ 1 α t J ¯ 1 α t H ¯ 0 α t = e α t π d d t 0 t e α τ t τ f t d τ δ t .
Equation (72) must possess the following special solution:
f t = e α t α erf α t + e α t π t .
By combining Equations (62) and (70), we can conclude that
p 2 α 2 p η t = e α t p e α t α 1 + 2 π α t e α t 2 π Γ 3 2 , α t + e α t π t δ t = 1 π e α t d d t 0 t e α τ t τ α 1 + 2 π α t e α t 2 π Γ 3 2 , α t + e α τ π τ d τ δ t .
By combining Equations (66) and (74), we can conclude that
e α t p e α t α 1 + 2 π α t e α t 2 π Γ 3 2 , α t + e α t π t δ t = 1 2 2 α t J ¯ 0 α t J ¯ 1 α t + π α t J ¯ 0 α t H ¯ 1 α t J ¯ 1 α t H ¯ 0 α t = e α t π d d t 0 t e α τ t τ α 1 + 2 π α τ e α τ 2 π Γ 3 2 , α τ + e α τ π τ d τ δ t .
Equation (75) shows that there is a fractional-order correlation between the Gamma function and modified Struve function. Therefore, we give the following proposition.
Proposition 6.
The correlation between the modified Struve function and the Gamma function can be established through a 1/2-order fractional differential operator p .
Equation (75) can also be regarded as a 1/2-order fractional differential equation:
e α t p f t δ t = 1 2 2 α t J ¯ 0 α t J ¯ 1 α t + π α t J ¯ 0 α t H ¯ 1 α t J ¯ 1 α t H ¯ 0 α t = e α t π d d t 0 t f t t τ d τ δ t ,
and the above equation must possess the following special solution:
f t = e α t α 1 + 2 π α t e α t 2 π Γ 3 2 , α t + e α t π t .
It is important to note that the Struve function has the following recursive properties:
H n 1 t H n + 1 t = 2 H n t 1 π Γ n + 3 2 t 2 n ,
H n 1 t + H n + 1 t = 2 n t H n t + 1 π Γ n + 3 2 t 2 n .
The modified Struve function H ¯ n t has similar recursive properties, which will not be further elaborated here. Equations (77) and (78) show that the correlation between the Struve function and Gamma function is diverse. As mentioned above, p in Equations (71) and (75) is a 1/2-order operator, which characterizes the fractional-order correlations. In Equation (77), H n t is the first derivative, which characterizes the integer-order correlations. There is no derivative in Equation (78), which characterizes the algebraic correlations.

9. Fractional-Order Correlation between Error Function and MeijerG Function

When studying the integral properties of the error function, we found that when the input signal f ( τ ) takes the logarithmic function ln τ , the following equation holds (proof can be found in Appendix C):
0 t ln α τ erf α τ d τ = 1 2 α 1 ln 4 γ e + 2 1 ln α t 1 π Γ 3 2 , α t + α t erf α t + 1 π α G 2 , 3 3 , 0 1 , 1 0 , 0 , 3 2 α t .
Usually, the Euler constant is denoted by γ e , γ e = lim n 1 + 1 2 + 1 3 + + 1 n ln n , and G is the MeijerG function, also known as the hyperbolic integral function, which is a common special function [27].
Equation (79) shows that the seemingly unrelated error function, Gamma function, and MeijerG function have a profound intrinsic correlation.

10. Fractional-Order Correlation between Error Function and Generalized Hypergeometric Function

In Section 5, we provide the expression for the n-th moment of the characteristic function when f ( τ ) takes any n-degree polynomial function (Equation (40)). In fact, when f ( τ ) takes τ n , the n-th moment of the characteristic function can also be expressed in the following form (the proof process can be found in Appendix D):
0 t τ n e α τ t τ d τ = 2 t n + 1 / 2 C n n + 1 / 2 F 1 1 n + 1 ; n + 3 / 2 ; α t .
In Equation (80), F 1 1 n + 1 ; n + 3 / 2 ; α t is the generalized hypergeometric function and C n n + 1 / 2 is the binomial coefficient. In mathematics, hypergeometric functions are defined using hypergeometric series, and many special functions are their special cases or limits [27]. The n-th moment at the left end of Equation (80) is related to the fractional-order operator. Therefore, the generalized hypergeometric function at the right end characterizes the fractional-order effect.
We can rewrite Equation (80) as
0 t τ n e α τ t τ d τ = π α e α t t n erf α t 1 π k = 0 n 1 p k γ n + 1 / 2 k , α t α n k ,
we can conclude that
0 t τ n e α τ t τ d τ = 2 t n + 1 / 2 C n n + 1 / 2 F 1 1 n + 1 ; n + 3 / 2 ; α t = π α e α t t n erf α t 1 π k = 0 n 1 p k γ n + 1 / 2 k , α t α n k .
Therefore, we give the following proposition.
Proposition 7.
There is an integer-order integral correlation between the hypergeometric series, error function, and Gamma function.
This proposition generally holds, but there are also special cases, such as those that follow.
Taking n = 0 in Equation (81), we can obtain the following special case:
0 t e α τ t τ d τ = 2 t C 1 / 2 , 0 F 1 1 1 ; 3 / 2 ; α t = π α e α t erf α t .
The second equation of Equation (82) is transformed into
erf α t = 2 α t π C 0 1 / 2 e α t F 1 1 1 ; 3 / 2 ; α t .
Equation (83) indicates that the error function can be decomposed into the product of an exponential function and a hypergeometric function. This indicates that there is an algebraic identity between the hypergeometric function and error function. It is important to note that hypergeometric functions constitute a vast set. Equation (83) shows that an element in this set, i.e., F 1 1 1 ; 3 / 2 ; α t , can be represented by the error function.
Taking n = 1 in Equation (81), we can obtain the following special case:
0 t τ e α τ t τ d τ = 2 t 3 / 2 C 1 3 / 2 F 1 1 2 ; 5 / 2 ; α t = π α e α t t erf α t 1 π γ 3 / 2 , α t α .
The second equation of Equation (84) indicates that there is also a simple algebraic transformation relationship between the error function, Gamma function, and hypergeometric function.
By combining Equations (83) and (84), we can conclude that
0 t τ e α τ t τ d τ = 2 t 3 / 2 C 1 3 / 2 F 1 1 2 ; 5 / 2 ; α t = π α e α t 2 t α t π C 0 1 / 2 e α t F 1 1 1 ; 3 / 2 ; α t 1 π γ 3 / 2 , α t α .
Equation (85) is organized as follows:
γ 3 / 2 , α t = 2 α t 3 / 2 e α t 1 C 0 1 / 2 F 1 1 1 ; 3 / 2 ; α t 1 C 1 3 / 2 F 1 1 2 ; 5 / 2 ; α t .
Equation (86) indicates that the Gamma function can be expressed as the product of a power function, an exponential function, and a hypergeometric function.
Similarly, Tremblay and Choi et al. [28,29] obtained a series of transformation formulas for hypergeometric functions through a fractional-order operator, providing a new perspective for establishing the correlations between special functions.
Thus, this article uses bone fractal operators as the backdrop, the error function as the link, and combines the symmetry evolution laws of functional fractal trees in the physical fractal space to reveal the fractional-order correlation between different special functions. Moreover, in the study of the calculus properties of special functions, a notable phenomenon of symmetry was also discovered. The network of interrelationships between special functions, with the error function as its core, is shown in Figure 2.

11. Conclusions

This article confirms that the integral
0 t e α τ t τ d τ
is an intrinsic form in the convolutional kernel function of bone fractal operators. It demonstrates that the output signal is still related to the error function under varying input signals, such as external loads or stimuli. By centering on the error function, it reveals fractional-order correlations among various special functions, including the error function, Bessel function, Struve function, Gamma function, Meijer G function, and generalized hypergeometric function. These correlations suggest that many special functions are interdependent rather than independent. Biological fractal operators in the physical fractal space offer a novel perspective for understanding these unique correlations.

Author Contributions

Conceptualization, Z.J., G.P., C.L., T.Z. and Y.Y.; methodology, Z.J., G.P., C.L., T.Z. and Y.Y.; software, Z.J. and Y.Y.; validation, Z.J., C.L. and Y.Y.; formal analysis, Z.J., C.L. and Y.Y.; investigation, Z.J., T.Z. and Y.Y.; writing—original draft preparation, Y.Y. and Z.J.; writing—review and editing, Y.Y. and Z.J. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China, grant number 12050001.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Calculation of t n p erf α t

The last term of Equation (40) can calculate convolution using the Laplace transform. For the Laplace transform of general convolution t n f t
L t n f t = L t n L f t ,
wherein
L t n = 0 t n e s t d t = t n + 1 n + 1 e s t 0 0 t n + 1 n + 1 d e s t = s n + 1 0 t n + 1 e s t d t = s n + 1 L t n + 1 .
Therefore, there is a recursive formula:
L t n + 1 = n + 1 s L t n , L t n = n s L t n 1 .
Note:
L t 0 = L 1 = 1 s .
Therefore, the following general formula can be written as
L t n = n ! s n + 1 .
Therefore, Equation (A1) can be transformed into
L t n f t = n ! s n + 1 L f t .
We can perform Laplace inverse transformation on Equation (A5):
t n f t = n ! p n + 1 f t .
If taking f t = p erf α t , we can obtain
t n p erf α t = n ! p n + 1 p erf α t = n ! p n erf α t .

Appendix B. Prove the Algebraic Correlation between Error Function and Gamma Function

We can prove that the following algebraic transformation holds:
erf α t + 2 π Γ 3 2 , α t = 1 + 2 π α t e α t .
Referring to the Laplace transformation table, we have
L erf α t = α s s + α , L Γ 3 2 , α t = Γ 3 2 s 1 1 1 + s / α 3 2 ,
L 1 = 1 s , L t e α t = L t s + α = Γ 3 2 s + α 3 2 .
Therefore, the Laplace transform for Equation (A8) can be written as
α s s + α + 2 π s Γ 3 2 1 1 + s / α 3 2 = 1 s + 2 α π Γ 3 2 s + α 3 2 .
Note:
Γ 3 2 = π 2 ,
therefore, the left end of Equation (A9) is further organized as follows:
LHS = α s s + α + 2 π s Γ 3 2 1 1 + s / α 3 2 = α s s + α + 1 s s 1 1 + s / α 3 2 = α s s + α α 3 2 s 1 1 + s / α 3 2 + 1 s = α s s + α 1 α s + α + 1 s = RHS .
Thus, Equation (62) is confirmed.

Appendix C. Prove the Fractional-Order Correlation between Error Function and MeijerG Function

According to the differential properties of error function,
p erf α t = α π t e α t ,
therefore, the left end of Equation (79) can be rewritten as
LHS = 1 p ln α t erf α t = 1 p 2 α π t e α t ln α t + 1 t erf α t .
Then, we can obtain the Laplace transform of Equation (A11)
L LHS s = 1 s 2 L α π t e α t ln α t + 1 t erf α t = α s 2 L ln α t π t s + α + 1 s 2 s + L erf α t p d p = α s 2 s + α ln 4 s α + 1 + γ e + 2 s 2 sinh 1 α s .
It is important to note the right end (RHS) of Equation (79). With the help of Equation (A8), the second term of RHS can be transformed into
2 π Γ 3 2 , α t + 2 α t erf α t = 1 + 2 α t π e α t + 2 α t 1 erf α t .
We can multiply both ends of Equation (A13) by 1 ln α t and perform a Laplace transform:
L 1 ln α t 2 π Γ 3 2 , α t + 2 α t erf α t = L 1 ln α t s + L 2 1 ln α t α t π s + α + L 2 α t 1 erf α t + L ln α t erf α t = 1 + γ e + ln s α s 2 α s s + α + 2 α α s 2 s + α ln 4 s α + 1 + γ + 2 s 4 α s 2 sinh 1 α s .
The general form of the MeijerG function satisfies the Laplace transform:
L t η G p , q m , n a p b q α t s = s η 1 G p + 1 , q m , n + 1 η , a p b q α s ,
we can derive
L 2 π G 2 , 3 3 , 0 1 , 1 0 , 0 , 3 2 α t = 2 π s G 3 , 3 3 , 1 0 , 1 , 1 0 , 0 , 3 2 α s = 2 s sinh 1 α s 2 α s s + α + 2 s ln s α + 2 ln 4 s .
Finally, considering all parts of RHS, we can calculate its Laplace transform:
L RHS = 1 λ ln 4 2 α s 1 2 α L 1 ln α t 2 π Γ 3 2 , α t + 2 α t erf α t s + 1 2 α L 2 π G 2 , 3 3 , 0 1 , 1 0 , 0 , 3 2 α t = α s 2 s + α ln 4 s α + 1 + γ e + 2 s 2 sinh 1 α s = L LHS .
Thus, Equation (79) is confirmed.

Appendix D. Prove the Fractional-Order Correlation between Error Function and Generalized Hypergeometric Function

Equivalent to proving Equation (80) holds:
0 t τ n e α τ t τ d τ = 2 t n + 1 / 2 C n n + 1 / 2 F 1 1 n + 1 ; n + 3 / 2 ; α t .
Equation (A17) can be proved using the Laplace transform or more directly by expanding it using a power series. The left end can be written in the convolutional form
0 t τ n e α τ t τ d τ = t n e α t 1 t .
Introducing the expansion of the exponential function,
e α t = m = 0 + α t m Γ m + 1 ,
therefore, the convolution at the right end of Equation (A18) can be written as
t n e α t 1 t = m = 0 + α m t m + n Γ m + 1 1 t .
According to the convolutional properties of power function:
t m Γ m + 1 t n Γ n + 1 = t m + n + 1 Γ m + n + 1 ,
therefore, we have
t m + n Γ m + 1 1 t = Γ m + n + 1 Γ 1 / 2 + 1 Γ m + 1 t m + n Γ m + n + 1 t 1 / 2 Γ 1 / 2 + 1 = Γ m + n + 1 Γ 1 / 2 Γ m + 1 t m + n + 1 / 2 Γ m + n + 3 / 2 = π Γ m + n + 1 Γ m + 1 t m + n + 1 / 2 Γ m + n + 3 / 2 .
Further export
t n e α t 1 t = = π m = 0 + α m t m + n + 1 / 2 Γ m + n + 3 / 2 Γ m + n + 1 Γ m + 1 = π t n + 1 / 2 m = 0 + α t m Γ m + n + 3 / 2 Γ m + n + 1 Γ m + 1 = π t n + 1 / 2 m = 0 + β m α t m ,
finally, we have
0 t τ n e α τ t τ d τ = π t n + 1 / 2 m = 0 + β m α t m ,
wherein β m represents the coefficient of the m-th item:
β m = Γ m + n + 1 Γ m + n + 3 / 2 Γ m + 1 ,
the ratio of adjacent coefficients is
β m + 1 β m = Γ m + n + 3 / 2 Γ m + 1 Γ m + n + 2 Γ m + n + 5 / 2 Γ m + 2 Γ m + n + 1 = m + n + 1 m + n + 3 / 2 m + 1 .
Equation (A25) satisfies the definition of the Generalized Hypergeometric function F q p a 1 , , a p ; b 1 , , b q ; z , wherein the parameter values in the Hypergeometric function are as follows:
p = q = 1 , a 1 = n + 1 , b 1 = n + 3 / 2 , z = α t ,
therefore, Equation (A23) can be written as
0 t τ n e α τ t τ d τ = π β 0 t n + 1 / 2 F 1 1 n + 1 ; n + 3 / 2 ; α t = 2 t n + 1 / 2 C n n + 1 / 2 F 1 1 n + 1 ; n + 3 / 2 ; α t .
Thus, Equation (80) is confirmed.

References

  1. Guo, J.Q.; Yin, Y.J.; Ren, G.X. Abstraction and operator characterization of fractal ladder viscoelastic hyper-cell for ligaments and tendons. Appl. Math. Mech.-Engl. Ed. 2019, 40, 1429–1448. [Google Scholar] [CrossRef]
  2. Guo, J.; Yin, Y.; Hu, X.; Ren, G. Self-similar network model for fractional-order neuronal spiking: Implications of dendritic spine functions. Nonlin. Dyn. 2020, 100, 921–935. [Google Scholar] [CrossRef]
  3. Peng, G.; Guo, J.Q.; Yin, Y.J. Self-Similar functional circuit models of arteries and deterministic fractal operators: Theoretical revelation for biomimetic materials. Int. J. Mol. Sci. 2021, 22, 12897. [Google Scholar] [CrossRef] [PubMed]
  4. Jacobs, B.; Driscoll, L.; Schall, M. Life-span dendritic and spine changes in areas 10 and 18 of human cortex: A quantitative golgi study. J. Comp. Neurol. 1997, 386, 661–680. [Google Scholar] [CrossRef]
  5. Duan, H.; Wearne, S.L.; Rocher, A.B.; Macedo, A.; Morrison, J.H.; Hof, P.R. Age-related dendritic and spine changes in corticocortically projecting neurons in macaque monkeys. Cereb. Cortex. 2003, 13, 950–961. [Google Scholar] [CrossRef] [PubMed]
  6. Jian, Z.; Peng, G.; Luo, C.; Zhou, T.; Yin, Y. Convolution Kernel Function and Its Invariance Properties of Bone Fractal Operators. Fractal Fract. 2024, 8, 151. [Google Scholar] [CrossRef]
  7. Gao, H.; Ji, B.; Jäger, I.L.; Arzt, E.; Fratzl, P. Materials become insensitive to flaws at nanoscale: Lessons from nature. Proc. Natl. Acad. Sci. USA 2003, 100, 5597–5600. [Google Scholar] [CrossRef]
  8. Shao, Y.; Zhao, H.P.; Feng, X.Q.; Gao, H. Discontinuous crack-bridging model for fracture toughness analysis of nacre. J. Mech. Phys. Solids. 2012, 60, 1400–1419. [Google Scholar] [CrossRef]
  9. Jäger, I.; Fratzl, P. Mineralized collagen Mbrils: A mechanical model with a staggered arrangement of mineral particles. Biophys. J. 2000, 79, 1737–1746. [Google Scholar] [CrossRef]
  10. Norman, T.L.; Vashishth, D.; Burr, D.B. Fracture toughness of human bone under tension. J. Biomech. 1995, 28, 309–320. [Google Scholar] [CrossRef]
  11. Rho, J.Y.; Kuhn-Spearing, L.; Zioupos, P. Mechanical properties and the hierarchical structure of bone. Med. Eng. Phys. 1998, 20, 92–102. [Google Scholar] [CrossRef] [PubMed]
  12. Fratzl, P.; Weinkamer, R. Nature’s hierarchical materials. Prog. Mater. Sci. 2007, 52, 1263–1334. [Google Scholar] [CrossRef]
  13. Reznikov, N.; Bilton, M.; Lari, L.; Stevens, M.M.; Kröger, R. Fractal-like hierarchical organization of bone begins at the nanoscale. Science 2018, 360, eaao2189. [Google Scholar] [CrossRef] [PubMed]
  14. Ji, B.H.; Gao, H.J. Mechanical properties of nanostructure of biological materials. J. Mech. Phys. Solids 2004, 52, 1963–1990. [Google Scholar] [CrossRef]
  15. Ji, B.H.; Gao, H.J. A study of fracture mechanisms in biological nano-composites via the virtual internal bond model. Mater. Sci. Eng. A. 2004, 366, 96–103. [Google Scholar] [CrossRef]
  16. Gao, H.J. Application of fracture mechanics concepts to hierarchical biomechanics of bone and bone-like materials. Int. J. Fracture. 2006, 138, 101–137. [Google Scholar] [CrossRef]
  17. Abdulnabi, F.F.; FAl-Janaby, H.; Ghanim, F.; Lupaș, A.A. Geometric Features of the Hurwitz–Lerch Zeta Type Function Based on Differential Subordination Method. Symmetry 2024, 16, 784. [Google Scholar] [CrossRef]
  18. Mondal, S.R.; Giri, M.K.; Kondooru, R. Sufficient Conditions for Linear Operators Related to Confluent Hypergeometric Function and Generalized Bessel Function of the First Kind to Belong to a Certain Class of Analytic Functions. Symmetry 2024, 16, 662. [Google Scholar] [CrossRef]
  19. Olteanu, O. Functional and Operatorial Equations Defined Implicitly and Moment Problems. Symmetry 2024, 16, 152. [Google Scholar] [CrossRef]
  20. Courant, R.; Hilbert, D. Methods of Mathematical Physics; Interscience Publishers: New York, NY, USA, 1962. [Google Scholar]
  21. Mikusinski, J. Operational Calculus, 2nd ed.; Pergamon Press: Oxford, UK, 1983. [Google Scholar]
  22. Watson, G.H. A Treatise on the Theory of Bessel Functions; Cambridge University Press: Cambridge, UK, 1952. [Google Scholar]
  23. Luke, Y.L. The Special Functions and Their Applications; Academic Press: New York, NY, USA; London, UK, 1969. [Google Scholar]
  24. Mohammed, P.O.; Sarikaya, M.Z.; Baleanu, D. On the Generalized Hermite–Hadamard Inequalities via the Tempered Fractional Integrals. Symmetry 2020, 12, 595. [Google Scholar] [CrossRef]
  25. Moreira, W.L.; Neves, A.A.R.; Garbos, M.K.; Euser, T.G.; Cesar, C.L. Expansion of arbitrary electromagnetic fields in terms of vector spherical wave functions. Opt. Express. 2016, 24, 2370–2382. [Google Scholar] [CrossRef]
  26. Zhou, T.; Yin, Y.; Peng, G.; Luo, C.; Jian, Z. Fractal Operators Abstracted from Arterial Blood Flow. Fractal Fract. 2024, 8, 420. [Google Scholar] [CrossRef]
  27. Erdélyi, A. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  28. Tremblay, R. Using the Well-Poised Fractional Calculus Operator to obtain transformations of the Gauss hypergeometric function with higher level arguments. Montes Taurus J. Pure Appl. Math. 2021, 3, 260–283. [Google Scholar]
  29. Choi, J.; Rathie, A.K. General summation formulas for the KAMPEDE FERIET function. Montes Taurus J. Pure Appl. Math. 2019, 1, 107–128. [Google Scholar]
Figure 1. Symmetry evolution laws of functional fractal trees in the physical fractal space: (a) fractal cell abstracted from nerve fibers [2]; (b) fractal cell abstracted from neuronal dendrites [2]; (c) fractal cell abstracted from compact bone fibers [6].
Figure 1. Symmetry evolution laws of functional fractal trees in the physical fractal space: (a) fractal cell abstracted from nerve fibers [2]; (b) fractal cell abstracted from neuronal dendrites [2]; (c) fractal cell abstracted from compact bone fibers [6].
Symmetry 16 01279 g001
Figure 2. The network of interrelationships between special functions.
Figure 2. The network of interrelationships between special functions.
Symmetry 16 01279 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jian, Z.; Luo, C.; Zhou, T.; Peng, G.; Yin, Y. Fractional-Order Correlation between Special Functions Inspired by Bone Fractal Operators. Symmetry 2024, 16, 1279. https://doi.org/10.3390/sym16101279

AMA Style

Jian Z, Luo C, Zhou T, Peng G, Yin Y. Fractional-Order Correlation between Special Functions Inspired by Bone Fractal Operators. Symmetry. 2024; 16(10):1279. https://doi.org/10.3390/sym16101279

Chicago/Turabian Style

Jian, Zhimo, Chaoqian Luo, Tianyi Zhou, Gang Peng, and Yajun Yin. 2024. "Fractional-Order Correlation between Special Functions Inspired by Bone Fractal Operators" Symmetry 16, no. 10: 1279. https://doi.org/10.3390/sym16101279

APA Style

Jian, Z., Luo, C., Zhou, T., Peng, G., & Yin, Y. (2024). Fractional-Order Correlation between Special Functions Inspired by Bone Fractal Operators. Symmetry, 16(10), 1279. https://doi.org/10.3390/sym16101279

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop