Dual Numbers and Operational Umbral Methods

TL;DR

This paper introduces a novel perspective on dual numbers by embedding them within an operational umbral calculus framework, enhancing their application in numerical computations and finite difference calculus.

Contribution

It presents a new formalism that integrates dual numbers with umbral calculus, offering fresh insights and potential computational advantages.

Findings

New algebraic formalism for dual numbers within umbral calculus

Enhanced understanding of dual numbers in finite difference calculus

Potential for improved numerical computation techniques

Abstract

Dual numbers and their higher order version are important tools for numerical computations, and in particular for finite difference calculus. Based upon the relevant algebraic rules and matrix realizations of dual numbers, we will present a novel point of view, embedding dual numbers within a formalism reminiscent of operational umbral calculus.

Figures

Equations

$$z=x+ϵy,(x,y)\in \mathbb{R}$$

$${ϵ}^2=0\text{mbox}andϵ\ne 0.$$

$$z+w=x+u+ϵ(y+\nu )$$

$$z\cdot w=xu+ϵ(x\nu +yu)$$

$$\displaystyle z^{-1}=\frac{1}{x}\Big{(}1-\epsilon\,\frac{y}{x}\Big{)}\quad(x\not=0)$$

$$\displaystyle z^{n}=x^{n}\Big{(}1+n\epsilon\frac{y}{x}\Big{)}\quad(n\in\mathbb{Z}_{\geq 0}\,,\;x\neq 0)$$

$$f(x+y){\approx }_k\sum _{m=0}^k\frac{y^m}{m!}{f}^{(m)}(x).$$

$${\left(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{\pm }^{}\right)}_{i,j}:={\delta }_{j,i\pm 1},{\left(\text{tensor}\ast {[}_k]\hat{1}\right)}_{i,j}:={\delta }_{i,j},{\left(\text{tensor}\ast {[}_k]\hat{0}\right)}_{i,j}:=0,$$

$${\left(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{\pm }^{\ell }\right)}_{i,j}={\delta }_{j,i\pm \ell }\Rightarrow {\left(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{\pm }^{}\right)}^k=0.$$

$$f\left(x\text{tensor}\ast {[}_k]\hat{1}+y\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{\pm }^{}\right)=\sum _{m=0}^{k-1}\frac{1}{m!}{y}^m{f}^{(m)}(x)\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{\pm }^m.$$

$$\begin{gathered} f\left(x\text{tensor}\ast {[}_2]\hat{1}+y\text{tensor}\ast {[}_2]\widehat{ϵ}{}_{+}^{}\right)=f(x)\text{tensor}\ast {[}_2]\hat{1}+y{f}^{\prime }(x)\text{tensor}\ast {[}_2]\widehat{ϵ}{}_{+}^{} \\ =\left(\begin{array}{cc}f(x)& y{f}^{\prime }(x)\\ 0& f(x)\end{array}\right). \end{gathered}$$

$$⟨{}_k{e}_1|:=(1,0,\dots ,0),|{}_{k}1⟩:=\left(\begin{array}{c}1\\ ⋮\\ 1\end{array}\right),$$

$$⟨F(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{+}^{}){⟩}_k:=⟨{}_k{e}_1|F(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{+}^{})|{}_{k}1⟩.$$

$${⟨f\left(x\text{tensor}\ast {[}_k]\hat{1}+y\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{+}^{}\right)⟩}_k=\sum _{m=0}^{k-1}\frac{y^m}{m!}{f}^{(m)}(x).$$

$$\exp \left(\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{+}^{}x\right)=\sum _{r=0}^{k-1}\frac{x^r}{r!}\text{tensor}\ast {[}_k]\widehat{ϵ}{}_{+}^r,$$

$$e_n(x):=\sum _{r=0}^n\frac{x^r}{r!},$$

$$\left(\begin{array}{c}{e}_n(x)\\ e_{n-1}(x)\\ ⋮\\ e_1(x)\\ 1\end{array}\right):=\exp \left(\text{tensor}\ast {[}_{n+1}]\widehat{ϵ}{}_{+}^{}x\right)|{}_{n+1}1⟩.$$

$$e_n^{\prime }(x)=e_{n-1}(x).$$

$$F\,\rightsquigarrow_{k}\,G\quad:\Leftrightarrow\quad G=\left(F\big{|}_{\epsilon^{k+1}\to 0}\right)\big{|}_{\epsilon\to 1}$$

$$G=⟨F(\text{tensor}\ast {[}_{k+1}]\widehat{ϵ}{}_{+}^{}){⟩}_{k+1}.$$

$$\hat{z}\equiv \hat{z}(a,b):=a+ϵb.$$

$$f(x)=e^{-\alpha x^2+\hat{z}(a,b)x}.$$

$$f(x){⇝}_2{e}^{-\alpha x^2+ax}\left[1+bx+\frac{1}{2}(bx{)}^2\right],$$

$${\int }_{-\infty }^{+\infty }f(x)dx=\sqrt{\frac{\pi }{\alpha }}{e}^{\frac{\hat{z}(a,b{)}^2}{4\alpha }}=\sqrt{\frac{\pi }{\alpha }}{e}^{\frac{a^2}{4\alpha }+\frac{ab}{2\alpha }ϵ+\frac{b^2}{4\alpha }{ϵ}^2}.$$

$$\sum _{n=0}^{\infty }\frac{t^n}{n!}{H}_n(x,y)$$

$$H_n(x,y)$$

$$\sqrt{\frac{\pi }{\alpha }}{e}^{\frac{a^2}{4\alpha }+\frac{ab}{2\alpha }ϵ+\frac{b^2}{4\alpha }{ϵ}^2}=\sqrt{\frac{\pi }{\alpha }}{e}^{\frac{a^2}{4\alpha }}\sum _{m\ge 0}\frac{{ϵ}^m}{m!}{H}_m\left(\frac{ab}{2\alpha },\frac{b^2}{4\alpha }\right){⇝}_k\sqrt{\frac{\pi }{\alpha }}{e}^{\frac{a^2}{4\alpha }}{}_H{e}_k\left(\frac{ab}{2\alpha },\frac{b^2}{4\alpha }\right).$$

$${}_H{e}_k(x,y):=\sum _{r=0}^k\frac{1}{r!}{H}_r\left(x,y\right).$$

$$g(x):=e^{-\hat{z}(a,b)x^2}$$

$$\begin{gathered} {\int }_{-\infty }^{+\infty }g(x)dx=\sqrt{\frac{\pi }{\hat{z}}}=\sqrt{\frac{\pi }{a+ϵb}} \\ {⇝}_k\sqrt{\frac{\pi }{a}}\sum _{r=0}^k\left(\genfrac{}{}{0ex}{}{-\frac{1}{2}}{r}\right){\left(\frac{b}{a}\right)}^r. \end{gathered}$$

$${\partial }_x^n{e}^{\alpha x^2}=H_n(2\alpha x,\alpha )e^{\alpha x^2}$$

$${\partial }_x^n{e}^{-\hat{z}{x}^2}\stackrel{\text{eqref}eq:HermiteIdA}{=}{H}_n(-2\hat{z}x,-\hat{z})e^{-\hat{z}{x}^2}\stackrel{\text{eqref}eq:defHermite}{=}n!\sum _{r=0}^{\lfloor \frac{n}{2}\rfloor }\sum _{s\ge 0}\frac{(-1{)}^{n-r+s}{2}^{n-2r}{x}^{n-2(r-s)}}{(n-2r)!r!s!}{\hat{z}}^{n-r+s}.$$

Full text

\history\Title

Dual Numbers and Operational Umbral Methods \AuthorNicolas Behr *1,**, Giuseppe Dattoli 2, Ambra Lattanzi 2,3 and Silvia Licciardi 2 \AuthorNamesNicolas Behr, Giuseppe Dattoli, Ambra Lattanzi, Silvia Licciardi

\corresCorrespondence: nicolas.behr@irif.fr

1 Introduction

The dual numbers (DNs), introduced during the second half of the XIXth century Clifford (1871); Grünwald (1906); Segre (1912); Yaglom (2012); Klawitter (2014), can be viewed as abstract entities much like the ordinary complex numbers, and are defined as

[TABLE]

where the corresponding “imaginary” unit or dual number unit (DNU) $\epsilon$ is a nilpotent number,

[TABLE]

The dual numbers were originally introduced within the context of geometrical studies, and later exploited to deal with problems in pure and applied mechanics Kotelnikov (1985); Study (1903). For instance, it has been demonstrated in Rooney (1975); Hsia and Yang (1981); Martínez and Duffy (1993) how to formulate the equations of rigid body motion in terms of just three “dual” equations instead of their six “real” counterparts (thereby realizing an equivalence between spherical and spatial kinematics). More recently, as further discussed in the present paper, their importance has been recognized in numerical analysis to reduce round-off errors Fike and Alonso (2011). We believe that the use of dual numbers in the applied sciences is not as widespread as it could be, and that many new fields of research would benefit from their relevant introduction. An important domain in which they may bring significant novelties is that of the perturbative techniques in classical and quantum mechanics.

The main contribution of this paper consists in fixing the underlying algebraic rules of the dual numbers in the wider context of umbral and operational calculus. The paper is organized as follows: Section 2 delivers a basic mathematical introduction to dual numbers. Section 3 is devoted to the description of the computational procedure based upon dual numbers and umbral calculus. In Section 4, we will provide insight into how this powerful method can be applied to deal with problems arising in different contexts. For illustration, we will consider the Schrödinger and the heat equation, cornerstones in the respective fields of Physics. Section 5 provides a conclusion with further considerations for future works.

2 Higher order dual numbers

The DN algebraic rules Harkin and Harkin (2004); Özdemir (2018), summarized below, are a straightforward consequence of the previous identity (2) (with $z=x+\epsilon y$ and $w=u+\epsilon v$):

[TABLE]

While the addition operation is entirely analogous to the component-wise addition operation on two-dimensional vectors, the last three operations (product, inverse and power) characterize the distinguishing special algebraic properties of dual numbers (DNs). The multiplication is commutative, associative and distributive, thus the DNs form a two-dimensional associative and commutative algebra over the real numbers.

We will now extend this traditional dual number formalism as motivated by the following type of problem. Consider the Taylor expansion up to some order $k$ (denoted $\approx_{k}$) of an at least $k$-fold continuously differentiable function $f$ around a point $x$,

[TABLE]

Following the automatic differentiation paradigm Rall (1981); Rall and Corliss (1996); Fike (2016), since in practice the function $f$ will be implemented in some algorithmic from, it may be advantageous to formulate truncations such as (4) in terms of generalized (or higher order) dual numbers. To this end, let us introduce the families of square matrices $\tensor*[_{k}]{\hat{\epsilon}}{{}^{\,}_{\pm}}$, $\tensor*[_{k}]{\hat{1}}{}$ and $\tensor*[_{k}]{\hat{0}}{}$ with entries (for $i,j=1,\dotsc,k$)

[TABLE]

where $\delta_{i,j}$ denotes the Kronecker symbol. It is straightforward to verify that for all $k\geq 2$ and $\ell\geq 0$

[TABLE]

Then, under the assumptions (4), endowing the function $f(x)$ suitably with a component-wise action on square matrices, we find (for $k\geq 2$)

[TABLE]

For example, setting $k=2$, reproduces the well-known dual number identity Fike (2016)

[TABLE]

It may be verified that e.g. for the choice “+” in (5), the first row of the resulting matrices in (7) contains the terms of the Taylor expansion up to order $k-1$. More explicitly, introducing the auxiliary notations for the row vector $\left\langle{}_{k}e_{1}\right|$ and the column vector $\left|{}_{k}\mathbf{1}\right\rangle$ of length $k\geq 2$,

[TABLE]

let us define the order $k$ evaluation operation acting on some function $F(\tensor*[_{k}]{\hat{\epsilon}}{{}^{\,}_{+}})$ depending on a generalized dual number as

[TABLE]

We thus find that

[TABLE]

Recently, expansions such as (7) have received considerable interest in the field of numerical analysis Berland (2006). Referring to Fike (2016) for an overview, various alternative types of "numbers" have been studied for the purpose of finding optimized numerical schemes for computing $k$-th order derivatives of functions. For example, it has been demonstrated that the use of so-called hyper-dual numbers results in first and second derivative calculations that are exact, regardless of the step size Fike and Alonso (2011).

For later convenience, motivated by the identity (for $k\geq 2$)

[TABLE]

we may introduce the so-called truncated exponential polynomials Dattoli et al. (2003) $e_{n}(x)$ defined through the series

[TABLE]

which may be expressed in terms of generalized dual numbers as

[TABLE]

One may thus easily verify the property

[TABLE]

Having provided a matrix-based extension of ordinary to $k$-th order dual numbers of arbitrary order $k\geq 2$, we will now proceed to develop a computational procedure embedding dual numbers with other techniques inspired by the operational umbral formalism.

3 Umbral-type methods and Dual Numbers

Starting from this section, we will employ the notational simplification of writing $\epsilon$ for the dual number unit (DNU) $\tensor*[_{k}]{\hat{\epsilon}}{{}^{\,}_{\pm}}$ of generalized dual numbers (cf. Eq. (5)), making the order $k\geq 2$ of the DN explicit only via the analogue of the notation (10), and masking the matrix nature of $\tensor*[_{k}]{\hat{\epsilon}}{{}^{\,}_{\pm}}$. Thus for some function $F\equiv F(\epsilon)$, we write

[TABLE]

for the truncation of $F$ via setting $\epsilon^{k+1}=0$ and afterwards $\epsilon=1$. It is straightforward to verify that this formal definition may be implemented in terms of the matrix representations introduced in Section 2 via use of (10) as

[TABLE]

Consider then the dual complex parameter

[TABLE]

Following the principles of umbral calculus111Albeit the term umbral calculus has been introduced in the seminal papers by Roman and Rota Roman and Rota (1978), in the following we will make reference to the formalism developed in Licciardi (2018) which enriches the original formalism with the wealth of techniques derived from the operational calculus.Roman and Rota (1978); Licciardi (2018), we will treat the dual complex parameter $\hat{z}$ as an ordinary algebraic quantity in calculations of integrals, derivatives and other operations, delaying the evaluation of $\hat{z}$ via performing the operation $\,\rightsquigarrow_{k}\,$ to the very end of the computations. We will now illustrate the computational benefits of this approach via a number of examples.

3.1 Dual shifted Gaussians

We first consider a Gaussian-type function explicitly containing in its argument the dual complex parameter (18), whence the dual-shifted Gaussian function

[TABLE]

Assuming for instance third order dual numbers (i.e. $\epsilon^{3}=0$), we may write the above function in more conventional terms as

[TABLE]

which is easily recognized as the product of a shifted Gaussian with a second degree polynomial.

In full analogy to the umbral operational methods of  Licciardi (2018), it is then straightforward to calculate the following integral of the function $f$ of (19) via the standard Gaussian integral formula

[TABLE]

The term on the right has in fact a definite meaning, since the use of the generating function of the two variable Hermite polynomials Appel and De Feriet (1926)

[TABLE]

permits to cast the r.h.s. of (21) into the form

[TABLE]

Here, ${}_{H}e_{k}(x,y)$ denotes the Hermite-based truncated exponential polynomial Dattoli (1999, 2000); Dattoli et al. (2008) defined as

[TABLE]

3.2 Another form of dual Gaussian

Let us consider as a further example

[TABLE]

and the following infinite integral (for $\mbox{Re}(a)>0$)

[TABLE]

Here, by invoking the operation $\,\rightsquigarrow_{k}\,$, we obtain a finite series, thus obviating the need to impose any condition on the relevant convergence range.

3.3 Examples from symbolic calculus

The calculus of higher order dual numbers may be further refined via combining it with the wealth of techniques available from the theory of special functions and symbolic calculus as put forward in Licciardi (2018); Babusci et al. (2014); Babusci and Dattoli (2011); Górska et al. (2012); Dattoli et al. (2017); Babusci et al. (2017). Consider for illustration the following identity, known from the theory of two-variable Hermite polynomials Dattoli and Migliorati (2006),

[TABLE]

which allows to simplify the task of calculating successive derivatives of the dual Gaussian introduced in (25), such as in the computation

[TABLE]

Another interesting type of calculus concerns infinite integrals involving rational functions such as

[TABLE]

For example, the infinite integral

[TABLE]

may be easily transformed into truncated form in full analogy to the calculation summarized in (26).

3.4 Umbral image type techniques

Referring to Behr et al. (2019) for the precise technical details (compare also Dattoli and Migliorati (2006)), suffice it here to provide the following definition for the action of the formal integration operator $\hat{\mathbb{I}}$ on the formal variable $v$ (for $\alpha\in\mathbb{C}$):

[TABLE]

Then an interesting variant of the example presented in (28) may be obtained as

[TABLE]

In summary, the combination of the concept of higher order dual numbers with techniques from symbolic and umbral-image type calculus appears to offer a large potential in view of novel tools of computation. To corroborate this claim, we will now present some first high-level results in this direction.

4 Dual numbers and solution of heat- and Schrödinger-type equations

Before entering the main topic of this section, let us recall a few useful “operational rules”, starting with the Glaisher identity Dattoli et al. (2008); Crofton (1880)

[TABLE]

which can also be understood as the solution of the heat equation with a Gaussian as initial function. It will prove particularly useful in the following to note that according to the definition of the Hermite polynomials $H_{n}(x,y)$ as given in (22b), an alternative interpretation of (32) is provided in terms of the double lacunary exponential generating function $\mathcal{H}_{2,0}(\lambda;x,y)$ of the polynomials $H_{n}(x,y)$, where we employ notations as in Behr et al. (2018)

[TABLE]

By specializing eq. (32) to $\alpha=\hat{z}$ (with $\hat{z}=a+\epsilon b$ the dual complex parameter (18)), we obtain the operational identity

[TABLE]

Via the simple factorizations

[TABLE]

we may transform the identity (34) as

[TABLE]

By re-inserting the definition of the first double-lacunary EGF, using the Glaisher-identity (32) for the second one and finally truncating to order $k$, we eventually arrive at the compact result

[TABLE]

For example, by evaluating the above expression for second order dual numbers, one finds

[TABLE]

The above result may be interpreted as the solution of the heat-type equation

[TABLE]

An analogous problem has been addressed in Dattoli and Migliorati (2006) within the framework of a different method. The techniques we have envisaged may be further exploited to treat the paraxial propagation of the so-called flattened distributions, introduced in Gori (1994) to study the laser field evolution in optical cavities employing super-Gaussian mirrors Siegman (1986). These cavities shape beams whose transverse distribution is not reproduced by a simple Gaussian, but by a function exhibiting a quasi-constant flat-top, expressible through a function of the type

[TABLE]

The paraxial propagation of these beams has less obvious properties than, say, Laguerre or Hermite Gauss modes Siegman (1986). In order to overcome this drawback, Gori introduced the so-called flattened beams Gori (1994) which permit a fairly natural expansion in terms of Gauss Laguerre/Hermite modes, thus providing a straightforward solution to the corresponding paraxial wave equation.

Invoking our formalism as developed so far, we may approximate the aforementioned Gori beams in the form

[TABLE]

Here, $e_{m}(x)$ denotes the truncated exponential polynomials introduced in (13), and both parameters $\alpha$ and $m$ depend on $p$ (see Dattoli and Migliorati (2006) for further details). Recalling from (18) the definition $\hat{z}(a,b):=a+b\epsilon$ of the dual complex parameter, the r.h.s. of (41) may be equivalently expressed as

[TABLE]

whence as an instance of a dual Gaussian as described in Section 3.2. The problem of the relevant propagation can accordingly be reduced to that of an ordinary Gaussian mode, namely to the solution of the Schrödinger type equation

[TABLE]

Consequently, by invoking the operational identity (34), the paraxial evolution of a flattened beam may be expressed in the form

[TABLE]

which reproduces indeed the known solution of our problem (compare Dattoli and Migliorati (2006)).

In a forthcoming paper we will discuss this specific application in further detail by applying the method to the problem of designing super-Gaussian optical systems.

5 Weyl formula and modified Hermite polynomials

The wide flexibility of the method we propose is corroborated by the following further example, relevant to the use of operational ordering tools. Let us consider an evolution equation of the form

[TABLE]

The relevant procedure for combining differential calculus with the umbral formalism is described in  Babusci and Dattoli (2011). Following this approach, the solution of (45) can be expressed as

[TABLE]

In order to evaluate the solution (46) explicitly, we need to suitably “factorize” the exponential operator. This so-called disentanglement operation may be implemented via the Weyl formula Dattoli et al. (1997)

[TABLE]

which is applicable whenever the identities $[\hat{X},[\hat{X},\hat{Y}]]=[[\hat{X},\hat{Y}],\hat{Y}]]=0$ hold. Applying the Weyl formula (47) to our solution (46), we obtain

[TABLE]

Thus the solution at any desired truncation order $k$ may be obtained by invoking the dual number evaluation operation $\,\rightsquigarrow_{k}\,$ of (16).

As already mentioned above, the Weyl formula applies in the example presented because the algebraic structure of the argument of the exponential in (46) satisfies a special property: the commutators of the associated generators reduce to a constant after the first commutation bracket. A more interesting extension is given by the case in which the generators are embedded into a solvable Lie algebra. In this case, the combined use of the dual number formalism and of the Wei-Norman ordering method Wei and Norman (1963) leads to new and interesting results. They deserve a separate treatment that will be reported in a forthcoming paper.

As a final example, we define modified Hermite polynomials $H_{n}(x,\hat{z})$, whence ordinary two-variable Hermite polynomials $H_{n}(x,y)$ as introduced in (22b) evaluated at $y=\hat{z}$, with $\hat{z}\equiv\hat{z}(a,b)$ the dual complex parameter of (18),

[TABLE]

It is straightforward to verify that these modified polynomials inherit all the relevant properties from the polynomials $H_{n}(x,y)$, such as the recurrences

[TABLE]

and we find that they satisfy the second order differential equation

[TABLE]

The explicit form of these truncated polynomials is easily obtained. For example, by using third order dual numbers, which implies

[TABLE]

we find the explicit formula

[TABLE]

where we have invoked the well-known identity

[TABLE]

6 Final Comments

The method we have outlined in this paper offers many computational advantages to treat problems where truncated expansions (not necessarily of Taylor type) of functions are involved. At its core, the umbral formalism and the notion of higher order dual numbers allow to delay the explicit expansions to later stages in a given calculation, thus opening the possibility to exploit numerous efficient computation strategies from the theory of operational calculus and special functions.

The technique we have introduced in this paper is amenable for new applications in various different fields. We have presented herein the solution of parabolic equations in transport problems, and within such a context a fairly important example has been provided by treating the propagation of flattened beams Dattoli and Migliorati (2006); Gori (1994) in optics. For brevity, we have just outlined the procedure in terms of a 1-dimensional computation. The relevant extension to the 3-dimensional case does not require any particular conceptual effort, but only a consistent numerical implementation. In a forthcoming investigation, we will further extend the method and study its potential for treating perturbative problems in classical and quantum mechanics.

\authorcontributions

conceptualization, G.D.; methodology, N.B., G.D.; validation, N.B., G.D., S.L.; formal analysis, N.B., G.D., A.L., S.L.; writing – original draft preparation, N.B., G.D., A.L.; writing – review and editing, N.B., A.L., S.L.

\funding

The work of N.B. is supported by funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 753750. A.L. was supported by the NCN research project OPUS 12 no. UMO-2016/23/B/ST3/01714 and by the NAWA project: Program im. Iwanowskiej PPN/IWA/2018/1/00098. S.L. was supported by a Enea-Research Center Individual Fellowship.

Acknowledgements.

N.B. would like to thank the LPTMC (Paris 06) and ENEA Frascati for warm hospitality. \conflictsofinterestThe authors declare no conflicts of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results. \reftitleReferences