\section{Introduction and preliminaries}

Throughout this paper $\mathbb{N}$ denotes the set of the positive integers, and $\mathbb{Z}, \mathbb{Q}$, and
$\mathbb{R}$ have the usual meaning.

The aim of this work is to prove characterization theorems on derivations as well as on linear functions.
Therefore, firstly we have to recall some definitions and auxiliary results.

A function $f:\mathbb{R}\rightarrow\mathbb{R}$ is called an \emph{additive} function
if,
\[
f(x+y)=f(x)+f(y)
\]
holds for all $x, y\in\mathbb{R}$.

We say that an additive
function $f:\mathbb{R}\rightarrow\mathbb{R}$ is a \emph{derivation} if
\[
f(xy)=xf(y)+yf(x)
\]
is fulfilled for all $x, y\in\mathbb{R}$.

The additive function $f\colon\mathbb{R}\to\mathbb{R}$ is termed to be a \emph{linear function} if
$f$ is of the form
\[
 f(x)=f(1)x
\qquad
\left(x\in\mathbb{R}\right).
\]

It is easy to see from the above definition that every derivation
$f:\mathbb{R}\rightarrow\mathbb{R}$ satisfies equation
\[
\tag{$\ast$} f(x^{k})=kx^{k-1}f(x)
\quad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
for arbitrarily fixed $k\in\mathbb{Z}\setminus\left\{0\right\}$.
Furthermore, the converse is also true, in the following sense:
if $k\in\mathbb{Z}\setminus\left\{0, 1\right\}$ is fixed and
an additive function $f:\mathbb{R}\rightarrow\mathbb{R}$ satisfies
$\left(\ast\right)$, then $f$ is a derivation, see e.g.,
Jurkat \cite{Jur65}, Kurepa \cite{Kur64}, and
Kannappan--Kurepa \cite{KK70}.

Concerning linear functions, Jurkat \cite{Jur65} and,
independently, Kurepa \cite{Kur64} proved that every additive function
$f:\mathbb{R}\rightarrow\mathbb{R}$ satisfying
\[
f\left(\frac{1}{x}\right)=\frac{1}{x^{2}}f(x)
\quad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
has to be linear.

In \cite{NH68} A.~Nishiyama and S.~Horinouchi investigated
additive functions $f:\mathbb{R}\rightarrow\mathbb{R}$
satisfying the additional equation
\[
f(x^{n})=cx^{k}f(x^{m})
\quad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right),
\]
where $c\in\mathbb{R}$ and $n, m, k\in\mathbb{Z}$ are arbitrarily
fixed.

Furthermore, the above problem was generalized by Pl.~Kannappan and S.~Kurepa
in \cite{KK70} by proving the following statement.

\begin{thm}\label{T1}
Let $f, g:\mathbb{R}\rightarrow\mathbb{R}$ be additive functions and
$n, m\in\mathbb{Z}\setminus\left\{0\right\}$, $n\neq m$.
Suppose that
\[
f(x^{n})=x^{n-m}g(x^{m})
\]
holds for all $x\in\mathbb{R}\setminus\left\{0\right\}$.
Then the functions $F, G:\mathbb{R}\rightarrow\mathbb{R}$ defined by
\[
F(x)=f(x)-f(1)x
\quad
\text{and}
\quad
G(x)=g(x)-g(1)x
\quad
\left(x\in\mathbb{R}\right)
\]
are derivations and $nF(x)=mG(x)$ is fulfilled for all $x\in\mathbb{R}$.
\end{thm}

Here we remark that the 'approximate' variant of the above theorem was dealt with in \cite{BG10}.

Furthermore, in a series of papers  (see \cite{HKR99, HKR00, HK00, HK00a, HKR01}) F.~Halter--Koch and
L.~Reich proved several characterization theorems concerning derivations as well as field
homomorphisms.

Let $\mathbb{K}$ be a field containing $\mathbb{Q}$, $n\in\mathbb{Z}\setminus\left\{0\right\}$,
$\left(\begin{array}{cc}
a&b\\
c&d
\end{array}
\right)\in\mathbf{GL}_{2}(\mathbb{Q})$ and let
$f, g\colon\mathbb{K}\to\mathbb{K}$ be additive function so that
\[
\tag{$\star$} f\left(\frac{ax^{n}+b}{cx^{n}+d}\right)=\frac{x^{n-1}g(x)}{(cx^{n}+d)^{2}},
\]
respectively,
\[
\tag{$\blacktriangle$} f\left(\frac{ax^{n}+b}{cx^{n}+d}\right)=\frac{ag(x)^{n}+b}{cg(x)^{n}+d}
\]
holds for all possible values of $x$. In \cite{HK00} it is proved that
equation $\left(\star\right)$ (under a mild condition) implies for the function $g$ that the function
$G\colon\mathbb{K}\to\mathbb{K}$ defined by
\[
G(x)=g(x)-g(1)x
\qquad
\left(x\in\mathbb{K}\right)
\]
is a derivation.
Furthermore, in \cite{HKR01} the authors succeed to prove that equation $\left(\blacktriangle\right)$
furnishes that the mapping $g(1)^{-1}\cdot g\colon\mathbb{K}\to\mathbb{K}$ is a field automorphism.

In this work we will extend the above mentioned results, and (among others) the following problem will be dealt with.
Let $n\in\mathbb{Z}$, $f, g\colon\mathbb{R}\to\mathbb{R}$ be additive functions,
$\left(\begin{array}{cc}
a&b\\
c&d
\end{array}
\right)\in\mathbf{GL}_{2}(\mathbb{Q})$ be arbitrarily fixed,
and let us assume that the mapping
\[
 \phi(x)=g\left(\frac{ax^{n}+b}{cx^{n}+d}\right)-\frac{x^{n-1}f(x)}{(cx^{n}+d)^{2}}
\quad
\left(x\in\mathbb{R}, cx^{n}+d\neq 0\right)
\]
satisfies some regularity on its domain (e.g. (locally) boundedness, continuity, measurability).
Is it true that in this case the above functions can be represented as a sum of a derivation and a linear function?

During the second section the concept of multi-additive functions will be used. Therefore, at this part we will list some definitions and
statements that we will use subsequently.

Let $G, H$ be abelian groups, let $h\in G$ be arbitrary and consider a function $f:G\rightarrow H$.
The \emph{difference operator} $\Delta_{h}$ with the span $h$ of the function $f$ is defined by
\[
 \Delta_{h}f(x)=f(x+h)-f(x)
\qquad
\left(x\in G\right).
\]
The iterates $\Delta^{n}_{h}$ of $\Delta_{h}$, $n=0, 1, \ldots$ are defined by the recurrence
\[
 \Delta^{0}_{h}f=f, \qquad
\Delta^{n+1}_{h}f=\Delta_{h}\left(\Delta^{n}_{h}f\right)
\qquad
\left(n=0, 1, \ldots\right)
\]
Furthermore, the superposition of several difference operators will be denoted shortly
\[
 \Delta_{h_{1} \ldots h_{n}}f=\Delta_{h_{1}}\ldots \Delta_{h_{n}}f,
\]
where $n\in\mathbb{N}$ and $h_{1}, \ldots, h_{n}\in G$.

Let $n\in\mathbb{N}$ and $G, H$ be abelian groups.
A function $F:G^{n}\rightarrow H$ is
called \emph{$n$--additive} if, for every
$ i \in \{\, 1 \,,\, 2 \,,\, \dots \,,\, n \,\} $
and for every
$ x_{1}, \ldots, x_{n}, y_{i}\in G \,$,
\begin{multline*}
F\left(x_{1}, \ldots, x_{i-1}, x_{i} + y_{i}, x_{i+1}, \ldots, x_{n}\right)
\\
=F\left(x_{1}, \ldots, x_{i-1}, x_{i}, x_{i+1}, \ldots, x_{n}\right)
+F\left(x_{1}, \ldots, x_{i-1}, y_{i}, x_{i+1}, \ldots, x_{n}\right),
\end{multline*}
i.e., $F$ is additive in each of its variables
$x_{i}\in G$, $i=1, \ldots, n$.
For the sake of brevity we use the notation $G^{0}=G$ and we call
constant functions from $G$ to $H$ $0$--additive functions.
Let $F:G^{n}\rightarrow H$ be an arbitrary function.
By the \emph{diagonalization (or trace)} of $F$ we understand the
function
$f:G\rightarrow H$ arising from
$F$ by putting all the variables (from $G$)
equal:
\[
f(x)=F(x, \ldots, x)
\qquad
\left(x\in G\right).
\]
It can be proved by induction that for any symmetric, $n$--additive function
$F:G^{n}\rightarrow H$ the equality
\[
 \Delta_{y_{1}, \ldots, y_{k}}f(x)=
\left\{
\begin{array}{lcl}
 n!F(y_{1}, \ldots, y_{n})& \text{for}& k=n \\
0 & \text{for} & k>n
\end{array}
\right.
\]
holds, whenever $x, y_{1}, \ldots, y_{n}\in G$, where $f:G\rightarrow H$ denotes the trace of the
function $F$. This means that a symmetric, $n$--additive function is uniquely determined by its trace.

The function $f:G\rightarrow H$ is called a \emph{polynomial function} of degree at most $n$, where
$n$ is a nonnegative integer, if
\[
 \Delta_{y_{1}, \ldots, y_{n+1}}f(x)=0
\]
is satisfied for all $x, y_{1}, \ldots, y_{n+1}\in G$.

\begin{thm}[Kuczma \cite{Kuc09}, Sz\'{e}kelyhidi \cite{Sze91}]\label{T1.1}
The function $p:G\rightarrow H$ is a polynomial at degree at most $n$ if and only if there
exist symmetric, $k$--additive functions $F_{k}:G^{k}\rightarrow H$, $k=0, 1, \ldots, n$ such that
\[
 p(x)=\sum_{k=0}^{n}f_{k}(x)
\qquad
\left(x\in G\right),
\]
where $f_{k}$ denotes the trace of the function $F_{k}$, $k=0, 1, \ldots, n$.
Furthermore, this expression for the function $p$ is unique in the sense that the
functions $F_{k}$, which are not identically zero,  are uniquely determined.
\end{thm}

During the proof of our main theorem, the basic idea is to apply two statements from
\cite{Sze85, Sze91} concerning polynomial functions.
Despite the fact that the following results will be used only in case the topological group is
$\mathbb{R}$ and the topological linear space is also $\mathbb{R}$, here we present the general case.
The reason for it is, that the problem, we investigate, can be formulated in more general circumstances.
Perhaps in the general case the two cited theorems of L.~Sz\'{e}kelyhidi would also play an important role.

\begin{thm}[Sz\'{e}kelyhidi \cite{Sze85, Sze91}]\label{T1.2}
 Let $G$ be an abelian group and let $X$ be a locally convex topological
linear space. If a polynomial $p:G\rightarrow X$ is bounded on $G$, then it is
constant.
\end{thm}

\begin{thm}[Sz\'{e}kelyhidi \cite{Sze85, Sze91}]\label{T1.3}
 Let $G$ be a topological abelian group which is generated by any neighbourhood of the
zero, and let $X$ be a topological linear space, and $p:G\rightarrow X$ be a polynomial function.
Then the following statements hold.
\begin{enumerate}[(i)]
\item If $p:G\rightarrow X$ is continuous at a point, then it is
continuous on $G$.
\item Assume that $G$ is locally compact and $X$ is locally convex.
If $p:G\rightarrow X$ is bounded on a measurable set of positive measure,
then it is continuous.
\item Suppose that $G$ is locally compact and $X$ is locally convex and locally bounded.
If $p:G\rightarrow X$ is measurable on a measurable set of positive measure,
then it is continuous.
\end{enumerate}
\end{thm}

\section{Main results}

\subsection*{Preparatory statements}

In order to avoid superfluous repetitions, henceforth we will say that the
function in question is \emph{locally regular} on its domain, if at least
one of the following statements are fulfilled.
\begin{enumerate}[(i)]
 \item bounded on a measurable set of positive measure;
\item continuous at a point;
\item there exists a set of positive Lebesgue measure so that the restriction of the function in question
is measurable in the sense of Lebesgue.
\end{enumerate}

Furthermore, a function will be called \emph{globally regular}, if instead of (ii), \\
(ii)' \; continuous on its domain \\
holds.

Firstly we prove a simple lemma.

\begin{lem}\label{L1}
 Let $\alpha\in\mathbb{R}$ be arbitrarily fixed and let us assume that for the
function $\phi\colon ]0, +\infty[\to\mathbb{R}$ the following statements are valid.
\begin{enumerate}[(a)]
 \item the function $\phi$ is $\mathbb{Q}$-homogeneous of order $\alpha$, that is,
\[
 \phi(rx)=r^{\alpha}\phi(x)
\qquad
\left(x\in ]0, +\infty[, r\in\mathbb{Q}\cap ]0, +\infty[\right).
\]
\item the function $\phi$ is continuous at a point.
\end{enumerate}
Then $\phi$ is continuous everywhere.
\end{lem}
\begin{proof}
 Let us assume that the function $\phi$ is continuous at the point $x_{0}\in ]0, +\infty[$ and let
$\tilde{x}\in ]0, +\infty[$ be arbitrary.
Then, there exists a sequence of positive rational numbers $(r_{n})_{n\in\mathbb{N}}$ so that
$\lim_{n\to\infty} r_{n}= \dfrac{\tilde{x}}{x_{0}}$. In this case
the sequence $\left(\frac{1}{r_{n}}\right)_{n\in\mathbb{N}}$ is also a sequence of positive rational numbers and it
converges to $\dfrac{x_{0}}{\tilde{x}}$.
Due to property (a),
\[
\frac{1}{r_{n}^{\alpha}}\phi(\tilde{x})=\phi\left(\frac{\tilde{x}}{r_{n}}\right)
\qquad
\left(n\in\mathbb{N}\right).
\]
Taking the limit $n\to \infty$, the left hand side converges to
$\left(\frac{x_{0}}{\tilde{x}}\right)^{\alpha}\phi(\tilde{x})$. Furthermore,
the sequence $\left(\frac{\tilde{x}}{r_{n}}\right)_{n\in\mathbb{N}}$ converges to
$x_{0}$, therefore the continuity of the function $\phi$ implies that
the right hand side tends to $\phi(x_{0})$ as $n\to \infty$.
This implies that
\[
 \left(\frac{x_{0}}{\tilde{x}}\right)^{\alpha}\phi(\tilde{x})=\phi\left(\frac{x_{0}}{\tilde{x}}\tilde{x}\right)
\]
is fulfilled.
Since $\tilde{x}\in ]0, +\infty[$ was arbitrary, we get that
\[
 \phi(\lambda x)=\lambda^{\alpha}\phi(x)
\qquad
\left(\lambda, x\in ]0, +\infty[\right),
\]
which obviously implies the (everywhere) continuity of the function $\phi$.

\end{proof}

Furthermore, it is important pointing out the following fact.
Fix $\alpha\in\mathbb{R}$ and let $\phi\colon ]0, +\infty[\to\mathbb{R}$
be a $\mathbb{Q}$-homogeneous function of order $\alpha$.
Suppose that $\phi$ fulfils property (i) or (iii) on the set of positive Lebesgue measure
$D$.
Then, the $\mathbb{Q}$-homogeneity of the function $\phi$ implies that
(i), respectively (iii) holds for the function $\phi$ on the set
$r\cdot D$ for all $r\in\mathbb{Q}$.

\begin{thm}\label{Thm2.1}
 Let $n, m\in\mathbb{Z}\setminus\left\{0\right\}$, $n\neq m$ so that
$n=-m$ or $\mathrm{sign}(n)=\mathrm{sign}(m)$, let further
$f, g\colon\mathbb{R}\to\mathbb{R}$ be additive functions.
Define the function $\phi\colon \mathbb{R}\setminus\left\{0\right\}\to\mathbb{R}$ by the
formula
\[
 \phi(x)=f\left(x^{n}\right)-x^{n-m}g\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right),
\]
and assume that $\phi$ is locally regular.
Then, the functions $F, G\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x \quad
\text{and}
\quad
G(x)=g(x)-g(1)x
\qquad
\left(x\in\mathbb{R}\right)
\]
are derivations and
\[
 nF(x)=mG(x)
\]
holds for arbitrary $x\in\mathbb{R}$.
\end{thm}
\begin{proof}
 Concerning the values of $n$ and $m$ we have to distinguish three cases.
At first, let us assume that $n, m>0$.
There is no loss of generality in assuming $n>m$.
Define the function $\Phi$ on $\mathbb{R}^{n}$ by
\begin{multline*}
 \Phi(x_{1}, \ldots, x_{n})=
f(x_{1}\cdots x_{n}) -
\frac{1}{\binom{n}{m}}
\sum_{\mathrm{card}(I)=m}
\left( \prod_{j \in \{1,2,\dots,n\} \setminus I} x_j \right)
g \left( \prod_{i \in I} x_i \right)
\\
\left(x_{1}, \ldots, x_{n}\in\mathbb{R}\right),
\end{multline*}
where the summation is considered for all subsets $I$ of cardinality $m$ of the index set
$\left\{1, 2,\ldots, n\right\}$.
Due to the additivity of the functions $f$ and $g$, the function
$\Phi$ is a symmetric and $n$--additive function. Furthermore, its trace, that is,
\[
 \Phi(x, \ldots, x)=\phi(x)=f\left(x^{n}\right)-x^{n-m}g\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\right)
\]
is a polynomial function. On the other hand $\phi$ is a locally regular function.
In view of Theorem \ref{T1.3}, this means that $\phi$ is a continuous polynomial function.
Therefore, there exists $c\in\mathbb{R}$ such that
\[
 \Phi(x_{1}, \ldots, x_{n})=cx_{1}\cdots x_{n}
\quad
\left(x_{1}, \ldots, x_{n}\in\mathbb{R}\right),
\]
that is,
\[
 \phi(x)=cx^{n}
\quad
\left(x\in\mathbb{R}\right).
\]
With the substitution $x=1$, we get $\phi(1)=c$.
On the other hand, the definition of the function $\phi$ yields that
$\phi(1)=f(1)-g(1)$.
Thus,
\[
 f(x^{n})-x^{n-m}g(x^{m})=\left[f(1)-g(1)\right]x^{n}
\qquad
\left(x\in\mathbb{R}\right).
\]
Define the functions $F, G\colon \mathbb{R}\to\mathbb{R}$ by
\[
 F(x)=f(x)-f(1)x
\quad
\text{and}
\quad
G(x)=g(x)-g(1)x
\qquad
\left(x\in\mathbb{R}\right).
\]
Then the above identity yields that
\[
 F\left(x^{n}\right)=x^{n-m}G(x^{m})
\qquad
\left(x\in\mathbb{R}\right).
\]
The statement of the theorem follows now from Theorem \ref{T1}.

Secondly, let us assume that $n, m<0$. In this case we get that the function
\[
 \phi(x)=f\left(x^{n}\right)-x^{n-m}g\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is locally regular on its domain. Let $u\in\mathbb{R}\setminus\left\{0\right\}$, with the substitution $x=\frac{1}{u}$ this
yields that
\[
 \phi\left(\frac{1}{u}\right)=
f\left(u^{-n}\right)-u^{-(n-m)}g\left(u^{-m}\right)
\qquad
\left(u\in\mathbb{R}\setminus\left\{0\right\}\right).
\]
Since $-n, -m>0$, the results of the previous case can be applied for the function
\[
 \psi(u)=\phi\left(\frac{1}{u}\right)
\qquad
\left(u\in \mathbb{R}\setminus \left\{0\right\}\right),
\]
which is, due to the local regularity of $\phi$, also locally regular.

Finally, let us assume that $n=-m$. Without the loss of generality $m>0$ can be assumed.
In this case
\[
 \phi(x)=f\left(x^{-m}\right)-x^{-2m}g\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is locally regular, or equivalently, the mapping
\begin{equation}\label{Eq1}
\psi(x)= \phi\left(\sqrt[m]{x}\right)=f\left(\frac{1}{x}\right)-\frac{1}{x^{2}}g(x)
\qquad
\left(x>0\right)
\end{equation}
has the local regularity property.
Let $u>0$ be arbitrary and  let us substitute $u(u+1)$ in place of $x$. Then
\begin{multline*}
\psi(u(u+1))= \phi\left(\sqrt[m]{u(u+1)}\right)
\\
=f\left(\frac{1}{u(u+1)}\right)-\frac{1}{u^{2}(u+1)^{2}}g\left(u(u+1)\right)
\\
\left(u>0\right).
\end{multline*}
Using the additivity of the function $f$,
\begin{multline*}
\psi(u(u+1))=\phi\left(\sqrt[m]{u(u+1)}\right)
\\
=f\left(\frac{1}{u}\right)-f\left(\frac{1}{u+1}\right) -\frac{1}{u^{2}(u+1)^{2}}g\left(u(u+1)\right)
\\
\left(u>0\right).
\end{multline*}
On the other hand,
\[
\psi(u)=\phi\left(\sqrt[m]{u}\right)=f\left(\frac{1}{u}\right)-\frac{1}{u^{2}}g(u)
\qquad
\left(u>0\right)
\]
and
\[
\psi(u+1)=\phi\left(\sqrt[m]{u+1}\right)=f\left(\frac{1}{u+1}\right)-\frac{1}{(u+1)^{2}}g(u+1)
\qquad
\left(u>0\right).
\]
Therefore,
\begin{multline*}
\psi(u(u+1))-\psi(u)+\psi(u+1)
\\
=
\phi\left(\sqrt[m]{u(u+1)}\right)-\phi(\sqrt[m]{u})+\phi(\sqrt[m]{u+1})
\\
=
f\left(\frac{1}{u}\right)-f\left(\frac{1}{u+1}\right) -\frac{1}{u^{2}(u+1)^{2}}g\left(u(u+1)\right)
\\
-f\left(\frac{1}{u}\right)+\frac{1}{u^{2}}g(u)+
f\left(\frac{1}{u+1}\right)-\frac{1}{(u+1)^{2}}g(u+1)
\qquad
\left(u>0\right)
\end{multline*}
Making use of the additivity of the function g, after rearrangement, we obtain that
\[
 \chi(u)=2ug(u)-g\left(u^{2}\right)
\qquad
\left(u>0\right),
\]
where
\[
 \chi(u)=u^{2}(u+1)^{2}\left[\psi\left(u(u+1)\right)-\psi(u)+\psi(u+1)\right]+u^{2}g(1)
\qquad
(u>0).
\]
By our assumptions, the function $\phi$  is locally regular on
$\mathbb{R}\setminus\left\{0\right\}$ and due to the additivity of $f$ and $g$,
it is $\mathbb{Q}$-homogeneous of order $n$. Thus, by Lemma \ref{L1},
$\phi$ is globally regular on $\mathbb{R}\setminus \left\{0\right\}$.
This implies that $\psi$ is globally regular on $]0, +\infty[$, which means that
$\chi$ is locally regular.
Due to the results of the first case this yields that the function
$G\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 G(x)=g(x)-g(1)x
\qquad
\left(x\in\mathbb{R}\right)
\]
is a derivation.
In view of \eqref{Eq1}, this implies that
\[
\phi\left(\sqrt[m]{x}\right)=f\left(\frac{1}{x}\right)-\frac{1}{x^{2}}\left[G(x)+g(1)x\right]
\qquad
\left(x>0\right),
\]
that is ,
\[
\phi\left(\sqrt[m]{x}\right)=f\left(\frac{1}{x}\right)+G\left(\frac{1}{x}\right)+g(1)\frac{1}{x}
\qquad
\left(x>0\right),
\]
since $G$ is a derivation.
Let $u>0$, with the substitution $x=\frac{1}{u}$ we get that
\[
\psi(u)=\phi\left(\sqrt[m]{\frac{1}{u}}\right)=f(u)+G(u)+g(1)u
\qquad
\left(u>0\right).
\]
Let us observe that the right hand side of this identity is an additive function, being the sum of additive functions.
Moreover, the left hand side is locally regular, due to the local regularity of $ \phi$.
Thus $\psi$ is a regular additive function, which means that there exists $c\in\mathbb{R}$ so that
\[
 f(u)+G(u)+g(1)u=cu
\qquad
\left(u\in\mathbb{R}\right).
\]
With $u=1$, $c=f(1)+g(1)$ can be obtained, therefore,
\[
 f(u)=\left[f(1)+g(1)\right]u-g(1)u+G(u)
\qquad
\left(u\in\mathbb{R}\right),
\]
i.e.,
\[
 f(u)=-G(u)+f(1)u
\qquad
\left(u\in\mathbb{R}\right).
\]
This means that the function $F\colon \mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x
\qquad
\left(x\in\mathbb{R}\right)
\]
is a derivation and
\[
 F(x)=-G(x)
\qquad
\left(x\in\mathbb{R}\right)
\]
holds.
\end{proof}

\begin{lem}\label{L2.1}
 Let $\kappa\in\mathbb{R}$, $n, m\in\mathbb{Z}$, $n \neq m$ and assume that
$f\colon\mathbb{R}\to\mathbb{R}$ is an additive function.
Define the function $\phi\colon\mathbb{R}\setminus\left\{0\right\}\to\mathbb{R}$ by
\[
 \phi(x)=f\left(x^{n}\right)-\kappa x^{n-m}f\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\setminus \left\{0\right\}\right)
\]
and assume that $\phi$ is locally regular.
Then, the function $F\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x \quad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is a derivation so that for any $x\in\mathbb{R}$
\[
\left(n-\kappa m\right)F(x)=0.
\]
\end{lem}
\begin{proof}
 In view of the previous theorem, it is enough to deal with the case $\mathrm{sign}(n)\neq\mathrm{sign}(m)$ and $n\neq -m$.
Due to the definition of the function $\phi$
\[
 \phi\left(x^{n}\right)=f\left(x^{n^{2}}\right)-\kappa x^{n(n-m)}f\left(x^{nm}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
and
\[
 \phi\left(x^{m}\right)=f\left(x^{nm}\right)-\kappa x^{m(n-m)}f\left(x^{m^{2}}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right),
\]
therefore
\[
 \phi\left(x^{n}\right)+\kappa x^{n(n-m)}\phi\left(x^{m}\right)=
f\left(x^{n^{2}}\right)-\kappa^{2} x^{n^{2}-m^{2}}f\left(x^{m^{2}}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right).
\]
By our assumptions, $\phi$ is a locally regular function on
$\mathbb{R}\setminus\left\{0\right\}$.
However, the additivity of $f$ implies that
\[
 \phi(rx)=r^{n}\phi(x)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}, r\in\mathbb{Q}\setminus \left\{0\right\}\right).
\]
Using Lemma \ref{L1}, we get that $\phi$ is globally regular.
Therefore, the function
\[
 \psi(x)=\phi\left(x^{n}\right)+\kappa x^{n(n-m)}\phi\left(x^{m}\right)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is locally regular.
Since $n^{2}, m^{2}>0$ and $n^{2}\neq m^{2}$, the results of the previous theorem can be applied
(with the choice $\psi(x)=\phi\left(x^{n}\right)+\kappa x^{n(n-m)}\phi\left(x^{m}\right)$ and $g(x)=\kappa^{2}f(x)$)
to obtain that
\[
 f(x)=F(x)+f(1)x
\qquad
\left(x\in\mathbb{R}\right),
\]
where $F\colon \mathbb{R}\to\mathbb{R}$ is a derivation and
\[
 nF(x)=m\kappa F(x)
\]
is also fulfilled for all $x\in\mathbb{R}$.

\end{proof}

From this lemma, the following corollary can be concluded immediately.

\begin{cor}\label{C2.1}
 Let $r\in\mathbb{Q}\setminus\left\{0, 1\right\}$ be arbitrarily fixed and
$f\colon\mathbb{R}\to\mathbb{R}$ be an additive function and define the function
 by
\[
 \phi(x)=f\left(x^{r}\right)-rx^{r-1}f\left(x\right)
\qquad
\left(x\in\mathbb{R},\, x>0\right),
\]
and assume that $\phi$ is locally regular.
Then, the function $F\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x \quad
\left(x\in\mathbb{R}\right)
\]
is a derivation.
\end{cor}

\subsection*{Derivations along rational functions}

In view of the results of the previous subsection, we are able to prove the following

\begin{thm}\label{T2.2}
 Let $n\in\mathbb{Z}\setminus\left\{0\right\}$ and
$\left(\begin{array}{cc}
a&b\\
c&d
\end{array}
\right)\in\mathbf{GL}_{2}(\mathbb{Q})$ be such that
\begin{enumerate}[--]
 \item if $c=0$, then $n\neq 1$;
\item if $d=0$, then $n\neq -1$.
\end{enumerate}
Let further $f, g\colon\mathbb{R}\to\mathbb{R}$ be additive functions and define the function
$\phi$ by
\[
 \phi(x)=f\left(\frac{ax^{n}+b}{cx^{n}+d}\right)-\frac{x^{n-1}g(x)}{\left(cx^{n}+d\right)^{2}}
\qquad
\left(x\in\mathbb{R}, \,
cx^{n}+d\neq 0\right).
\]
Let us assume $\phi$ to be globally regular.
Then, the functions $F, G\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x \quad \text{and}
\quad
G(x)=g(x)-g(1)x
\quad
\left(x\in\mathbb{R}\right)
\]
are derivations.
\end{thm}

\begin{proof}
 Firstly, let us suppose that  $c=0$. This means that the function
\[
 \phi(x)=f\left(\frac{a}{d}x^{n}+\frac{b}{d}\right)-\frac{1}{d^{2}}x^{n-1}g(x)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is globally regular.
In this case, the statement immediately follows from Theorem \ref{Thm2.1}.

Similarly, if $d=0$, then
\[
 \phi(x)=f\left(\frac{a}{c}+\frac{b}{c}x^{-n}\right)-x^{-n-1}g(x)
\qquad
\left(x\in\mathbb{R}\setminus\left\{0\right\}\right)
\]
is globally regular. Therefore, due to Theorem \ref{Thm2.1}, we obtain that the functions
\[
 F(x)=f(x)-f(1)x \quad
\text{and}
\quad
G(x)=g(x)-g(1)x
\qquad
\left(x\in\mathbb{R}\right)
\]
are derivations.

Thus, henceforth $cd\neq 0$ can be assumed. Furthermore, due to the $\mathbb{Q}$-homogeneity of the functions
$f$ and $g$, $c=1$ can be supposed. That is,
\[
 \phi(x)=f\left(\frac{ax^{n}+b}{x^{n}+d}\right)-\frac{x^{n-1}g(x)}{\left(x^{n}+d\right)^{2}}
\qquad
\left(x\in\mathbb{R}, \,
x^{n}+d\neq 0\right).
\]
Since the function $f$ is additive,
\[
 f\left(\frac{ax^{n}+b}{x^{n}+d}\right)=f(a)-f\left(\frac{D}{x^{n}+d}\right),
\qquad
\left(x\in\mathbb{R}, \,
x^{n}+d\neq 0\right),
\]
therefore,
\begin{equation}\label{Eq2}
 \phi(x)=f(a)-f\left(\frac{D}{x^{n}+d}\right)-\frac{x^{n-1}g(x)}{(x^{n}+d)^{2}}
\qquad
\left(x\in\mathbb{R}, \,
x^{n}+d\neq 0\right),
\end{equation}
where $D=\det\left(\begin{array}{cc}
                    a&b\\
1&d
                   \end{array}
\right)$.
Let us observe that
\[
 \frac{D}{x^{n}+d}=\frac{D}{d}-\frac{D}{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d}
\]
holds for all $x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0.$
Using this identity, we receive
\begin{multline}\label{Eq3}
\phi(x)=
f(a)-f\left(\frac{D}{d}\right)+f\left(\frac{D}{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d}\right)
-\frac{x^{n-1}g(x)}{\left(x^{n}+d\right)^{2}}
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right),
\end{multline}
where the additivity of the function $f$ was also used.
Let us replace $x$ by $\sqrt[n]{d^{2}}\dfrac{1}{x}$ in \eqref{Eq2} to acquire
\begin{multline*}
 \phi\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)=
f(a)-f\left(\frac{D}{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d}\right)-
\frac{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n-1}g\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)}
{\left(\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d\right)^{2}}
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right).
\end{multline*}
Since
\[
 \left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d= d\frac{1}{x^{n}}\left(x^{n}+d\right),
\]
the above identity yields that
\begin{multline*}
\phi\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)=
f(a)-f\left(\frac{D}{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d}\right)-
\frac{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n-1}g\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)}
{\left(d\frac{1}{x^{n}}\right)^{2}\left(\left(x^{n}+d\right)\right)^{2}}
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right).
\end{multline*}
After some rearrangement, we arrive at
\begin{multline}\label{Eq4}
 \phi\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)=
f(a)-f\left(\frac{D}{\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)^{n}+d}\right)-
\frac{x^{n-1}\frac{1}{\sqrt[n]{d^{2}}}x^{2}g\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)}{\left(x^{n}+d\right)^{2}}
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right).
\end{multline}
In case we add \eqref{Eq3} and \eqref{Eq4} together,
\begin{multline*}
 \phi(x)+\phi\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)=
f\left(2a-\frac{D}{d}\right)-
\frac{x^{n-1}}{\left(x^{n}+d\right)^{2}}
\left[
g(x)+x^{2}\frac{1}{\sqrt[n]{d^{2}}}g\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)
\right]
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right).
\end{multline*}
Let us define the functions
\[
 h(x)=\frac{1}{\sqrt[n]{d^{2}}}g\left(\sqrt[n]{d^{2}}x\right)
\qquad
\left(x\in\mathbb{R}\right)
\]
and
\begin{multline*}
 \psi(x)=
-\frac{1}{x^{2}}\frac{\left(x^{n}+d\right)^{2}}{x^{n-1}}
\left[
\phi(x)+\phi\left(\sqrt[n]{d^{2}}\frac{1}{x}\right)-f\left(2a-\frac{D}{d}\right)
\right]
\\
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right).
\end{multline*}
In this case
\[
 \psi(x)=h\left(\frac{1}{x}\right)+\frac{1}{x^{2}}g(x)
\quad
\left(x\in\mathbb{R}, \, x\neq 0, x^{n}+d\neq 0\right)
\]
holds.
By our assumptions $\phi$ is a globally regular mapping, therefore the function
$\psi$ has the local regularity property.
Due to Theorem \ref{Thm2.1}, this gives that the functions
$F, G\colon\mathbb{R}\to\mathbb{R}$ defined by
\[
 F(x)=f(x)-f(1)x
\quad
\text{and}
\quad
G(x)=g(x)-g(1)x
\quad
\left(x\in\mathbb{R}\right)
\]
are derivations.
\end{proof}

\subsection*{A characterization of linearity}

Finally, in the last part of this paper we will present a characterization of linearity.
Just as in the proof of Theorem \ref{T2.2}, Theorem \ref{Thm2.1} will again play an important role.

\begin{thm}
 Let $n\in\mathbb{N}, \, n\neq 1$ and $f\colon\mathbb{R}\to\mathbb{R}$ be an additive function.
Define $\phi$ on $\mathbb{R}$ by
\[
 \phi(x)=f\left(x^{n}\right)-f(x)^{n}
\qquad
\left(x\in\mathbb{R}\right).
\]
Let us assume that $\phi$ is locally regular.
Then the function $f$ is linear, that is,
\[
 f(x)=f(1)x
\]
holds for all $x\in\mathbb{R}$.
\end{thm}
\begin{proof}
 Let us define the function $\Phi\colon\mathbb{R}^{n}\to\mathbb{R}$ by
\[
 \Phi(x_{1}, \ldots, x_{n})=f\left(x_{1}\cdots x_{n}\right)-f(x_{1})\cdots f(x_{n})
\qquad
\left(x_{1}, \ldots, x_{n}\in\mathbb{R}\right).
\]
Due to the additivity of $f$, the function $\Phi$ is a symmetric, $n$-additive function. Furthermore,
\[
 \Phi(x, \ldots, x)=\phi(x)=f\left(x^{n}\right)-f(x)^{n}
\qquad
\left(x\in\mathbb{R}\right).
\]
From the local regularity of the function $\phi$  we immediately deduce that
$\phi$ is a continuous polynomial function.
Consequently,
\begin{equation}\label{Eq5}
 \Phi(x_{1}, \ldots, x_{n})=cx_{1}\cdots x_{n}
\qquad
\left(x_{1}, \ldots, x_{n}\in\mathbb{R}\right)
\end{equation}
holds with a certain $c\in\mathbb{R}$. Due to the definition of the
function $\phi$, we have $\phi(1)=f(1)-f(1)^{n}$. On the other hand
\[
 \phi(1)=\Phi(1, \ldots, 1)=c.
\]
Hence $c=f(1)-f(1)^{n}$.
Let $u\in\mathbb{R}$, with the substitution
\[
 x_{1}=u, \qquad x_{i}=1 \, \text{for $i\geq 2$},
\]
equation \eqref{Eq5} yields that
\[
 f(u)-f(u)f(1)^{n-1}=\left(f(1)-f(1)^{n}\right)u
\qquad
\left(u\in\mathbb{R}\right).
\]
In case $f(1)^{n-1}\neq 1$, this furnishes
\[
 f(u)=f(1)u
\qquad
\left(u\in\mathbb{R}\right).
\]
If $f(1)^{n-1}=1$, then
\[
 c=f(1)-f(1)^{n}=f(1)\left[1-f(1)^{n-1}\right]=0.
\]
Therefore equation \eqref{Eq5} with the substitutions
\[
 x_{1}=u, \quad x_{2}=v, \quad \text{and} \quad x_{i}=1 \quad \text{for $i\geq 3$}
\qquad
\left(u, v\in\mathbb{R}\right)
\]
yields that
\[
 f(uv)=f(u)f(v)f(1)^{n-2}
\qquad
\left(u, v\in\mathbb{R}\right),
\]
that is, $f(1)^{n-2}\cdot f$ is a non identically zero real homomorphism.
In view of Kuczma \cite[Theorem 14.4.1.]{Kuc09} this implies that
\[
 f(1)^{n-2}f(u)=u
\qquad
\left(u\in\mathbb{R}\right).
\]
Since $f(1)^{n-2}\cdot f(1)=f(1)^{n-1}=1$,
\[
 \frac{f(u)}{f(1)}=u
\]
holds for all $u\in\mathbb{R}$, that is, $f$ is a linear function, indeed.
\end{proof}

\subsection*{Acknowledgements}
This paper is dedicated to the 60\textsuperscript{\textrm{t}h} birthday of
Professor L\'{a}szl\'{o} Sz\'{e}kelyhidi (University of Debrecen, Hungary).
Furthermore, the author wishes to express her gratitude to the two anonymous referees
for their work.


\section{Introduction}
\label{}

Scafetta has argued that there is a strong coherence between the global temperature series over the last 160 years and the movement data of the sun around the barycenter (center of mass) of the solar system. This claim is based on coincidence between spectral peaks of the two time series. Such a coherence would, if proven to be correct, go  a long way towards showing that there is a celestial origin of the climate oscillations on the 10-62 year scale.

The claimed coherence is presented in \citet{Scafetta2010} and it depends strongly on the result of maximum entropy spectral analysis of the climate series and the identification of specific peaks in the spectra. In that paper a total of 11 peaks with periods ranging from 6 years and above are identified and approximately matched with a planetary or moon--earth oscillation.

In a follow-up work, only the five main periodicities from period 9 years and above are discussed \citep{Scafetta2012shared}. We will not discuss the 9-year periodicity claimed to correlate with the lunar orbit here, so the four periodicities which are claimed to have an origin in the solar orbit are at 10-10.5, 20-21, 30 and 60-62 years.

The underlying hypothesis of these papers is that there are periodicities in the climate series which correspond to periodicities in the solar system where the movement of the sun around the barycenter is used as a proxy. The logical consequence is drawn from this that the better the climate series periodicities can be understood and in particular the better their origin can be understood, the more they can be be used for predicting future climate development. This last point is elaborated in the papers \citep{Scafetta2012b,solheim2012long} which will not be discussed here.

The purpose of this paper is to reanalyze the data for coherence. The main tools are time--frequency spectral analysis and the magnitude squared coherence.

\section{Analysis Methods}

\subsection{The maximum entropy method}

Maximum entropy spectral analysis  is a high resolution method which fits an all-pole spectral model to the data. It is able to bring out spectral peaks that many other methods cannot show. One of the the most used numerical implementations is due to \citet{burg1967maximum}. During the 70's the Burg method was investigated for its properties and many applications were explored. One example application which we have been involved in is ocean wave spectral analysis \citep{Holm1979}.

When the maximum entropy method (MEM) is applied one needs to be aware of its properties and in particular some pitfalls. One of them is that if a too high order is used, the spectral peaks will become sharper and sharper and eventually spurious peaks will appear. \citet{ghil2002advanced}, the source for the analysis methods of \citet{Scafetta2010}, also warn against using a too high order saying that an upper bound is generally taken as half the number of samples. An even stricter criterion is the recommendation that the maximum order should be between one-third and one-half of the number of samples \citep{kay1979effects}.

A second property is that too high order increases the chance that spectral peaks may split into two \citep{fougere1976spontaneous}. Third, there can be bias in the frequency location of peaks especially for frequencies where only a few periods are contained in the time series. The bias may depend on initial phase \citep{chen1974experiments}.

Fourth, \citet{ghil2002advanced} caution against applying MEM uncritically to nonstationary time series or data which are not well approximated by an AR process and advice that cross testing with other spectral analysis methods then is particularly important. This is particularly relevant here as the periodicities will turn out to be time-varying in the climate series data.

Although MEM is a high resolution spectral analysis method, high resolution is not the same as high accuracy. Resolution shows how well two lines can be distinguished from each other. This is not so much of an issue for the data analyzed here. Therefore we will rely to a large extent on the periodogram and following \citet{ghil2002advanced} will validate MEM against the windowed periodogram, although MEM will have finer widths and seemingly better accuracy.

\subsection{The windowed periodogram}

The periodogram has the desirable property that for a single sinusoid or for well isolated sinusoids in noise, it is the maximum likelihood estimator of frequency \citep{kay1988modern}. But one needs to be cautious in applying the periodogram as both the length of the window and the choice of a window function may be crucial for good performance. This is especially the case here where there may be very strong spectral peaks that can leak into and obscure other interesting parts of the spectrum. Therefore a Kaiser window is used here since its performance can be adjusted by a single parameter, $\beta$. Note that Kaiser's original definition uses the $\alpha=\beta/\pi$ parameter \citep{kaiser1974nonrecursive}, while we follow the convention in Matlab (\textregistered Mathworks Inc., Natick, MA) which characterizes the window by $\beta$. The larger the value, the lower the sidelobes and the poorer the frequency resolution will be. Prior to windowing, the data is detrended by removing any linear trend. The main purpose of that is to remove any zero-frequency components that may leak into the interesting high period part of the spectrum.

Another consideration is that nonlinear processing of the amplitude of the spectral estimator may make features more visible. \citet{kay1984high} show how one spectral estimator with what seems to be low resolution may be transformed into another which in the display seems to have higher resolution. The transform is a monotonic transform. Inspired by this, we have used a power law function to transform the periodogram into a visually as attractive plot as possible using the transformation:
\begin{equation}
T(P_{xx}(f)) = P_{xx}(f)^{a},
\label{eq:fft}
\end{equation}
where $P_{xx}$ is the power spectrum estimate and $0<a\le1$. Only the periodogram output is transformed in this way, the Burg MEM result is unprocessed ($a=1$) and converted to dB.

\subsection{Time-Frequency Analysis}

It will be clear later that the global temperature data in particular is not time-invariant. As spectral estimation methods that assume stationarity will produce many spurious peaks when applied to nonstationary data, this is the main criticism of the results of \citet{Scafetta2010}.

In the time--frequency analysis, we have used Burg's Maximum Entropy Method and the windowed periodogram with an offset of 1 year between each analysis window. The length of the data window must be chosen and the longer the window, the better the frequency resolution. On the other hand, the shorter the window, the better the ability to track time variations in the data. For the examples shown here, a window length of 60 years was found to be a reasonable compromise.

Such a compromise between time and frequency resolution is the background for the development of multi-resolution methods such as wavelets. They often have poorer frequency resolution and better time resolution at high frequencies compared to the low frequencies. We have chosen not to use wavelets here, primarily because it is not used in \citet{Scafetta2010,Scafetta2012shared}, but also because it may often be hard to understand and evaluate the performance and accuracy of these methods, such as the time- and frequency-resolution in the various frequency bands. This may also make it hard to separate features of the data from properties of the analysis method.

The order used for MEM analysis is found by trial and error and comparison with the periodogram as advised in  \citet{ghil2002advanced}, keeping in mind that it should not exceed between $1/3$ and $1/2$ of the number of samples and that the window length should not smooth out nonstationarities. Based on a number of analyses where Figs.\ 3 and 6A of \citet{Scafetta2010} were reproduced, we ended up choosing order 300 for a time series of length 60 years. This is $300/(12\cdot 60) \approx 42\%$ of the data length which is somewhat more conservative than the order 1000 in \citet{Scafetta2010} which is equivalent to $1000/(12 \cdot 160) \approx 52\%$. Such a high order in combination with a window that spans the entire data record of 160 years contributes to additional peaks in the spectra shown in Figs.\ 3 and 6A of \citet{Scafetta2010}. Although in a later paper Scafetta has realized that half the length of the record is the highest allowed pole order for MEM,  \citet{scafetta2012multi}, there is still no discussion of the smearing effect of analyzing the whole record in one analysis window.

\subsection{Magnitude squared coherence estimation}
The magnitude squared  coherence (MSC) is a well-established correlation factor in the frequency domain. It has been used for decades as the best method for identifying coupling between signals \citep{bendat1980}. The MSC is defined as the normalized squared cross spectrum between two time series $x(t)$ and $y(t)$

\begin{equation}
C_{xy}(f) = \frac{|P_{xy}(f)|^2}{P_{xx}(f)\cdot P_{yy}(f)},
\label{Eq:msc}
\end{equation}
where $P_{xx}(f)$ and $P_{yy}(f)$ are the power spectra of the two time series and $P_{xy}(f)$ is the complex cross spectrum. $C_{xy}(f)$ will always be between 0 and 1. In order to estimate it, both time series have to be divided into several weighted and partly overlapping segments which are Fourier transformed and averaged.

The estimator finds relationships between data based on linear systems theory where the effect, $y(t)$, is a filtered version of the cause, $x(t)$, with noise added. The smaller the contribution of the noise at a particular frequency, the closer the MSC will be to 1. If on the other hand there is a chaotic system coupling between input and output, the MSC will not be able to discriminate between that and additive noise. In such cases it is not so useful. However, there is nothing in \citet{Scafetta2010} that indicates chaotic coupling, so the MSC should be an appropriate measure to use.

Statistics of the estimator is given in \citet{carter1973statistics} which says that the maximum bias, the maximum standard deviation, and the maximum rms error only decay as $N_d^{-1/2}$, where $N_d$ is the number of disjoint segments or independent averages. If only one segment is used to cover the whole length of the time series, the estimator will always give $C_{xy}(f)=1$ independent of the data and thus the estimate is useless. The maximum bias is shown by \citet{carter1973statistics} to be approximately bounded by $0.5 \sqrt{\pi/N_d}$ and the largest bias occurs for values of MSC close to 0. The number of independent averages should therefore be large in order to get a good estimate, on the order of $N_d=10$ or more.

Since a Kaiser window which tapers the ends heavily down is used, 75\% overlap is used between the $N$ segments. The number of segments with this overlap is:

\begin{equation}
N = 1 + \lfloor{\frac{N_{tot}-N_{seg}}{0.25 N_{seg}}}\rfloor
\label{Eq:N}
\end{equation}
where $N_{tot}$ is the total number of monthly samples and $N_{seg}$ is the analysis interval, e.\ g.\ $N_{seg}=40\cdot12$ for a 40 year interval.

The  actual value for the number of independent estimates, $N_d$, will be less than $N$ and is estimated from the correlation between the data windows at the 50\% and 75\% points ($c(0.5)$ and $c(0.75)$) as given by \citet{harris1978use}:

\begin{equation}
N_d \approx N/[1+2 c^2(0.75) + 2c^2(0.5)]
\label{Eq:Nd}
\end{equation}

Confidence intervals for the MSC have been estimated using the method of \citet{wang2004exact}. The independence threshold has also  been estimated using the method of \citet{carter1987coherence} and \citet{wang2004optimising}. Coherence estimates below this value indicates that the two processes are independent. It is:

\begin{equation}
IT = 1-(1-P_d)^{1/(N_d-1)}
\label{Eq:Indep}
\end{equation}
where $P_d$ is the confidence level, $0.95$ or $0.99$ in our case.

It should be mentioned that the coherence measure used here is different from the measure mentioned in  \citet{Scafetta2010,Scafetta2012shared}. In the first paper, although a term which is called coherence is listed and discussed, there is little documentation for how it is estimated beyond a comparison of spectral peaks. The coherence measure in \citet{Scafetta2012shared} is based on the joint power statistics or the cross-spectrum. Therefore it is close to the numerator  of Eq.\ (\ref{Eq:msc}). However, without the normalization by the power spectra in the denominator, this measure may indicate large peaks where there is large energy, regardless of whether there is spectral correlation or not. Thus it is not a good indicator of coherence.

\subsection{Implementation}
Only standard functions of Matlab\textregistered have been used. The MEM-method is \texttt{pburg.m}, the periodogram is estimated by using \texttt{periodogram.m}, the window function is \texttt{kaiser.m}, and detrending is performed by using \texttt{detrend.m}. The function \texttt{mscohere.m} is used for magnitude squared coherence estimation.

The number of samples in the frequency domain was set to a very large number ($2^{16}= 65536$) in order to interpolate the spectral plots and avoid sampling artifacts. With monthly sampling, the maximum frequency supported is 6 year$^{-1}$, while the periods of interest between 60 and 6 years are in the frequency range from $60^{-1} \approx 0.017$ to $6^{-1} \approx 0.17$ year$^{-1}$. Thus only the $6^{-1}/6 \approx 2.8 \%$ lowest part of the spectrum is really of interest. Despite the interpolation there will be only about $65536 \cdot 2.8\%/2 \approx 918$ samples in this interval, which is a reasonable number.

\section{Data Analysis}

\subsection{Temperature data}

The data is the monthly sampled global surface temperature (HadCRUT3) series from 1850 to 2010.75 ($N_{tot}=1929$ samples). The dataset is a combination of land and marine temperature anomalies on a $5^{\circ}$ by $5^{\circ}$ grid-box basis, see plot in Fig.\ 1 of \citet{Scafetta2010} and description in \citet{brohan2006uncertainty}.

The time--frequency spectrogram of the temperature data is shown in Figs.\ \ref{fig:Temp-spectrogram60} and \ref{fig:Temp-spectrogramMEM60}. The periodogram is transformed with $a = 0.4$ in Eq.\ (\ref{eq:fft}) and the MEM estimate uses order 300. Both use a window length of 60 years which is a length that does not seem to smear the nonstationarity of the data too much. It should be noted that there is more temporal variation in these data than in e.g.\ sunspot data, so a shorter data window than that used for MEM analysis of sunspot data \citep{currie1973fine} must be used. A Kaiser window with $\beta=3.0$ was used for the periodogram. This is a mild form of windowing with first sidelobe at $-23.8$ dB. Both figures show in general the same result, although the MEM plot has sharper and sometimes more continuous lines.

The 11 periodicities claimed by \citet{Scafetta2010} are denoted by arrows in the margin. From the time--frequency spectrogram only three of these can be recognized with some confidence. The strongest one is at about 60 years and is close to time-invariant. It is in agreement with the 65-70 year cycle reported by \citet{schlesinger1994oscillation}. There is also a line in the 15-20 year range which varies with time as well as one around 9 years which comes and  goes and varies in frequency.

\begin{figure}[tb]
\begin{center}
	\includegraphics[width=.9\columnwidth]{Temp_spectrogram_60yrs_bar}
	\caption{Periodogram-based temperature spectrogram, 60 year window and Kaiser window with $\beta=3$ and a shift of 1 year between windows, $a=0.4$. Period in years as estimated by \citet{Scafetta2010} shown with arrows on right-hand side.}
	\label{fig:Temp-spectrogram60}
\end{center}
\end{figure}

\begin{figure}[tb]
\begin{center}
	\includegraphics[width=.9\columnwidth]{Temp_spectrogramMEM_60yrs_bar}
	\caption{MEM-based temperature spectrogram, 60 year window, order 300, shift 1 year between windows. Period in years as estimated by \citet{Scafetta2010} shown with arrows on right-hand side.}
	\label{fig:Temp-spectrogramMEM60}
\end{center}
\end{figure}

\subsection{Solar movement data}

The solar orbit data has been taken from the Horizons system of JPL  by generating equal spaced monthly data between 15 January 1750 and 16 December 2100. Prior to this analysis the data set was confirmed to reproduce Fig.\ 4 and the spectrum of Fig.\ 6A of \citet{Scafetta2010}.

The MEM spectrogram for speed of the center of the mass of the solar system (SCMSS) data is shown in Fig.\ \ref{fig:SCMSS-spectrogramMEM60} with a window length of 60 years. The 10 periodicities of \citet{Scafetta2010} are denoted by arrows in the margin. Due to the near time-invariance of this dataset, several of these are also visible in the time--frequency spectrogram.

The strongest stationary component is that at about 20 years which can be explained by the period for which the position of the Sun, Jupiter, and Saturn are re-aligned (the synodic period of Jupiter and Saturn). There is also a fairly consistent periodicity at 11.9 years which is the revolution period of Jupiter. There are weaker almost constant-frequency lines at 9.8, 7.5, and 5.9 years also.

\begin{figure}[bt]
\begin{center}

	\includegraphics[width=.9\columnwidth]{SCMSS_Burg_60yrs_Correct-Fig3_}
	\caption{MEM-based spectrogram of the SCMSS (speed of the center of mass of the solar system), 60 year window, order 300, shift 1 year between windows. Period in years as estimated by \citet{Scafetta2010} shown with arrows on right-hand side.}
	\label{fig:SCMSS-spectrogramMEM60}
\end{center}
\end{figure}

Although there is some overlap between the frequencies of the main lines of Figs. \ref{fig:Temp-spectrogramMEM60} and \ref{fig:SCMSS-spectrogramMEM60}, the time-varying nature of the lines in the temperature data is quite different from the almost constant frequencies in the SCMSS data. The importance of these differences will be even clearer when the magnitude squared coherence is estimated.

\subsection{Coherence Analysis}
The magnitude squared coherence (MSC) is estimated by averaging the Fourier transforms of windowed, overlapped segments of the data. As the SCMSS data has strong periodicities it is necessary to use a Kaiser window with $\beta=7$, i.e.\ with a first sidelobe as low as $-51$ dB in order to avoid spectral leakage. This can also be seen from  the large dynamic range of Fig.\ \ref{fig:SCMSS-spectrogramMEM60} compared to that of Fig.\ \ref{fig:Temp-spectrogramMEM60}.  Due to the strong tapering effect, 75\% overlap is used for the MSC analysis.

First an analysis is done over windows of length 20 years. Over an interval of 160 years this gives $N=29$ analysis windows, Eq.\ (\ref{Eq:N}). Using the properties of the window and the overlap, Eq.\ (\ref{Eq:Nd}) gives as a result that there are $N_d=15.9$ degrees of freedom in the estimate of MSC with this particular window. This gives an independence threshold of $IT=0.18$ at the 95\% confidence level and $IT=0.27$ at the 99\% level. The coherence between the solar movement (SCMSS) and the global temperature data in Fig.\ \ref{fig:MSC-SCMSS-Temp} is not above 0.15 for any frequency. This is a very low value which is below the independence thresholds.

The analysis was then repeated with segment lengths of 30, 40, and 50 years. The resulting decrease in the number of independent estimates, $N_d$, the position and the peak value of the magnitude square coherence and the independence threshold at the 95  and 90\% confidence levels are shown in Table \ref{table:MSC}.

\begin{table}[t]
\footnotesize
\centering
\begin{tabular}{c | c c c c c c}
$N_{seg}/12$  & $N$ &$N_d$ & Peak    & Peak & IT & IT\\
(years)     &        &     & (years)&  MSC & 95\%        & 99\% \\
\hline
20  & 29 & 15.9 & 11.7 & \bf{0.14} & 0.18 & 0.27  \\
30  & 18& 9.9   & 14.9 & \bf{0.35} & 0.29 & 0.40  \\
40  &13 &7.1   & 16.8 & \bf{0.59} & 0.39 & 0.53 \\
\hline
\end{tabular}
\caption{Number of averages $N$, number of independent averages $N_d$, MSC peak location and value, and independence thresholds (IT) at the 95\% and 99\% confidence levels  as a function of segment length.}
\label{table:MSC}
\end{table}

The maximum MSC-value increases with window length, but so does the uncertainty as the number of independent averages, $N_d$, decreases. When 30 year windows are used with the same data window and the same overlap there is only $N_d=9.9$ degrees of freedom and the independence threshold increases to $0.29$ at the 95\% confidence level and $0.4$ at the 99\% level, see Fig.\ \ref{fig:MSC-SCMSS-Temp-30}. In this case there is a peak in the coherence at $14.9$ years with a value of $0.35$ which is slightly above the 95\% confidence level, but below the 99\% confidence level.

\begin{figure}[tb]
\begin{center}
	\includegraphics[width=.9\columnwidth]{MSC-SCMSS-Temp-20yrs-years}
	\caption{Magnitude squared coherence between speed of center of mass of the solar system (SCMSS) and global temperature data (solid line) with segment length 20 years. The 95\% confidence intervals (dash-dot lines) and the 95\% independence threshold of 0.18  (dashed line) are also shown. Peaks are at 5.3 and 11.7 years.}
	\label{fig:MSC-SCMSS-Temp}
\end{center}
\end{figure}

\begin{figure}[tb]
\begin{center}
	\includegraphics[width=.9\columnwidth]{MSC-SCMSS-Temp-30yrs-years}
	\caption{Magnitude squared coherence between speed of center of mass of the solar system (SCMSS) and global temperature data (solid line) with segment length 30 years. The 95\% confidence intervals (dash-dot lines) and the 95\% independence threshold of 0.29 (dashed line) are also shown. Peaks are at 6.3 and 14.9 years.}
	\label{fig:MSC-SCMSS-Temp-30}
\end{center}
\end{figure}

The value increases when 40 year windows are used and there is a peak in the coherence at $16.8$ years which is above the independence thresholds, see Table \ref{table:MSC}. However, as the number of independent averages is only $N_d=7.1$ the reliability of this result is less than those for the lower window lengths.

\section{Discussion and Conclusion}
\citet{Scafetta2010} claimed the global temperature series for the last 160 years to have  spectral lines at 21, 30 and 62 years. Time--frequency analysis shows that the lines are time-varying (Figs.\ \ref{fig:Temp-spectrogram60} and \ref{fig:Temp-spectrogramMEM60}) and very different from the nearly constant lines in the time--frequency plot for the speed of the center of mass of the solar system (SCMSS) (Fig.\ \ref{fig:SCMSS-spectrogramMEM60}).  The supposed periodicity around 30 years in \citet{Scafetta2010} is not really present in the climate series at all and could be an artifact due to a combination of model overfitting and smearing due to the time-invariance assumption which has been forced on the data. The claimed spectral peaks by \citet{Scafetta2010} for the global temperature series are therefore not reproducible if proper consideration is taken of the time-varying nature of the data.
The only significant coherence between the climate series and the sun's movement that was possible to find was at 15-17 years (Table \ref{table:MSC}). However, both the low number of independent averages that it builds on as well as the lack of a physical explanation for this coherence, makes us hesitate to claim that it is significant.

\section{Acknowledgement}

I would like to thank Prof.\ ShouYan Wang, Suzhou Institute of Biomedical Engineering and Technology for stimulating discussions and for sharing his code for confidence interval estimation for the MSC. I also thank  Prof.\ F.\ Albregtsen and Dr.\ A.\ Austeng, University of Oslo, as well as Prof.\ K.\ Rypdal and Dr.\ M.\ Rypdal, University of Troms\o, Norway for helpful discussions during the work with this paper.

\bibliographystyle{elsarticle-harv}


\section{Introduction}
We  consider the  Fokker-Planck
 equation
 \begin{equation}\label{1}
 \frac{\partial{f}}{\partial{t}}=\Delta_{v}{f}-v\cdot \bigtriangledown_{x}{f}+
 {\mathrm{div}}_{v}\big(vf(1+kf)\big),\end{equation}
 where  $t{\in}{\mathbb R}$ stands for time,   $x=(x_{1},x_{2},...,x_{n}){\in}\Omega$ stands for
 space coordinates,  $v=(v_{1},v_{2},...,v_{n})\, {\in} \, {\mathbb R}^{n}$ is velocity
 and $\, {\mathbb R} \,$ denotes the real axis. This non-linear Fokker-Planck equation serves  as a kinetic model for bosons $(k>0)$ and fermions $(k<0)$.

The important notion of Kullback-Leibler distance \cite{KL} is essentially generalized in our paper and new
conditional extremal problems, which appear in this way, are solved. The solutions $f(t,x,v)$
of the Fokker-Planck equation are studied in the bounded domain
$\Omega$ of the $x$-space. Such an approach essentially changes the usual situation, that is,
the total energy depends on $t$ and
the notion of distance (between
a stationary solution and an  arbitrary solution of the Fokker-Planck equation)
 includes the $x$-space.  Thus,  the notion of distance remains well-defined also in the spatially inhomogeneous case.
 Recall that the Kullback-Leibler distance, which has many applications (see, e.g., \cite{Hab, SH, Vi2} and references therein),
 is  defined only in the spatially homogeneous case.

 In our previous paper \cite{ASLS} we studied a model case of the one dimensional $x$-space.
 Here the case $\dim \Omega \geq 1$ is dealt with.
Furthermore, using our generalization of the Kullback-Leibler distance,  we show
local stability of global Maxwellians
 in the spatially inhomogeneous case.

 The comparison of the energy and entropy in the classical and quantum cases is an important subject
(see  \cite{Caret, LAS1, LAS2, LAS3, LASPLet, LAS4, We} and references therein).
 Here, we compare these energy and entropy for the situation described by the Fokker-Planck
 equation. It is especially interesting for the applications that the fermion and  boson cases are essentially different.

Our definition  of the quantum entropy $S_k$ ($k\not=0$) is slightly different from the previous definitions (see \cite{Dol, Lu2}).
We  show that the natural requirement
\begin{equation}\label{q1}
S_{k}{\to}S_{c},\quad k{\to}0 \quad (S_c=S_0 \,\,\mathrm{is} \,\, \mathrm{the} \,\,
\mathrm{classical} \,\, \mathrm{entropy})
\end{equation}
is not fulfilled in the case of the old definition, however \eqref{q1}  holds for   our modified definition (see Section \ref{Se3}).
Some necessary definitions are given in Subsection \ref{Prel}.
An important functional, which attains its maximum at the
function $M_{k}$ is introduced there.
The distance between solutions and the corresponding extremal problem
are studied in Sections \ref{Se3} and \ref{Lya}.
Our results on the Fokker-Planck equation are mostly related to the corresponding results
 (from \cite{LASPLet, LAS4})
on the Boltzmann equation but the theorems from Section \ref{Lya}
have no analogs in the Boltzmann case.

We use the standard notation $|v|=\sqrt{v_1^2+\ldots + v_n^2}$ and $C_{0}^{1}$  denotes the class
of differentiable functions $f(x,v)$, which tend to zero sufficiently rapidly
when $v$ tends to infinity.
\section{The extremal problem} \label{Se3}
 \subsection{Preliminaries}\label{Prel}
Here, we  present some  well-known  notions and results
 connected with   equation \eqref{1}.
It is required that the distribution function $f(t,x,v)$ satisfies the inequalities
\begin{equation}\label{3}
f(t,x,v) \geq 0, \quad 1+kf(t,x,v) \geq 0 \quad (k \in {\mathbb R}),
\end{equation}
and we set $f\log f=0$ for the case that $f=0$ and $(1+kf)\log(1+kf)=0$ for $1+kf=0$.
Then, the mapping
\begin{align}\label{5'}&
\Phi(f):=f\log{f}\quad {\mathrm{for}} \quad k=0, \\
\label{5}&
 \Phi(f):=f\log{f}-\frac{1}{k}(1+kf)\log(1+kf)+f \quad {\mathrm{for}} \quad k\not=0,
 \end{align}
 is well-defined. Now, the entropy is given by the equality
 \begin{equation}\label{4}
 S(t,f)=S(t)=-\int_{\Omega}\int_{{\mathbb R}^n}\Phi(f)d{v}dx.
\end{equation}
The notions of density $\rho(t,x)$,  total density $\rho(t)$, mean velocity $u(t,x)$,
energy $E(t,x)$, and total  energy $E(t)$ are introduced via formulas:
 \begin{align}\label{7}&
 \rho(t,x)=\int_{{\mathbb R}^n} f(t,x,v)dv , \quad \rho(t)=\int_{\Omega}\rho(t,x){dx},\\
 \label{8}&
 u(t,x)=\big(1/\rho(x,t)\big)\int_{{\mathbb R}^n} v f (t,x,v){d}v,
 \\ \label{9}&
 E(t,x)=\int_{{\mathbb R}^n} \frac{|v|^{2}}{2}f(t,x,v)dv , \quad
E(t)=\int_{\Omega}\int_{{\mathbb R}^n} \frac{|v|^{2}}{2}f(t,x,v)dv dx.
\end{align}
We assume that the domain $\Omega$ is bounded, and so its volume
is also bounded:
\begin{equation}\label{11}
\mathrm{Vol}(\Omega)=V_{\Omega}<\infty.
\end{equation}
 \subsection{The free energy functional and extremal problem}\label{Free}
We  introduce the "free energy" functional
\begin{equation} \label{12}
F(f)=F\big(f(t)\big)=S(t) - E(t),
\end{equation}
 where  $S(t)$ and $E(t)$ are defined by formulas \eqref{4} and \eqref{9}, respectively.
Next,   we use the calculus
of variations (see \cite{Ha}) and find the function $f_{max}$ which maximizes
the functional \eqref{12}, where the parameters  $\, t\,$ and $\rho(t)=\widetilde \rho>0$ are fixed.
 The corresponding Euler's equation takes the form
 \begin{equation} \label{14}
 -\frac{|v|^{2}}{2}- \log{f}+\log(1+kf)+
{\mu}=0.
\end{equation}
From the last relation we obtain
 \begin{equation} \label{15}
 f/(1+kf)=C\exp\left\{-{|v|^{2}}/({2})\right\}, \quad C:=\mathrm{e}^{-\mu}.
 \end{equation}
 Formula \eqref{15} implies  that
 \begin{equation} \label{16}  f=M_{k}=\frac{C\exp\left\{-{|v|^{2}}/2\right\}}
{1-kC\exp\left\{-{|v|^{2}}/2\right\}},
\end{equation}
that is, $f$ coincides with the global Maxwellian $M_k$.

In view of the requirement $\rho(t)=\widetilde \rho$, the constant $C$ in the equality \eqref{16} is derived from the relation
 \begin{equation} \label{16'}
 V_{\Omega}\int_{{\mathbb R}^n}\frac{C\exp\left\{-{|v|^{2}}/2\right\}}
{1-kC\exp\left\{-{|v|^{2}}/2\right\}}dv=\widetilde \rho.
\end{equation}
The function  $f$ given by \eqref{15} (or, equivalently, by \eqref{16}) is nonsingular and satisfies conditions \eqref{3} and $\rho(t)=\widetilde \rho>0$
if and only if
\begin{equation}\label{op1}
C>0, \quad 1-kC>0.
\end{equation}
In Subsection \ref{comp} we prove that there is a unique value $C$ satisfying  relations
\eqref{16'} and \eqref{op1}.  Let us show that $F$ indeed attains  its maximum on the global Maxwellian $M_k$
corresponding to such $C$.
According to \eqref{4}, \eqref{9} and \eqref{12}  the "free energy" $F$ admits the representation
\begin{equation}\label{op2}
F=\int_{\Omega}\int_{{\mathbb R}^n}\Psi(f)dvdx, \quad \Psi:=-\frac{|v|^2}{2}f-\Phi.
\end{equation}
Taking into account \eqref{5'}, \eqref{5}, \eqref{op1} and \eqref{op2}, we have the inequality
 \begin{equation} \label{17}  \frac{\delta^{2}}{\delta{f}^{2}}\Psi=-\frac{1}{(1+kf)f}<0 ,
 \end{equation}
and the next proposition follows.
\begin{Pn}\label{Pn1} Under condition \eqref{op1}, the functional  $F$ given by \eqref{12}
attains its maximum on the function $M_k$ of the form \eqref{16}
$($where $C$ is defined in \eqref{16'}$)$, that is,
 \begin{equation} \label{18}
 G(f)=F(M_{k})-F(f)>0\quad (f{\ne}M_{k}).
\end{equation}
\end{Pn}
\begin{Rk} In this subsection we introduced the important functional $F$, which attains its maximum
on the global Maxwellians and the generalization $G$ $($see \eqref{18}$)$ of the Kullback-Leibler distance.
As opposed to the Kullback-Leibler distance, which is defined in the $x$-homogeneous case, the distance $G$ is well-defined for the functions $f$,
which depend on $x$.

Later in this section we will consider the energy $E$, entropy $S$ and free energy $F$ of the
global Maxwellians. The following sections are dedicated to the study of the general solutions
of the Fokker-Planck equation.
\end{Rk}
\subsection{Comparison of the classical and quantum characteristics}\label{comp}
 Let us calculate the integral on the left-hand side of \eqref{16'}. Using  spherical coordinates, we have
 \begin{equation} \label{3.9}
C \int_{{\mathbb R}^n}\frac{\exp\left\{-{|v|^{2}}/2\right\}}
{1-kC\exp\left\{-{|v|^{2}}/2\right\}}dv=\omega_{n-1}C\int_{0}^{\infty}
\frac{r^{n-1}\exp\left\{-{r^{2}}/2\right\}}
{1-kC\exp\left\{-{r^{2}}/2\right\}}dr,
\end{equation}
where the surface area of the $(n-1)$-sphere of radius $1$ is
 \begin{equation} \label{3.10}
 \omega_{n-1}=\frac{2{\pi}^{n/2}}{\Gamma(n/2)},
 \end{equation}
and $\Gamma(z)$ is the gamma function. Euler's integral representation of the gamma function
easily yields
 \begin{equation} \label{3.11}
 \int_{0}^{\infty}\mathrm{e}^{-ar^{2}}r^{n-1}dr=\frac{1}{2}a^{-n/2}\Gamma(n/2) \quad (a>0).
 \end{equation}
 \begin{Nn}\label{NnC} According to \eqref{16} and \eqref{16'}, the value $C$  corresponding to $M_k$
 depends on $k$.  We denote this value by $C_k$.
 \end{Nn}
 Taking into account \eqref{16'} and \eqref{3.9}--\eqref{3.11},
we obtain
 \begin{equation} \label{3.12}
 (2{\pi})^{n/2}V_{\Omega}C_kL_{n/2}(kC_k)=\widetilde \rho,
\end{equation}
where
 \begin{equation} \label{3.13}
  L_{n/2}(z)=\frac{2^{1-(n/2)}}{\Gamma(n/2)}\int_{0}^{\infty}\frac{\mathrm{e}^{-\frac{r^{2}}{2}}}
{1-z\mathrm{e}^{-\frac{r^{2}}{2}}}r^{n-1}dr=\sum_{m=1}^{\infty}\frac{z^{m-1}}{m^{n/2}}.
\end{equation}
We note that the series representation $L_{n/2}(z)=\sum_{m=1}^{\infty}\frac{z^{m-1}}{m^{n/2}}$
does not hold for $|z|>1$, and we use only the first equality in \eqref{3.13} for the case that $z< -1$.
Using the first equality in  \eqref{3.13}, we derive the following statement.
\begin{Pn}\label{Proposition 3.1.}
The function $L_{n/2}(z)$  increases strictly monotonically in the interval $-\infty \leq z<1$ and
 \begin{align} \label{3.13'} &
L_{n/2}(0)=1;
 \\ & \label{3.14}
 L_{1/2}(1)=L_{1}(1)=\infty;\quad L_{n/2}(1)<\infty,\,\, n>2.
 \end{align}
 \end{Pn}

Proposition \ref{Proposition 3.1.} implies the next two corollaries.

\begin{Cy}\label{CyB}
If $k>0$ $($boson case$)$ and either $n=1$ or $n=2$, then
equation \eqref{3.12} has one and only one solution $C_k$ such that $C_k>0$ and $kC_k<1$.
\end{Cy}

\begin{Cy}
\label{Cy 2.5}
If $k>0$ $($boson case$)$, $n>2$ and
 \begin{equation} \label{3.15}
 (2{\pi})^{n/2}V_{\Omega}L_{n/2}(1)>k\widetilde \rho,
 \end{equation}
then
equation \eqref{3.12} has one and only one solution $C_k$ such that $C_k>0$ and $kC_k<1$.
\end{Cy}

\begin{Rk}
The function $L_{n/2}(z)$ belongs to the class of the $L$-functions \cite{LaGa} and
is connected with the famous $($see, e.g., \cite{T}$)$ Riemann zeta-function
 \begin{equation} \label{3.16}
 \zeta(z)=\sum_{p=1}^{\infty}\frac{1}{p^{z}};\quad \Re{z}>1
 \end{equation}
 by the relation
 \begin{equation} \label{3.17}
 L_{n/2}(1)=\zeta(n/2).
 \end{equation}
Hence, we have the equalities
 \begin{equation} \label{3.18}
 L_{3/2}(1)=\zeta(3/2)=2.612,\quad L_{2}(1)=\zeta(2)=1.645,
 \end{equation}
 \begin{equation} \label{3.19}
 L_{5/2}(1)=\zeta(5/2)=1.341,\quad L_{3}(1)=\zeta(3)=1.202.
\end{equation}
\end{Rk}

Let us study the fermion case $k<0$.  Taking into account Notation \ref{NnC},
we may later  (differently from the constants $C_k$) consider $C$ as a variable.
In view of \eqref{3.13}, the next proposition is valid.

\begin{Pn}
\label{Proposition 3.2}
Assume that $k<0$. Then, $CL_{n/2}(kC)$ increases strictly monotonically with   respect to $C$ $\,(0 \leq C<\infty)$,
 and $CL_{n/2}(kC){\to}\infty$ for $C{\to}\infty$.
 \end{Pn}

\begin{Cy}\label{CyF}
If $k<0$ $($fermion case$)$, then
equation \eqref{3.12} has one and only one solution $C_k$ such that $C_k>0$.
\end{Cy}
The second inequality in \eqref{op1}  holds in the fermion case (i.e., in the case $\,\,C=C_k$ and $k<0$) automatically.

Finally,  we consider in this section
 the   energy of the global Maxwellian:
 \begin{align} \nonumber
 E_k=E(M_k)&=\frac{1}{2}\int_{\Omega} \int_{{\mathbb R}^n}|v|^2\frac{C_k\exp\left\{-{|v|^{2}}/2\right\}}
{1-kC_k\exp\left\{-{|v|^{2}}/2\right\}}dvdx
\\ \label{3.20}
&=\frac{1}{2}\omega_{n-1}V_{\Omega}C_k\int_{0}^{\infty}r^{n+1}
\frac{\exp\left\{-{r^{2}}/2\right\}}
{1-kC_k\exp\left\{-{r^{2}}/2\right\}}dr.
\end{align}
It is immediate from \eqref{3.12} that
 \begin{equation} \label{op3}
\widetilde \rho/ L_{n/2}(kC_k)=(2 \pi )^{n/2}V_{\Omega}C_k,
\end{equation}
and so, using \eqref{3.10} and \eqref{3.13}, we rewrite \eqref{3.20} in the form
 \begin{equation} \label{3.21}
 E_{k}=\left(\frac{n{\widetilde \rho}}{2}\right)\frac{L_{(n/2)+1}(kC_k)}{L_{n/2}(kC_k)}.
\end{equation}
We note that the corresponding classical energy $E_c$ (i.e., the energy for the case $k=0$) is given by the formula
 \begin{equation} \label{3.22}
 E_{c}=\frac{n{\widetilde \rho}}{2}.
 \end{equation}
The points $kC=\pm 1$ are called the critical points in boson and fermion theories.
(Recall that the series representation of $L_{n/2}(z)$ in \eqref{3.13}
does not hold for $|z|>1$.)
 \begin{Pn}\label{Proposition 3.3.}
{Let the  condition}
 \begin{equation} \label{3.23}
-1\leq  kC_k{<}1
 \end{equation}
be fulfilled. Then, we have the inequalities
 \begin{equation} \label{3.24}
E_{q,B}<E_{c} <  E_{q,F}  \quad {\mathrm{for}} \quad n \geq 1,
 \end{equation}
 where $E_{q}$ denotes the energy in the quantum case $($i.e., the case $k\not=0)$, $E_{q,B}$
 stands for the energy in the boson case $k>0$ and  $E_{q,F}$ stands for the energy in the fermion case $k<0$.
 \end{Pn}
\begin{proof}
Taking into account the second equality in \eqref{3.13}, we obtain:
 \begin{equation} \label{3.25}
  \frac{L_{(n/2)+1}(kC)}{L_{n/2}(kC)}<1\quad \mathrm{for}\quad k>0, \quad kC<1.
\end{equation}
Moreover, we will show that
 \begin{equation} \label{3.25'}
\frac{L_{(n/2)+1}(kC)}{L_{n/2}(kC)}>1\quad\mathrm{for} \quad k<0, \quad kC \geq -1.
\end{equation}
For this purpose, we compare sums of two consequent terms (with numbers $2p$ and $2p+1$) in the Taylor series representations
\eqref{3.13} of $L_{(n/2)+1}(z)$ and $L_{n/2}(z)$, and derive that
\begin{align}\nonumber &
 |z|^{2p-1}
 \left(
 \left(-\frac{1}{(2p)^{l+1}}+\frac{|z|}{(2p+1)^{l+1}}\right)
 -\left(-\frac{1}{(2p)^{l}}+\frac{|z|}{(2p+1)^{l}}\right)
 \right)
 \\ & \label{op4}
 ={|z|^{2p-1}}\left(\frac{2p-1}{(2p)^{l+1}}-\frac{2p|z|}{(2p+1)^{l+1}}\right).
\end{align}
 Furthermore, it is immediate that
\begin{align}\label{op5} &
 {|z|^{2p-1}}\left(\frac{2p-1}{(2p)^{l+1}}-\frac{2p|z|}{(2p+1)^{l+1}}\right)=
 \frac{|z|^{2p-1}}{(2p)^l}\left(\frac{2p-1}{2p}-|z|\left(\frac{2p}{2p+1}\right)^{l+1}\right).
\end{align}
Finally, it is easy to see that, for $-1\leq z<0$, $l=n/2$, $n \geq 2$, we have
\begin{align}\label{op6} &
 \frac{2p-1}{2p}-|z|\left(\frac{2p}{2p+1}\right)^{l+1}\geq  \frac{2p-1}{2p}-\left(\frac{2p}{2p+1}\right)^{2}>0,
\end{align}
and relations \eqref{op4}--\eqref{op6} imply \eqref{3.25'}. It remains to prove \eqref{3.25'} for the case $n=1$.
We easily calculate directly that
\begin{align}\label{opp1} &
 \frac{2p-1}{2p}-\frac{8}{10}\left(\frac{2p}{2p+1}\right)^{3/2}>0,
\end{align}
which yields
\begin{align}\label{opp2} &
 \frac{2p-1}{2p}-|z|\left(\frac{2p}{2p+1}\right)^{3/2}>0 \quad {\mathrm{for}} \quad -0.8 \leq z <0.
\end{align}
Formulas \eqref{op4}, \eqref{op5} and \eqref{opp2} show that
\begin{align}\label{opp3} &
L_{3/2}(kC)>L_{1/2}(kC) \quad {\mathrm{for}} \quad -0.8 \leq kC <0.
\end{align}
We use connections between Lerch zeta functions and Riemann zeta functions and
take into account the estimate \cite[sequence A078434]{oeis} in order
to calculate that up to the first two decimal places we have
\begin{align}\label{opp4} &
L_{3/2}(-1) =0.76 \, .
\end{align}
Taking into account that (for $-1\leq z<0$) the series in \eqref{3.13} is an alternating series
satisfying the Leibniz criterion, we obtain the inequalities
\begin{align}\label{opp5} &
0.65<L_{1/2}(-0.8)<0.6589 \, .
\end{align}
In view of Proposition \ref{Proposition 3.1.}, the functions $L_{1/2}$ and $L_{3/2}$ increase
monotonically on the interval $[-1, \,-0.8]$. Hence, relations \eqref{opp4} and \eqref{opp5} imply that
\begin{align}\label{opp6} &
L_{3/2}(kC)>L_{1/2}(kC) \quad {\mathrm{for}} \quad -1 \leq kC \leq -0.8 \, .
\end{align}
Inequalities \eqref{opp3} and \eqref{opp6} prove that \eqref{3.25'} holds also for $n=1$.
Thus, it is proved that \eqref{3.25'} is valid for all $n\geq 1$.
Inequalities \eqref{3.24} follow directly from \eqref{3.21}, \eqref{3.22}, \eqref{3.25} and \eqref{3.25'}.
\end{proof}
\begin{Rk}  The proof of formula \eqref{3.25'},  for the case that  $n=1$,
shows  that Conjecture 6.1 from \cite{LASPLet} is valid.
\end{Rk}
Formulas \eqref{3.12} and \eqref{3.13'} yield
 \begin{equation} \label{op7}
 C_0=(2\pi)^{-n/2}\widetilde \rho /V_{\Omega}.
 \end{equation}
\begin{La}\label{La2.10}
The following inequalities are valid:
 \begin{align} \label{3.30}&
 C_k>C_{0} \quad {\mathrm{for}} \quad k<0; \qquad  C_k<2C_{0} \quad {\mathrm{for}} \,\, -1<2kC_0<0;
 \\ \label{3.30'}& C_k<C_{0} \quad {\mathrm{for}} \quad 0<kC_0< 1.
 \end{align}
 \end{La}
 \begin{proof}
 Let $k<0$. Then, according to Proposition \ref{Proposition 3.1.}, we have
\begin{equation} \label{op8}
 L_{n/2}(kC_0)<L_{n/2}(0)=1.
 \end{equation}
 Using relations \eqref{3.12}, \eqref{op7} and \eqref{op8}, we obtain
 \begin{equation} \label{op9}
 C_kL_{n/2}(kC_k)=C_0>C_0L_{n/2}(kC_0).
 \end{equation}
 Hence, Proposition \ref{Proposition 3.2} implies that the first inequality in \eqref{3.30} holds.

 Next, let $-1<2kC_0<0$. Again using Proposition \ref{Proposition 3.1.}, we see that
  \begin{equation} \label{op10}
2 L_{n/2}(2kC_0)>2L_{n/2}(-1).
 \end{equation}
Taking into account \eqref{3.13} and \eqref{3.13'}, from \eqref{op10} we derive
 \begin{equation} \label{op11}
2 L_{n/2}(2kC_0)>L_{n/2}(0)=1.
 \end{equation}
 Thus, we obtain
 \begin{equation} \label{op12}
 2 C_0 L_{n/2}(2kC_0)>C_0=(2\pi)^{-n/2}\widetilde \rho /V_{\Omega}.
 \end{equation}
 In view of  \eqref{3.12} and \eqref{op12},  we have $2 C_0 L_{n/2}(2kC_0)>C_k L_{n/2}(kC_k)$.
 Then, as in the proof of the first inequality in \eqref{3.30}, we apply
 Proposition \ref{Proposition 3.2} and see
  that  the second inequality in \eqref{3.30} holds.

 Finally, let $0<kC_0< 1$. Since $L_{n/2}(z)$ is increasing (see Proposition~\ref{Proposition 3.1.}),
 we have
 \begin{equation} \label{op8'}
 L_{n/2}(kC_0)>L_{n/2}(0)=1,
 \end{equation}
 and, moreover, $C L_{n/2}(kC)$ also increases strictly monotonically. According to \eqref{3.12}, \eqref{op7}
 and \eqref{op8'}, the relations
 \begin{equation} \label{op13}
 C_0L_{n/2}(kC_0)>C_0=C_kL_{n/2}(kC_k)
 \end{equation}
 hold. Therefore, since $C L_{n/2}(kC)$ is monotonic, we see that  $C_k<C_0$.
 \end{proof}
 It follows from Lemma \ref{La2.10} that $C_k$ is bounded  in the  neighborhood of $k~=~0$.
 The behavior of $E_q$, $F_q=F(M_k)$ and the entropy $S_q=S(M_k)$ in the punctured neighborhood
 of  $k=0$ is given in the next proposition.
 \begin{Pn}\label{Pn3.4}
For $k\to 0$, we have  the following asymptotic  relations:
 \begin{align} \label{3.31}&
 E_{q}-E_{c}=-\frac{n{\widetilde \rho}kC_{0}}{{4}(2^{n/2})}+O(k^{2}) \quad (E_c=E(M_0)),
\\ & \label{3.32}
 S_{q}-S_{c}=-\frac{(n-2){\widetilde \rho}kC_{0}}{{4}(2^{n/2})}+O(k^{2}) \quad (S_c=S(M_0)),
\\ &  \label{3.33}
 F_{q}-F_{c}=\frac{{\widetilde \rho}kC_{0}}{{2}(2^{n/2})}+O(k^{2}) \quad (F_c=F(M_0)).
\end{align}
\end{Pn}
\begin{proof}
Step 1.
In order to calculate the entropy $S(M_{k})$ we recall definitions \eqref{4} and \eqref{16} of $S$ and $M_k$, respectively,
 and use the equalities
  \begin{equation} \label{n7}
M_{k}=g/(1-k g), \quad 1+k M_{k}=(1-k g)^{-1}, \quad g:=C_k\mathrm{e}^{-|v|^2/2},
\end{equation}
which simplify the expressions  $\Phi(M_k)$ for $\Phi$ given by \eqref{5'} and \eqref{5}:
\begin{align} \label{n5}&
\Phi(M_k)=M_{k}(1+\log g)+(1/k)\log(1-k g) \quad \mathrm{for} \quad k\not=0, \\
\label{n8}&
\Phi(M_0)=g\log g \quad \mathrm{for}  \quad M_0=g.
\end{align}
Recall definitions \eqref{7} and \eqref{9} of $\rho$ and $E$ and recall that for $f=M_k$ we have
$\rho(t) \equiv {\mathrm{const}}=\widetilde \rho$.
Substituting $\log g=\log C_k -(1/2)|v|^2$ into \eqref{n5}  and then substituting  \eqref{n5}
 into
\eqref{4}, we obtain
 \begin{align} \label{n9}&
 S(M_{k})=E_{q}-(1+\log{C_k})\widetilde \rho -\frac{1}{k}V_{\Omega}\int_{{\mathbb R}^n}
 \log(1-k g)dv \quad  \mathrm{for} \quad k\not=0.
 \end{align}
Substituting $\log g=\log C_k -(1/2)|v|^2$ into \eqref{n8} and then substituting   \eqref{n8}
 into
\eqref{4}, we obtain
\begin{align}
  \label{3.26}& S(M_{0})=S_c= E_{c}-{\widetilde \rho}\log{C_{0}} ,
\end{align}
where $E_c=E(M_0)$ and $C_0$ is the value of $C_k$ for the case that $k=0$
(recall Notation \ref{NnC}). Taking into account the definition \eqref{9} of energy and
using spherical coordinates and integration by parts,  we rewrite \eqref{n9}:
 \begin{equation} \label{3.27}
 S(M_{k})=E_{q}-(1+\log{C_k})\widetilde \rho
+\frac{2}{n}{E_{q}} =\left(1+\frac{2}{n}\right){E_{q}}-(1+\log{C_k})\widetilde \rho \quad  \mathrm{for} \quad k\not=0.
\end{equation}

According to \eqref{3.22}, we have $\widetilde \rho=2E_c/n$. Therefore, for $k \not=0$ formulas \eqref{3.26} and \eqref{3.27} imply that
 \begin{equation} \label{3.28}
 S(M_k)-S(M_0)=S_{q}-S_{c}=\frac{n+2}{n}(E_{q}-E_{c})-{\widetilde \rho}\log({C_k}/{C_{0}}),
\end{equation}
Hence, taking into account  \eqref{12} and \eqref{3.28} we derive
 \begin{equation} \label{3.29}  F_{q}-F_{c}=\frac{2}{n}(E_{q}-E_{c})-{\widetilde \rho}\log({C_k}/{C_{0}}).
\end{equation}

Step 2. The equalities in  \eqref{op9} and \eqref{op13} yield
 \begin{equation} \label{op15}
\frac{kC_0}{L_{n/2}(kC_k)}=kC_k.
\end{equation}
In view of formula \eqref{op7} and Lemma \ref{La2.10}, we see that the values
$|kC_k|$ and $\sup_{|z|\leq \varepsilon}\left|\frac{d}{dz}\left(\frac{kC_0}{L_{n/2}(z)}\right)\right|$
are small
for small values of $|k|$.
Thus, we can apply the iteration method to the equation $z=\frac{kC_0}{L_{n/2}(z)}$
in order to derive
 \begin{equation} \label{op16}
C_k=C_0+O(k), \quad k \to 0.
\end{equation}
Taking into account the series expansion in \eqref{3.13} and formulas \eqref{op15} and \eqref{op16},
we obtain
 \begin{equation} \label{op17}
C_k/C_0=1/{L_{n/2}(kC_k)}=1-kC_k/2^{n/2}+O(k^2).
\end{equation}
Furthermore, from \eqref{op16} and \eqref{op17} we have
 \begin{equation} \label{op18}
\log\big(C_k/C_0\big)=-kC_0/2^{n/2}+O(k^2).
\end{equation}
Using formulas \eqref{3.21}, \eqref{3.22}, \eqref{op16} and the series expansion in \eqref{3.13},
we see that  \eqref{3.31} holds. According to \eqref{3.31} and \eqref{op18} we may rewrite \eqref{3.28}
in the form \eqref{3.32}. Finally, in view of \eqref{3.31} and \eqref{op18}, we rewrite \eqref{3.29}
in the form \eqref{3.33}.
\end{proof}

\begin{Cy}
The following inequalities hold for small values of $k:$
 \begin{equation} \label{3.34}
E_{q,B}<E_{c}<E_{q,F}; \quad S_{q,B}<S_{c}<S_{q,F} \quad {\mathrm{for}} \quad n>2; \quad F_{q,F}<F_{c}<F_{q,B}.
 \end{equation}
 \end{Cy}

\section{General-type solutions }\label{sec4}
\subsection{Dissipative and conservative solutions }
\paragraph{1.}
In this section we study  general solutions $f(t,x,v)$ (satisfying \eqref{3}) of the Focker-Planck equation \eqref{1}.
The total energy flux through the surface $\partial\Omega$ per unit time is given by the equalities
 \begin{align}\nonumber
A(f,\Omega): &= \int_{\Omega}\int_{{\mathbb R}^n}(|v|^{2}/2)v{\cdot}{\bigtriangledown_{x}}f(t,x,v)d{v}dx
\\  \label{4.1}  &= \int_{\partial\Omega}\int_{{\mathbb R}^n}(|v|^{2}/2)\big(v{\cdot}n(y)\big)f(t,y,v)d{v}d\sigma,
\end{align}
where $\partial\Omega$ is the boundary of the $\Omega$, and the integral $\int_{\partial\Omega}g d\sigma$ is the surface integral with $n(y)$ being the outward
unit normal to that surface, $y{\in}\partial\Omega.$ The second equality in \eqref{4.1} is immediate from Gauss-Ostrogradsky divergence formula.

 The total density flux through the surface $\partial\Omega$ per unit time has the form
\begin{align}&
  \label{4.2}  B(f,\Omega):=\int_{\Omega}\int_{{\mathbb R}^n}v{\cdot}{\bigtriangledown_{x}}f(t,x,v)d{v}dx= \int_{\partial\Omega}\int_{{\mathbb R}^n}\big(v{\cdot}n(y)\big)f(t,y,v)d{v}d\sigma.
 \end{align}

\begin{Dn} By  $D(\Omega)$, we denote the class of
functions $f(t,x,v)$ satisfying the Fokker-Planck equation \eqref{1}, inequalities \eqref{3}
and the  condition $A(f,\Omega){\geq}0$ for all $t$ $($i.e., the class of the dissipative solutions $f$$)$.
\end{Dn}

\begin{Dn} By  $C(\Omega)$, we denote the class of
functions $f(t,x,v)$ satisfying the Fokker-Planck equation \eqref{1}, inequalities \eqref{3}
and the  condition $A(f,\Omega){=}0$ for all $t$ $($i.e., the class of the conservative solutions $f$$)$.
\end{Dn}

It is obvious that $C(\Omega){\subset}D(\Omega).$

\begin{Pn}
Let  $f(t,x,v)$ satisfy \eqref{1} and \eqref{3}, and assume that for all $y\in \partial\Omega$ the equality
 \begin{equation} \label{4.3}
 f(t,y,v)=f(t,y,-v)
 \end{equation}
holds.  Then, $f(t,x,v)\in C(\Omega)$.
 \end{Pn}

\begin{proof}
Taking into account  \eqref{4.3}, we derive
 \begin{equation} \label{4.4}
 \int_{{\mathbb R}^n}(|v|^{2}/2)\big(v{\cdot}n(y)\big)f(t,x,v){v}d{v}=0.
 \end{equation}
It is immediate from \eqref{4.1} and \eqref{4.4} that  $A(f,\Omega){=}0$.
\end{proof}

\begin{Cy}
Maxwellians  $M_{k}$ given by
\eqref{16} are conservative, that is, $M_k\in C(\Omega)$.
\end{Cy}

The \textbf{bounce-back condition} \eqref{4.3} means that particles arriving with a certain velocity
to the boundary $\partial\Omega$ will bounce back with an opposite velocity (see \cite[p.16]{Vi1}).

\subsection{Boundedness}
Introducing the function
\begin{equation}
\label{5.1}
s(r)=r\log{r}-\frac{1}{k}(1+{k}r) \log(1+{k}r),\end{equation}
we see that
\begin{equation}\label{5.2} s^{\prime}(r)=\log\big(r/(1+{k}r)\big),\quad
s^{\prime\prime}(r)=\big(r(1+{k}r)\big)^{-1},\end{equation}
and obtain the proposition below.
\begin{Pn}\label{Pn5.1}
The function $s(r)$ is convex on the semi-axis $r \geq 0$ for the case that $k>0$ and on the interval $0{\leq}r<1/|k|$ for the case that $k<0$.
\end{Pn}
We consider functions $g$ such that
\begin{equation}\label{5.5}g\geq 0, \quad 1+kg>0, \quad \frac{1}{2}\int_{{\mathbb R}^n}|v|^{2}g(v)dv<\infty.
\end{equation}
For $k<1$ and $Z$ given by
\begin{equation}\label{5.3}Z=1/(\mathrm{e}^{|v|^{2}/2}-k),
\end{equation}
the convexity of $s$ implies that
\begin{equation}\label{5.4}s(g)-s(Z){\geq}s^{\prime}(Z)(g-Z).\end{equation}
Using \eqref{5.2}, \eqref{5.3}, \eqref{5.4}  and the equality
\begin{equation}\label{5.6}Z/(1+{k}Z)=\mathrm{e}^{-|v|^{2}/2},\end{equation}
we easily derive
\begin{equation}\label{5.7}-s(g){\leq}\frac{|v|^{2}}{2}g-s(Z)-\frac{|v|^{2}}{2}Z.
\end{equation}
Taking into account the equality
\begin{equation}\label{5.8}1+{k}Z=1/(1-{k}\mathrm{e}^{-|v|^{2}/2}), \end{equation}
we rewrite \eqref{5.7} in the form
\begin{equation}\label{5.9}-s(g){\leq}\frac{|v|^{2}}{2}g
-\frac{1}{k}\log(1-{k}e^{-|v|^{2}/2}).
\end{equation}
The  proposition below follows from \eqref{5.9} and the inequality $\log(1+a){\leq}a$ for $a{\geq}0$.
\begin{Pn}\label{Pr5.2} Assume that $k<0$ and $g$ satisfies  \eqref{5.5}.
Then we have
\begin{equation}
\label{5.11}
-{|v|^{2}}g/2-s(g){\leq}
\mathrm{e}^{-|v|^{2}/2}.
\end{equation}
\end{Pn}
\begin{Rk}\label{Rk5.1.} For the case
$k=-1$, the
inequality \eqref{5.11} was  derived  in \cite{CPR}.
\end{Rk}
Recall definitions  \eqref{4}, \eqref{9} and \eqref{12} of $S$, $E$ and $F$, respectively,
and note that $S$ is expressed via $\Phi$.
Compare the expression \eqref{5} for $\Phi$
 with the expression for $s$.
Thus, we see that  formula \eqref{5.11} yields the following corollary.

\begin{Cy}\label{Cy5.1} Let the conditions of Proposition \ref{Pr5.2} be fulfilled.
Then there exists a positive constant $C$ such that
\begin{equation}\label{5.12}F(g){\leq}C.\end{equation}
\end{Cy}
In a similar way, from \eqref{5.9} and series representation of $\log(1+a)$ we obtain the next corollary.
\begin{Cy}\label{Cy5.2} Let $0<k<1$ and assume that $g$ satisfies  \eqref{5.5}.
Then there exists a positive constant $C$ such that
\begin{equation}\label{5.13} F(g){\leq}C.\end{equation}
\end{Cy}
If either conditions of Corollary \ref{Cy5.1} or Corollary \ref{Cy5.2} hold, we put
\begin{equation}\label{5.13'} \widetilde C=\inf F(g).\end{equation}
Passing to the limit and using Corollaries \ref{CyB}, \ref{Cy 2.5} and \ref{CyF}, we could show
(under rather general conditions) the existence of the global Maxwellian
$M_{\widetilde C}$ such that
\begin{align}&
\label{7.9} F(M_{\widetilde C})=\widetilde C.
\end{align}
Later we assume that $M_{\widetilde C}$ exists.
\section{Lyapunov functional}
\label{Lya}
The Lyapunov functional for equation \eqref{1} has the form
 \begin{equation}\label{7.1} \widetilde G(f)=\widetilde C-F(f),\end{equation}
 where
 $F(f)$ is defined by \eqref{12}. Recall that, for the generalization $G(f)$ (given by \eqref{18}) of the Kullback-Leibler distance,
 we always assumed that $\rho(f)$ is fixed at the moment $t$ (i.e., $\rho(f,t) \equiv \widetilde \rho$).
 We do not assume this in the present section. In other words, we consider $\widetilde G$ on a wider set
 of solutions.
 Clearly, under conditions of Corollaries \ref{Cy5.1} or \ref{Cy5.2}, the inequality
\begin{equation}
\label{7.2} \widetilde G(f){\geq}0
\end{equation}
 holds. Thus, if a stationary solution $M_{\widetilde C}$ of \eqref{1} satisfies the equality \eqref{7.9},
the value $\widetilde G(f)=F(M_{\widetilde C})-F(f)$ at $t$ may be considered as a distance  between $M_{\widetilde C}$ and $f$ at the time $t$.

In this section we substitute condition \eqref{3} by a stronger condition with the strict inequalities
\begin{equation}\label{3'}
f(t,x,v) > 0, \quad 1+kf(t,x,v) > 0 \quad (k \in {\mathbb R}).
\end{equation}
 Hence, we can consider equation \eqref{1} in the form
  \begin{equation}\label{7.3}
 \frac{\partial{f}}{\partial{t}}=-v \cdot \bigtriangledown_{x}{f}+
 {\mathrm{div}}_{v}\left(f(1+kf){\bigtriangledown}_{v}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)\right).\end{equation}
 We multiply \eqref{7.3} { by $\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)$ and integrate over ${\mathbb R}^n$ (applying also integration by parts
 with respect to the  variables $v_i$)
 and $\Omega.$
Under natural assumptions on the decay of $f$ and $\partial f/\partial v_i$ at infinity (so that the values at $\pm \infty$ disappear in the formulas
for the integration by parts),  we obtain
\begin{align}& \int_{\Omega}\int_{{\mathbb R}^n}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)\frac{\partial{f}}{\partial{t}}dvdx {\nonumber} \\
 &=- \int_{\Omega}\int_{{\mathbb R}^n}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)v \cdot \bigtriangledown_{x}{f}dvdx \label{7.4}\\ & \quad
 {\nonumber} -\int_{\Omega}\int_{{\mathbb R}^n}f(1+kf)\left|{\bigtriangledown}_{v}\left(\log\frac{f}{1+kf}+1+v^{2}/2\right)\right|^{2}dvdx.
 \end{align}
 Let us introduce the total flow of the entropy across the boundary $\Omega$:
 \begin{equation}\label{7.5}U(f, \Omega)=-\int_{\Omega}\int_{{\mathbb R}^n}(\log\frac{f}{1+kf}+1)v \cdot \bigtriangledown_{x}{f}dvdx=
-\int_{\Omega}\int_{{\mathbb R}^n}v\cdot \bigtriangledown_{x}\Phi(f)dvdx,
\end{equation}
where $\Phi$ is defined by \eqref{5'} and \eqref{5}.
 Using Gauss-Ostrogradsky formula, we also have
\begin{equation}\label{7.5'}U(f, \Omega)=-
\int_{\partial\Omega} \int_{{\mathbb R}^n}\big(v{\cdot}n(y)\big)\Phi(f)dvd\sigma .\end{equation}
It follows from \eqref{7.4} and \eqref{7.5}  that
\begin{equation}\label{7.6}\frac{d\widetilde G}{dt}{\leq}U(f,\Omega)-A(f,\Omega),\end{equation}
where the function $A(f,\Omega)$ is defined by relation \eqref{4.1}.
Indeed, relations \eqref{4.1} and \eqref{7.5}  yield the equality
\begin{align}& U(f,\Omega)-A(f,\Omega)=- \int_{\Omega}\int_{{\mathbb R}^n}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)v \cdot \bigtriangledown_{x}{f}dvdx,
\label{7.6'}
\end{align}
whereas formulas \eqref{3}, \eqref{7.4} and \eqref{7.6'} imply that
\begin{align}&  \label{7.6''} U(f,\Omega)-A(f,\Omega)
\geq \int_{\Omega}\int_{{\mathbb R}^n}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)\frac{\partial{f}}{\partial{t}}dvdx.
\end{align}
Furthermore, according to \eqref{4}, \eqref{9}, \eqref{12} and \eqref{7.1}, we have
\begin{equation}\nonumber
\frac{\partial\widetilde G}{\partial t}=\frac{\partial E}{\partial t}-\frac{\partial S}{\partial t}=\int_{\Omega}\int_{{\mathbb R}^n}\left(\log\frac{f}{1+kf}+1+|v|^{2}/2\right)\frac{\partial{f}}{\partial{t}}dvdx,
\end{equation}
and so \eqref{7.6} is immediate from \eqref{7.6''}.

Using inequality \eqref{7.6}, we derive the following assertion.
\begin{Tm}\label{Tm7.1}
Assume that $f{\in}C_{0}^{1}$ is a dissipative solution of  \eqref{1}, which satisfies
\eqref{3},  and  that the inequality
\begin{equation}\label{7.8} U(f,\Omega){\leq}0 \end{equation}
holds.  Then the inequality
 $( \partial{\widetilde G}/{\partial t}){\leq}0$ is valid.
\end{Tm}
\begin{Cy}
Assume that the  conditions of Theorem \ref{Tm7.1} are fulfilled and
$\widetilde G(f,t_0)<\delta$. Then, the inequality $\widetilde G(f,t)<\delta$
holds for all $ t>t_{0}$.
\end{Cy}
Thus, the distance $\widetilde G$ between a Maxwellian $M_{\widetilde C}$ satisfying \eqref{7.9} and $f$
satisfying \eqref{7.8} decreases. Taking into account the definition of the Lyapunov stability,
we proved the following important result.
\begin{Tm}\label{Tm7.2} The stationary solution $M_{\widetilde C}$ is locally  stable $($i.e., Lyapunov stable$)$
in the class of the dissipative solutions $f$ satisfying inequalities \eqref{3} and
\eqref{7.8}. \end{Tm}
We note that $M_{\widetilde C}$ does not depend on $x$, and so definitions \eqref{4.1} and \eqref{7.5} imply
that $A(M_{\widetilde C}, \, \Omega)=U(M_{\widetilde C}, \, \Omega)=0$.

\begin{Rk}The earlier results on the local stability for the Fokker-Planck equation
$($see \cite{CPR}$)$ were obtained  only for   the spatially homogeneous case. Theorems \ref{Tm7.1} and \ref{Tm7.2} have no analogs also
in the case of
the Boltzmann equation.\end{Rk}

{\bf Concluding remark.}
Following \cite{LAS1, LAS2, LAS3, LASPLet, LAS4}, we compare classical and quantum results, that is,
determinate and probabilistic cases.
As in the case of the Boltzmann equation \cite{LASPLet, LAS4},
we use a special extremal principle, which is based on the ideas of  game theory. (We note that extremal principles  remain central in modern physics.)
The "players" in the game described by the Fokker-Planck equation
are the total energy $E$ and entropy $S$, the "gain" in the game is the functional $F$ and the strategy in the game  is determinate in the classical case
and probabilistic in the quantum case. It is of interest that the inequalities $E_{q,B}<E_{c}<E_{q,F}$, $S_{q,B}<S_{c}<S_{q,F}$ (for $n>2$), and $F_{q,F}<F_{c}<F_{q,B}$
(see \eqref{3.34}) hold for small $k$ in our game.

\noindent{\bf Acknowledgments.}
 {The research of A.L. Sakhnovich was supported by the
Austrian Science Fund (FWF) under Grant  No. P24301.}