A. M. Mathai Centre for Mathematical and Statistical Sciences: A Brief History of the Centre and Prof. Dr. A. M. Mathai’s Research and Education Programs at the Occasion of His 85th Anniversary

Hans J. Haubold

Office for Outer Space Affairs, United Nations Vienna International Centre, Vienna, Austria.

DOI: 10.4236/ce.2020.113028                

Abstract

A brief history of the Centre for Mathematical and Statistical Sciences, Kerala, India, is given and an overview of Mathai’s research and education programs in the following topics is outlined: Fractional Calculus; Functions of Matrix Argument—M-transforms, M-convolutions; Krätzel integrals; Pathway Models; Geometrical Probabilities; Astrophysics—reaction rate theory, solar neutrinos; Special Functions—G and H-functions; Multivariate Analysis; Algorithms for Non-linear Least Squares; Characterizations—characterizations of densities, information measure, axiomatic definitions, pseudo analytic functions of matrix argument and characterization of the normal probability law; Mathai’s Entropy—entropy optimization; Analysis of Variance; Population Problems and Social Sciences; Quadratic and Bilinear Forms; Linear Algebra; Probability and Statistics.

Keywords

Fractional Calculus, Functions of Matrix Argument, Special Functions of Mathematical Physics, Geometrical Probabilities, Multivariate Analysis, Characterization, Algorithms, Quadratic and Bilinear Forms, Entropy, Reaction, Diffusion, Solar Neutrinos

The Centre for Mathematical Sciences (CMS) was established in 1977 and registered in Trivandrum, Kerala, India, as a non-profit scientific society under the Travancore-Cochin Literary, Scientific and Charitable Societies Registration Act XII of 1955, and as a research and education centre covering all aspects of mathematics, statistics, mathematical physics, computer and information sciences. Applicants for incorporation were K.T. Chandy, A. Abraham, J. Zachariah, V.I. Subramoniam, K. Sankara Rao, Y. Sitaraman, and A. George, with Professor George taking the lead.

Research and Education in Mathematical and Statistical Sciences

Initial activities of research concentrated on population studies and applied statistics, the areas of interest of Professor George. Research and project activities started at CMS and at one time there were 35 project staff. CMS has executed a large number of research and education projects for central and state governmental agencies.

By the end of 1984 Professor George passed away and then governing council of CMS requested A. M. Mathai to take over and build it up into an international centre of excellence. Professor Mathai agreed to volunteer for five years, from 1985 to 1990, and build it up into an international centre of excellence, provided support came from the State Government in the form of basic building and ground. All the senior scientists in mathematical sciences in Kerala gave moral support to CMS and volunteered in running research and education programs at CMS (Figure 1).

In 2002, CMS has shifted its office and library into St Thomas College Pala compound into a one-story building donated by the Diocese of Palai, Kerala, India. From 2006 CMS became a Department of Science and Technology (DST), Government of India, Centre for Mathematical Sciences with sufficient financial support from DST. Then CMS could recruit up to 18 PhD students and a core faculty of five. Full-fledged activities of research, education at the undergraduate and research level started, conferences, workshops and lecture series by leading researchers around the world initiated. By this time CMS had become a research

Figure 1. A. M. Mathai’s Research Program in mathematics, statistics, and physics developed at the Department of Mathematics and Statistics, McGill University, Montreal, Canada, in the 1970’s.

centre of Banaras Hindu University (BHU) in Varanasi, India, Anna University, Chennai, India, and Mahatma Gandhi University, Kottayam, India.

Since 2006, the CMS has grown into a leading research centre of international repute in aspects of mathematical sciences, particularly those concerning research in astrophysics, special functions of mathematical physics, statistics, and fractional calculus. CMS had a publication series (books, proceedings, collections of research papers, lecture notes), a newsletter of two issues per year, a mathematics modules series (self-study books on basic topics) and a mathematical sciences series for the general public (documents archive available at http://neutrino.aquaphoenix.com/cmsintl.org/).

CMS has been making efforts to keep pace with developments all around the world and making contributions to bring India into the frontline research in mathematical sciences and allied disciplines. In the publications series CMS brought out around 40 publications, two of which were taken as such and published by Springer, New York, one in 2008 ( Mathai & Haubold, 2008; Figure 8) and one in 2010 ( Mathai, Saxena, & Haubold, 2010; Figure 9). Newsletter series continued until 2015. Ten modules were brought out in the modules series. Three of these in the area of linear algebra were combined and published by De Gruyter, Berlin and Boston in 2017 ( Mathai & Haubold, 2017a; Figure 2) and another three in the area of probability and statistics were combined and published by De Gruyter, Berlin and Boston in 2017 ( Mathai & Haubold, 2017b; Figure 3). These two volumes also served as education material for the UN-affiliated Regional Centres for Space Science and Technology Education in

Figure 2. Mathai & Haubold (2017a): LinearAlgebra for Physicists and Engineers (see also Mathai & Haubold, 2018a).

Figure 3. Mathai & Haubold (2017b): Probability and Statistics for Physicists and Engineers (see also Mathai & Haubold, 2018a).

Asia and the Pacific, Latin America and the Caribbean, Africa, and Western Asia (Mathai & Haubold, 2018a).

CMS has been functioning as a national centre for mathematical sciences. The main aim of the centre was to encourage youngsters to do research by organizing activities in mathematical sciences. The research groups at CMS cover major disciplines in mathematical sciences. The following are the research groups at CMS: Astrophysics Research Group, Fractional Calculus Research Group, Special Functions Research Group, Statistical Distribution Theory Research Group, Geometric Probability Research Group, Discrete Mathematics Research Group, Algebra and Analysis Research Group. The research groups at CMS are active in interdisciplinary research.

Astrophysics Research Group

The seniors in the group are HJH and A. M. Mathai ( Mathai & Haubold, 1988, 2018a; Figure 4). R. K. Saxena also joined this group. Senior research fellows in this group are D. Kumar, D. P. Joseph, S. S. Nair, and N. Sebastian. The research group focused on the solar neutrino problem and related issues of thermonuclear reaction rates, closed-form solutions for differential equations governing the internal structure of the Sun, and solar neutrino physics. Experimental data coming from solar neutrino experiments were analysed by Fourier analysis, wavelet analysis and diffusion entropy analysis. Application of special functions of mathematical physics has been encouraged in all research topics.

Fractional Calculus Research Group

The seniors in the group are R.K. Saxena, A. M. Mathai, and HJH ( Mathai & Haubold, 2017c, 2018b, 2018c; Figures 5-7). Senior research fellows working in this group are Seema S. Nair, Anitha Kattuveettil and Nicy Sebastian. Dilip

Figure 4. Mathai & Haubold (2018d): Erdélyi-Kober fractional calculus from a statistical perspective, inspired by solar neutrino physics.

Figure 5. Mathai & Haubold (2018c): An introduction to fractional calculus.

Kumar also has work in fractional calculus. A national level workshop in this area was held at CMS in November 2009. As an emerging field, CMS conducted several short courses and included fractional calculus in the syllabus of many SERC Schools at CMS. Differentiation and integration are usually regarded as

Figure 6. Mathai & Haubold (2018b): Matrix methods and fractional calculus.

Figure 7. Mathai & Haubold (2017c): Fractional and multivariable calculus: Model building and optimization problems.

discrete operations, in the sense that one differentiates or integrates a function once, twice, or any whole number of times. However, in some circumstances it is useful to evaluate a fractional derivative. In a letter to LHospital in 1695, Leibniz raised the possibility of generalizing the operation of differentiation to non-integer orders, and LHospital asked what would be the result of half-differentiating x. Leibniz replied It leads to a paradox, from which one day useful consequences will be drawn. The paradoxical aspects are due to the fact that there are several different ways of generalizing the differentiation operator to non-integer powers, leading to inequivalent results. Major contributions in this area from CMS are that fractional calculus was extended to functions of matrix argument in the real and complex domains and extended fractional calculus to the complex domain. Further, research focused on the well established fact that fractional time derivatives do not appear in any established fundamental theory of physics such as statistical mechanics, electrodynamics, or quantum mechanics. Particular attention was placed on the correlation between the diffusion equation and Schroedinger equation and their fractional counterparts.

Special Functions Research Group

The seniors in the group are R.K. Saxena, P. N. Rathie, S. B. Provost, P. Moschopoulos, HJH, and A. M. Mathai. “The H-function: Theory and Applications” (Figure 8) and “Special Functions for Applied Scientists” (Figure 9) are the latest books of this group, which were out from Springer, New York, in 2010 and 2008, respectively. Senior research fellows working in this area are Dilip Kumar, Nicy Sebastian, Dhannya P. Joseph, Seema S. Nair, Anitha Kattuveettil, and Naiju M. Thomas. There are mathematical theories in which some classes of special functions appear naturally. A familiar classification is by increasing complexity, starting with polynomials and algebraic functions and progressing through the elementary transcendental functions etc to the higher transcendental functions. Functions of hypergeometric type can be ordered by the behavior of singular points of the differential equations representing them, or by a group-theoretical analysis of their symmetries. But all these classifications are incomplete. For example, Mathieu functions fall outside the hypergeometric class, and gamma and zeta functions are not the solutions of simple differential equations. For rate and diffusion equations and their fractional generalizations, generalized hypergeometric functions like Meijer’s G-function and Fox’ H-function need to be made available analytically and numerically. They were of central interest to research and teaching of this research group.

Statistical Distribution Theory Research Group

The seniors in the group are A. M. Mathai, S. B. Provost, W. J. Anderson, HJH, and P. N. Rathie. Generalizations of Dirichlet integrals and Dirichlet densities are introduced by this group. Seemon Thomas, St Thomas College Palai, who is an associate of CMS, has received his PhD in this area. The pathway model introduced by Mathai in 2005 is popular in statistical distribution theory and non-extensive statistical mechanics. The first PhD on the pathway model was awarded to Shanoja S. Pai in 2010. Other researchers in this group are Nicy Sebastian, Seema S. Nair, Dhannya P. Joseph, Naiju M. Thomas, Prajitha P. and Princy T. CMS has conducted short-term courses in this area by eminent faculties of national and international standing. Major contributions in this area from

Figure 8. Mathai, Saxena, and Haubold (2010): The h-function: Theory and applications.

Figure 9. Mathai and Haubold (2008): Special functions for applied scientists.

CMS is that statistical distribution theory is shown to be a powerful tool in interpreting and developing Mellin convolutions of products and ratios, in establishing a connection to fractional integrals, in giving physical interpretations to the concept of M-convolutions introduced in 1970s (Figure 10).

Figure 10. Mathai, Provost, and Hayakawa (1995): Bilinear forms and zonal polynomials.

Geometrical Probability Research Group

A. M. Mathai is the senior in this group and Seemon Thomas has joined and is working in the group. Major contribution in this area is the replacements of differential and integral geometry techniques with Jacobians of matrix transformations techniques in the study of certain random geometrical configurations and extending the theory to functions of matrix argument ( Mathai, 1999a; Figure 12).

Discrete Mathematics Research Group

The seniors in this group are B.D. Acharya (former Advisor to Government of India) who joined as a Visiting Full Professor at CMS, R. Natarajan of Lakehead University, Canada, who was appointed as Reader at CMS, K.A. Germina of the Hill Area Campus of CMS. Two junior research fellows working in this group are Miss Alphy Joseph and Miss Sona Jose. The activities of this group also include building up a strong computer science base at CMS. A national level workshop in this area was held at CMS in February 2010 and another one was conducted in August 2010 in connection with the visit of Professor T. Zaslavsky (USA). Short-term courses are conducted for the enhancement of research in this area.

Algebra and Analysis Research Group

The seniors in this group are S.C. Mathew and P.G. Romeo. The junior research fellows working in this group are Miss Diana Mary George and Miss Ginu Varghese. Short-term courses are conducted to motivate the students in this area.

Research Recognitions

Figure 11. Mathai (1997): jacobians of matrix transfromations and functions of matrix argument.

Figure 12. Mathai (1999a): An introduction to geometrical probability: Distributional aspects with applications.

CMS is an approved research centre of three universities, Banaras Hindu University (BHU), Varanasi, Anna University, Coimbatore, and Mahatma Gandhi University (MGU), Kottayam. Students joining CMS for research can register themselves in any of the above universities for PhD and fulfill the PhD requirements of the respective universities. Some collaborative research is being done with other foreign and national institutions also. Six students from CMS got their PhDs in the area of statistics from BHU. Three PhD scholars from CMS received their PhDs from Anna University in the areas of applied analysis, graph theory and fuzzy mathematics. Four PhD scholars from CMS received their PhDs from MGU in the areas of fuzzy mathematics, graph theory and mathematical physics.

Library Development

Back volumes of more than 400 journals (with over 17,000 volumes) and more than 6000 books were available in mathematical sciences for reference purposes at the CMS library. Collections of reprints of individual articles are also arranged in CMS library for ready reference. CMS had built up a good library with books and journals collected by Professor Mathai from Canada and USA.

Portraits Collection

CMS library had a collection of the portraits of over 300 world mathematicians. The enlarged and laminated color portraits were placed on the library walls at CMS for motivating the school/colleges students to do work in mathematical sciences. Students of mathematical sciences have attended lectures about masters in mathematical sciences, they have learnt about the discoveries of these masters and the students had the opportunity to see their portraits and learn more about these great men and women by studying at CMS library.

The summary of achievements of CMS (later CMS was renamed CMSS (Centre for Mathematical and Statistical Sciences)) from 2007 to 2014, when it was a DST Centre, is the following:

Achievements of CMSS 2007-2014

Number of national/international level awards won by CMSS scholars = 19.

Number of research papers published in international refereed journals & proceedings = 170.

Number of PhDs produced by CMSS = 13 [6 from BHU, 4 from MGU, and 3 from Anna University, Chennai].

Number of research level students educated through all-India SERC Schools > 285 (Figure 13).

Number of undergraduate students educated through Undergraduate Mathematics Camps > 1020.

Number of paper presentations abroad by CMSS research scholars = 8.

Number of paper presentations in India by CMSS scholars > 65.

Number of invited and keynote addresses by the faculty of CMSS, including a lead keynote address and 3 keynote addresses at a United Nations’ workshops > 63.

Number of conferences sponsored or co-sponsored by CMSS = 7.

Number of distinguished visitors to CMSS from abroad > 32.

Number of distinguished national visitors to CMSS > 63.

Number of international research collaborators with CMSS faculty and scholars = 8.

Figure 13. Haubold & Mathai (2010): Proceedings of the third UN/ESA/NASA Workshop on the international heliophysical year 2007 and basic space science.

Number of national level research collaborators with CMSS faculty and scholars = 7.

The Honorary Director

Professor A. M. Mathai was honored by several institutions and national societies in India and thrice by the United Nations, once in Tokyo, Japan, once in Daejeon, Korea, and once in Quito, Ecuador and he gave keynote addresses also. He was the President of the Indian Mathematical Society for one term (the oldest professional society in India), current President of the Society for Special Functions and their Applications, Chairman of the Kerala State Statistical Commission (2014-2018), he was given Lifetime Achievement Award and Fellowship by the Indian Society for Probability and Statistics, and he is a Member of the Loka Kerala Sabha.

Achievements in Research and Education by A. M. Mathai

In this article Arak M. Mathai’s research accomplishments in mathematics and statistics are discussed which he achieved and shared with collaborators, postdocs, and students at McGill University, at the Centre for Mathematical and Statistical Sciences, and at annual United Nations workshops over a period of time of more than 60 years (Figure 14, Figure 15).

In the 1970s Mathai and collaborators developed an axiomatic theory providing axiomatic definitions for basic concepts in information theory and statistics. Focus was placed on functional equations, dispersion, principles of statistical

Figure 14. A.M. Mathai.

Figure 15. Centre for mathematical and statistical sciences.

inference, and inequalities among information theory and statistical measures. Conditional distributions, characterizations through independence of linear and quadratic statistics and characterizations with the help of differential equations and functional equations were developed in detail. Special functions, persistent in research and education in mathematics, statistics, and natural sciences, today available in the Digital Library of Mathematical Functions and accessible as a free Web-based collection of formulas (http://dlmf.nist.gov), cross-linked and with live graphics that can be magnified and rotated and downloaded became increasingly important in Mathai’s research. New results for theory and applications of generalized hypergeometric functions, among them Meijer’s G-function and Fox’ H-function extended to real scalar functions of matrix argument when the argument matrix is real or in the complex domain, became part of Mathai’s research. Eventually the full programme for research and education in these fields of mathematics and statistics with applications to natural sciences was published in a triology of monographs ( Mathai & Rathie, 1975; Mathai & Pederzoli, 1977; Mathai & Saxena, 1978; Figure 1). This research and education program set the stage for CMS (Figure 15).

A.M. Mathai was born on 28 April 1935 in Arakulam, near Palai, in the Idukki district of Kerala, India as the eldest son of Aley and Arakaparampil Mathai. After completing his high school education in 1953 from St. Thomas High School, Palai, he joined St. Thomas College, Palai with record marks and obtained his B.Sc. degree in Mathematics in 1957. In 1959, he completed his Master’s Degree in Statistics from the University of Kerala, Thiruvananthapuram, Kerala, India; he achieved a first class degree, the first rank and a gold medal. Then, he joined St. Thomas College, Palai, University of Kerala, as a Lecturer in Mathematics and served there until 1961. He obtained a Canadian Commonwealth scholarship in 1961 and went to the University of Toronto, Canada, to complete his M. A. degree in Mathematics in 1962. He was awarded a PhD from the University of Toronto, Canada, in 1964. Then, he joined McGill University, Canada as an Assistant Professor in 1964. From 1968 to 1978, he was an Associate Professor there. He became a Full Professor of McGill in 1979 and served the Department of Mathematics and Statistics until he took early retirement in 2000. He is also the founder of the Canadian Journal of Statistics and the Statistical Science Association of Canada which became the Statistical Society of Canada. As of this date, A. M. Mathai is an Emeritus Professor of Mathematics and Statistics at McGill University, Canada (http://www.math.mcgill.ca/people/mathai), and Honorary Director of the Centre for Mathematical and Statistical Sciences, India. He has published over 300 research papers and more than 25 books on topics in mathematics, statistics, physics, astrophysics, and biology (selected titles Figure 1 to Figure 12). He is a Fellow of the Institute of Mathematical Statistics, National Academy of Sciences of India, Indian Society for Probability and Statistics, President of the Mathematical Society of India (2015-16), and an Elected Member of the International Statistical Institute. He was the Chairman of the Kerala State Statistical Commission (2014-2018) and President of the Society for Special Functions and their Applications (2018-). He was the principal organizer of the elements of Physics and Mathematics in the implementation of the annual workshops of the United Nations Basic Space Science Initiative (1991-2012; Figure 13), organized by ESA, NASA, and JAXA and hosted by member States of the United Nations for the benefit all Member States of the United Nations (Mathai & Haubold, 2018a; Pyenson, Mathai, & Haubold, 2019).

1. Fractional Calculus: Reaction and Diffusion

Fractional integrals, fractional derivatives and fractional differential equations were available only for real scalar variables. The most popular fractional integrals in the literature are Riemann-Liouville fractional integrals given by the following:

${}_{a}D{}_x^{-\alpha }f=\frac{1}{\Gamma \left(\alpha \right)}{\displaystyle {\int }_a^x}{\left(x-t\right)}^{\alpha -1}f\left(t\right)\text{d}t\mathrm{,}\Re \left(\alpha \right)>\mathrm{0,}$ (1.1)

where $\Re (\cdot )$ denotes the real part of $(\cdot )$ .

${}_{x}D{}_b^{-\alpha }f=\frac{1}{\Gamma \left(\alpha \right)}{\displaystyle {\int }_x^b}{\left(t-x\right)}^{\alpha -1}f\left(t\right)\text{d}t\mathrm{,}\Re \left(\alpha \right)>0.$ (1.2)

Here $-\alpha$ in the exponent of D indicates an integral. The D with positive exponent ${}_{a}D{}_x^{\alpha }f\mathrm{,}{}_{x}D{}_b^{\alpha }f$ is used to denote the corresponding fractional derivatives. Here (1.1) is called Riemann-Liouville left-sided or first kind fractional integral of order $\alpha$ and (2.2) is called Riemann-Liouville fractional integral of order $\alpha$ of the second kind or right-sided. If $a=-\infty$ and $b=\infty$ then (1.1) and (1.2) are called Weyl fractional integrals of order $\alpha$ and of the first kind and second kind, respectively, or the left-sided and right-sided ones. Mathai was trying to find an interpretation or connection of fractional integrals in terms of statistical densities and random variables. In Mathai (2009), an interpretation is given for Weyl fractional integrals as densities of sum (first kind) and difference (second kind) of independently distributed real positive random variables having special types of densities. Fractional integrals were also given interpretations as fractions of total integrals coming from gamma and type-1 beta random variables. Also Weyl fractional integrals were extended to real matrix-variate cases there.

1.1. Mellin Convolutions of Products and Ratios

Then while working on Mellin convolutions of products and ratios, Mathai found that a fusion of fractional calculus and statistical distribution theory was possible which also opened up ways of extending fractional calculus to real scalar functions of matrix argument, when the argument matrix is real or in the complex domain. Let us consider real scalar variables first. The Mellin convolution of a product of two functions $f_1\left(x_1\right)$ and $f_2\left(x_2\right)$ says the following: Consider the integral

$g_2\left(u_2\right)={\displaystyle {\int }_v}\frac{1}{v}{f}_1\left(\frac{u}{v}\right)f_2\left(v\right)\text{d}v\mathrm{.}$ (1.3)

Then the Mellin transform of $g_2\left(u_2\right)$ , with Mellin parameter s, is the product of the Mellin transforms of $f_1$ and $f_2$ . That is

$M_{g_2}\left(s\right)=M_{f_1}\left(s\right)M_{f_2}\left(s\right),$ (1.4)

where

$M_{f_1}\left(s\right)={\displaystyle {\int }_0^{\infty }}{x}_1^{s-1}{f}_1\left(x_1\right)\text{d}{x}_1\text{and}{\displaystyle {\int }_0^{\infty }}{x}_2^{s-1}{f}_2\left(x_2\right)\text{d}{x}_2=M_{f_2}\left(s\right)\mathrm{,}$

whenever they exist. If $g_2\left(u_2\right)$ is written as

$g_2\left(u_2\right)={\displaystyle {\int }_v}\frac{1}{v}{f}_1\left(v\right)f_2\left(\frac{u}{v}\right)\text{d}v$ (1.5)

then also the formula in (1.4) holds. Thus, the Mellin convolution of a product has the two integral forms in (1.3) and (1.5). But, Mellin convolution of a ratio will have four different representations. Two of these are the following:

$M_{g_1}\left(u_1\right)=M_{f_1}\left(s\right)M_{f_2}\left(2-s\right)$ (1.6)

where

$g_1\left(u_1\right)={\displaystyle {\int }_v}v{f}_1\left(uv\right)f_2\left(v\right)\text{d}v$ (1.7)

and

$M_{g_1}\left(s\right)=M_{f_1}\left(2-s\right)M_{f_2}\left(s\right)$ (1.8)

where

$g_1\left(u_1\right)={\displaystyle {\int }_v}\frac{v}{u_1^2}{f}_1\left(\frac{v}{u_1}\right)f_2\left(v\right)\text{d}v\mathrm{.}$ (1.9)

1.2. Statistical Interpretations of Mellin Convolutions

Let $x_1$ and $x_2$ be real scalar positive random variables, independently distributed, with densities $f_1\left(x_1\right)$ and $f_2\left(x_2\right)$ , respectively. Let $u_2=x_1{x}_2$ and $u_1=\frac{x_2}{x_1},v=x_2$ . Then the Jacobians are $\frac{1}{v}$ and $-\frac{v}{u_1^2}$ respectively or

$\text{d}{x}_1\wedge \text{d}{x}_2=\frac{1}{v}\text{d}{u}_2\wedge \text{d}v$

and

$\text{d}{x}_1\wedge \text{d}{x}_2=-\frac{v}{u_1^2}\text{d}{u}_1\wedge \text{d}v\mathrm{.}$

The joint density of $x_1$ and $x_2$ is $f_1\left(x_1\right)f_2\left(x_2\right)$ due to statistical independence and then the marginal densities of $u_2$ and $u_1$ , denoted by and are the following:

(1.10)

and

(1.11)

In, if is taken as v then the roles of and change in (1.10). If in is taken as v then we get (1.7) with the roles of and interchanged. Hence gives two forms and gives two forms for Mellin convolution of ratios. The Mellin convolutions in (1.4) and (1.8) can be interpreted in terms of random variables. gives due to independence where denotes the expected value. That is,

Similarly

This means, which is (1.4) the Mellin convolution of a product. Now, consider. Then due to statistical independence. This means

and

In other words, , which is (1.8), one form of Mellin convolution of a ratio. Mellin convolutions of products and ratios make direct connection to product and ratio of real positive random variables.

1.3. Mellin Convolutions, Statistical Densities and Fractional Integrals

Let and be statistical densities as in Section 1.2. Let have a type-1 beta density with parameters or

and zero elsewhere. In statistical problems, the parameters are real but the integrals hold for complex parameters and hence the conditions are given for complex parameters. Let an arbitrary density. Then the Mellin convolution of a product is the following:

(1.12)

where

(1.13)

is an Erdélyi-Kober fractional integral of order and of the second kind with parameter. This is a direct connection among Erdélyi-Kober fractional integral of the second kind, Mellin convolution of a product and statistical density of product of two independently distributed real scalar positive random variables where one has a type-1 beta density with parameters and the other has an arbitrary density.

1.4. General Definition for Fractional Integrals

Motivated by this observation, Mathai has given a new definition for fractional integrals of the first and second kinds of order. Let

and elsewhere, and where and are pre-fixed functions and is an arbitrary function, where and need not be statistical densities. Consider the Mellin convolution of a product. Then using the same notation as there is

(1.14)

(1.15)

But (1.15) gives a Weyl fractional integral of the second kind of order. If v is bounded above by a constant b then (1.15) is a Riemann-Liouville fractional integral of the second kind of order. Thus, by specifying and it can be seen that all the various definitions of fractional integrals of order of the second kind can be obtained from (1.14). Evidently when and then one has (1.13) or Erdélyi-Kober fractional integral of order of the second kind with parameter.

1.5. First Kind Fractional Integrals

Let and be as given above in Section 1.4. Consider (1.9), the Mellin convolution of a ratio. Then

(1.16)

Let and. Then (1.16) reduces to the following form:

(1.17)

Then

(1.18)

is a statistical density, when f ia a statistical density. This is an Erdélyi-Kober fractional integral operator of the first kind of order and parameter, denoted by. From (1.16), by specializing and one can get all the various definitions of fractional integrals of order of the first kind in the real scalar case. This formal definition of fractional integrals as Mellin convolutions of ratio and product was introduced formally in Mathai (2013). A geometrical interpretation of fractional integrals as fractions of integral over a simplex in n-space is given in Mathai (2014).

1.6. Extension of Fractional Integrals to Real Matrix-Variate Case

Mathai (2009) introduced fractional integrals in the real matrix-variate case but they could not be given any physical interpretations. In Mathai (2013) there are interpretations in terms of statistical distribution problem and M-convolutions introduced in Mathai (1997) (Figure 11). A meaningful interpretation of M-convolutions is given as the densities of product and ratio of matrix-variate random variables.

All the matrices appearing here are real positive definite matrices, unless specified otherwise. The notation means the real symmetric matrix, , is positive definite. means the positive definite square root of the positive definite matrix. If is then the wedge product of differentials will be denoted by, that is,

and it is if and. Also, means the integral over all of the real-valued scalar function of X, such that

where A and B are positive definite constant matrices. Here Jacobians of matrix transformations are needed. These will be given as lemmas, without proofs. For proofs and for other such Jacobians see Mathai (1997).

Lemma 1.1. Let be matrix of distinct real scalar variables’s. Let A be and B be nonsingular constant matrices. Then

(1.19)

where denotes the determinant of.

Lemma 1.2. Let be. Let A be a nonsingular constant matrix. Then

(1.20)

Lemma 1.3. Let X be a nonsingular matrix. Let. Then

(1.21)

Lemma 1.4. Let be. Let be a lower triangular matrix with positive diagonal elements, that is, . Then

(1.22)

With the help of (1.22) one can evaluate a matrix-variate gamma integral and write the result as

(1.23)

where the real matrix-variate gamma integral is

(1.24)

with denoting the trace of. Apply Lemma 1.4 to X in the integrand of (1.24). Then the integral splits into integrals over’s and’s, for, and both types of integrals can be evaluated by using a real scalar variable gamma integral. Then the final result is that of (1.23). Combining (1.24) and (1.20) one can define a matrix-variate gamma density as

(1.25)

Since the total integral is 1, the identity follows from (1.25),

(1.26)

This identity will be used to establish fractional derivatives in a class of matrix-variate functions. The real matrix-variate type-1 beta density is defined as

(1.27)

There is a corresponding type-2 beta density, which is of the form

(1.28)

1.7. Fractional Integrals for the Real Matrix-Variate Case

With the preliminaries in Section 1.6 one can define fractional integrals in the real matrix-variate case. Let and be real matrix-variate random variables, independently distributed. Let be defined as the symmetric product of and and let be defined as the symmetric ratio of over. Let. Then with the help of the above lemmas one can show that, ignoring sign,

(1.29)

Denoting the densities of and as and, one can compute these by using the densities and through transformation of variables and they will be the following:

(1.30)

and

(1.31)

Let be a type-1 beta density of the type in (1.27) with parameters. Note that in (1.27) the parameters are. Then will be of the following form:

(1.32)

This in (1.32) for is an Erdélyi-Kober fractional integral of the second kind of order and parameter and hence Mathai called the integral as Erdélyi-Kober fractional integral of order and parameter of the second kind in the real matrix-variate case. In a similar fashion, Erdélyi-Kober fractional integral of order of the first kind with parameter, available from (1.31) by taking as a real matrix-variate type-1 beta with parameters is the following:

(1.33)

The density of, again denoted by, is given by

(1.34)

where the first kind Erdélyi-Kober fractional integral in the matrix-variate case is given in(1.33). The above notations as well as a unified notation for fractional integrals and fractional derivatives were introduced by Mathai (2013, 2014, 2015).

The above results in the real matrix-variate case are extended to complex matrix-variate cases, see Mathai (2013), to many matrix-variate cases, see Mathai (2014) and also the corresponding fractional derivatives in the matrix-variate case are worked out in Mathai (2015). The matrix differential operator introduced in Mathai (2015) is not a universal one, even though it works on some wide classes of functions. The matrix differential operator is introduced through the following symbolic representation. Let D be a differential operator defined for real matrix-variate case. Then and represent αth order fractional derivative and fractional integral respectively. Then

This is the αth order fractional derivative in Riemann-Liouville sense. Consider

is the αth order fractional derivative in the Caputo sense. In the Caputo case, operates on f first and then the fractional integral is taken, whereas in the Riemann-Liouville sense, the th order fractional integral is taken first and then operates on this. A universal differential operator D in the real as well as complex matrix-variate case is still an open problem for further research.

2. Krätzel Integral: Thermonuclear Functions

Let x be a real scalar positive variable. Consider the integrals

(2.1)

and

(2.2)

Structures such as the ones in (2.1) and (2.2) appear in many different areas. This (2.2) for is the basic Krätzel integral, see Krätzel (1979). For and general is the generalized Krätzel integral. An integral transform of the form

(2.3)

where is arbitrary so that exists, is known as Krätzel transform. Mathai has investigated various aspects of (2.1) and (2.2) in detail and he has also introduced a statistical density in terms of Krätzel integral. The structures in (2.2) and (2.1) can be generated as Melin convolutions of product and ratio. Consider the real scalar variables and and the corresponding functions and. Then it is seen from (1.10) that the Mellin convolution of a product is given by

(2.4)

or, and the Mellin convolution of a ratio, from (1.7), as

(2.5)

or. Let and be generalized gamma functions of the form

(2.6)

Then of (2.5) reduces to the form

(2.7)

This is the form in (2.1). Now, consider Mellin convolution of a product when and are generalized gamma functions in (2.6). Then reduces to the form

(2.8)

This is the form in (2.2). Hence (2.1) and (2.2) can be treated as Mellin convolutions of ratio and product when and are generalized gamma functions.

Note that if and are multiplied by the corresponding normalizing constants and then and become statistical densities. Let and be independently distributed real scalar positive random variables. Let. Then the densities of and are given by (2.4) and (2.5) multiplied by the appropriate constants and reduce to the forms in (2.8) and (2.7), multiplied by appropriate constants. In other words, (2.1) and (2.2), multiplied by appropriate constants, can be looked upon as the density of a ratio and product respectively.

The integrand in (2.2) for and normalized is the inverse Gaussian density available in stochastic processes. The integral in (2.2) for is the basic reaction-rate probability integral, which will be considered later. Mathai (2012) has introduced a Krätzel density associated with (2.1) and (2.2) and it is shown that one has general Bayesian structures in (2.1) and (2.2). For example, let us consider a conditional density of y, given x, in the form

(2.9)

and can act as the normalizing constant. In other words, the conditional density is a generalized gamma density. Let the marginal density of x be given by, and can act as a normalizing constant, a generalized gamma density. Then the joint density of y and x is given by

(2.10)

Then the unconditional density of y, , is available by integrating out x from this joint density. That is,

(2.11)

Now, compare (2.2) and (2.11). They are of one and the same forms. Hence (2.2), multiplied by an appropriate constant, can be considered as an unconditional density in a Bayesian structure.

For in (2.1) and (2.2) one can extend the integrals to the real and complex matrix-variate cases. Mathai has also looked into this problem of Krätzel integrals in the matrix-variate cases. There will be difficulty with the Jacobians if one considers general parameters and in the matrix-variate case. The type of difficulties that can arise is described in Mathai (1997) by considering the transformation when. In the real matrix-variate case the scalar quantity is replaced by the determinant and exponent is replaced by for or if a is also replaced by a positive definite constant matrix. Mathai has also extended Baysian structures, densities of product and ratio, inverse Gaussian density, Krätzel integral and Krätzel density, to matrix-variate cases. When the matrix is in the complex domain is replaced by = absolute value of the determinant of, where is a matrix in the complex domain.

3. Pathway Model: Entropy, Probability, Dynamics

In a physical system the stable solution may be exponential or power function or Gaussian. This is the idealized situation. But in reality the solution may be somewhere nearby the ideal or the stable situation. In order to capture the ideal situation as well as the neighboring unstable situations, a model with a switching mechanism was introduced by Mathai (2005). A form of this was proposed in the 1970’s by Mathai in connection with population studies. This was a scalar variable case. Then the ideas were extended to matrix-variate cases and brought out in 2005. For the real scalar positive variable situation, the model is the following:

(3.1)

If (3.1) is to be used as a statistical density then is the normalizing constant there. Otherwise is a constant, may be and then (3.1) will be a mathematical model. For we can write and then (3.1) becomes

(3.2)

When then and go to

(3.3)

Note that in (3.1) is in the family of generalized type-1 beta family of functions, whereas is in the family of generalized type-2 beta family of functions and belongs to the generalized gamma family of functions. Thus, when the pathway parameter q, goes from to 1 we have one family of functions, when q is from 1 to we have another family of functions and when we have a third family of functions. Thus, all the three cases are contained in (3.1), which is the pathway model for the real positive scalar variable case. Replace x by, , to extend the families over the real line.

When are statistical densities, then (3.1) to (3.3) give a distributional pathway. Mathai has also established a parallel pathway in terms of entropy optimization and in terms of differential equations. These give entropic and differential pathways as well. For example, consider the optimization of Mathai’s entropy, namely

(3.4)

where is a density function of x, and x can be real scalar or vector or matrix variable. A density means that for all x and. If one takes the limit when then (3.4), for real scalar x, reduces to

(3.5)

where is Shannon’s entropy or measure of “uncertainty’’ or the complement of “information’’. In (3.5), is taken as zero when. Consider the optimization of (3.4) subject to the conditions (a): and (b):. For, condition (b) becomes since the total probability is 1. For, (a) means that the first moment is fixed. This can correspond to the physical law of conservation of energy when dealing with energy distribution. If one uses calculus of variation to optimize (3.4) then the Euler equation is

(3.6)

where and are Lagrangian multipliers. Note that (3.6) gives the structure

for some and, which means

(3.7)

for some and. For and one has the model in (3.1) with replaced by. Thus, for one has an entropic pathway. Similarly one can consider the corresponding differential equations to obtain a differential pathway.

The original paper Mathai (2005) deals with rectangular matrix-variate case. Let be and of rank m be a matrix of distinct real scalar variables’s. Let A be and B be constant positive definite matrices. Consider the function

(3.8)

where are scalars, I is a identity matrix and is a constant. If (3.8) is to be taken as a density then is the normalizing constant there. For, goes to

(3.9)

and when, and go to

(3.10)

If a location parameter matrix is to be introduced then replace X by where M is a constant matrix.

Note that the structure is the structure of the volume content of a parallelotope in Euclidean n-space. Look at the m rows of X. These are vectors. These can be taken as m points in n-dimensional Euclidean space. These m vectors, , are linearly independent when the rank of X is m. These taken in a given order can form a convex hull and a m-parallelotope. The volume of this m-parallelotope is the determinant. Hence is the volume content of a generalized m-parallelotope.

Also is a generalized quadratic form. For and it is a quadratic form in the vector variable. Thus the theory of quadratic form and generalized quadratic form can be extended to a wider class represented by the pathway model (3.8). The current theory of quadratic form and bilinear form in random variables is confined to samples coming from a Gaussian population, see Mathai and Provost (1992), Mathai, Provost and Hayakawa (1995) (Figure 10). The results on quadratic and bilinear forms can now be extended to the wider class of pathway models. One problem in this direction is discussed in Mathai (2007). The matrix-variate pathway model in Mathai (2005) is extended to complex domain in Mathai and Provost (2005, 2006). Some works in the scalar complex variable case, associated with normal or Gaussian population, are available in the literature with applications in sonar, radar, communication and engineering problems. Some applications of hermitian forms, corresponding to the quadratic forms in Mathai and Provost (1992), in light scattering and quantum mechanics are also available in the literature.

Note that (3.8) for is a matrix-variate type-1 beta density or is a type-1 beta matrix. This is the exact form of the matrix appearing in the generalized analysis of variance and design of experiments areas, in the likelihood ratio test involving one or more multivariate normal or Gaussian populations etc, a summary of the contributions of Mathai and his co-workers is available from Mathai and Saxena (1973). The theory available there is based on Gaussian populations. Now, generalized analysis of variance can be examined in a wider pathway family so that the limiting form corresponding to (3.10), will be the Gaussian case.

While exploring a reliability problem, Mathai (2003) came across a multivariate family of densities, which could be taken as a generalization of type-1 Dirichlet family of densities. Then Mathai and his co-workers introduced several generalizations of type-1 and type-2 Dirichlet densities, see for example Thomas and Mathai (2009). For the different generalizations of type-1 and type-2 Dirichlet family, a number of characterization results are established showing that these models could also be generated by products of statistically independently distributed real scalar random variables. This is exactly the same structure available in the likelihood ratio criteria in the null cases of testing hypotheses on the parameters of one or more Gaussian populations as well as in the determinant or in the model (3.8) for. Thus, it is already shown that these three areas are connected.

In (3.1) if one puts then one gets Tsallis statistics in non-extensive statistical mechanics. Also, (3.2) for as well as for some general is superstatistics (Beck, 2004; Cohen, 2004). (3.2) and its limiting form (3.3) are covered in superstatistics but (3.1) is not covered because superstatistics considerations deal with a conditional density of generalized gamma form as well as the marginal density a generalized gamma form then the unconditional density, which is superstatistics in statistical terms from a Bayesian point of view, can only produce a type-2 beta form, namely (3.2) form and not (3.1) form. Thus, superstatistics is also a special case of the pathway model in the real scalar positive variable case.

In the pathway idea itself there is an open area which is not yet explored. The scalar version of the pathway model in (3.1) to (3.3) can be looked upon as the behavior of a hypergeometric series (binomial series) going to (exponential series). That is,

(3.11)

(3.12)

From the point of view of a hypergeometric series, the process (3.11) to (3.12) is the process of a binomial series going to an exponential series. But a Bessel series can also be sent to an exponential series. For example, consider the Bessel series

(3.13)

Therefore, a generalized form, covering the path towards the exponential form

, is also a Bessel form. This path of a Bessel formgoing to an exponential form can produce a large variety of results. This area has open problems for further research.

4. Special Functions of Matrix Argument

A multivariate function usually means a function of many scalar variables. This is different from a matrix-variate function or a function of matrix argument. Functions of matrix argument are real-valued scalar functions where X is a square or rectangular matrix. For example, for a matrix X, = determinant of X, = trace of X are real-valued scalar functions when X is real. Even for a square matrix X, the square root cannot be uniquelydetermined unless further conditions are imposed on X. If one uses the definition, then the square root of A, one can have manycandidates for B. For example, for a simple matrix like a identity matrix, are square roots:

If one restricts A and to be positive definite matrices then is the

only candidate here. Hence, if X is real positive definite or Hermitian positive definite then can be uniquely defined. Therefore, functions of matrix argument are developed mainly when the argument matrix is either real positive definite or Hermitian positive definite. There are three approaches available in the literature for functions of matrix argument, that is, real-valued scalar functions of matrix argument X. For convenience, all the matrices appearing in this section are positive definite denoted by, real or Hermitian, unless stated otherwise. One definition is through Laplace and inverse Laplace transforms. This development is due to Herz (1955) and others. Here the basic assumption of functional commutativity is used, that is, even if. For example, determinant and trace will satisfy this property. When X is real symmetric then there exists an orthonormal matrix Q such that where are the eigenvalues of X. Then

(4.1)

or, which is a function of real variables’s, when and real, has become a function of D which is of p real variables, under this assumption of functional commutativity. If and then

(4.2)

Therefore

(4.3)

the Laplace transform of because (4.3) is not consistent with the definition of multivariate Laplace transform. In (4.2) the non-diagonal terms appear twice. In the multivariate Laplace transform, the variables and the corresponding parameters must appear only once each. If one considers a modified parameter matrix for all i and j, and

(4.4)

then

(4.5)

is the Laplace transform in the real symmetric positive definite matrix-variate case, where is the parameter matrix and stands for the wedge product of the differentials’s or

(4.6)

Under this approach, a hypergeometric function of matrix argument, denoted by

where and are scalar parameters and X is a real positive definite matrix, is defined by a Laplace and inverse Laplace pair. Under this definition, explicit forms are available only for and. Details of the definition and properties may be seen from Herz (1955) and from Mathai (1997).

The second approach is through zonal polynomials, developed by James (1961)and Constantine (1963). Here also functional commutativity is implicitly assumed, though not stated explicitly. Under this definition, a hypergeometric series is defined as follows:

(4.7)

where are zonal polynomials of order k, and

(4.8)

Here is the Pochhammer symbol and is the generalized Pochhammer symbol. All terms of the series in (4.7) are explicitly available but since zonal polynomials are complicated to compute, only the first few terms up to are computed. Details of zonal polynomials may be found, for example from the book Mathai, Provost and Hayakawa (1995). The definition through (4.5) and its inverse Laplace form and the definition through (4.7) are not very powerful in extending results in the univariate case to the corresponding matrix-variate case. When (4.7) is used to extend univariate results to matrix-variate cases the following two basic results will be essential. These will be stated here as lemmas without proofs.

Lemma 4.1.

(4.9)

where

(4.10)

with defined in (4.8).

Lemma 4.2.

(4.11)

Starting from 1970, Mathai developed functions of matrix argument through M-transforms and M-convolutions. Under M-transform definition, a hypergeometric function with matrix argument is defined as that class of functions satisfying functional commutativity and the integral equation

(4.12)

where

For example

(4.13)

Since the left side in (4.12) is a function of only one parameter, one cannot normally recover a function of p scalar variables. It is conjectured that when is analytic in the cone of positive definite matrices, one has uniquely recovered from the right side of (4.12). This is not established yet and also an explicit form of an inverse or through the right-side of (4.12) is an open problem. But (4.12) is the most convenient form to extend univariate results on hypergeometric functions to the corresponding class of matrix-variate cases. In general, when one goes from a univariate case, such as a univariate function, to a multivariate case, there is nothing called a unique multivariate analogue. Whatever be the properties of the univariate function that one wishes to preserve in the multivariate analogue, one may be able to come up with different functions as multivariate analogues of a univariate function. Hence, a class of multivariate analogues is more appropriate than a single multivariate analogue. Properties of M-transforms and properties of hypergeometric family coming from (4.12) are available in the book Mathai (1997). When Mathai introduced M-convolutions and M-transforms, details in Mathai (1997), no physical meaning could be found. Now, a physical interpretation is available for M-convolutions as densities of products and ratios of matrix random variables, as illustrated in Sections 1.7.

5. Geometric Probabilities: Probability Density Function

The work until 1999 is summarized in Mathai (1999a). The work started as an off-shoot of the work in multivariate statistical analysis. Mathai noted that the moment structure for many types of random geometric configurations was that of product of independently distributed type-1 beta, type-2 beta or gamma random variables. Such structures were already handled by Mathai and his co-workers in connection with problems in multivariate statistical analysis. Earlier contributions of Mathai in this area are available from Chapter 4 of the book Mathai (1999a). Then Ruben, a colleague of Mathai at McGill University, one day gave a copy of his paper showing a conjecture in geometrical probabilities, called Miles’ conjecture about a re-scaled, relocated random volume, generated by uniformly distributed random points in n-space, as asymptotically normal when. The proof was very roundabout. Mathai noted that it could be proved easily with the help of the asymptotic expansions of gamma functions. This paper was published in Mathai (1982). Then Mathai formulated and proved parallel conjectures regarding type-1 beta distributed, type-2 beta distributed points and gamma distributed random points and published a series of papers. Then Mathai noted that many European researchers were working on distances between random points, and random areas when the random points are in particular shapes such as triangles, parallelograms, squares, rhombuses etc. As generalizations of all these classes of problems, Mathai generalized Buffon’s clean-tile problem, the starting point of geometrical probabilities. He considered placing a ball at random in a pyramid with polygonal base, defining “at random’’ in terms of kinematic measure, Mathai (1999c). When mixing geometry with probability or measure theory, or in the area of stochastic geometry, the basic axioms of probability are not sufficient, as pointed out by Bertrand’s or Russell’s paradoxes. We need an additional axiom of invariance under Euclidean motion. Another contribution of Mathai in this area is Mathai (1999b) where he has shown that the usual complicated procedures coming from integral geometry and differential geometry are not necessary for handling certain types of random volumes but only the simple properties of functions of matrix argument and Jacobians of matrix transformation are sufficient. The procedure is illustrated in the distribution of volume content of parallelotope generated by random points in Euclidean n-space. The work on geometrical probabilities is currently progressing in the areas of random sets, image processing etc. The book, Mathai (1999a), only deals with distributional aspects of random geometric configurations.

As an application of geometrical probabilities, Mathai and his co-authors looked into a geography problem of city designs of rectangular grid cities, as in North America, versus circular cities as in Europe, with reference to travel distance, and the associated expense and loss of time, from suburbs to city core, see Mathai (1998), Mathai and Moschopoulos (1999a).

6. Astrophysics: Solar Neutrinos

After publishing the books Mathai and Saxena (1973, 1978) physicists were using results in special functions in their physics problems. HJH came to Montreal, Canada in 1982 with open problems on reaction-rate theory, solar models, solar neutrinos, and gravitational instability. The idea was to get exact analytical results and analytical models where computations and computer models were available. Mathai figured out that the problems connected with reaction-rate theory and solar neutrinos could be tackled once the following integral was evaluated explicitly (Critchfield, 1972; Fowler, 1984):

(6.1)

The corresponding general integral is

(6.2)

In 1982 Mathai could not find any mathematical technique of handling (6.2) or its particular case (6.1). He noted that (6.2) could be written as a product of two integrable functions and thereby as statistical densities by multiplying with appropriate normalizing constants. Then the structure in (6.2) could be converted into the form

(6.3)

and the right side of (6.3) is the density of a product of two real scalar positive independently distributed random variables with densities and respectively with. Take

where are normalizing constants. When the density of u, denoted by, is given by

(6.4)

Now, it is only a matter of evaluating the density by using some other method. Note that where and are independently distributed means

where denotes the expected value of. Then for (6.4)

which also shows that since at is 1. Evaluations of and are not necessary for our procedure to hold.

Therefore

Hence is available from the inverse Mellin transform. That is,

(6.5)

Comparing (6.4) with (6.5) the required integral is given by the following:

The right side of (6.6) is a H-function.

For the reaction-rate probability integral, and. In this case, the H-function in (6.6) reduces to a G-function and explicit computable series forms are also given by Mathai and his co-workers. Problems considered were resonant reactions, non-resonant reactions, depleted case, and high energy tail cut off. A summary of the work until 1988 is available in the research monograph Mathai and Haubold (1988). After publishing papers in physics by using statistical techniques it was realized that the density of a product of independently distributed real positive random variables was nothing but the Mellin convolution of a product. Hence, one could have applied Mellin and inverse Mellin transform techniques there. The work in this area of reaction-rate also resulted in two encyclopedia articles, see Haubold and Mathai (1997, 1998).

Analytic Solar Models

Another attempt was to replace the current computer model for the Sun with analytic models. The idea was to assume a basic model for the matter density distribution in the Sun or in main sequence stars which could be treated as a sphere in hydrostatic equilibrium. Let r be an arbitrary distance from the center of the Sun and let be the radius of the Sun. Let so that. The model for the matter density distribution is taken as

(6.7)

where c is the central core density when. The parameters and are selected to agree with observational data. Then by using (6.7), expressions for mass, pressure, temperature and luminosity are computed by using physical laws. Then by using known observations, or comparing with known data on mass, pressure etc the best values for and are estimated so that close agreement is there with observational values of mass, pressure etc. Some of the results until 1988 are given in the monograph Mathai and Haubold (1988).

Another area that was looked into was the gravitational instability problem concerning the evolution of large scale structure in the Universe. The problem was formulated in the form of differential equations. Mathai tried to change the operator to. Then the differential equation got simplified. Then he changed the dependent variable and found that the differential equation became a particular case of G-function differential equation. This resulted in the first paper of Mathai in integer order differential equations and it was published in Mathai (1989). Then the results were applied to gravitational instability problem (Haubold & Mathai, 1988).

Another area looked into was the solar neutrino problem (Davis Jr., 2003; Sakurai, 2014). HJH and Mathai tried to come up with appropriate models to model the solar neutrino data. Mathai had noted that the graph of the time series data looked similar to the pattern that he had seen when working on modeling of the chemical called Melatonin in human body. Usually what is observed is the residual part of what is produced minus what is consumed or converted or lost. Hence the basic model should be an input-output type model. The necessary theory is available in Mathai (1993c). The simplest input-output model is an exponential type input and an exponential type output so that the residual part. When and are independently and identically exponentially distributed then u has a Laplace distribution. One model HJH and Mathai tried was Laplace type random variables over time so that the graph will look like blips at equal or random points on a horizontal line. If the time-lag is shortened then the blips will start joining together. If exponential models of different intensities, that is, in the input-output model , if is different for different blips then the pattern can be brought to the pattern seen in nature or the pattern seen from the data.

7. Special Functions of Mathematical Physics

Mathai and his co-workers are credited with popularizing special functions, especially G and H-functions, in statistics and physics. Major part of the special function work was done with co-worker Saxena. They thought that they were the first one to use G and H-functions in statistical literature. But D. G. Kabe pointed out that he had expressed a statistical density in terms of a G-function in 1958. This may be the first paper in statistics where G or H-function was used. Most probably the use of G and H-function in physics an engineering areas started after the publication of the books Mathai and Saxena (1973, 1978). The first work on the fusion of statistical distribution theory and special functions started by creating statistical densities by using generalized special functions. In this connection the most general such density is based on a product of two H-functions, which appeared in Mathai and Saxena (1969). Another area that was looked into was Bayesian structures. The unconditional density in Bayesian analysis is of the form

(7.1)

What are the general families of functions for and so that the integral in (7.1) can be evaluated? One can construct some general mixing families of and.

Another family of problems that was looked into were the null and non-null distributions of the likelihood ratio criterion or -criterion for testing hypotheses on the parameters of one or more multinormal populations. Consider the vector having the density

(7.2)

where is a constant vector, known as the mean-value vector here. For if ‘s are independently distributed with the same density in (7.2) then we say that we have a simple random sample of size N from the p-variate normal or Gaussian population (7.2). Suppose that we want to test a hypothesis H_o:V = is diagonal. This is called the test for independence in the Gaussian case. Then the -criterion can be shown to have the structure:

(7.3)

where and are independently distributed matrix-variate gamma variables of (1.25) with the same B. Then the structure in (7.3) is distributed as a product of independently distributed type-1 beta random variables,. Then the density of u can be written as a G-function of the type. The density of will go in terms of a H-function. The H-function is more or less the most generalized special function in real scalar variable case and it is defined by the following Mellin-Barnes integral and the following standard notation is used:

(7.4)

(7.5)

where

(7.6)

where are real and positive numbers,’s and’s are complex numbers. L is a contour separating the poles of to one side and those of to the other side. Existence of the contours and convergence conditions are available from the books Mathai (1993a), Mathai and Saxena (1973, 1978), Mathai and Haubold (2008), Mathai, Saxena and Haubold (2010). When then the H-function reduces to a G-function denoted as

(7.7)

Explicit computable series forms of and for the general, were given by Mathai in a series of papers. The first three forms correspond to product of independently distributed gamma variables, type-1 beta variables and type-2 beta variables respectively. The details of the computable representations are available in the book Mathai (1993a). This is achieved by developing an operator which can handle poles of all orders. This operator may be seen from Mathai & Rathie (1972) and its use from Mathai (1993a). This is a modification of a procedure developed in Mathai and Rathie (1971) to handle generalized partial fractions. Let

(7.8)

where for and the coefficients’s are to be evaluated. The technique developed in Mathai and Rathie (1971) enables one to compute’s explicitly.

The G and H-functions are also established by Mathai for the real matrix-variate cases through M-transforms, along with extensions of all special functions of scalar variables to the matrix-variate cases. Also Mathai extended multivariate functions such as Apple functions, Lauricella functions, Kampé de Fériet functions etc to many matrix-variate cases. Some details may be seen from Mathai (1993a, 1997).

By making use of the explicit series forms, MAPLE and MATHEMATICA have produced computer programs for numerical computations of G-functions and MATHEMATICA has a computer program for the evaluation of H-function also. Solutions of fractional differential equations usually end up in terms of Mittag-Leffler function, its generalization as Wright’s function and its generalization as H-function. In connection with fractional differential equations for reaction, diffusion, reaction-diffusion problems HJH, Mathai and Saxena have given solutions for a large number of situations, which may be seen from the joint works of Haubold, Mathai and Saxena (2011), see also Mathai and Haubold (2018c). In all these solutions, either Mittag-Leffler function or Wright’s function or H-function appears. Also many other physicists, mathematicians and engineers have tried other fractional partial differential equations where also the solutions are available in terms of H-functions.

A Pseudo Dirichlet Integral

A type-1 Dirichlet integral is over a simplex , namely

(7.9)

But one can construct a multivariate integral over a hypercube giving rise to the same where the integrand is different from type-1 Dirichlet format of (7.9). Mathai (2018) constructed such a function which he called it the pseudo Dirichlet function. Consider the following integral:

(7.10)

The method of proving this result is to expand by using a binomial expansion, integrate out variables one by one and then use the properties of Gauss’ hypergeometric function of argument 1 to obtain the result in (7.10). Mathai also extended the integral (7.10) to the real matrix-variate case. In the real matrix-variate case the corresponding integral gives a constant multiple of the form where

where is the real matrix-variate gamma given by the following:

(7.11)

The integral is the following:

where A is a symmetric product of matrices given by

(7.13)

with denoting the positive definite square root of the real positive

definite matrix. The structure in (7.10) gives the same gamma product in (7.9) with replaced by k.

8. Multivariate Statistical Analysis and Statistical Distribution Theory: Fractional Reaction and Diffusion

In the area of multivariate analysis, almost all exact null distributions in the most general cases and a large number of non-null distributions of -criteria for testing hypotheses on one or more multivariate Gaussian populations and exponential populations were given by Mathai and his co-workers. The -criterion is explained in (7.3). Null distributions mean the distributions when the null hypotheses are assumed to hold and non-null distributions mean without the restrictions imposed by the hypotheses. In the non-null situations some of the cases are still open problems. In the null cases, u, a one-to-one function of the -criterion, has usually the following representations:

(8.1)

(8.2)

(8.3)

where are independently distributed real scalar type-1 beta random variables, are the same type of type-2 beta random variables and are the same type of gamma random variables. The density of u in (8.1) can be written in terms of a, that of (8.2) as a and that of (8.3) as a. Computational aspects of these forms are already discussed in Section 7 above. In geometrical probabilities also the squares of the volume content of a p-parallelotope can be written as (8.1) when the random points are type-1 beta distributed, as (8.2) when the random points are type-2 beta distributed and as (8.3) when the random points are gamma distributed. There also densities can be evaluated in terms of the three types of G-functions, as explained above.

Also, Mathai and his co-workers have established a connection between -criterion in testing of statistical hypotheses, connected with multivariate normal populations, and certain generalizations of type-1 Dirichlet models. Various generalizations of type-1 and type-2 Dirichlet models were introduced by Mathai and his co-workers starting with Mathai (2003). In this area also G and H-functions appear. The forms and, coming from products of scalar variables of type-1 and type-2 beta, appear in this area of generalized Dirichlet models.

Exact 11-digit accurate percentage points connected with the null distributions of the -criteria were developed by Mathai and Katiyar starting with the Biometrika paper Mathai and Katiyar (1979). As a byproduct, an algorithm for non-linear least squares was also developed by them, see Mathai and Katiyar (1993a). Mathai has contributions in integer programming and optimization also, see Kounias and Mathai (1988).

Mittag-Leffler Function and Mittag-Leffler Density

HJH, Mathai and Saxena have solved fractional differential equations, starting from 2000, where the solutions invariably come in terms of Mittag-Leffler function, Wright’s function or H-function. Exponential type solutions of integer order differential equations automatically change to Mittag-Leffler functions when we go from integer order to fractional order differential equations. There is also a Mittag-Leffler stochastic process based on a Mittag-Leffler density, which is a non-Gaussian stochastic process. Work in this area is summarized in Haubold, Mathai and Saxena (2011). Mathai has also introduced a generalized Mittag-Leffler density and has shown that it is attracted to heavy-tailed models such as Lévy and Linnik densities, rather than to Gaussian models.

9. Characterization Problems: Gauss and Beyond

In this area, two basic books are Mathai and Rathie (1975) and Mathai and Pederzoli (1977). Characterization is the unique determination through some given properties. Characterization of a density means to show that certain property or properties uniquely determine that density. Unique determination of a concept means to give an axiomatic definition to that concept. That is, to show that the proposed axioms will uniquely determine the concept. The techniques used in this area, to go from the given properties to the density or from the given axioms to the concept such as “uncertainty’’ or its complement “information’’ etc, are functional equations, differential equations, Laplace, Mellin, and Fourier transforms. For example, look at the distribution of error. The error may be the error in measurement in an experiment, the error between observed and predicted values etc. If the factors contributing to the error are known then the experimenter will try to control these factors. Very often the error is contributed by infinitely many unknown factors each factor contributing infinitesimal quantities towards. Put some conditions on this. Let

(i)

where are assumed to be independently distributed. Assume that each or -a with equal probabilities. That is,

(ii)

Assume that the total variance of is finite or

(iii)

Check the consequence of these three assumptions., where a is fixed and finite. For large n one may take. The moment generating function of is

Hence

That is

Therefore, as,

which is the moment generating function of a normal density with mean value zero and variance or the density is

(9.1)

This is the derivation of the Gaussian or normal density given by Gauss, and hence it is also called the Gaussian density or error curve. Mathai and Pederzoli (1977) contains such characterizations of the normal probability law by using structural properties, regression properties etc. One fundamental idea was introduced in this area by Gordon and Mathai (1972). They tried to come up with pseudo-analytic functions of matrix argument involving rectangular matrices and by using this, characterization theorems were established for a multivariate normal density.

In Mathai and Rathie (1975), axiomatic definitions of information theory measures and basic statistical concepts are given. This is the first book giving axiomatic definitions of information measures. The techniques used are mainly from functional equations by using the proposed axioms create a functional equation and obtain its unique solution by imposing more conditions, if necessary, thus coming up with a unique definition or characterization of the concept. One such measure there is the one introduced as Havrda-Charvát measure, which for the continuous case is the following:

(9.2)

where is a density of the real scalar variable x. There is a corresponding discrete analogue, which is given by

(9.3)

where. A modified form of (9.2) and (9.3) is Tsallis entropy given by

(9.4)

for the continuous case, with a corresponding discrete analogue. Optimization of (9.4) under the condition that the total energy is preserved or the first moment is fixed, leads to Tsallis statistics of non-extensive statistical mechanics. Tsallis statistics is of the following form:

(9.5)

which is also a power law in the sense. Note that a direct optimization of (9.4), under the assumption that the first moment in is fixed, does not yield (9.5) directly. One has to go through an escort density

and then assume that the first moment is fixed in the escort density, to get the form in (9.5). Mathai’s entropy

(9.6)

when optimized under the condition of first moment in being fixed leads to Tsallis statistics directly. Also the optimization of (9.6) under two moment-type conditions leads to the pathway model, discussed in Section 3, where (9.5) will be a particular case.

10. Biological Modeling: Formation of Pattern

The most significant contribution in this area is the proposal of a theory of growth and form in nature and the explanation of the emergence of beautiful patterns in sunflower, along with explanation for the appearance of Fibonacci sequence and golden ratio there. The mathematical reconstruction of the sunflower head, with all the features that are seen in nature, is still the cover design of the journal of Mathematical Biosciences. The paper of Mathai and Davis appeared in that journal in 1974 and in 1976 the journal adopted the mathematically reconstructed sunflower head of Mathai and Davis (1974) as the journal’s cover design with acknowledgement to the authors. When this paper was sent for publication to this journal, the editor wrote back saying: “enthusiastically accepted for publication’’ because this was the first time all natural features were explained in full. As per Davis and Mathai, the programming of the sunflower head is like a point moving along an Archimedes’ spiral at a constant speed so that when the point makes an angle, a second point starts and moves at the same speed. When the second point comes to, a third point starts, and so on. The rule governing the movement is or

where k is a constant, giving Archimedes’ spiral. When or, one obtains sunflower, coconut tree crown, certain cactus heads and so on. Such a movement can be generated by a viscous fluid flowing up through a capillary with valves so that when a certain pressure is built up in one chamber the liquid moves up to the next chamber. The continuous flow is made pulse-like at the end. The upward motion can be effected by an evaporation process in the leaves, and there is no need for a heart-like mechanism in trees, pumping the fluid up. Mathai and Davis (1973) showed that the arrangments of leaves on a coconut tree crown is ideal from many mathematical points of view.

11. Design of Experiments and Analysis of Variance

The first paper of Mathai (Mathai, 1965), was on an approximate analysis of variance. It was on the analysis of a two-way classification with multiple observations per cell. Here the orthogonality is lost, and when estimating the main effects, one ends up in a singular system of linear equations of the form

(11.1)

where for all i and j, is called the incidence matrix and the sum of the elements in each row is equal to 1. Thus is singular and hence one cannot write it as where A and b are known and the vector is unknown and is to be estimated. Mathai noted that one could profitably use the conditions in the design and make a nonsingular matrix. One condition in the design is that where are the elements in. Let C be a matrix where all elements in the i-th row of C are the median of the i-th row elements in A, namely the median of for. Then evidently (null). Then

(11.2)

where and is the median of the i-th row elements in A. Then a norm of B is. But since the mean deviation from the median is the least, is the minimum under the circumstances. Therefore, not only that is nonsingular but the series is the fastest converging series for the problem at hand. Then

A good approximation for is available as. This is found to be sufficient for all practical purposes of testing of statistical hypotheses on the components of.

12. Population Problems and Social Sciences

A problem that was looked into was how to come up with a measure of “distance’’ or “closeness’’ or “affinity’’ between two sociological groups or how to say that one community is close to another community with respect to a given characteristic. Mathai introduced the concepts of “directed divergence’’, “affinity’’ etc from information theory to social statistics. Let and be two discrete populations. Consider the representation of P and Q as points on a hypersphere of radius 1,. Then the points are and. Consider where is the angle between these vectors or points on the hypersphere. Note that the lengths of the vectors are

and hence

(12.1)

This is a measure of angular dispersion and it is usually called “affinity’’ between P and Q or Matusita’s measure of affinity between two discrete distributions. George and Mathai computed “affinity” between communities with reference to the characteristic of production of children and found that the politicians’ statements did not match with the realities. Thus, some politicians’ claims of certain communities producing more children, was repudiated in a scientific way in George and Mathai (1974). They also studied the most important variable responsible for population growth, namely the interval between two live births in woman of child-bearing age group and proposed a model, George and Mathai (1975). They also gave a formula for estimating an event from information supplied by different agencies, replacing the popular Deming formula in this regard.

13. Quadratic and Bilinear Forms

Major contributions in these areas are summarized in the books Mathai and Provost (1992), and Mathai, Provost and Hayakawa (1995). There is a very important concept in quadratic forms in Gaussian random variables called chisquaredness of quadratic forms. That is, if and only if A is idempotent and of rank r, where X the vector having the standard normal distribution, that is, and is a chisquare random variable with r degrees of freedom. Is there a corresponding concept when dealing with bilinear forms? When the samples come from a bivariate Gaussian or normal population it is not difficult to work out the density of the sample correlation coefficient. But what about the density of the sample covariance, without the scaling factors of the standard deviations? Both these questions were answered by Mathai (1993c) where Mathai introduced a concept called Laplacianness of bilinear forms (Mathai, 1993b) and also worked out the density of the covariance structures observing that covariance structure is a bilinear form. The necessary and sufficient conditions for a bilinear form to be noncentral generalized Laplacian are given in Corollary 2.5.2 of Mathai, Provost and Hayakawa (1995). For a noncentral generalized Laplacian the moment generating function is of the form

where is the non-centrality parameter.

14. Reliability Analysis: Extension to Pathway Model and Matrix-Variate Case

In the area of reliability analysis, the basic concepts are survival function, hazard function, cumulative hazard, system reliability, reliability in the presence of other variables such as covariates etc. In a series of papers, Mathai and Princy in 2016 introduced the pathway model into the area so that the desired shapes for hazard function and the desired reliability for systems with components in series and parallel architecture could be obtained by selecting appropriate models from the pathway family of functions. Then, these ideas were extended to situations where the input variable or the variable under consideration is a rectangular matrix. As a byproduct, Maxwell-Boltzmann distribution, Raleigh distribution, Dirichlet averages etc were extended to matrix-variate cases, see for example Mathai and Princy (2017a, 2017b).

15. Mellin Convolutions of Products and Ratios and M-Convolutions

Mellin convolutions of products and ratios involving two functions are available in the literature. Mathai illustrated how these concepts are connected to statistical distribution theory and fractional calculus. In fact, a general definition for fractional integrals is given by Mathai using Mellin convolutions of products and ratios involving two functions. Corresponding M-convolutions involving two functions of matrix argument is Mathai’s contribution. He has also given physical interpretations for M-convolutions as densities of symmetric products and symmetric ratios of matrices. Mathai also extended Mellin convolutions and M-convolutions to three or more functions, see Mathai (2018). When three functions are involved, one can obtain several integral representations for the same Mellin convolution and M-convolution. For example, consider the symmetric product of three real symmetric positive definite matrices. One can take symmetric products as, , etc. When the original densities are assumed to be functionally symmetric then all such symmetric forms will produce the same densities whereas each symmetric product produces an integral representation which will be all different. Thus, one gets a large number of different integral representations for the same density or M-convolution of a product.

One can also obtain further representations by taking for example, , and consider the original symmetric product of three matrices as

symmetric products of two matrices each, which will produce several more integral representations. In the real scalar case the Mellin convolutions can be evaluated in terms of generalized special functions, thus producing integral representations for these special functions. For example, let be real scalar positive random variables, independently distributed and let the product.

where will be gamma products when’s have densities belonging to the pathway family of functions, namely type-1 beta, type-2 beta and gamma. Then the inverse of is a G-function. Then this G-function has several different types of integral representations. A series of papers are written by Mathai by using these ideas connecting statistical distributions, fractional calculus, Mellin convolutions and M-convolutions, see for example Mathai (2017, 2018, 2019).

16. Above Topics Itemized: New Concepts and Procedures

The following are the new concepts, new ideas and new procedures introduced by Mathai (Figure 16):

- Developed “dispersion theory’’ in 1967;

- Developed a generalized partial fraction technique, with Rathie, in 1971;

- Developed an operator to evaluate residues when poles of all types of orders occur (1971);

- Introduced the phrases “statistical sciences’’ in 1971 thereby the phrase “mathematical sciences’’ came into existence;

- Proposed a theory of growth and forms in nature (1974), the theory still standed, mathematically reconstructed a sunflower head;

- Introduced the concepts of “affinity’’, “distance’’ etc. in social sciences and created a procedure to compare sociological groups (1974);

- Introduced functions of matrix argument through M-transforms and M-convolutions;

- Introduced a non-linear least square algorithm (1993);

- Solved Miles’ conjecture in geometrical probabilities, created and solved parallel conjectures (1982);

- Introduced Jacobians of matrix transformations in solving problems of random volumes, replacing the complicated integral and differential geometry procedures (1982);

- Now meaningful physical interpretations are given for M-convolutions; Unique recovery of from its M-transform is still a conjecture;

Figure 16. H.-J. Treder.

- Extended Jacobians of matrix transformations from the real case to complex matrix-variate cases in a large number of situations;

- Introduced the concept of Laplacianness of bilinear forms and established the density of covariance structures (1993);

- Introduced pathway model and pathway idea (2005);

- Extended fractional calculus to real matrix-variate cases (2007);

- Established a connection between fractional calculus and statistical distribution theory (2007);

- Introduced Mathai’s (2007) entropy;

- Geometrical interpretation and a general definition for fractional integrals were given (2013-2015);

- Extended fractional calculus to complex matrix-variate case and complex domain in general (2013);

- Extended fractional calculus to many matrix-variate cases (real and complex) (2014);

- Developed a fractional differential operator in the matrix-variate case (2015);

- Extended reliability analysis concepts to rectangular matrix-variate case (2017).

Acknowledgements

HJH expresses his deep appreciation for a lifelong support from and cooperation with Prof. Dr. A. M. Mathai, Department of Mathematics and Statistics, McGill University, Montreal, Canada, and Director of the Centre for Mathematical and Statistical Sciences, Peechi, Kerala, India. HJH also takes the opportunity to place on record his gratefulness for the encouragment of research by AkM Prof. Dr. Dr. e.h. mult. H.-J. Treder, Director of the Einstein Laboratory for Theoretical Physics, Caputh, Germany (see Schulz-Fieguth, 2018). Treder was the director of the Central Institute for Astrophysics of the Academy of Sciences (Berlin, GDR). In 1965 he was the principal organizer of the widely respected Einstein Symposium at the 50th anniversary of the invention of general relativity theory. For the Berlin international celebrations of Einstein’s 100th birthday, 1979, he managed to secure the summer house of Einstein in Caputh, Brandenburg, as the Einstein Laboratory of Theoretical Physics in consultation with the administrators of the estate of Otto Nathan and Einstein. In 1981 he hosted the Michelson Colloquium at Potsdam to celebrate and recall the first Michelson experiment performed in 1881 at the Astrophysical Observatory in Potsdam. Treder was able to secure space and time for intense research work in his professional environment ranging from the solar neutrino problem (Treder, 1974) to fractional calculus (Treder, 1989). He supported actively the United Nations efforts to make available education and research in science to nations worldwide.

Conflicts of Interest

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

References

[1]	Beck, C. (2004). Superstatistics: Theory and Applications. Continuum Mechanics and Thermodynamics, 16, 293-304. https://doi.org/10.1007/s00161-003-0145-1
[2]	Cohen, E. G. D. (2004). Superstatistics. Physica D: Nonlinear Phenomena, 193, 35-52. https://doi.org/10.1016/j.physd.2004.01.007
[3]	Constantine, A. G. (1963). Some Noncentral Distribution Problems in Multivariate Analysis. Annals of Mathematical Statistics 34, 1270-1285. https://doi.org/10.1214/aoms/1177703863
[4]	Critchfield, C. L. (1972). Analytic Forms of the Thermonuclear Function. In F. Reines (Ed.), Cosmology, Fusion & Other Matters: George Gamow Memorial Volume (pp. 186-191). Colorado: Colorado Associated University Press.
[5]	Davis Jr., R. (2003). A Half Century with Solar Neutrinos. Reviews of Modern Physics, 75, 985-994. https://doi.org/10.1103/RevModPhys.75.985
[6]	Fowler, A. W. (1984). Experimental and Theoretical Nuclear Astrophysics: The Quest for the Origin of the Elements. Reviews of Modern Physics, 56, 149-179. https://doi.org/10.1103/RevModPhys.56.149
[7]	George, A., & Mathai, A. M. (1974). Applications of the Concepts of Affinity and Distance to Population Problems. Journal of Biosocial Sciences (London), 6, 347-356. https://doi.org/10.1017/S0021932000009706
[8]	George, A., & Mathai, A. M. (1975). Distribution of Inter-Live-Birth Intervals. Sankhya Series B, 37, 332-342.
[9]	Gordon, F. S., & Mathai, A. M. (1972). Characterization of Multivariate Normal Distribution using Regression Properties. Annals of Mathematical Statistics, 43, 53-65. https://doi.org/10.1214/aoms/1177692714
[10]	Haubold, H. J., & Mathai, A. M. (1988). Gravitational Instability in a Multicomponent Cosmological Medium. Journal of Mathematical Physics, 29, 2069-2077. https://doi.org/10.1063/1.527864
[11]	Haubold, H. J., & Mathai, A. M. (1997). Sun. In J. H. Shirley, & R. W. Fairbridge (Eds.), Encyclopaedia of Planetary Sciences (pp. 786-794). London: Chapman and Hall. https://doi.org/10.1007/1-4020-4520-4_390
[12]	Haubold, H. J., & Mathai, A. M. (1998). Structure of the Universe,. In G. L. Trigg (Ed.), Encyclopedia of Applied Physics (pp. 47-51).
[13]	Haubold, H. J., & Mathai, A. M. (Eds.) (2010). Proceedings of the Third UN/ESA/NASA Workshop on the International Heliophysical Year 2007 and Basic Space Science. Berlin and Heidelberg: Springer. https://doi.org/10.1007/978-3-642-03325-4
[14]	Haubold, H. J., Mathai, A. M., & Saxena, R. K. (2011). Mittag-Leffler Functions and Their Applications. Journal of Applied Mathematics, 2011, Article ID: 298628. https://doi.org/10.1155/2011/298628
[15]	Herz, C. S. (1955). Bessel Functions of Matrix Argument. Annals of Mathematics, 63, 474-523.
[16]	James, A. T. (1961). Zonal Polynomials of the Real Positive Definite Symmetric Matrices. Annals of Mathematics, 74, 456-469. https://doi.org/10.2307/1970291
[17]	Kounias, S., & Mathai, A. M. (1988). Maximizing the Sum of Integers When Their Sum of Squares Is Fixed. Optimization, 19, 123-131. https://doi.org/10.1080/02331938808843325
[18]	Kratzel, E. (1979). Integral Transforms the Bessel Type. In Generalized Functions & Operational Calculus, Proc. Conf. Verna, 1975 (pp. 148-165). Sofia: Bulg. Acad. Sci.
[19]	Mathai, A. M. (1965). An Approximate Analysis of a Two-Way Layout. Biometrics, 21, 376-385. https://doi.org/10.2307/2528097
[20]	Mathai, A. M. (1982). On a Conjecture in Geometric Probability Regarding Asymptotic Normality of a Random Simplex. Annals of Probability, 10, 247-251. https://doi.org/10.1214/aop/1176993929
[21]	Mathai, A. M. (1989). On a System of Differential Equations Connected with the Gravitational Instability in a Multi-Component Medium in Newtonian Cosmology. Studies in Applied Mathematics, 80, 75-93. https://doi.org/10.1002/sapm198980175
[22]	Mathai, A. M. (1993a). A Handbook of Generalized Special Function for Statistical and Physical Sciences. Oxford: Oxford University Press.
[23]	Mathai, A. M. (1993b). On Generalized Laplacianness of Bilinear Forms in Normal Variables. Journal of Multivariate Analysis, 45, 239-246. https://doi.org/10.1006/jmva.1993.1036
[24]	Mathai, A. M. (1993c). The Residual Effect of a Growth-Decay Mechanism and the Distribution of Covariance Structures. Canadian Journal of Statistics, 21, 277-283. https://doi.org/10.2307/3315753
[25]	Mathai, A. M. (1997). Jacobians of Matrix Transformations and Functions of Matrix Argument. New York: World Scientific Publishing. https://doi.org/10.1142/3438
[26]	Mathai, A. M. (1998). Pollution by Vehicular Travels from the Suburbs to the City Core. Environmetrics, 9, 617-628. https://doi.org/10.1002/(SICI)1099-095X(199811/12)9:6<617::AID-ENV326>3.0.CO;2-E
[27]	Mathai, A. M. (1999a). An Introduction to Geometrical Probability: Distributional Aspects with Applications. Newark: Gordon and Breach.
[28]	Mathai, A. M. (1999b). Random p-Parallelotope in Euclidean n-Space. Advances in Applied Probability, 31, 343-354. https://doi.org/10.1239/aap/1029955138
[29]	Mathai, A. M. (1999c). Buffon Type Problem in the Euclidean Space R3. Rendiconti del Circolo Matematico di Palermo, XLVIII, 487-506. https://doi.org/10.1007/BF02844338
[30]	Mathai, A. M. (2003). Order Statistics from a Logistic Distribution and Applications to Survival & Reliability Analysis. IEEE Transactions on Reliability, 52, 200-206. https://doi.org/10.1109/TR.2003.813432
[31]	Mathai, A. M. (2005). A Pathway to Matrix-Variate Gamma and Normal Densities. Linear Algebra and Its Applications, 396, 317-328. https://doi.org/10.1016/j.laa.2004.09.022
[32]	Mathai, A. M. (2007). Random Volumes under a General Matrix-Variate Model. Linear Algebra and its Applications, 425, 162-170. https://doi.org/10.1016/j.laa.2007.03.022
[33]	Mathai, A. M. (2009). Fractional Integrals in the Matrix-Variate Cases and Connection to Statistical Distributions. Integral Transforms and Special Functions, 20, 871-882. https://doi.org/10.1080/10652460902867791
[34]	Mathai, A. M. (2012). Generalized Kratzel Integral and Associated Statistical Densities. International Journal of Mathematical Analysis, 6, 2501-2510.
[35]	Mathai, A. M. (2013). Fractional Integral Operators in the Complex Matrix-Variate Case. Linear Algebra and its Applications, 439, 2901-2913. https://doi.org/10.1016/j.laa.2013.08.023
[36]	Mathai, A. M. (2014). Fractional Integral Operators Involving Many Matrix Variables. Linear Algebra and its Applications, 446, 196-215. https://doi.org/10.1016/j.laa.2014.01.003
[37]	Mathai, A. M. (2015). Fractional Differential Operators in the Complex Matrix-Variate Case. Linear Algebra and its Applications, 478, 200-217. https://doi.org/10.1016/j.laa.2015.03.038
[38]	Mathai, A. M. (2017). A Variant of Type-1 Beta and Dirichlet Integrals. Journal of the Ramanujan Society of Mathematics and Mathematical Sciences, 6, 25-40.
[39]	Mathai, A. M. (2018). Mellin Convolutions, Statistical Distributions and Farctional Calculus. Fractional Calculus & Applied Analysis, 21, 376-398. https://doi.org/10.1515/fca-2018-0022
[40]	Mathai, A. M. (2019). Extended Fractional Integrals, Products and Ratios of Matrices and Statistical Distributions. 67th Annual Conference of Bharat Ganita Parishad, Lucknow, India, November 2019.
[41]	Mathai, A. M., & Davis, T. A. (1973). A Mathematical Explanation of the Emergence of Foliar Spirals in Coconut and Other Palm Trees. Proceedings of the National Science Academy (India), 39, 194-202.
[42]	Mathai, A. M., & Davis, T. A. (1974). Constructing a Sunflower Head. Mathematical Biosciences, 20, 117-133. https://doi.org/10.1016/0025-5564(74)90072-8
[43]	Mathai, A. M., & Haubold, H. J. (1988). Modern Problems in Nuclear and Neutrino Astrophysics. Berlin: Akademie-Verlag.
[44]	Mathai, A. M., & Haubold, H. J. (2008). Special Functions for Applied Scientists. New York: Springer. https://doi.org/10.1007/978-0-387-75894-7
[45]	Mathai, A. M., & Haubold, H. J. (2017a). Linear Algebra for Physicists and Engineers. Berlin: De Gruyter. https://doi.org/10.1515/9783110562507
[46]	Mathai, A. M., & Haubold, H. J. (2017b). Probability and Statistics for Physicists and Engineers. Berlin: De Gruyter. https://doi.org/10.1515/9783110562545
[47]	Mathai, A. M., & Haubold, H. J. (2017c). Fractional and Multivariable Calculus: Model Building and Optimization Problems. Switzerland: Springer. https://doi.org/10.1007/978-3-319-59993-9
[48]	Mathai, A. M., & Haubold, H. J. (2018a). United Nations Basic Space Science Initiative (UNBSSI) 1991-2012 and beyond. Creative Education, 9, 192-248. https://doi.org/10.4236/ce.2018.92015
[49]	Mathai, A. M., & Haubold, H. J. (2018b). Matrix Methods and Fractional Calculus. London: World Scientific Publishing. https://doi.org/10.1142/10639
[50]	Mathai, A. M., & Haubold, H. J. (2018c). An Introduction to Fractional Calculus. New York: Nova Science Publishers. https://doi.org/10.1142/10639
[51]	Mathai, A. M., & Haubold, H. J. (2018d). Erdélyi-Kober Fractional Calculus: From a Statistical Perspective, Inspired by Solar Neutrino Physics. In Springer Briefs in Mathematical Physics Volume 31. Singapore: Springer. https://doi.org/10.1007/978-981-13-1159-8
[52]	Mathai, A. M., & Katiyar, R. S. (1979). Exact Percentage Points for Testing Independence. Biometrika, 66, 353-356. https://doi.org/10.1093/biomet/66.2.353
[53]	Mathai, A. M., & Katiyar, R. S. (1993). A New Algorithm for Nonlinear Least Squares. Researches in Mathematical Statistics (Russian), 207, 143-157. (Translation by the American Mathematical Society)
[54]	Mathai, A. M., & Moschopoulos, P. G. (1999). Pollution by Vehicular Travels from Satellite Townships to the City. Environmetrics, 10, 791-802. https://doi.org/10.1002/(SICI)1099-095X(199911/12)10:6<791::AID-ENV396>3.0.CO;2-R
[55]	Mathai, A. M., & Pederzoli, G. (1977). Characterization of the Normal Probability Law. New Delhi: Wiley Eastern; and New York: Wiley Halsted.
[56]	Mathai, A. M., & Princy, T. (2017a). Maxwell-Boltzmann and Rayleigh Distributions for the Rectangular Matrix-Variate Case. Physica A, 532, 287-311.
[57]	Mathai, A. M., & Princy, T. (2017b). Analogues of Reliability Analysis for the Matrix-Variate Case. Linear Algebra and Its Applications, 532, 287-311. https://doi.org/10.1016/j.laa.2017.06.022
[58]	Mathai, A. M., & Provost, S. B. (1992). Quadratic Forms in Random Variables: Theory and Applications. New York: Marcel Dekker.
[59]	Mathai, A. M., & Provost, S. B. (2005). Some Complex Matrix-Variate Statistical Distributions on Rectangular Matrices. Linear Algebra and its Applications, 410, 198-216. https://doi.org/10.1016/j.laa.2005.07.016
[60]	Mathai, A. M., & Provost, S. B. (2006). On q-Logistic and Related Distributions. IEEE Transactions on Reliability, 55, 237-244. https://doi.org/10.1109/TR.2006.874927
[61]	Mathai, A. M., & Rathie, P. N. (1971). Exact Distribution of Wilks’ Criterion. Annals of Mathematical Statistics, 42, 1010-1019. https://doi.org/10.1214/aoms/1177693329
[62]	Mathai, A. M., & Rathie, P. N. (1972). Enumeration of Almost Cubic Maps. Journal of Combinatorial Theory, 13, 83-90. https://doi.org/10.1016/0095-8956(72)90012-3
[63]	Mathai, A. M., & Rathie, P. N. (1975). Basic Concepts in Information Theory and Statistics: Axiomatic Foundations and Application. New York: Wiley Halsted; and New Delhi: Wiley Eastern.
[64]	Mathai, A. M., & Saxena, R. K. (1969). Distribution of a Product and Structural Setup of Densities. Annals of Mathematical Statistics, 40, 1439-1448. https://doi.org/10.1214/aoms/1177697515
[65]	Mathai, A. M., & Saxena, R. K. (1973). Generalized Hypergeometric Functions with Applications in Statistics and Physical Sciences. Heidelberg: Springer. https://doi.org/10.1007/BFb0060468
[66]	Mathai, A. M., & Saxena, R. K. (1978). The H-Function with Applications in Statistics and Other Disciplines. New York: Wiley Halsted; and New Delhi: Wiley Eastern.
[67]	Mathai, A. M., Provost, S. B., & Hayakawa, T. (1995). Bilinear Forms and Zonal Polynomials. In Lecture Notes in Statistics 102. New York: Springer. https://doi.org/10.1007/978-1-4612-4242-0
[68]	Mathai, A. M., Saxena, R. K., & Haubold, H. J. (2010). The H-Function: Theory and Applications. New York: Springer. https://doi.org/10.1007/978-1-4419-0916-9
[69]	Pyenson, L., Mathai, A. M., & Haubold, H. J. (2019). United Nations Education Program in Space Science and Technology 1988-2018. Creative Education, 10, 2219-2231. https://doi.org/10.4236/ce.2019.1010160
[70]	Sakurai, K. (2014). Solar Neutrino Problems—How They Were Solved. Tokyo: Terrapub.
[71]	Schulz-Fieguth, M. (2018). Hans-Jürgen Treder—Ein Portrat. Potsdam, Germany: Christian and Cornelius Rüss.
[72]	Thomas, S., & Mathai, A. M. (2009). P-Content of p-Parallelotope and Its Connection to Likelihood Ratio Statistics. Sankhya Series A, 71, 49-63.
[73]	Treder, H.-J. (1974). Gravitation und weitreichende schwache Wechselwirkungen bei Neutrino-Feldern (Gedanken zu einer Theorie der solaren Neutrinos). Astronomische Nachrichten, 295, 169-184. https://doi.org/10.1002/asna.19742950403
[74]	Treder, H.-J. (1989). Fraktale und gebrochene Ableitungen. Annalen der Physik, 46, 74-75. https://doi.org/10.1002/andp.19895010110