Advances in Mathematical Physics
Volume 2018 (2018), Article ID 9575626, 9 pages
https://doi.org/10.1155/2018/9575626

Research Article

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

Pawel Blasiak,¹ Gérard H. E. Duchamp,² and Karol A. Penson³

¹Institute of Nuclear Physics Polish Academy of Sciences, 31342 Kraków, Poland
²Université Paris 13, Sorbonne Paris Cité, LIPN, CNRS UMR 7030, 93430 Villetaneuse, France
³Université Pierre et Marie Curie (Paris 06), Sorbonne Universités, 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

Academic Editor: Andrei D. Mironov

Copyright © 2018 Pawel Blasiak et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We give a purely combinatorial proof of the Glaisher-Crofton identity which is derived from the analysis of discrete structures generated by the iterated action of the second derivative. The argument illustrates the utility of symbolic and generating function methodology of modern enumerative combinatorics. The paper is meant for nonspecialists as a gentle introduction to the field of graphical calculus and its applications in computational problems.

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

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

where functions and are specified by the following equations

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 . 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 [8] and the Glaisher-Crofton identity [9–11]

and the Glaisher-Crofton identity [9–11]

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 ; 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 which is defined as a denumerable collection of objects built of atoms represented by according to some well specified procedure. A typical combinatorial problem consists in enumeration of objects in according to the size which is usually the number of atoms. In other words, one seeks the sequence which counts the number of objects comprised of exactly atoms. This sequence can be encoded in a generating function

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 acting on some well-defined class which consists in

 “selecting in all possible ways a single atom of type in each element of and replacing it with an atom of a new type .”

In other words, one may think of the new class as formed of all structures taken from in which one of the atoms gets “repainted” into a new colour . Since each structure in built of atoms of type gives rise to new ones with atoms and a single , then the generating function enumerating objects in the new class is given by

where counts objects according to the number of ’s and ’s, respectively. This is substantiated in the standard transfer rule:

Now, let us consider the -th derivative acting on . Combinatorially it means

 “select in all possible ways an unordered collection of atoms of type and replace (repaint) them by atoms of type .”

Clearly, for each structure of size in we have possible choices, and hence the generating function of the new class evaluates to

In consequence we get the following transfer rule:

which gives the combinatorial interpretation of the -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 as a generating function of some combinatorial class . Then differentiation yields a new generating function which enumerates objects in the new class comprised of structures taken from in which some of the atoms were replaced with ’s (how many are replaced depends on the order of the derivative). Hence, the derivative of a function can be understood as a well-defined combinatorial transformation of the associated combinatorial class 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

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 applied to some function which can be considered as the generating function of some class of objects built of atoms . Then from our previous discussion the exponential corresponds to

 “selecting in all possible ways an arbitrary number of ’s and replacing them by ’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 either with atom (which makes no real effect) or with atom (which means the replacement). Hence we have the following combinatorial equivalence

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)):

which applies to any combinatorial class and its generating function . 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 acting on consists in

 “selecting in all possible ways an unordered pair of atoms and replacing the chosen ’s by atoms of type .”

We will call such (unordered) pair a doubleton. Clearly, for an object composed of atoms this can be done in ways, which agrees with the algebraic identity .

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

which consists in

 “selecting in all possible ways a set of unordered pairs (doubletons) of atoms and replacing them by atoms.”

Note that, as in (11), we deem the order in the sequence irrelevant by introducing factor 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

coincides with the number of possible ways of choosing a set of unordered pairs from the set of 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:

whose combinatorial meaning comes down to

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

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.,

For it specialises to the Hermite polynomial . More generally, we may also write

which stems from the fact that is an eigenvector of the derivative operator to eigenvalue . Again, for it is the exponential generating function of Hermite polynomials (cf. (5)). For the purpose at hand we will leave variable 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 as the exponential generating function of the labelled class of sets (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 carry integer labels , and additionally to each an atom (or weight) of type is attached; see Figure 1 on the left. We have the following translation rule (see (A.5)):

(Clearly, for given one can built one such set and its weight is .) Now, following the combinatorial description (17) action of on an individual set in 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 , while each atom forming the doubleton (arising from nontrivial action of the second derivative) carries the label . Yet another way of seeing the resulting class of objects is to understand them simply as a set of singletons and doubletons . See Figure 1 for illustration. Formally, one writes the following sequence of combinatorial equivalences:

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

Figure 1: On the left, there is an instance of a structure in class , i.e., a set of labelled atoms which are additionally weighted by an atom of type . In the middle, example of a structure arising from the action of on which is obtained by selecting a set of dubletons of atoms and replacing them by atoms of type . Selected pairs are depicted by wavy lines. On the right, decomposition of such a structure into a product of two sets comprising, respectively, singletons (untouched by the derivatives) and doubletons (selected by the derivatives).

As a consequence of this discussion we get a simple combinatorial model of Hermite polynomials. Namely, coefficients of count the number of possible ways of

 “selecting a subset of doubletons (each weighted by ) out of a set of distinguishable objects.”

It is exactly the coefficient of in (18) for . Another way to see it directly from our combinatorial description of the exponential of the second derivative is to interpret as the generating function of a single element set; i.e., we have . 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-Kampé de Fériet 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 [3] (see also the appendix).

Our goal is to calculate the explicit form of the expression (cf. the left hand side of (6) with and )

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

We know from the discussion of (17) that on the combinatorial level exponential consists in selecting in all possible ways unordered pairs of ’s (not necessarily attached to the same ) and replacing each in the chosen pairs by . Figure 2 in the middle illustrates a generic structure arising in this procedure. If we denote the resulting class of structures by , then we may write

where is the exponential generating function enumerating structures in . Note that this gives a precise combinatorial meaning to (22).

Figure 2: On the left, there is an instance of a structure in , i.e., a set of labelled atoms which have additional finer structure comprised of two atoms of type (say the left and the right one). In the middle, example of a structure arising from the action of on which is obtained by selecting a set of dubletons of atoms and replacing them by atoms of type . Selected pairs are depicted by wavy lines. On the right, decomposition of such a structure into a product of two sets comprising, respectively, closed () and open () chains.

Since we have reduced our goal to finding the exponential generating function one needs to come up with a systematic specification of structures in . For this purpose let us observe that doubletons, initially detached from one another in , are now tied together to form either open or closed chains. Then, we may group them together splitting each structure in 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)):

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

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

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

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

The class of open chains 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 with two additional blocks of types and attached at the ends (left and right, respectively). Here as well, each block occurs in two possible configurations arising from two possible choices of by the derivative in the same doubleton. See Figure 4 for illustration. This gives the following combinatorial specification:

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

Figure 4: If embedded in the line, each open chain in forms a sequence whose inner part (in the grey box) is built out of blocks of type with two blocks of type 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).

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

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 [16], where polynomial sequences can be considered as generated by the action of differential operators (see also monomiality principle [17]). Such description is attainable for a large family of Sheffer-type polynomials [18] (including binomial-type and Appell sequences [19]) 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 Kampé de Fériet 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 [3]. 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 which consists of a denumerable collection of objects built of the labelled atoms (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 , where which is conveniently encoded in the exponential generating function (e.g.f.)

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 class, which comprises a single element of size 1 and has e.g.f. , and the neutral class, which consists of a single element of size 0 and has e.g.f. . 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

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

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

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 . The corresponding generating functions are given by the following dictionary (assuming , i.e., ):

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].

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, Generatingfunctionology, 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,” Séminaire 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,” Séminaire Lotharingien de Combinatoire, vol. 61A, 2009.

V. Strehl, “Lacunary Laguerre series from a combinatorial perspective,” Séminaire Lotharingien de Combinatoire, vol. 76, 2017.