# Combinatorics of Second Derivative: Graphical Proof of Glaisher-Crofton Identity.

1. Introduction

Many computational problems involve action of complex expressions in derivatives on functions. A typical example is the exponential of a Hamiltonian acting on some initial condition which is a formal solution to the evolution equation. Applications of the latter range from classical heat and diffusion theory, financial mathematics, and economy to quantum field theory, hence practical interest in operational formulas enabling explicit evaluation of such expressions. Methods used to this effect usually involve operator and special function techniques, integral transforms, umbral calculus methods, etc. See comprehensive review of the subject [1, 2]. In this paper we develop another approach based on modern combinatorial methods of the analysis and enumeration of structures via generating functions [3-5].

The most known operational identities that involve the exponential of the first derivative are formulas for the shift and dilation operators

exp ([lambda] d/dx]) F (x) = F (x + [lambda])

exp ([lambda]x d/dx) F (x) = F ([e.sup.[lambda]]x) (1)

where F(x) is an arbitrary function. (Here, we leave subtle problems of convergence aside and consider F(x) as a formal power series in one variable x.) They are a special case of the general closed-form operational expression

exp [[lambda] (q(x) d/dx + v(x))] F(x) = g([lambda], x) x F(T ([lambda], x)), (2)

where functions T([lambda], x) and g([lambda], x) are specified by the following equations

[partial derivative]T ([lambda], x)/[partial derivative][lambda] = q(T ([lambda], x)), T(0, x) = x, (3)

[partial derivative]g ([lambda], x)/[partial derivative][lambda] = v(T ([lambda], x)) x g([lambda], x), g(0, x) = 1. (4)

See [1, 2, 6] for the proof based on operator techniques and [7, Sect. 6] for a recently developed combinatorial approach.

Note that formulas (1) and (2) are valid for any function F(x). However, this is a very unique situation which holds only for the exponential of an expression linear in the first derivative. For the second derivative the closed-form formulas are not known and the best that can be done is evaluation on specific functions. There are only a few explicit examples which include the formula for the exponential generating function of Hermite polynomials  and the Glaisher-Crofton identity [9-11]

[mathematical expression not reproducible] (5)

and the Glaisher-Crofton identity [9-11]

exp([alpha] [d.sup.2]/[dx.sup.2]) exp (-[x.sup.2]) = 1/[square root of (1 + 4[alpha])] exp (-[x.sup.2]/1 + 4[alpha]). (6)

Formulas of this type are usually derived using integral representations in the complex domain. (For example, to derive (6) one may use the integral representation exp(-[x.sup.2]) = [[integral].sup.+[infinity].sub.-[infinity]] exp(-[[xi].sup.2] + 2i[xi]x)d[xi]; see [10, 11].) However, in this paper we demonstrate that it is also possible to prove these identities on the basic algebraic level by analysing combinatorial structures generated by the iterated action of the second derivative. In this approach functions are treated as generating functions enumerating simple combinatorial objects (like sets of subsets, cycles, sequences, etc.) and consequently expressions in the derivatives transform objects in the initial class into richer structures whose generating functions can be quickly identified with the methods of symbolic combinatorics [3,4]. We remark that methodology exposed in this note treads in the steps of combinatorial approach to algebraic identities developed by D. Foata in [12, 13].

Our goal in this paper is to develop and promote general combinatorial methodology for solving computational problems which, in many cases, provides better insight into algebraic and analytic manipulations. We illustrate this approach by explaining combinatorial meaning of the exponential of the second derivative and use this interpretation to derive equations (5) and (6).

2. Combinatorics of Derivatives

Action of derivatives on a function can be seen as a transformation of combinatorial structures. This viewpoint comes from interpreting the function as a generating function enumerating objects in some combinatorial class. In the following we formalise this intuition by describing the relevant constructions and develop a broader picture which includes higher derivatives and their exponentials. This framework will be illustrated by a simple combinatorial proof and interpretation of Taylor's formula in Section 2.1 and equations (5) and (6) in Sections 2.2 and 3.

2.1. Generating Functions, First Derivative, and Taylor's Formula. Let us consider a combinatorial class F which is defined as a denumerable collection of objects built of atoms represented by X according to some well specified procedure. A typical combinatorial problem consists in enumeration of objects in F according to the size [absolute value of (x)] : F [right arrow] N which is usually the number of atoms. In other words, one seeks the sequence [f.sub.n] = #{[phi] [member of] F : [absolute value of ([phi])] = n} which counts the number of objects comprised of exactly n atoms. This sequence can be encoded in a generating function

[mathematical expression not reproducible] (7)

which is a convenient tool for enumeration of complex structures via the so-called transfer rules. The latter translate combinatorial constructions into algebraic manipulations of the corresponding generating functions (see [3-5] for a comprehensive treatment of the subject and the appendix for a quick extract of a few transfer rules used in this paper). In the following we will briefly review combinatorics of the first derivative and recall a simple combinatorial interpretation of Taylor's formula.

Here we will be concerned with the derivative operation [yD.sub.x] acting on some well-defined class F which consists in

"selecting in all possible ways a single atom of type X in each element of F and replacing it with an atom of a new type Y."

In other words, one may think of the new class [yD.sub.x]F as formed of all structures taken from F in which one of the atoms X gets "repainted" into a new colour Y. Since each structure in F built of n atoms of type X gives rise to n new ones with (n-1) atoms X and a single Y, then the generating function enumerating objects in the new class yDxF is given by

[mathematical expression not reproducible] (8)

where [g.sub.n,l] counts objects according to the number of X's and Y's, respectively. This is substantiated in the standard transfer rule:

Selecting singleton:

G = [yD.sub.x]F [??]

G(x, y) = y d/dx F(x). (9)

Now, let us consider the k-th derivative (1/k!)[([yD.sub.x]).sup.k] acting on F. Combinatorially it means

"select in all possible ways an unordered collection of k atoms of type X and replace (repaint) them by atoms of type Y."

Clearly, for each structure of size n in F we have ([mathematical expression not reproducible]) possible choices, and hence the generating function of the new class (1/k!)[([yD.sub.x]).sup.k] F evaluates to

[mathematical expression not reproducible] (10)

In consequence we get the following transfer rule:

Selecting k - subset:

G = 1/k! [([yD.sub.x]).sup.k] F [??]

G(x, y) = 1/k! [y.sup.k] [d.sup.k]/[dx.sup.k] F(x), (11)

which gives the combinatorial interpretation of the k-th derivative on the level of combinatorial structures. This brings an interesting perspective on the derivative of a function which we will develop throughout the paper. Namely, one may think of a function F(x) as a generating function of some combinatorial class F. Then differentiation yields a new generating function which enumerates objects in the new class comprised of structures taken from F in which some of the atoms X were replaced with Y's (how many are replaced depends on the order of the derivative). Hence, the derivative of a function F(x) can be understood as a well-defined combinatorial transformation of the associated combinatorial class F in a sense that on the level of generating functions it corresponds to simple differentiation (cf. (9) and (11)).

For illustration of this viewpoint let us recall the usual Taylor's formula

[mathematical expression not reproducible] (12)

Surprisingly, it admits a transparent combinatorial interpretation (see [3, Note III.31] or [7, Note 3]). To see this we observe that the l.h.s. is the sum of derivatives exp([yD.sub.x]) [equivalent to] [[summation].sub.k[greater than or equal to]0](l/k!)[([yD.sub.x]).sup.k] applied to some function F(x) which can be considered as the generating function of some class of objects built of atoms X. Then from our previous discussion the exponential exp([yD.sub.x]) corresponds to

"selecting in all possible ways an arbitrary number of X's and replacing them by Y's."

(since the sum contains derivatives of arbitrary order we may choose subsets of arbitrary cardinality). On the other hand, this is the same as substituting each atom X either with atom X (which makes no real effect) or with atom Y (which means the replacement). Hence we have the following combinatorial equivalence

exp ([yD.sub.x]) F = F [??] (X + Y), (13)

which on the level of generating functions, by virtue of (11) and the transfer rule for substitutions (see the appendix, (A.4)), directly translates into (12). Hence from the combinatorial point of view Taylor's formula is a simple manifestation of the following transfer rule (cf. (13)):

Selecting arbitrary subset:

G = exp ([yD.sub.x]) F [??]

G (x, y) = F (x + y), (14)

which applies to any combinatorial class F and its generating function F(x). This is a typical example of combinatorial methodology which draws on the fact that in many cases the same combinatorial structure allows different specifications.

2.2. Second Derivative and Hermite Polynomials. Combinatorial interpretation of the second derivative can be developed along similar lines. From the above we know that the 2-nd derivative (1/2)[([yD.sub.x]).sup.2] acting on F consists in

"selecting in all possible ways an unordered pair of X atoms and replacing the chosen X's by atoms of type Y."

We will call such (unordered) pair a doubleton. Clearly, for an object composed of n atoms this can be done in n(n - l)/2 ways, which agrees with the algebraic identity (1/2)[([yD.sub.x]).sup.2][x.sup.n] = (n(n - 1)/2)[x.sup.n-2][y.sup.2].

More generally, by iterating k times one picks out a sequence of doubletons in the original structure. Hence, we define the following construction:

Selecting k - subset of doubletons:

G = 1/k! [(1/2 [([yD.sub.x]).sup.2]).sup.k] F [??]

1/[2.sup.k]k! [y.sup.2k] [d.sup.2k]/[dx.sup.2k] F(x), (15)

which consists in

"selecting in all possible ways a set of k unordered pairs (doubletons) of X atoms and replacing them by Y atoms."

Note that, as in (11), we deem the order in the sequence irrelevant by introducing factor 1/k! in front of the iterated derivative (hence the "set" and not the "sequence" in the description). For a quick check of this specification we observe that the coefficient on the r.h.s. of the identity

[mathematical expression not reproducible] (16)

coincides with the number of possible ways of choosing a set of k unordered pairs from the set of n objects, and hence by linearity we establish correctness of the description and transfer rule (15).

Now, we are in position to give interpretation of the exponential of second derivative:

Selecting arbitrary subset of doubletons:

[mathematical expression not reproducible] (17)

whose combinatorial meaning comes down to

"selecting in all possible ways an arbitrary subset of (unordered) pairs of X atoms and replacing each chosen X by atom of type Y."

This is the sum of the derivative operations of the type (15) which on the level of generating functions is the sum of the derivatives. Unfortunately it does not come close to any neat expression like (12) or (14). Indeed, Taylor's formula does not generalise in a straightforward manner. Innocuous as it may seem, selecting pairs instead of singletons introduces considerable complexity into the picture and requires careful analysis which quickly gets intractable. However, in particular cases of simple combinatorial structures (and their generating functions) it is possible to carry all the calculus through. An example that we will consider in detail is the Glaisher-Crofton formula (6) which evaluates action of the exponential of the second derivative on the gaussian. Before we proceed to this result, discussed in Section 3, we will illustrate our combinatorial methodology on a simpler case of the action on monomials and provide a link with combinatorial model of Hermite polynomials (cf. [4, 12, 13]).

Let us start with the explicit expression which is obtained from expanding the exponential and differentiating the monomial, i.e.,

[mathematical expression not reproducible] (18)

For y = i it specialises to the Hermite polynomial [He.sub.n](x). More generally, we may also write

exp (1/2 [([yD.sub.x]).sup.2]) exp (xt) = exp (1/2 [(yt).sup.2]) x exp (xt) = exp (xt + 1/2 [(yt).sup.2]), (19)

which stems from the fact that exp(xt) is an eigenvector of the derivative operator [D.sub.x] to eigenvalue t. Again, for y = i it is the exponential generating function of Hermite polynomials (cf. (5)). For the purpose at hand we will leave variable y unspecified so as to deal with positive integers only (cf. coefficients in (18)). (It is a typical combinatorial trick to introduce additional labels (or weights) which often allows getting rid of negative or noninteger factors entering multiplicatively in the expressions. Then enumeration of structures proceeds also with respect to this additional label (or weight) which, if needed, can be specified to the required value at the end.) Our aim is to understand these formulas in terms of enumeration of structures.

In order to use combinatorial description of (17) we interpret exp(xt) as the exponential generating function of the labelled class of sets Set(XT) (see [3, Sect. II] for a precise definition and discussion of labelled classes and their relation with exponential generating functions). It is comprised of sets whose atoms T carry integer labels 1, 2,3, ..., and additionally to each T an atom (or weight) of type X is attached; see Figure 1 on the left. We have the following translation rule (see (A.5)):

[mathematical expression not reproducible] (20)

(Clearly, for given n one can built one such set and its weight is [x.sup.n]) Now, following the combinatorial description (17) action of exp((1/2)[([yD.sub.x]).sup.2]) on an individual set in Set(XT) consists in selecting in all possible ways a subset of doubletons. This amounts to splitting of each original set into products of two subsets: one comprising singletons and the other doubletons. Additionally, these subsets differ in that the atoms in the singletons (untouched by the derivatives) carry the weight X, while each atom forming the doubleton (arising from nontrivial action of the second derivative) carries the label Y. Yet another way of seeing the resulting class of objects is to understand them simply as a set of singletons XT and doubletons (1/2)[(YT).sup.2]. See Figure 1 for illustration. Formally, one writes the following sequence of combinatorial equivalences:

[mathematical expression not reproducible] (21)

which on the level of generating functions readily transform (cf. the appendix) into a sequence of algebraic equalities providing a combinatorial proof of (19).

As a consequence of this discussion we get a simple combinatorial model of Hermite polynomials. Namely, coefficients of [He.sub.n](x) = [[summation].sup.n.sub.k=0][h.sub.n,k][x.sup.k] count the number of possible ways of

"selecting a subset of k doubletons (each weighted by -1) out of a set of n distinguishable objects."

It is exactly the coefficient of [x.sup.n-2k] in (18) for y = i. Another way to see it directly from our combinatorial description of the exponential of the second derivative is to interpret [x.sup.n] as the generating function of a single n element set; i.e., we have F = [Set.sub.n](X) [??] F(x) = [x.sup.n]. Then combinatorial model in terms of the choice of doubletons is a simple consequence of the specification of (17). We note that this model is a rephrasing of the interpretation of Hermite polynomials in terms of weighted involutions [4, Sect. 2.3.]. Moreover, it can be straightforwardly extended to provide a combinatorial interpretation of a larger class of multivariate Hermite-Kampe de Feriet polynomials [1, 2].

3. Proof of the Glaisher-Crofton Identity

Here we prove the identity (6) by a purely combinatorial argument by analysing structures generated by the second derivative discussed above. We will proceed in a step by step manner explaining the details of combinatorial constructions and structures that appear along the way. Although most of them are standard in combinatorial community we take a rather explicit and methodological route that may be of help for an unaccustomed reader. In the following we adopt the standard notation from the book  (see also the appendix).

Our goal is to calculate the explicit form of the expression (cf. the left hand side of (6) with [alpha] = [y.sup.2]/2 and t = -1)

[mathematical expression not reproducible] (22)

Following our combinatorial strategy we will treat exp([x.sup.2]t) as an exponential generating function enumerating some combinatorial objects. Let us define them as a labelled class of sets Set([X.sup.2]T). This class is comprised of sets of labelled atoms T (i.e., each atom carries the integer label) which are weighted with two atoms of type X. We will depict them as sets of doubletons of unlabelled atoms X such that each doubleton carries the labelled marker T. See Figure 2 on the left for illustration. Clearly, we have the following transfer rule (cf. (A.5)):

[mathematical expression not reproducible] (23)

We know from the discussion of (17) that on the combinatorial level exponential exp((1/2)[([yD.sub.x]).sup.2]) consists in selecting in all possible ways unordered pairs of X's (not necessarily attached to the same T) and replacing each X in the chosen pairs by Y. Figure 2 in the middle illustrates a generic structure arising in this procedure. If we denote the resulting class of structures by R, then we may write

R = exp [(1/2([yD.sub.x]).sup.2]) C [??]

exp (1/2 [([yD.sub.x]).sup.2]) exp([x.sup.2]t) = R(x, y, t), (24)

where R(x, y, t) = [[summation].sub.k,l,n[greater than or equal to]0] [R.sub.k,l,n][x.sup.k][y.sup.l] ([t.sup.n]/n!) is the exponential generating function enumerating structures in R. Note that this gives a precise combinatorial meaning to (22).

Since we have reduced our goal to finding the exponential generating function R(x, y, z) one needs to come up with a systematic specification of structures in R. For this purpose let us observe that doubletons, initially detached from one another in C, are now tied together to form either open or closed chains. Then, we may group them together splitting each structure in R into a product of two sets: one containing only closed and the second open chains. This entails the following combinatorial decomposition and its translation to the exponential generating functions (cf. (A.3) and (A.5)):

R = Set (A) x Set (B) [??]

R (x, y, t) = exp (A (x, y, t)) x exp (B (x, y, t)), (25)

where A is the class of open chains and B is the class of closed chains; see Figure 2 on the right. Hence the problem comes down to finding both exponential generating functions A(x, y, t) and B(x, y, t) enumerating, respectively, objects of types A and B.

Let us start with the class of closed chains A. Its elements embedded in the plane can be seen as cycles whose building blocks have a finer structure of type [Y.sup.2]T which occurs in two possible configurations arising from two possible choices of X by the derivative in the same doubleton; see Figure 3 for pictorial explanation. Therefore the whole class is specified as follows:

A = CYC (2[Y.sup.2]T), (26)

and by the standard transfer rule for labelled classes, (A.6), we get

A(x, y, t) = log 1/1 - 2[y.sup.2]t, (27)

The class of open chains B can be described in a similar manner. First we observe that each such chain can be embedded in the line in two possible ways. Then it forms a sequence which can be decomposed into the inner part which is a sequence of blocks of type YTY with two additional blocks of types XYT and YTX attached at the ends (left and right, respectively). Here as well, each block occurs in two possible configurations arising from two possible choices of X by the derivative in the same doubleton. See Figure 4 for illustration. This gives the following combinatorial specification:

B = XTX + 1/2 (2XYT) x Seq (2YTY) x (2YTX). (28)

Note that the coefficient 1/2 stems from the double counting due to embedding in a line and the term XTX makes up for a single structure left out by the above description. Having specified B we obtain exponential generating function by means of the standard transfer rules, (A.7), which yield

B (x, y, t) = [x.sup.2]t + 1/2xyt x 1/1 - 2[y.sup.2]t x 2xyt = [x.sup.2]t/1 - 2[y.sup.2]t. (29)

Now, by substituting (27) and (29) to (25), we get

R(x, y, t) = 1/1 - 2[y.sup.2]t exp ([x.sup.2]t/1 - 2[y.sup.2]t). (24)

This completes the proof of the Glaisher-Crofton identity, (6), which is readily obtained from (22). In conclusion, let us remark that we benefit from the proof by a deeper combinatorial insight into the nature of both factors on the r.h.s. of (6). The latter can be interpreted as exponential generating functions enumerating sets of closed and open chains, respectively, formed by derivatives acting on the gaussian.

4. Discussion and Outlook

Many computational problems require keeping track and skilful rearrangement of terms involved in algebraic expressions. It often comes down to the analysis of their structural properties and counting terms grouped with respect to some relevant characteristics. This is a natural domain of application for modern combinatorics which has developed a large array of tools for systematic treatment of such problems. In this paper, we have considered a few examples where it can be effectively used for evaluation of the action of the exponential in the derivatives on a function. Fundamental in this approach is treatment of the function as a generating function enumerating some simple combinatorial objects. This shift in perspective allows interpreting the derivatives (and their exponentials) as combinatorial constructors which produce a new class of objects which often can be enumerated with combinatorial flair. We have illustrated this approach by showing simple combinatorial proofs and interpretation of Taylor's formula, connecting exponential in second derivative with a model of Hermite polynomials, and deriving the Glaisher-Crofton identity. It is worth emphasising that the proofs are purely combinatorial and thus do not require any arguments involving integral representations or analyticity. We note that our exposition builds up on the seminal papers of D. Foata [12-14] who gave combinatorial insight into the related Mehler formula. In this paper, the focus is shifted towards graphical calculus based on symbolic methods in enumeration of combinatorial structures in [3, 4].

Crucial to our development was combinatorial understanding of the exponential of the second derivative which consists in selecting a collection of unordered pairs in a structure it acts on. It can be also seen as superposition of the doubleton structure on the other one which is connected with the Hadamard product of generating functions considered in various combinatorial contexts (see, e.g., [3, 4, 15]). We should also mention a natural link with a rich framework of umbral calculus , where polynomial sequences can be considered as generated by the action of differential operators (see also monomiality principle ). Such description is attainable for a large family of Sheffer-type polynomials  (including binomial-type and Appell sequences ) and therefore admits combinatorial interpretation along the lines considered in the present paper. This theme will be the subject of subsequent publication.

We observe that our discussion is not limited only to the first and second derivatives. It can be straightforwardly extended to derivatives of higher order which correspond to selecting subsets of higher cardinality. Moreover, one can generalise this framework to partial derivatives in several variables and multivariate polynomials (e.g., Hermite or Kampe de Feriet polynomials [1, 2]) by interpreting them as enumerating combinatorial structures built of atoms of several kinds.

Finally, let us remark that combinatorial approach to derivatives also provides an interesting insights into operator identities. One example is a systematic treatment of the normal ordering problem [7, 20, 21]. Clearly, majority of operator identities admit combinatorial interpretation as they typically arise from algebraic manipulation of discrete structures [4, 22, 23]. As such it opens the whole field of application for combinatorial approach. In this paper we have illustrated this point only on a few examples which can be seen as instances of a broad class of operator identities amenable to combinatorial methodology (cf. Sack identity, Baker-Campbell-Hausdorff formula, Rodrigues-type formulas, Crofton identities, etc. [1, 2]).

Appendix

A. Combinatorial Constructions

Our primary reference for combinatorial analysis is the standard book Analytic Combinatorics by Ph. Flajolet and R. Sedgewick . Here, we briefly recall basic terminology and a few standard translation rules for labelled constructions used in this paper.

Suppose we are given a combinatorial class C which consists of a denumerable collection of objects built of the labelled atoms T (see [3, Ch. II] for precise definition of the labelled class). Usually size of an object is the number of atoms it is built of, and a typical problem is to count the number of structures of a given size. In other words, one seeks the sequence [c.sub.n] = #[C.sub.n], where [C.sub.n] = {[phi] [member of] C : [absolute value of ([phi])] = n} which is conveniently encoded in the exponential generating function (e.g.f.)

[mathematical expression not reproducible] (A.1)

Let us remark that the reason for the use of exponential generating functions, rather than ordinary generating functions (o.g.f.), is simplicity of transfer rules in the domain of labelled classes. (Ordinary generating functions are typically used for enumeration of unlabelled structures, cf. [3, Ch. I.].)

The point of combinatorial analysis of structures is construction of complex classes from simpler ones. The initial building blocks include the atomic classT, which comprises a single element of size 1 and has e.g.f. C(t) = t, and the neutral classE, which consists of a single element of size 0 and has e.g.f. C(t) = 1. Then complex structures are built by well defined set of theoretical constructions which provide a precise specification of the class. Remarkably, these constructions can be translated into algebraic equations for the corresponding generating functions which solve the enumeration problem. Below, we give a short list of such constructions and translation rules that we exploit in this paper.

The most basic one is the disjoint union, henceforth denoted by "+", which corresponds to

C = A + B [??]

C(t) = A(t) + B(t). (A.2)

Another one is the labelled product, denoted by "x", which forms a cartesian product of objects and relabels the atoms in order-consistent manner. We have the following translation rule:

C = A x B [??]

C(t) = A(t) x B(t). (A.3)

If objects of one structure B are substituted into atoms of another structure A and relabelled in the order-consistent way, then the e.g.f. of such constructed class is given by (assuming [B.sub.0] = 0, i.e., B(0) = 0)

C = A [??] B [??]

C(t) = A(B(t)). (A.4)

It is then possible to form the class of all (labelled) sequences, sets, and cycles (respectively denoted by Set, Seq, and Cyc) built from objects in A. The corresponding generating functions are given by the following dictionary (assuming [A.sub.0] = 0, i.e., A(0) = 0):

C = Set (A) [??]

C(t) = exp (A(t)), (A.5)

C = CYC (A) [??]

C (t) = log 1/1 - A(t), (A.6)

C = Seq (A) [??]

C(t) = 1/1 - A(t). (A.7)

This is a nonexhaustive selection of possible constructions which is used in the present paper. For a comprehensive survey of the methods of combinatorial enumeration via generating functions, we refer to the classic books on this subject [3-5].

https://doi.org/10.1155/2018/9575626

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper is dedicated to the memory of the late Philippe Flajolet and Allan I. Solomon. We would like to thank Andrzej Horzela for helpful discussions.

References

 G. Dattoli, P. L. Ottaviani, A. Torre, and L. Vazquez, "Evolution operator equations: integration with algebraic and finite-difference methods. Applications to physical problems in classical and quantum mechanics and quantum field theory," La Rivista del Nuovo Cimento della Societa Italiana di Fisica, vol. 20, no. 2, pp. 1-133, 1997.

 P. E. Ricci and I. Tavkhelidze, "An introduction to operational techniques and special polynomials," Journal of Mathematical Sciences, vol. 157, no. 1, pp. 161-189, 2009.

 P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009.

 F. Bergeron, G. Labelle, and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge University Press, 1998.

 H. S. Wilf, Generating functionology, A K Peters, Ltd., Wellesley, MA, USA, 3rd edition, 2006.

 P. Blasiak, A. Horzela, K. A. Penson, G. H. E. Duchamp, and A. I. Solomon, "Boson normal ordering via substitutions and Sheffer-type polynomials," Physics Letters A, vol. 338, no. 2, pp. 108-116, 2005.

 P. Blasiak and P. Flajolet, "Combinatorial models of creation-annihilation," Seminaire Lotharingien de Combinatoire, vol. 65, 2011.

 H. M. Srivastava and H. L. Manocha, A Treatise on Generating Functions, Ellis-Horwood Limited, 1984.

 M. W. Crofton, "Theorems in the calculus of operations," Quarterly Journal of Mathematics, vol. 16, pp. 323-352, 1879.

 G. Dattoli, S. Khan, and P. E. Ricci, "On Crofton-Glaisher type relations and derivation of generating functions for Hermite polynomials including the multi-index case," Integral Transforms and Special Functions, vol. 19, no. 1-2, pp. 1-9, 2008.

 G. Dattoli, "Generalized polynomials, operational identities and their applications," Journal of Computational and Applied Mathematics, vol. 118, no. 1-2, pp. 111-123, 2000.

 D. Foata, "A combinatorial proof of the Mehler formula," Journal of Combinatorial Theory, Series A, vol. 24, no. 3, pp. 367-376, 1978.

 D. Foata and A. M. Garsia, "A combinatorial proof to the Mehler formulas for Hermite polynomials," in Relations between Combinatorics and other parts of Mathematics, vol. 34 of AMS, pp. 136-179, 1979, Proc. Symposia in Pure Mathematics, Columbus 1978.

 D. Foata, "Some Hermite polynomial identities and their combinatorics," Advances in Applied Mathematics, vol. 2, no. 3, pp. 250-259, 1981.

 P. Blasiak, K. A. Penson, A. I. Solomon, A. Horzela, and G. H. Duchamp, "Some useful combinatorial formulas for bosonic operators," Journal of Mathematical Physics, vol. 46, 2005.

 S. Roman, The Umbral Calculus, Academic Press, 1984.

 P. Blasiak, G. Dattoli, A. Horzela, and K. A. Penson, "Representations of monomiality principle with Sheffer-type polynomials and boson normal ordering," Physics Letters A, vol. 352, no. 1-2, pp. 7-12, 2006.

 S. Khan, M. W. Al-Saad, and G. Yasmin, "Some properties of Hermite-based Sheffer polynomials," Applied Mathematics and Computation, vol. 217, no. 5, pp. 2169-2183, 2010.

 S. Khan, G. Yasmin, R. Khan, and N. A. Makboul Hassan, "Hermite-based Appell polynomials: properties and applications," Journal of Mathematical Analysis and Applications, vol. 351, no. 2, pp. 756-764, 2009.

 P. Blasiak, A. Horzela, K. A. Penson, A. I. Solomon, and G. H. Duchamp, "Combinatorics and Boson normal ordering: a gentle introduction," American Journal of Physics, vol. 75, no. 7, pp. 639-646, 2007.

 P. Blasiak, "On urn models, non-commutativity and operator normal forms," Physics Letters A, vol. 374, no. 48, pp. 4808-4813, 2010.

 G. Labelle and C. Lamathe, "General combinatorial differential operators," Seminaire Lotharingien de Combinatoire, vol. 61A, 2009.

 V. Strehl, "Lacunary Laguerre series from a combinatorial perspective," Seminaire Lotharingien de Combinatoire, vol. 76, 2017.

Pawel Blasiak (iD), (1) Gerard H. E. Duchamp, (2) and Karol A. Penson (3)

(1) Institute of Nuclear Physics Polish Academy of Sciences, 31342 Krakow, Poland

(2) Universite Paris 13, Sorbonne Paris Cite, LIPN, CNRS UMR 7030, 93430 Villetaneuse, France

(3) Universite Pierre et Marie Curie (Paris 06), Sorbonne Universites, LPTMC, CNRS UMR 7600, 75252 Paris Cedex 05, France

Correspondence should be addressed to Pawel Blasiak; pawel.blasiak@ifj.edu.pl

Received 30 November 2017; Accepted 26 September 2018; Published 22 October 2018

Caption: Figure 1: On the left, there is an instance of a structure in class F = Set(XT), i.e., a set of labelled atoms T which are additionally weighted by an atom of type X. In the middle, example of a structure arising from the action of exp((1/2)[([yD.sub.x]).sup.2]) on F which is obtained by selecting a set of dubletons of X atoms and replacing them by atoms of type Y. Selected pairs are depicted by wavy lines. On the right, decomposition of such a structure into a product of two sets comprising, respectively, singletons XT (untouched by the derivatives) and doubletons (1/2)[(YT).sup.2] (selected by the derivatives).

Caption: Figure 2: On the left, there is an instance of a structure in C = Set([X.sup.2]T), i.e., a set of labelled atoms T which have additional finer structure comprised of two atoms of type X (say the left and the right one). In the middle, example of a structure arising from the action of exp((1/2)[([yD.sub.x]).sup.2]) on C which is obtained by selecting a set of dubletons of X atoms and replacing them by atoms of type Y. Selected pairs are depicted by wavy lines. On the right, decomposition of such a structure into a product of two sets comprising, respectively, closed (A) and open (B) chains.

Caption: Figure 3: A generic closed chain in class A embedded in the plane can be seen as a cycle built out of two kinds of blocks, each of type [Y.sup.2]T, arising from different choices of X's attached to the same T by the derivative on either the left or right (hence possible crossings of the wavy lines).

Caption: Figure 4: If embedded in the line, each open chain in B forms a sequence whose inner part (in the grey box) is built out of blocks of type [Y.sup.2]T with two blocks of type XYT attached at the ends. Here again, each block can have two configurations due to different choices made by the derivatives (possible crossings of the wavy lines).
Title Annotation: Printer friendly Cite/link Email Feedback Research Article Blasiak, Pawel; Duchamp, Gerard H.E.; Penson, Karol A. Advances in Mathematical Physics Report 1USA Jan 1, 2018 5889 Completeness and Nonclassicality of Coherent States for Generalized Oscillator Algebras. Abundant Lump-Type Solutions and Interaction Solutions for a Nonlinear (3+1) Dimensional Model. Combinatorics Derivatives (Mathematics)