About the Use of Generalized Forms of Derivatives in the Study of Electromagnetic Problems

: The use of non-local operators, deﬁning Riemann–Liouville or Caputo derivatives, is a very useful tool to study problems involving non-conventional diffusion problems. The case of electric circuits, ruled by non-integer derivatives or capacitors with fractional dielectric permittivity, is a fairly natural frame of relevant applications. We use techniques, involving generalized exponential operators, to obtain suitable solutions for this type of problems and eventually discuss speciﬁc problems in applications.


Introduction
The use of non-integer and integral operators [1][2][3][4][5][6] has evolved into the Differ-Integral Calculus (DIC). The mathematical environment in which DIC has grown is that of the generalized transform. Most of the associated concepts are well framed within the context of appropriate integral transforms [2,7]. The techniques associated with the DIC, combined with other methods of operational nature, has provided significant progress within the context of applied mathematics in general and of applied classical electromagnetic field theory in particular and the relevant roots can be traced back to the original works of Heaviside [8].
The associated technicalities have played, in the past, important roles in the study of equations describing evolution problems in relativistic Quantum Mechanics and tracing back to the work of Klein Gordon and Dirac [9]. Within this context, an important step has been the understanding of the role played of DIC within the most general framework of pseudo and non-local operator theory [10][11][12][13]. The most natural context for the study of differential equations, be they of ordinary or partial nature, is the use of techniques tracing back to the evolution operator and to generalized integral transforms, including those of Lévy-type [14]. We will provide a description of these methods by reviewing a number of applications of what is defined as (in a broad sense) fractional calculus.
An application field where fractional calculus has received an increasing attention over the years is CMOS-based applications. According to the International Technology Roadmap for Semiconductors (ITRS), the expected cut-off frequency of CMOS transistors is 0.9 THz in 2021. At such high frequencies, the primary challenge in CMOS-based THz design is to develop accurate device models that allow one to take the loss from strong frequency-dependent dispersion and non quasi-static effects in THz into account. Traditionally, on-chip interconnects are characterized by the distributed integer order RLGC model. Drude's classical relaxation-effect model is used for the skin effect modeling [15]. In addition, the loss due to dielectric polarization and dipole rotation can be modeled by a dielectric-loss still using integer-order models. However, such an integer-order model is insufficient to describe the transmission line performance at the THz region because the loss term in transmission line theory is difficult to model the dispersion loss and non-quasistatic effects [16] which may cause a large deviation in the THz frequency region. Some valuable models of fractional-order models of transmission lines are presented in [17,18].
A generalization of the CMOS circuit driving a distributed fractional RLC load into the fractional-order domain has been presented in [19], making the on-chip interconnect structure easy to adjust by including the effect of the fractional orders. Hence, incorporation of fractional-order models becomes crucial to properly address both electromagnetic compatibility as well as signal integrity issues.
Another field where fractional-order models have become important is that of pulsed electric fields (PEFs). The main goal of this discipline is the study of the interaction between biological tissues and electromagnetic fields aiming at the investigation of the use of ultrashort pulses in biomedical applications [20]. It is known that the electric field distribution excited in biological media mainly depends on the electric properties of tissues. The dielectric properties of biological tissues result from the interaction of electromagnetic energy with the tissue constituents at the cellular and molecular level which is strongly affected by the bound water content. The dielectric properties of many biological materials exhibit an experimental dielectric response in the frequency domain that cannot be described by a simple exponential expression with a single relaxation time (Debye model). On the contrary, empirical fractional-order relationships including Cole-Cole (C-C), Cole-Davidson (C-D), and Havriliak--Negami (H-N) equations have been proposed in order to fit such dielectric spectra. The evaluation of the transient propagation of electromagnetic waves in these types of fractional-order dielectric media is not trivial since it contains the transformation of the fractional-order expressions from the frequency domain to the time-domain. It has been addressed in the context of the finite different time-domain (FDTD) method by resorting to the auxiliary differential equation (ADE) approach [21,22]. An alternative approach based on the application of the Riemann-Liouville theory to the analysis of C-C media has been proposed in [23,24]. The double fractional derivative operator relevant to the general H-N response is directly incorporated in the FDTD basic algorithm in [25].
What we would like to convey here is the strong interplay between DIC and integral transforms. They constitute the relevant backbone and probably the intimate mathematical flavor. The paper is therefore planned as follows.
In Section 2 we show how the concept of fractional integral emerges from the study of the Lamb-Bateman integral equation. The emergency of the Liouville-Riemann definition of fractional derivative is then shown to naturally emerge from the Laplace transform, which is then exploited to treat fractional ODE's.
The rest of the papers deal with specific applications of "fractional" electric circuits and extension to RC circuits with non-standard permittivity.
Finally, we go back to quantum-like models by treating the case of relativistic heat equations.

Lamb-Bateman Equation, Operational Methods, and Fractional Calculus
We give an idea of the underlying technicalities, also exploited in this article, by discussing the solution of the Lamb-Bateman equation [26][27][28].
Example 1. The Lamb-Bateman equation, originally proposed for the solution of a solitary wave, is written as where u(x) is a function to be determined in terms of f (x). The relevant solution is readily achieved by the use of operational methods [29], allowing one to recast it in the form which is just a consequence of the shift operator identity [29] e a∂ x f (x) = f (x + a).
If we take the freedom of treating the derivative operator, in the definition ofL, as an ordinary algebraic quantity, we can work out the Gaussian integral and end up with the explicit definition of the operatorL asL and eventually We accordingly expressed the solution of the Lamb-Bateman equation in terms of the derivative of order 1 2 of the function f (x).
The next step is that of stipulating the action of a fractional derivative on a given function.

Exercise 1.
We recall first the Laplace transform identity which, by using Equation (3), allows writing of the following Changing the integration variable and arranging the limits of integration, we eventually find The lowest integration limit, following from (8), is −∞; we put a generic a to be consistent with the current definition of the differintegral operator. If a = 0 the sub-index on the left of the derivative operator will be neglected. The use of the semigroup property of the fractional differintegral operator [1][2][3][4][5][6] finally leads to The apex in front of g(x) denotes first-order derivative.
In the previous Exercise, we adopted 0 as the lower limit. The relevant value should however be chosen in agreement with the specific problem under study. Further comments will be provided in the final comments. Remark 1. More in general, any derivative of order α can be written as which is the definition of the fractional derivative according to Caputo [30,31]. Equation (12) states, with the ordinary derivative, that the fractional derivative of a constant is vanishing. This feature is not contained in other forms of fractional derivatives (e.g., the Riemann-Liouville and Euler derivatives [1][2][3][4][5][6]).
In the following, we limit ourselves to the derivative order with 0 < α < 1 and just touch on the cases with α > 1 later in the article.
The Caputo derivative (as well as the others) is a non-local integral operator, not dissimilar from those entering the treatment of pseudo-differential operators in relativistic quantum mechanic [9][10][11]13,32,33].

Observation 1.
Regarding the use of this family of operators we note, e.g., that the solution of non-homogeneous differential equation where y(x) is the unknown, can be formally written as The use of the Laplace transform identity applied to the operator and of the exponential shift operator (3) yields By exploiting Observation 1 we can provide the solution to the following problem involving partial fractional derivatives.

Problem 1 (Heat-like Equation)
. The solution of evolution-like equations involving partial fractional derivatives as, e.g., can be viewed as a kind of heat equation [13,29,[34][35][36] and solved by the use of an appropriate transform. The use of the formalism of the evolution operator [29,37] yields The action of the evolution operatorÛ(t) on the initial function γ(x) can be obtained via the so-called Lévy transform method [13,14]. To this aim, we note that the following integral transform, concerning the stretching integral, is known as one-sided Lévy stable distribution (Equation (19) is also known as Kohlrausch-Williams-Watts function [38].) g γ (η) [14] and acts in the following way We can therefore manage to transform the exponential evolution operator in Equation (18) aŝ which eventually yields, for the solution of our problem, Other forms of fractional derivatives have appeared in the literature. (21), we consider the case with α = 1 2 , the Lévy distribution writes

Remark 2. If, in Equation
which, for an initial Gaussian distribution, yields the temporal behaviour reported in Figure 1 which displays the interesting feature that, while time increases, the initial gaussian distribution develops the long tail shape characteristic of a Lévy distribution. The same procedure can be applied when dealing with slightly more complicated problems, like in the following problem. (23) whose formal solution reads 1+∂ ν Explicitly, in terms of the Lévy functions g 1 2 (η), according to Equation (19) we find In Figure 2 we reported a comparison between the solutions of the fractional PDE in (17) and in (23) for ν = 1 2 . They are qualitatively similar; they exhibit the development of an asymmetric long tail for increasing time, with a significantly larger reduction in the peaks in the case (23). An analogous technique has been discussed in Ref. [13], where a generalized form of the heat equation (α = 2) has been considered in terms of a double transformation (Lévy & Gauss-Weierstrass).

Remark 3.
In order to obtain a more complete picture of the formalism we have foreseen, we consider the solution of Equation (23) for the cases with α = 1 3 and α = 2 3 . To this aim, we remind that [39,40] with K α (x) being the second kind modified Bessel function and where W (κ,µ) (x) is the Whittaker function and U(a, c; x) the confluent hyper-geometric function of the second kind. The comparison between the solution of Equation (23) with α = 1 3 , 2 3 is eventually reported in Figure 3. As a further example of application of the composite transform method, we consider the solution of the problem given below.
where the unknown function u(x) is expressed in terms of f (x) as Observation 2. Before concluding this introductory section, we note that for 0 < β < 1, we can obtain the solution of the Example 2 in alternative way. We consider the fractional differential equation (30) and recast it according to the following where the action of the fractional exponential operator on the function κ(x) can be specified as It is worth stressing that the use of Equation (12) eventually yields The numerical computation of Equation (31) is accordingly complicated by the necessity of evaluating three nested integrals or four integrals if we consider the alternative following calculation and finally The quadrupole integral on the right is a caveat on the slowness of the associated numerical procedure. Further comments on the subject will be however provided in the concluding section.

Generalized Exponential Operators and Fractional Electrical Circuit
In this section, we consider the extension to the fractional case of the differential equations characterizing an electric circuit.
In the previous section we fixed the few elements of the formalism we will use in the following. We privileged, in our definition of fractional derivative, the Caputo form. The reason for this choice stems from the fact that we prefer to avoid ambiguity due to the non-vanishing of the derivative of a constant.
The modeling of an electrical circuit in terms of non-integer differential equations does not appear, in general, well motivated and often the mathematical formalism is privileged with respect to the effective physical meaning. It should be however noted that the Ohm law is at the basis of the circuital phenomena. This law summarizes, from a macroscopic point of view, the microscopic dynamic of charge carriers inside a conductor, along with the relevant interactions defining the charge drift velocity. Classical Newtonian mechanics is used to provide a link between the applied electric field and charge velocities, whereas the diffusion process may be governed by other mechanisms not strictly associated with classical diffusion and more similar to Lévy type process which, according to the discussion of the previous section, requires fractional heat type equations.
We assume, for the moment, that the fractional derivative can model an actual physical circuit and discuss some aspects of the relevant study, in terms of the solutions of the associated differential equations.

Example 3 (RC Circuit
). An ordinary RC circuit is ruled by the first-order differential equation Equation (38), in non dimensional form, gives The non-homogeneous part of the solution can, according to the methods outlined in the previous section, be written as and a graphical example is reported in Figure 7a, for f (η) = cos(ση), σ = ωτ. The solution was checked numerically to be sure that it is consistent with the original Equation (38).
Before proceeding further it should, however, be noted that, by standard means, the solution of first-order non-homogeneous ODE's of the type (39) gives (41) or, after rearranging the integration variable, We denoted by Q * the initial condition (at η = 0); the first term is therefore the solution of the homogeneous part, while the integral accounts for the non-homogeneous term. If we rearrange the integration variable in Equation (40), it is evident that the two solutions coincide for Q * = 0 and η 0 = −∞. Observation 3. Whilst naïve, in view of the extension to a fractional RC circuit, we note that Equation (38) is just the Kirchhoff voltage law [41] and can also be written as where V C is the voltage across the capacitor. Since RC has the physical dimensions of a time, Equation (43) is dimensionally consistent. This is a point of non-secondary importance whenever one is trying its extension to the fractional case, which can be written as where the constant α 1−γ has been introduced to avoid problems with physical dimensions. Our fractional ordinary differential equation (Frac-ODE) can be written as We accordingly chose τ γ to have the dimensions of [T γ ]; therefore, it can eventually be cast in the d dη which yields If we assume that the source potential has a cosine form and apply the operator identity we can write the solution of our problem as In Figure 7b we provide a comparison between integer and fractional solutions. The non-integer solution differs either in terms of the amplitude and phase but there is no significant difference in the qualitative behavior.
If the potential is not given by an oscillating function, but writes in terms of a generic time-dependent function f (t), we can handle Equation (47) as follows and examples, for γ = 1 3 , with f (x) Gaussian and super-Gaussian, are given in Figure 5a,b, where we have plotted the solution given in (29) which, when β = 1, reduces to Equation (51). Example 4 (RL Circuit). The analytical tools are not different if we deal with an RL circuit. The problem is still ruled by a first-order differential equation. If we cast the Kirchhoff circuit law as where i L is the current flowing through the inductance. The relevant fractional version reads and, in non-dimensional form, eventually reads d dη It is evident that, from the mathematical point of view, there are no differences with respect to the RC example. In order to prove the reliability of the method, we consider a slightly different form of the external potential, namely The solution of our problem reads therefore In Figure 8 we have shown the associated behaviors. In the forthcoming section, we benchmark the results so far obtained by the use of a different procedure employing solutions of the Frac-ODE based on the Mittag-Leffler function.
Before going further we touch on the possibility of extending the above formalism to non-oscillating external potentials.

Remark 4.
The technique we have outlined is appropriate for potentials which can be expressed as an eigenfunction, or a combination of eigenfunctions, of the derivative operator. If f (η) has a Fourier transformf (k), the solution can be obtained in a similar way as shown below In the case of a decaying exponential potential (like a Gaussian or a super-Gaussian) the method, as we have outlined, yields the result reported in Figure 9.

Mittag-Leffler Functions and Fractional Electric Circuits
Before proceeding further it is necessary to fix at least two points: (a) To state the form of the general solution of a non-homogeneous fractional equation including the initial conditions; (b) To establish the consistency between the different solutions we have obtained.
Regarding the point (a), we note that the Caputo (namely for a fractional derivative defined according Caputo) fractional differential equation is solved as where is the Mittag-Leffler function [42][43][44][45] and its first-order derivative is Equation (59) reduces to the ordinary exponential for γ = 1. It is also easily checked that, in the same limit, the solution becomes that of a non-homogeneous first-order ODE.
Regarding point (b), in Figure 10, we have reported a comparison between the solutions of Equation (58) for y 0 = 0 and f (τ) a Gaussian, by using the Fourier transform (57), the Laplace/Lévy transform (51) and solution (59). The different solutions compare well and the uncertainties are due to a poor efficiency of the routine we have used in making multiple integrations. Regarding the integration with Mittag-Leffler, the truncation of the series gives rise to the discrepancy. The agreement can be improved by the increase of the series contributions, at the cost of larger computational times. In Figure 11 we have also plotted a comparison between two solutions of (58) with different γ values. We checked the correctness of the solution procedure by comparing the different integration procedures.

Fractional RLC Circuit
In the previous section, we dealt with first-order ODE and we managed to obtain the corresponding non-homogeneous contributions. Now we provide a natural further extension by studying an RLC circuit [46][47][48][49][50][51].

Problem 3 (RLC Circuit). Regarding a second-order ODE of the type
we can proceed by just extending the outlined procedure and writing with L being the circuit inductance. The solution of Equation (63) for vanishing initial conditions i 0 , d dτ i(τ) | τ=0 can be written as In order to make the previous formal solution effective, we use the identity e −s d 2 and eventually end up with which is easily integrated and the relevant solutions are reported in Figure 12. Alternatively, when the previous solution displays convergence problems, the use of the Fourier-Transform method yields In Figure 12 we report the current wave form for V (τ) ∝ e −τ 2 and different values of A = ω 0 τ L . Small A values determine a strongly damped behavior.
For future convenience, we consider a capacitor defined as an operator specified in terms of fractional derivatives aŝ whose physical meaning will be discussed in the forthcoming section and providing The effect of this new term on the current wave-form is provided in Figure 13.
The behavior of the current, generally oscillating, exhibits a critical dependence on the various parameters defining the operator function P (α,γ) (.), whose physical meaning will be discussed in the forthcoming section. In the previous sections we provided a solution to fractional ODE using a technique privileging the FT method. The procedure we have envisaged is useful when non-homogeneous terms are present and does not apply for homogeneous cases. We would like however to underline that the search for a solution of fractional ODE is not an easy task and a general theory is not yet available. Effective solution techniques can however be employed and, for example, the Laplace transform method offers an efficient tool, some time hampered by the difficulties encountered with the corresponding anti-transform. Before closing this section we like to quote the possibility of reducing a non-homogeneous ODE to a Volterra type integral equation [52] and apply to fractional ODE the wealth of techniques available for that problem.

Summary
Non Debye dielectrics are described by a fractional relative permittivity, which can be put in the form [25] where ε r ∞ and ε r s are the infinite and the static relative permittivity, τ is the principal relaxation time and α and β are positive adjustable parameters between zero and one. The parameters of α and β substantially determine the model used for the dielectric; it is clear that if both α and β are set to zero, we obtain the standard one-pole Debye model. Apart from the Debye model, we consider other three dielectric representations whose characteristics are summarized in Table 1. In the following, for all the proposed models we assume τ = 153 ps.

RC Circuit Tests Including Fractionary Dielectrics
The dispersive nature of dielectrics can be easily included in circuit models through the Fourier transform technique, which permits to work in the frequency domain; finally Inverse Fast Fourier Transform (IFFT) techniques can be employed to restore the timedomain behavior of the quantities under investigation. A possible very useful application for non Debye dielectric models is in the framework of the Partial Elements Equivalent Circuit (PEEC) method, an integral equation based method that permits performance of electromagnetic (EM) simulations of complex structures introducing a circuit representation of the EM phenomena. The PEEC method is very useful in Electromagnetic Compatibility (EMC) and Signal integrity (SI) areas, allowing prediction and solving of electromagnetic interference (EMI) issues.
Dielectric materials are a fundamental part of modern electric and electronic devices and for this reason it is important to accurately include their presence; this can be done in the PEEC method through the concept of excess capacitance [53]. In PEEC models is feasible to include in the excess capacitance also the fractional relative permittivity of Equation (70), in order to provide a more accurate representation of the EM phenomena that affect the structures under observation.
It is reported here the example of an RC series circuit fed by a Gaussian pulse voltage, with the capacitor supposed to consist of flat and parallel plates, filled by a dielectric characterized by a relative permittivity generally described by Equation (70); all the four dielectric laws presented are examined here.
The ordinary differential equation, describing an RC series circuit through the electric charge q(t) in the time-domain, is the following: where R is the resistance value of the RC circuit and V(t) = V 0 f (t) is the forcing function which represents the voltage source applied to the circuit; the capacitance C(t) is a timedependent function because of the dispersive nature of the dielectric inside the capacitor.
In the frequency domain, the capacitor can be described by its complex capacitance where S is the surface of the plates and d is the distance between them. Therefore, its time-domain counterpart is of the form is the time-domain fractional dielectric law described in Equation (71). We now define the parameter: τ c = RC 0 and introduce the dimensionless variable η = t τ c . Hence, the forcing function f (t) is also a f (τ c η) = f * (η) and the variable function q(t) is also expressible as a q(τ c η). The last positions permit the writing of Equation (72) as By considering that dt = τ c dη and the Equation (74), we write Equation (76) as which, using the relation τ c = RC 0 and expressing that ε r (t) = ε r (τ c η), becomes By introducing the parameter q 0 = C 0 V 0 , which has the dimensions of the electric charge, we can write and finally we can obtain the differential equation for the normalized chargeq(τ c η) in the where the relative permittivity can be expressed as A solution of the Equation (80) can be found first by expressing it in the frequency domain, where d dη = jk, k being the dimensionless pulsation whereQ(k) is the Fourier transform with respect to the variable η of the normalized chargẽ q(τ c η). Furthermore, it is necessary to express the fractional dielectric permittivity law in the dimensionless frequency domain as Finally, in the frequency domain, we have the following final expression for the normalized chargeQ where F * (k) is the Fourier transform of f * (η).
We considered a forcing function f * (η) = e −(η−10) 2 , represented in Figure 14, and the data for the test in Table 2. The functionq * (η) is restored by the IFFT algorithm, which is an approximation of the integral below. The results for the time-domain behavior of the normalized charge using the IFFT technique are depicted in Figure 15. It is clear, from Figure 15, that considering different frequency domain dielectric laws for the capacitor changes necessarily also the time-domain behavior of the overall circuit. In order to validate the IFFT technique, in Figure 16 is reported a comparison with the approximated integral computation, for the Havriliak-Negami model, showing a very good agreement between the two techniques. Moreover, it is possible to observe in Figure 17 the normalized charge behavior for a dielectric described by the Havriliak-Negami model, considering three different values of the parameter κ. The main advantage of the IFFT technique is its simplicity and ease of coding. On the other hand, its inclusion in a 3D time-stepping electromagnetic solver is not straightforward since it requires an a priori knowledge of the source term in the integral (15) which is not the case when the fractionary dielectric is part of a more complex 3D model. In this case, recursive schemes are recommended which are based on the knowledge of a rational approximation of the fractional permittivity.

Final Comments
This article has been devoted to the use of mathematical tools, employing operational calculus, for the solution of non-homogeneous fractional ODEs. We have noted that fractional calculus is an extension of ordinary calculus, having the merit of dealing with a new point of view to the differential calculus, expressed in terms of integral operators, interpreted as non-integer derivatives. Even though not explicitly stated, the distinguishing feature of such a formalism is the use of operators suited to deal with non-local problems and physical mechanisms employing memory effects as e.g., dissipation in various materials [3].
We started our discussion considering the search for a solution of the Lamb-Bateman equation, a genuine non-local problem, by the use of a purely operational technique and we have shown that, by keeping the freedom of treating differential operators as ordinary algebraic quantities, we have been able to provide a "natural" meaning for the derivatives of fractional order. In order to corroborate the point we have just raised, we discuss an example from Ref. [54] and mention that the Abel integral equation of the first kind is given by where g(s) is the unknown, b is an upper limit specified by the nature of the kernel k(x, s)(x − s) −α . By assuming that it is just a constant (the unit in our example) we easily obtain the solution, which is contained in the definition of Caputo derivative given in Equation (12); we find indeed which, once inverted yields The example is underscored to stress the interplay between the concept of fractional derivatives and integral operators.
Regarding the second point we have raised (non-locality), we would like to mention the spinless relativistic Schrodinger's equation [9], which, in terms of pseudo-operators, where the operator on the left realizes a pseudo fractional operator. According to what we have learned so far we can cast Equation (89) in the form [13,55] The exponential operator containing the second-order derivative is a diffusion operator, whose action on a given function is provided by the so-called Gauss-Weierstrass transform, namely which holds if the integral on the r.h.s. converges. Therefore, in conclusion, we find We accordingly transformed the original problem containing the square root of a differential operator into an integro-differential equation in which the non-locality is specified by a kind of diffusive kernel. The relevant solution is not easy to achieve using this form, which is however amenable for a perturbative expansion. We furthermore employed a solution procedure based on an evolution operator and Fourier transform method to solve Equation (89). By using a Gaussian as the initial function, we obtain the result reported in Figure 18, where we showed the square modulus of the wave function Ψ(η, τ) at different times. The relevant evolution is characterized by an asymmetric diffusion, including a shift of the peak and the development of a long tail. We also checked that this behavior is compatible (at short times) with the results obtained with the perturbative solution. The use of pseudo-operators, involving square roots of derivative operators in relativistic Quantum Mechanics is delicate either from the mathematical and physical points of view. However, it is scientifically sound and does not imply any doubtful steps. The operators on both sides of Equation (89) are indeed hermitian, and the associated evolution operator is unitary. This is the minimum requirement when a Quantum Mechanics problem is posed, since the physical consequences from the associated mathematical treatment is not flawed by intrinsic "fatal" errors.
Some authors "generalize" the Schroedinger equation by directly replacing (time or spatial) derivative with fractional counterparts. The development of a credible formulation of Fractional Quantum Mechanics (FQM), requires the change of the specific paradigms leading to the formulation of the Schroedinger equation itself. For this reason, in Ref. [56], the problem has been afforded starting from a path integral formulation. This study has led to the conclusion that standard Quantum Mechanics follows from the fractality of Brownian motion and FQM deepen its root into the fractality of Lévy flights. The formal basis of the (non-relativistic) FQM has been considered with adequate mathematical rigor. The Schroedinger equation derived in [56] reads where D a is a constant, adjusting the associated physical units and defining the link between energy and momentum in FQM [56]. The operator on the r.h.s. of Equation (93) is not Hermitian; this implies that the associated eigenvalues are not real and the corresponding evolution operator is not unitary. An important possibility within this respect is offered by the formalism adopted in dissipative Quantum Mechanics [57] whose mathematical (and physical as well) treatment requires the use of bi-orthogonal states [56,57] yielding a more appropriate definition of the operator self-adjointedness condition and a correct tool for the evaluation of matrix elements between different quantum states.
Regarding the Free Particle evolution, Equation (93) yields where we used a more convenient non-dimensional form. The solution of Equation (94) can be obtained using the discussed method. Using the definition of the Riesz fractional derivative [56,58], we obtain the behavior reported in Figure 19. It is qualitatively similar to that displayed in Figure 18, thus suggesting a non-local nature of the underlying dynamics. In order to complete this discussion, we would like to add as a further example of application of the discussion so far developed to the solution of propagation equations of D'Alembert type, namely