\documentclass[reqno]{amsart}
\usepackage{setspace,amssymb}%,amsaddr
%\doublespacing%\onehalfspacing%
\usepackage{ifpdf}
\ifpdf
 \usepackage[hyperindex,pagebackref]{hyperref}%
\else
 \expandafter\ifx\csname dvipdfm\endcsname\relax
 \usepackage[hypertex,hyperindex,pagebackref]{hyperref}
 \else
 \usepackage[dvipdfm,hyperindex,pagebackref]{hyperref}
 \fi
\fi
\allowdisplaybreaks[4]
\numberwithin{equation}{section}
\theoremstyle{plain}
\newtheorem{theorem}{Theorem}[section]
\theoremstyle{definition}
\newtheorem{definition}{Definition}[section]
\theoremstyle{remark}
\newtheorem{remark}{Remark}[section]
\DeclareMathOperator{\td}{d}
\newcommand{\cmdeg}[1]{\sideset{}{_\mathrm{cm}^{#1}}\deg}

\begin{document}

\title[Complete monotonicity and completely monotonic degree]
{Complete monotonicity, completely monotonic degree, integral representations, and an inequality related to the exponential, trigamma, and modified Bessel functions}

\author[F. Qi]{Feng Qi}
\address[Feng Qi]{College of Mathematics, Inner Mongolia University for Nationalities, Tongliao City, Inner Mongolia Autonomous Region, 028043, China}
\email{\href{mailto: F. Qi <qifeng618@gmail.com>}{qifeng618@gmail.com}, \href{mailto: F. Qi <qifeng618@hotmail.com>}{qifeng618@hotmail.com}, \href{mailto: F. Qi <qifeng618@qq.com>}{qifeng618@qq.com}} \urladdr{\url{http://qifeng618.wordpress.com}}

\author[S.-H. Wang]{Shu-Hong Wang}
\address[Shu-Hong Wang]{College of Mathematics, Inner Mongolia University for Nationalities, Tongliao City, Inner Mongolia Autonomous Region, 028043, China}
\email{\href{mailto: S.-H. Wang <shuhong7682@163.com>}{shuhong7682@163.com}}

\begin{abstract}
In the paper, the authors verify the complete monotonicity of the difference $e^{1/t}-\psi'(t)$ on $(0,\infty)$, compute the completely monotonic degree and establish integral representations of the remainder of the Laurent series expansion of $e^{1/z}$, and derive an inequality which gives a lower bound for the first order modified Bessel function of the first kind. These results show us some new properties and relations of the exponential, trigamma, the first kind modified Bessel functions and the hypergeometric series.
\end{abstract}

\keywords{property; connection; completely monotonic function; completely monotonic degree; integral representation; difference; exponential function; trigamma function; hypergeometric series; inequality; modified Bessel function}

\subjclass[2010]{Primary 26A48, 33B10, 33C10; Secondary 26A12, 26D07, 30E20, 33B15, 33C20, 44A10}

\thanks{Tel: 86 22 83956502}

\maketitle

\section{Introduction}

In~\cite[Lemma~2]{Yang-Fan-2008-Dec-simp.tex}, the inequality
\begin{equation}\label{e-1-t-1}
\psi'(t)<e^{1/t}-1
\end{equation}
on $(0,\infty)$ was obtained and applied, where $\psi(t)$ stands for the digamma function which may be defined by the logarithmic derivative
\begin{equation*}
\psi(t)=[\ln\Gamma(t)]'=\frac{\Gamma'(t)}{\Gamma(t)}
\end{equation*}
and $\Gamma(t)$ is the classical Euler gamma function which may be defined for $\Re z>0$ by
\begin{equation*}%\label{gamma-dfn}
\Gamma(z)=\int^\infty_0t^{z-1} e^{-t}\td t.
\end{equation*}
The derivatives $\psi'(z)$ and $\psi''(z)$ of $\psi(z)$ are respectively called the tri- and tetra-gamma functions. As a whole, the derivatives $\psi^{(k)}(z)$ for $k\in\{0\}\cup\mathbb{N}$ are called the polygamma functions.
\par
The first aim of this paper is to generalize the inequality~\eqref{e-1-t-1} to complete monotonicity of a difference between both sides of~\eqref{e-1-t-1}.
\par
Our first result can be formulated as Theorem~\ref{CM-Exp-thm} below.

\begin{theorem}\label{CM-Exp-thm}
The function
\begin{equation}\label{alpha-exp=psi-eq}
h(t)=e^{1/t}-\psi'(t)
\end{equation}
is completely monotonic, that is, $(-1)^{k-1}h^{(k-1)}(t)\ge0$ for $k\in\mathbb{N}$, on $(0,\infty)$ and
\begin{equation}\label{h(t)-limit=1}
\lim_{t\to\infty}h(t)=1.
\end{equation}
\end{theorem}

The notion ``completely monotonic degree'' was innovated in~\cite[Definition~1]{psi-proper-fraction-degree-two.tex} which may be regarded as a slight but essential modification of~\cite[Definition~1.5]{Koumandos-Pedersen-09-JMAA}. This may be further fixed as follows.

\begin{definition}\label{x-degree-dfn}
Let $f(x)$ be a completely monotonic function on $(0,\infty)$ and denote $f(\infty)=\lim_{x\to\infty}f(x)$. If for some $r\in\mathbb{R}$ the function $x^r[f(x)-f(\infty)]$ is completely monotonic on $(0,\infty)$ but $x^{r+\varepsilon}[f(x)-f(\infty)]$ is not for any positive number $\varepsilon>0$, then we say that the number $r$ is the completely monotonic degree of $f(x)$ with respect to $x\in(0,\infty)$; if for all $r\in\mathbb{R}$ each and every $x^r[f(x)-f(\infty)]$ is completely monotonic on $(0,\infty)$, then we say that the completely monotonic degree of $f(x)$ with respect to $x\in(0,\infty)$ is $\infty$.
\end{definition}

For convenience, the notation
$\cmdeg{x}[f(x)]$
was designed in~\cite[p.~9890]{psi-proper-fraction-degree-two.tex} to denote the completely monotonic degree $r$ of $f(x)$ with respect to $x\in(0,\infty)$.
\par
It is clear that the completely monotonic degree $\cmdeg{x}[f(x)]$ of any completely monotonic function $f(x)$ on $(0,\infty)$ is at leat $0$.
It may be proved that the completely monotonic degree $\cmdeg{x}[f(x)]$ equals $\infty$ if and only if $f(x)$ is nonnegative and identically constant.
\par
The second aim of this paper is to compute the completely monotonic degree and to establish integral representations of the remainder of the Laurent series expansion of the exponential function $e^{1/z}$.
\par
Our second result may be stated as the following theorem.

\begin{theorem}\label{exp=k=degree=k+1-int-thm}
For $k\in\{0\}\cup\mathbb{N}$ and $z\ne0$, let
\begin{equation}\label{exp=k=sum-eq-degree=k+1}
H_k(z)=e^{1/z}-\sum_{m=0}^k\frac{1}{m!}\frac1{z^m}.
\end{equation}
\begin{enumerate}
\item
The completely monotonic degree of $H_k(t)$ on $(0,\infty)$ meets
\begin{equation}\label{H-k(t)-degree}
\cmdeg{t}[H_k(t)]=k+1.
\end{equation}
\item
For $\Re z>0$, the function $H_k(z)$ has the integral representations
\begin{equation}\label{exp=k=degree=k+1-int}
H_k(z)=\frac1{k!(k+1)!}\int_0^\infty {}_1F_2(1;k+1,k+2;t)t^k e^{-zt}\td t
\end{equation}
and
\begin{equation}\label{exp=k=degree=k+1-int-bes}
H_k(z)=\frac1{z^{k+1}}\biggl[\frac1{(k+1)!}+\int_0^\infty \frac{I_{k+2} \bigl(2\sqrt{t}\,\bigr)}{t^{(k+2)/2}} e^{-zt}\td t\biggr],
\end{equation}
where the hypergeometric series
\begin{equation}
{}_pF_q(a_1,\dotsc,a_p;b_1,\dotsc,b_q;x)=\sum_{n=0}^\infty\frac{(a_1)_n\dotsm(a_p)_n} {(b_1)_n\dotsm(b_q)_n}\frac{x^n}{n!}
\end{equation}
for $b_i\notin\{0,-1,-2,\dotsc\}$, the shifted factorial $(a)_0=1$ and
\begin{equation}
(a)_n=a(a+1)\dotsm(a+n-1)
\end{equation}
for $n>0$ and any real or complex number $a$, and the modified Bessel function of the first kind
\begin{equation}\label{I=nu(z)-eq}
I_\nu(z)= \sum_{k=0}^\infty\frac1{k!\Gamma(\nu+k+1)}\biggl(\frac{z}2\biggr)^{2k+\nu}
\end{equation}
for $\nu\in\mathbb{R}$ and $z\in\mathbb{C}$.
\end{enumerate}
\end{theorem}

As an application of Theorems~\ref{CM-Exp-thm} and~\ref{exp=k=degree=k+1-int-thm}, the following inequality for the first order modified Bessel function of the first kind $I_1$ may be derived.

\begin{theorem}\label{Bessel-2-ineq-thm}
For $t>0$, we have
\begin{equation}\label{I=1-exp-ineq}
I_1(t)>\frac{(t/2)^3}{1-e^{-(t/2)^2}}.
\end{equation}
\end{theorem}

\section{Proofs of Theorem~\ref{CM-Exp-thm}}

We supply two proofs of Theorem~\ref{CM-Exp-thm}.

\begin{proof}[First proof]
From the well known formula
\begin{equation}\label{polygamma}
\psi^{(n)}(z)=(-1)^{n+1}\int_0^{\infty}\frac{u^n}{1-e^{-u}}e^{-zu}\td u
\end{equation}
for $\Re z>0$ and $n\in\mathbb{N}$, see~\cite[p.~260, 6.4.1]{abram}, it is ready that $\lim_{t\to\infty}\psi^{(n)}(t)=0$ for $n\in\mathbb{N}$. So, the limit~\eqref{h(t)-limit=1} may be deduced immediately and, by
\begin{equation}\label{exp-frac1x-expans}
\bigl(e^{1/t}\bigr)^{(i)}=(-1)^ie^{1/t}\frac1{t^{2i}}\sum_{k=0}^{i-1}a_{i,k}t^{k}
\end{equation}
for $i\in\mathbb{N}$ and $t\ne0$, where $a_{i,k}=\binom{i}{k}\binom{i-1}{k}{k!}$ for $0\le k\le i-1$, in~\cite[Theorem~2.1]{exp-reciprocal-cm-IJOPCM.tex},
\begin{equation}\label{h(t)-i-der-to0}
h^{(i)}(t)=\bigl(e^{1/t}\bigr)^{(i)}-\psi^{(i+1)}(t) =(-1)^ie^{1/t}\sum_{k=0}^{i-1}\frac{a_{i,k}}{t^{2i-k}}-\psi^{(i+1)}(t)\to0
\end{equation}
for $i\in\mathbb{N}$ as $t\to\infty$.
\par
Utilizing the recurrence formula
\begin{equation}\label{abram-6.4.7}
\psi^{(n)}(z+1)=\psi^{(n)}(z)+(-1)^n\frac{n!}{z^{n+1}}
\end{equation}
in~\cite[p.~260, 6.4.7]{abram} and calculating reveal
\begin{align*}
h(t+1)-h(t)&=e^{1/(t+1)}-e^{1/t}+\psi'(t)-\psi'(t+1)\\
&=e^{1/(t+1)}-e^{1/t}+\frac1{t^2}\\
&=\frac1{t^2}+\sum_{k=0}^\infty\frac1{(k+1)!}\biggl[\frac1{(t+1)^{k+1}}-\frac1{t^{k+1}}\biggr],\\
[h(t+1)-h(t)]^{(i)}&=(-1)^i\frac{(i+1)!}{t^{i+2}} +\sum_{k=0}^\infty\frac{(-1)^i(i+k)!}{(k+1)!k!} \biggl[\frac1{(t+1)^{i+k+1}}-\frac1{t^{i+k+1}}\biggr],
\end{align*}
and
\begin{align*}
(-1)^i[h(t+1)-h(t)]^{(i)}&=\frac{(i+1)!}{t^{i+2}} +\sum_{k=0}^\infty\frac{(i+k)!}{(k+1)!k!} \biggl[\frac1{(t+1)^{i+k+1}}-\frac1{t^{i+k+1}}\biggr]\\
&<\frac{(i+1)!}{t^{i+2}} +\sum_{k=0}^2\frac{(i+k)!}{(k+1)!k!} \biggl[\frac1{(t+1)^{i+k+1}}-\frac1{t^{i+k+1}}\biggr]\\
&=\frac{i!}{12t^{i+3}(t+1)^{i+3}}f_i(t),
\end{align*}
where
\begin{align*}
f_i(t)&=6(i+1)t(t+1)\bigl[(t+1)^{i+2}+t^{i+2}\bigr] -12 t^2(t+1)^2\bigl[(t+1)^{i+1}-t^{i+1}\bigr]\\
&\quad-(i+1)(i+2)\bigl[(t+1)^{i+3}-t^{i+3}\bigr]\\
&=6(i+1)t(t+1)\Biggl[\sum_{\ell=0}^{i+2}\binom{i+2}{\ell}t^{\ell}+t^{i+2}\Biggr] -12 t^2(t+1)^2\sum_{\ell=0}^{i}\binom{i+1}{\ell}t^{\ell}\\
&\quad-(i+1)(i+2)\sum_{\ell=0}^{i+2}\binom{i+3}{\ell}t^{\ell}\\
&=\frac{(i-1)(i+4)(i+5)}{2}\biggl[\frac{(2-i)(i+3)}{3}t-i\biggr]t^2-i(i+1)(i+5)t\\
&\quad-\sum_{\ell=4}^{i}\biggl[(i+1)(i+2)\binom{i+3}{\ell}-6(i+1)\binom{i+3}{\ell-1} +12\binom{i+3}{\ell-2}\biggr]t^\ell\\
&=\frac{(i-1)(i+4)(i+5)}{2}\biggl[\frac{(2-i)(i+3)}{3}t-i\biggr]t^2-i(i+1)(i+5)t\\
&\quad-(1+i)(2+i)-(i+4)(i+5)\sum_{\ell=4}^{i} \frac{(i-\ell+1)(i-\ell+2)}{\ell(i-\ell+5)}\binom{i+3}{\ell-1}t^\ell
\end{align*}
and an empty sum is understood to be nil. As a result, the function $f_i(t)$ is negative and
\begin{equation}\label{(-1)=i[h(t+1)-h(t)]=(i)}
(-1)^i[h(t+1)-h(t)]^{(i)}=(-1)^i[h(t+1)]^{(i)}-(-1)^i[h(t)]^{(i)}<0
\end{equation}
for all $i\ge0$ and $t\in(0,\infty)$. Hence, by consecutive recursion and~\eqref{h(t)-i-der-to0},
\begin{multline*}
(-1)^i[h(t)]^{(i)}\ge(-1)^i[h(t+1)]^{(i)}\ge(-1)^i[h(t+2)]^{(i)}\ge \dotsm\\*
\ge(-1)^i[h(t+k)]^{(i)}\ge(-1)^i\lim_{k\to\infty}[h(t+k)]^{(i)}=0
\end{multline*}
for $i\in\mathbb{N}$ and $t\in(0,\infty)$. This implies that the function $h(t)$ is decreasing on $(0,\infty)$. Combining this monotonicity with~\eqref{h(t)-limit=1} gives $h(t)>1$ on $(0,\infty)$. In conclusion, by definition, the function $h(t)$ is completely monotonic on $(0,\infty)$. The first proof of Theorem~\ref{CM-Exp-thm} is complete.
\end{proof}

\begin{proof}[Second proof]
To see the complete monotonicity of the function $h$, one writes
\begin{align*}
h(t)-h(t+1)&=e^{1/t}-e^{1/(t+1)}-\frac1{t^2} =\sum_{k=1}^3A_k(t)+\sum_{k=4}^\infty A_k(t)-\frac1{t^2}
\end{align*}
with
\begin{equation*}
A_k(t)=\frac1{k!}\frac{(t+1)^k-t^k}{t^k(t+1)^k} =\frac1{k!}\sum_{j=0}^{k-1}\frac1{(t+1)^kt^{k-j}}.
\end{equation*}
Now since
\begin{equation*}
\sum_{k=1}^3A_k(t)-\frac1{t^2}=\frac1{6t^3(t+1)^3},
\end{equation*}
the difference $h(t)-h(t+1)$ is a sum of completely monotonic functions and hence completely monotonic on $(0,\infty)$. Finally, the function
\begin{equation*}
h(t)=\lim_{n\to\infty}\sum_{k=0}^n[h(t+k)-h(t+k+1)]
\end{equation*}
is completely monotonic on $(0,\infty)$. The second proof of Theorem~\ref{CM-Exp-thm} is complete.
\end{proof}

\section{Proofs of Theorem~\ref{exp=k=degree=k+1-int-thm}}

It is general knowledge that the exponential function $e^{1/z}$ for $z\in\mathbb{C}$ with $z\ne0$ can be expanded into the Laurent series
\begin{equation}\label{exp-reciproc-series}
e^{1/z}=\sum_{m=0}^\infty\frac1{m!}\frac1{z^m},\quad z\ne0.
\end{equation}
Therefore, it is clear that
\begin{equation}\label{exp-residue-term}
H_k(z)=\sum_{m=k+1}^\infty\frac{1}{m!}\frac1{z^m},\quad z\ne0
\end{equation}
and $x^{k+1}H_k(x)$ is completely monotonic on $(0,\infty)$. That is,
\begin{equation}\label{H-k(t)-degree>k+1}
\cmdeg{t}[H_k(t)]\ge k+1.
\end{equation}
Since, for any $\varepsilon>0$, the function
\begin{equation*}
x^{k+1+\varepsilon}H_k(x)=x^\varepsilon\sum_{m=0}^\infty\frac{1}{(m+k+1)!}\frac1{x^m}
\end{equation*}
tends to $\infty$ as $x\to\infty$, we see that for any $\varepsilon>0$ the function $x^{k+1+\varepsilon}H_k(x)$ is not completely monotonic on $(0,\infty)$. That is,
\begin{equation}\label{H-k(t)-degree<k+1}
\cmdeg{t}[H_k(t)]\le k+1.
\end{equation}
Combining~\eqref{H-k(t)-degree>k+1} and~\eqref{H-k(t)-degree<k+1} leads to~\eqref{H-k(t)-degree}.
\par
For $\Re z>0$ and $\Re k>0$, it was listed in~\cite[p.~255, 6.1.1]{abram} that
\begin{equation*}%\label{Gamma(z)=k-z-int}
\Gamma(z)=k^z\int_0^\infty t^{z-1}e^{-kt}\td t.
\end{equation*}
This formula can be rearranged as
\begin{equation}\label{Gamma(z)=k-z-int-rearr}
\frac1{z^w}=\frac1{\Gamma(w)}\int_0^\infty t^{w-1}e^{-zt}\td t
\end{equation}
for $\Re z>0$ and $\Re w>0$. Substituting the formula~\eqref{Gamma(z)=k-z-int-rearr} into~\eqref{exp-residue-term} yields
\begin{align*}
H_k(z)&=\int_0^\infty \Biggl[\sum_{m=k+1}^\infty\frac{1}{m!}\frac1{\Gamma(m)}t^{m-1}\Biggr]e^{-zt}\td t \\
&=\frac1{k!(k+1)!}\int_0^\infty {}_1F_2(1;k+1,k+2;t)t^k e^{-zt}\td t.
\end{align*}
The integral representation~\eqref{exp=k=degree=k+1-int} follows.
\par
The function $H_k(z)$ can be rewritten as
\begin{align*}
H_k(z)&=\frac1{z^{k+1}}\sum_{m=0}^\infty\frac{1}{(k+1+m)!}\frac1{z^m}\\
&=\frac1{z^{k+1}}\Biggl\{\frac1{(k+1)!}+\int_0^\infty \Biggl[\sum_{m=1}^\infty\frac{1}{(k+1+m)!} \frac1{\Gamma(m)}t^{m-1}\Biggr] e^{-zt}\td t\Biggr\}\\
&=\frac1{z^{k+1}}\biggl[\frac1{(k+1)!}+\int_0^\infty \frac{I_{k+2} \bigl(2 \sqrt{t}\,\bigr)}{t^{(k+2)/2}} e^{-zt}\td t\biggr].
\end{align*}
The integral representation~\eqref{exp=k=degree=k+1-int-bes} follows.
Theorem~\ref{exp=k=degree=k+1-int-thm} is thus proved.

\section{Proof of Theorem~\ref{Bessel-2-ineq-thm}}

When $k=0$, the integral representations~\eqref{exp=k=degree=k+1-int} and~\eqref{exp=k=degree=k+1-int-bes} become
\begin{equation}\label{open-answer-1}
e^{1/z}=1+\int_0^\infty \frac{I_1\bigl(2\sqrt{t}\,\bigr)}{\sqrt{t}\,} e^{-zt}\td t
\end{equation}
and
\begin{equation}\label{open-answer-2}
e^{1/z}=1+\frac1{z}\biggl[1+\int_0^\infty \frac{I_{2} \bigl(2\sqrt{t}\,\bigr)}{t} e^{-zt}\td t\biggr]
\end{equation}
for $\Re z>0$. Hence, by~\eqref{polygamma} for $n=1$, the function $h(z)$ defined by~\eqref{alpha-exp=psi-eq} has the following integral representation
\begin{equation}\label{h(t)-int-rep-eq}
h(z)=1+\int_0^\infty \biggl[\frac{I_1\bigl(2\sqrt{u}\,\bigr)}{\sqrt{u}\,} -\frac{u}{1-e^{-u}}\biggr]e^{-zu}\td u
\end{equation}
and
\begin{equation}\label{h(t)-int-rep-eq-der}
(-1)^kh^{(k)}(t)=\int_0^\infty \biggl[\frac{I_1\bigl(2\sqrt{u}\,\bigr)}{\sqrt{u}\,} -\frac{u}{1-e^{-u}}\biggr]u^ke^{-tu}\td u
\end{equation}
for $k\ge1$ are completely monotonic on $(0,\infty)$.
The famous Hausdorff-Bernstein-Widder Theorem~\cite[p.~161, Theorem~12b]{widder} states that a necessary and sufficient condition that $f(x)$ should be completely monotonic for $0<x<\infty$ is that
\begin{equation} \label{berstein-1}
f(x)=\int_0^\infty e^{-xt}\td\alpha(t),
\end{equation}
where $\alpha(t)$ is non-decreasing and the integral converges for $0<x<\infty$. Consequently, the function in the bracket of~\eqref{h(t)-int-rep-eq-der} is not less than zero, that is,
\begin{equation}\label{lca-qi-wang-final-ineq}
\frac{I_1\bigl(2\sqrt{u}\,\bigr)}{\sqrt{u}\,} \ge\frac{u}{1-e^{-u}}
\end{equation}
in which replacing $2\sqrt{u}\,$ by $t$ yields the inequality~\eqref{I=1-exp-ineq}. The proof of Theorem~\ref{Bessel-2-ineq-thm} is complete.

\section{Remarks}

In this final section, we would like to remark something.

\begin{remark}
An alternative proof of Theorem~\ref{CM-Exp-thm} was given in~\cite[Theorem~3.1]{QiBerg.tex}.
\end{remark}

\begin{remark}
The result~\eqref{H-k(t)-degree} in Theorem~\ref{exp=k=degree=k+1-int-thm} can be generalized as follows.
Suppose that
\begin{equation*}
f(t)=\sum_{m=0}^\infty a_mt^m
\end{equation*}
has infinite radius of convergence and that $a_m\ge0$ for all $m$. Then
\begin{equation*}
\mathcal{H}_k(t)=f\biggl(\frac1t\biggr)-\sum_{m=0}^k\frac{a_m}{t^m} =\sum_{m=k+1}^\infty\frac{a_m}{t^m}
\end{equation*}
is completely monotonic as well as
\begin{equation*}
t^{k+1}\mathcal{H}_k(t)=\sum_{m=k+1}^\infty\frac{a_m}{t^{m-k-1}},
\end{equation*}
so the function $\mathcal{H}_k(t)$ has the completely monotonic degree $k+1$ if $a_{k+1}\ne0$ and, if $a_{k+1}=\dotsm=a_{k+p}=0$ but $a_{k+p+1}\ne0$, then $\mathcal{H}_k(t)$ has the completely monotonic degree $k+p+1$.
\end{remark}

\begin{remark}
The integral representations~\eqref{open-answer-1} and~\eqref{open-answer-2} supply affirmative answers to an open problem posed in~\cite[p.~127, Section~4]{exp-reciprocal-cm-IJOPCM.tex}. See also~\cite{Zhang-Xiao-Jing-Thesis.tex}.
\end{remark}

\begin{remark}
The inequality~\eqref{I=1-exp-ineq} in Theorem~\ref{Bessel-2-ineq-thm} or, equivalently, the inequality~\eqref{lca-qi-wang-final-ineq} has been further investigated in the subsequent papers~\cite{Bessel-ineq-Dgree-CM.tex, QiBerg.tex}.
\end{remark}

\begin{remark}
The formulas~\eqref{exp=k=degree=k+1-int} and~\eqref{open-answer-1} have been applied in~\cite{Lah-Number-Int-Property.tex, Bell-Stirling-HyperGeom-JIS.tex} to combinatorics and number theory.
\end{remark}

\begin{remark}
This paper is a corrected and extended version of the preprint~\cite{simp-exp-degree.tex}.
\end{remark}

\begin{thebibliography}{99}

\bibitem{abram}
M. Abramowitz and I. A. Stegun (Eds), \textit{Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables}, National Bureau of Standards, Applied Mathematics Series \textbf{55}, 9th printing, Washington, 1970.

\bibitem{psi-proper-fraction-degree-two.tex}
B.-N. Guo and F. Qi, \textit{A completely monotonic function involving the tri-gamma function and with degree one}, Appl. Math. Comput. \textbf{218} (2012), no.~19, 9890\nobreakdash--9897; Available online at \url{http://dx.doi.org/10.1016/j.amc.2012.03.075}.

\bibitem{Yang-Fan-2008-Dec-simp.tex}
B.-N. Guo and F. Qi, \textit{Refinements of lower bounds for polygamma functions}, Proc. Amer. Math. Soc. \textbf{141} (2013), no.~3, 1007\nobreakdash--1015; Available online at \url{http://dx.doi.org/10.1090/S0002-9939-2012-11387-5}.

\bibitem{Lah-Number-Int-Property.tex}
B.-N. Guo and F. Qi, \textit{Some integral representations and properties of Lah numbers}, available online at \url{http://arxiv.org/abs/1402.2367}.

\bibitem{Koumandos-Pedersen-09-JMAA}
S. Koumandos and H. L. Pedersen, \textit{Completely monotonic functions of positive order and asymptotic expansions of the logarithm of Barnes double gamma function and Euler's gamma function}, J. Math. Anal. Appl. \textbf{355} (2009), no.~1, 33\nobreakdash--40; Available online at \url{http://dx.doi.org/10.1016/j.jmaa.2009.01.042}.

\bibitem{Bell-Stirling-HyperGeom-JIS.tex}
F. Qi, \textit{An explicit formula for Bell numbers in terms of Stirling numbers and hypergeometric functions}, available online at \url{http://arxiv.org/abs/1402.2361}.

\bibitem{Bessel-ineq-Dgree-CM.tex}
F. Qi, \emph{Properties of modified Bessel functions and completely monotonic degrees of differences between exponential and trigamma functions}, available online at \url{http://arxiv.org/abs/1302.6731}.

\bibitem{QiBerg.tex}
F. Qi and C. Berg, \textit{Complete monotonicity of a difference between the exponential and trigamma functions and properties related to a modified Bessel function}, Mediterr. J. Math. \textbf{10} (2013), no.~4, 1685\nobreakdash--1696; Available online at \url{http://dx.doi.org/10.1007/s00009-013-0272-2}.

\bibitem{simp-exp-degree.tex}
F. Qi and S.-H. Wang, \textit{Complete monotonicity of a difference between the exponential and trigamma functions and completely monotonic degree of the exponential function}, available online at \url{http://arxiv.org/abs/1210.2012}.

\bibitem{widder}
D. V. Widder, \textit{The Laplace Transform}, Princeton University Press, Princeton, 1946.

\bibitem{Zhang-Xiao-Jing-Thesis.tex}
X.-J. Zhang, \textit{Integral Representations, Properties, and Applications of Three Classes of Functions}, Thesis supervised by Professor Feng Qi and submitted for the Degree of Master of Science at Tianjin Polytechnic University in January 2013. (Chinese)

\bibitem{exp-reciprocal-cm-IJOPCM.tex}
X.-J. Zhang, F. Qi, and W.-H. Li, \textit{Properties of three functions relating to the exponential function and the existence of partitions of unity}, Int. J. Open Probl. Comput. Sci. Math. \textbf{5} (2012), no.~3, 122\nobreakdash--127.

\end{thebibliography}
\end{document}
