\section{Introduction}

This paper is devoted to positional determinacy in turn-based infinite-duration games over \emph{finite} arenas. An arena is a finite directed graph whose edges are colored into elements of some finite set of colors $C$ and whose nodes are partitioned between two players called Alice and Bob. They play by traveling over the nodes of the arena. In each turn, one of the players chooses an edge from the current node, and the players move towards the endpoint of this edge. Whether it is an Alice's or a Bob's turn to choose depends on whether the current node is an Alice's node or a Bob's node. This continues for infinitely many turns. As a result, the players obtain an infinite word over $C$ (by concatenating colors of edges that appear in the play). A \emph{winning condition} $W$, which is a set of infinite words over $C$, defines the aims of the players. Alice wants to obtain an infinite word which belongs to $W$, while Bob wants it to be outside $W$.

A vast amount of literature in this area is devoted to \emph{positional strategies}. A strategy of Alice or Bob is positional if it never makes two different moves from the same node. Implementing such strategies is easy because we only have to specify one edge for each node of the corresponding player. This makes these strategies relevant for such areas as controller synthesis~\cite{bloem2018graph}, where an implementation of a controller can be seen as its strategy against an environment.

Correspondingly, of great interest are winning conditions for which positional strategies are always sufficient to play optimally (for one of the players or even for both of them). This area has the following terminology.
A winning condition $W$ is \emph{half-positional} if for every arena the following holds: either Alice has a positional winning strategy w.r.t.~$W$ or Bob has a winning strategy w.r.t.~$W$ (not necessarily positional). A winning condition $W$ is \emph{bi-positional} if additionally the same requirement as for Alice holds for Bob (or, in other words, if both $W$ and its complement are half-positional).

The most famous example of a bi-positional winning condition is the parity condition. It has an abundance of applications in logic, including decidability of logical theories~\cite{zielonka1998infinite} and modal $\mu$-calculus~\cite{gradel2002automata}.

Gimbert and Zielonka performed a general study of bi-positional winning conditions in a series of two papers~\cite{gimbert2004can,gimbert2005games}. In~\cite{gimbert2004can}, they gave a powerful sufficient condition for bi-positionality. It is suitable for almost all known bi-positional winning conditions. Moreover, in~\cite{gimbert2005games} they gave a sufficient and necessary condition for bi-positionality. Unfortunately, it is far more complex, and hence is less convenient for applications. Nevertheless, it has a corollary called \emph{1-to-2-player lifting}, which is of great interest in practice. It states that as long as Alice can play positionally without the presence of Bob and Bob can play positionally with the presence of Alice, they both can play positionally with the presence of each other.

In turn, no sufficient and necessary condition for half-positionality is known. This problem is of great importance in controller synthesis because half-positional winning conditions are more relevant there than bi-positional ones  (we are interested in positional strategies only for the player representing a controller). One example of half-positional winning conditions that are not bi-positional are Rabin conditions~\cite{klarlund1991rabin}. These are finite unions of parity conditions.

 In~\cite{kopczynski2006half}, Kopczy\'{n}ski undertook a study of half-positional winning conditions that are \emph{prefix-independent} (which means closed under removing or adding finite prefixes). He defined several properties of prefix-independent winning  conditions that imply half-positionality (concavity, monotonicity, etc.). He noticed that all these properties are preserved under finite unions.  This led him to the following conjecture.

\begin{conjecture}[Kopczy\'{n}ski, \cite{kopczynski2006half}]
\label{conj}
Prefix-independent half-positional winning conditions are closed under finite unions.
\end{conjecture}
He refuted this conjecture for uncountable unions. However, it remained open even for countable unions. In this paper, we refute Conjecture \ref{conj}. Moreover, we show that there are two \emph{bi-positional} prefix-independent winning conditions whose union is not \emph{half-positional}.

To this end, we first introduce a new class of bi-positional winning conditions called \emph{energy conditions over totally ordered groups}, or ETOG conditions for short. They are defined as follows (see more details in Section \ref{sec:def}). We consider elements of some totally ordered group (we stress that it should be bi-ordered) as colors of edges. Given an infinite sequence of these elements, we arrange them into a formal series. Alice wants the sequence of its partial sums to have an infinite decreasing subsequence. This generalizes canonical energy\footnote{In case of $\mathbb{Z}$, a color of an edge is interpreted as an amount of energy needed to pass through this edge (and negative edges are edges where one can recharge). In this sense, Bob wins if there is a finite amount of initial energy allowing him to never run out of it.} conditions~\cite{chakrabarti2003resource} that are defined over $\mathbb{Z}$ with the standard ordering.

In Section \ref{sec:pos}, we establish bi-positionality of energy conditions over totally ordered group using a sufficient condition of Gimbert and Zielonka. Next, we refute Conjecture \ref{conj} in Section \ref{sec:refute}. A key factor allowing us to do this is that free groups can be totally ordered. We construct two energy conditions over a free group with 2 generators whose union is not half-positional. We also observe in Section \ref{sec:refute} that energy conditions over free groups are non-permuting, and that they can be used to refute 1-to-2-player lifting for half-positionality.

We believe that the class of energy conditions over totally ordered groups is interesting on its own. Namely, we find this class suitable for the following conjecture.
\begin{conjecture}
\label{our_conj}
Every bi-positional prefix-independent winning condition coincides on periodic sequences with some energy condition over a totally ordered group.
\end{conjecture}
We cannot expect it to hold for all sequences, but periodic once are sufficient, say, for algorithmic applications.
If our conjecture is true, it gives an explicit description of the class of bi-positional prefix-independent winning condition. This would in line with an explicit description of the class of \emph{continuous} bi-positional payoffs from~\cite{kozachinskiy:LIPIcs.CONCUR.2021.10}.
We discuss our conjecture in more detail in Section \ref{sec:conj}, where we reduce it to a problem about free groups.

\medskip

\textbf{Open problems and related works.} Conjecture \ref{conj} is still open for \emph{infinite arenas} (that are not considered in this paper at all). Our approach does not work because energy conditions over totally ordered groups are not in general half-positional in infinite arenas. In a recent breakthrough, Ohlmann~\cite{pierrephd} obtained a characterization of  half-positionality in infinite arenas (with the use of so-called well-monotonic universal graphs). As far as we know, this did not yet resolve Conjecture \ref{conj} over infinite arenas (but there is a hope that it will).

Among winning conditions, of particular interest are $\omega$-regular ones, due to their connections to automata theory and logic. We do not yet have a good understanding which of them are half-positional, although this property is decidable~\cite{kopczynski2007omega}, given a deterministic parity automaton recognizing a condition. As far as we know, Conjecture \ref{conj} is open in restriction to $\omega$-regular conditions. In a recent preprint~\cite{https://doi.org/10.48550/arxiv.2205.01365}, Bouyer et al.~(among other results) prove Conjecture \ref{conj} for  $\omega$-regular conditions that can be recognized by deterministic B\"{u}chi automata.

More sufficient conditions for half-positionality over finite arenas can be found in~\cite{bianco2011exploring}.
\section{Preliminaries}

If $C$ is a set, we denote by $C^*$ (resp., by $C^\omega$) the set of all finite (resp., infinite) words over $C$. For $x\in C^*$, by $|x|$ we denote the length of $x$. Additionally, by $C^+$ we denote the set of all finite non-empty words over $C$. If $x\in C^+$, then by $x^\omega$ we denote an infinite word obtained by repeating $x$ infinitely many times. The free group over $C$ is denoted by $F_C$.

\medskip

An arena $\mathcal{A}$ over a non-empty finite set (of colors) $C$ is a tuple $\langle V_A, V_B, E\rangle$, where $V_A$ and $V_B$ are disjoint finite sets and $E\subseteq (V_A\cup V_E) \times C \times (V_A\cup V_B)$ is such that for every $s\in V_A\cup V_B$ there exist $c\in C$ and $t\in V_A\cup V_B$ for which $(s, c, t)\in E$. Elements of $V_A$ are called Alice's nodes, and elements of $V_B$ are called Bob's nodes. Elements of $E$ are called edges of $\mathcal{A}$. An edge $e = (s, c, t)\in E$ is represented as a $c$-colored arrow from $s$ to $t$. We use the notation $\mathsf{source}((s, c, t)) = s, \mathsf{col}((s, c, t)) = c$ and $\mathsf{target}((s, c, t)) = t$. Our definition guaranties that every node $v\in V_A\cup V_B$ has an out-going edge, that is, an edge $e$ such that $\mathsf{source}(e) = v$.

An infinite-duration game over $\mathcal{A}$ from a node $s\in V_A\cup V_B$ is played as follows. At the beginning, one of the players chooses an edge $e_1\in E$ with $\mathsf{source}(e_1) = s$. Namely, if $s\in V_A$, then Alice chooses $e_1$, and if $s\in V_B$, then Bob chooses $e_1$. More generally, in the first $n$ turns players choose $n$ edges $e_1, e_2, \ldots, e_n\in E$, one edge per turn. These edges always form a \emph{path} in $\mathcal{A}$, that is, we have $\mathsf{target}(e_1) = \mathsf{source}(e_2), \ldots, \mathsf{target}(e_{n-1}) = \mathsf{source}(e_n)$. Then the $(n + 1)$st turn is played as follows. Players consider the endpoint node of the current path, which is $\mathsf{target}(e_n)$. One of the players chooses an edge $e_{n + 1}$ with $\mathsf{source}(e_{n + 1}) = \mathsf{target}(e_n)$. Namely, if $\mathsf{target}(e_n) \in V_A$, then Alice chooses $e_{n + 1}$, and if $\mathsf{target}(e_n)\in V_B$, then Bob chooses $e_{n + 1}$. After infinitely many turns, players get an infinite sequence of edges $p = (e_1, e_2, e_3, \ldots)$ called a play (it forms an infinite path in $\mathcal{A}$).

A winning condition over a set of colors $C$ is a subset $W\subseteq C^\omega$. A strategy of Alice is winning from $s\in V_A\cup V_B$ w.r.t.~$W$ if any play $p = (e_1, e_2, e_3, \ldots)$ with this strategy in the infinite-duration game over $\mathcal{A}$ from $s$ is such that its sequence of colors $\mathsf{col}(e_1)\mathsf{col}(e_2)\mathsf{col}(e_3)\ldots$ belongs to $W$. Similarly, a strategy of Bob is winning from $s\in V_A\cup V_B$ w.r.t.~$W$ if any play $p = (e_1, e_2, e_3, \ldots)$ with this strategy in the infinite-duration game over $\mathcal{A}$ from $s$ is such that $\mathsf{col}(e_1)\mathsf{col}(e_2)\mathsf{col}(e_3)\ldots\notin W$.

A positional strategy of Alice is a function $\sigma\colon V_A\to E$ such that $\mathsf{source}(\sigma(u)) = u$ for any $u\in V_A$. It is interpreted as follows: for any $u\in V_A$, whenever Alice has to choose an edge from $u$, she chooses $\sigma(u)$. Similarly, a positional strategy of Bob is a function  $\tau\colon V_B\to E$ such that $\mathsf{source}(\tau(u)) = u$ for any $u\in V_B$. It is interpreted analogously.

 A winning condition $W\subseteq C^\omega$ is \emph{half-positional} if for every finite arena $\mathcal{A}$ over $C$ there exists a positional strategy $\sigma$ of Alice such that for every node $s$ of $\mathcal{A}$ the following holds: if $\sigma$ is not winning w.r.t.~$W$ from $s$, then Bob has a winning strategy w.r.t.~$W$ from $s$. A winning condition $W$ is \emph{bi-positional} if both $W$ and its complement $C^\omega\setminus W$ are half-positional.

A winning condition $W\subseteq C^\omega$ is \emph{prefix-independent} if for all $x\in C^*$ and $\alpha\in C^\omega$ we have $\alpha\in W\iff x\alpha \in W$.

We state the following sufficient condition for bi-positionality due to Gimbert and Zielonka.

\begin{definition}
\label{def:fairly_mixing}
Let $W\subseteq C^\omega$ be a winning condition over a finite set of colors $C$. We call $W$ \textbf{fairly mixing} if the following 3 conditions hold:
\begin{enumerate}[label=\Alph*)]
\item For every $x\in C^*$ and $\alpha, \beta\in C^\omega$ we have that
\[(x\alpha\notin W \land x\beta\in W) \implies (\alpha\notin W \land \beta\in W).\]
\item For every $S\in\{W, C^\omega\setminus W\}$, for every $x\in C^+$ and for every $\alpha\in C^\omega$ we have that
\[(x^\omega\in S, \alpha\in S)\implies (x\alpha\in S).\]
\item For every $S\in\{W, C^\omega\setminus W\}$ and for every infinite sequence $x_1, x_2, x_3, \ldots\in C^+$ it holds that:
\begin{align*}
\bigg[( x_1 x_3 x_5\ldots\in S\big) \land (x_2 x_4 x_6\ldots\in S) \land(\forall n \ge 1\,\, x_n^\omega\in S) \bigg] \implies x_1 x_2 x_3\ldots\in S.
\end{align*}
\end{enumerate}
\end{definition}

\begin{theorem}[\cite{gimbert2004can}]
\label{thm:gztheorem}
 Any fairly mixing winning condition is bi-positional.
\end{theorem}

\section{Definition of Energy Games over Totally Ordered Groups}
\label{sec:def}

Recall that a totally ordered group is a triple $(G, +, \le)$, where $(G, +)$ is a group and $\le$ is a total order on $G$ such that
\[a\le b \implies x + a + y \le x + b + y \qquad \mbox{for all } a, b, x, y\in G.\]

Consider any finite set $C$ of colors and any totally ordered group $(G, +,\le)$. By a \emph{valuation of colors} over $(G, +, \le)$ we mean any function $\mathbf{val}\colon C\to G$. It can be extended to a homomorphism $\mathbf{val}\colon C^*\to G$ by setting
\[\mathbf{val}(\mbox{empty word}) = 0, \qquad \mathbf{val}(c_1 c_2\ldots c_n) = \mathbf{val}(c_1) + \mathbf{val}(c_2) + \ldots + \mathbf{val}(c_n).\]
Additionally, for every infinite sequence of colors $c_1 c_2 c_3\ldots\in C^\omega$, we denote by $\mathbf{\overline{val}}(c_1 c_2 c_3\ldots)$ the sequences of valuations of its finite prefixes:
\[\mathbf{\overline{val}}(c_1 c_2 c_3\ldots) = \{\mathbf{val}(c_1\ldots c_n)\}_{n = 1}^\infty.\]
In other words,  $\mathbf{\overline{val}}(c_1 c_2 c_3\ldots)$ is the sequence of  partial sums of the formal series $\sum_{n = 1}^\infty \mathbf{val}(c_n)$.

\emph{An energy condition over  $(G, +, \le)$}, defined by a valuation of colors $\mathbf{val}\colon C\to G$, is the set $W\subseteq C^\omega$ of all $\alpha\in C^\omega$ such that $\mathbf{\overline{val}}(\alpha)$
has an infinite decreasing subsequence. It is immediate that any energy condition over a totally ordered group is prefix-independent.

As an illustration, we show that parity conditions fall into this definition. The parity condition over $d$ priorities is a winning condition $W_{par}^d\subseteq\{1,2,\ldots,d\}^\omega$,
\[W_{par}^d = \{ c_1 c_2 c_3\ldots\in \{1, 2, \ldots, d\}^\omega\mid \limsup_{n\to\infty} c_i \mbox{ is odd}\}.\]
Observe that $W_{par}^d$ is an energy condition over $\mathbb{Z}^d$ with the lexicographic ordering, defined by the following valuation:
\begin{align*}
    \mathbf{val}(d) &= ((-1)^d, 0, \ldots 0) \\
    \mathbf{val}(d - 1) &= (0, (-1)^{d - 1}, \ldots 0)\\
    &\vdots\\
     \mathbf{val}(1) &= (0, 0, \ldots, -1).
\end{align*}

As far as we know, the most general class of bi-positional prefix-independent winning conditions that were previously considered are \emph{priority mean payoff conditions}~\cite{gimbert2006deterministic}. They can also be defined as energy conditions over $\mathbb{Z}^d$. Moreover, to define them, it is sufficient to consider only valuations that map each color to a vector with at most 1 non-zero coordinate, as in the case of parity conditions.

\section{Bi-positionality of Energy Conditions over Totally Ordered Groups}
\label{sec:pos}
In this section, we establish

\begin{theorem}
\label{thm:pos}
Every ETOG condition  is bi-positional.
\end{theorem}

We derive it from the following technical result (which will also be useful in Section \ref{sec:conj}). If $C$ is a non-empty finite set and $W\subseteq C^\omega$, define $\mathsf{per}(W) = \{x\in C^+\mid x^\omega \in W\}$ to be the set of periods of periodic words from $W$.

\begin{proposition}
\label{prop:technical}
Let $C$ be a non-empty finite set. Consider any set $P\subseteq C^+$ such that both $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts. Define a winning condition $W_P\subseteq C^\omega$ as follows:
\[W_P = \{xy_1 y_2 y_3\ldots \mid x\in C^*, y_1, y_2, y_3,\ldots\in P\}.\]
Then $W_P$ is a prefix-independent fairly mixing winning condition with $P = \mathsf{per}(W_P)$.
\end{proposition}

Let us start with a derivation of Theorem \ref{thm:pos}.
\begin{proof}[Proof of Theorem \ref{thm:pos} (modulo Proposition \ref{prop:technical})]
Assume that $W\subseteq C^\omega$ is an energy condition over a totally ordered group $(G, +, \le)$, defined by a valuation of colors $\mathbf{val}\colon C \to G$. Set $P = \{y\in C^+ \mid \mathbf{val}(y) < 0\}$. We claim that $W = W_P$. Indeed, $W$ consists of all $\alpha = c_1 c_2 c_3\ldots\in C^\omega$ such that
\[\mathbf{\overline{val}}(\alpha) = (\mathbf{val}(c_1), \mathbf{val}(c_1 c_2), \mathbf{val}(c_1 c_2 c_3), \ldots)\]
has an infinite decreasing subsequence. Consider any $i < j$. Observe that the $j$th element of $\mathbf{\overline{val}}(\alpha)$ is smaller than the $i$th element of $\mathbf{\overline{val}}(\alpha)$ if and only if
\[-\mathbf{val}(c_1\ldots c_i) + \mathbf{val}(c_1 \ldots c_j)  = \mathbf{val}(c_{i + 1}\ldots c_j) < 0.\]
In other words, $\mathbf{\overline{val}}(\alpha)$ has an infinite decreasing subsequence if and only if $\alpha = c_1 c_2 c_3\ldots$ can be represented, except for some finite prefix, as a as a sequence of words with negative valuations. This means that $W = W_P$.

We now show that both $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts. By Proposition \ref{prop:technical}, this would imply that $W = W_P$ is fairly mixing. In turn, by Theorem \ref{thm:gztheorem}, this implies that $W$ is bi-positional.

Consider any two words $x, y\in C^+$. Obviously:
\begin{align*}
    \mathbf{val}(x) < 0, \mathbf{val}(y) < 0 &\implies \mathbf{val}(xy) = \mathbf{val}(x) + \mathbf{val}(y) < 0,\\
    \mathbf{val}(x) \ge 0, \mathbf{val}(y) \ge 0 &\implies \mathbf{val}(xy) = \mathbf{val}(x) + \mathbf{val}(y) \ge 0.
\end{align*}
This demonstrates that both $P$ and $C^+\setminus P$ are closed under concatenations. Now, we claim that $\mathbf{val}(c_1 c_2\ldots c_n) < 0 \iff \mathbf{val}(c_2\ldots c_n c_1) < 0$ for any word $c_1c_2\ldots c_n\in C^+$ (this implies that both $P$ and $C^+\setminus P$ are closed under cyclic shifts). Indeed,
\begin{align*}
     &\mathbf{val}(c_1) + \mathbf{val}(c_2) +  \ldots + \mathbf{val}(c_n) < 0 \\
    &\iff -\mathbf{val}(c_1) + (\mathbf{val}(c_1) + \mathbf{val}(c_2) +  \ldots + \mathbf{val}(c_n)) + \mathbf{val}(c_1) < -\mathbf{val}(c_1) + 0 + \mathbf{val}(c_1)\\
    &\iff \mathbf{val}(c_2) + \ldots + \mathbf{val}(c_n) + \mathbf{val}(c_1) < 0.
\end{align*}
\end{proof}

\begin{proof}[Proof of Proposition \ref{prop:technical}]
Prefix-independence of $W_P$ is immediate. We now show that $P = \mathsf{per}(W_P)$. We have $z^\omega\in W_P$ for any $z\in P$ by definition. Hence, $P\subseteq \mathsf{per}(W_P)$. Now, take any $z\in \mathsf{per}(W_P)$. We show that $z\in P$. By definition of $\mathsf{per}(W_P)$, we have $z^\omega = xy_1 y_2 y_3\ldots$ for some $x\in C^*$ and $y_1,y_2,y_3\ldots \in P$. There exist $i < j$ such that $|xy_1\ldots y_i|$ and $|xy_1\ldots y_j|$ are equal modulo $|z|$. This means that $y_{j+1} \ldots y_j$ must be a multiple of some cyclic shift of $z$. We have that $y_{j+1} \ldots y_j\in P$ because $P$ is closed under concatenations. This means that this cyclic shift of $z$ also belongs to $P$. Indeed, otherwise we could write  $y_{j+1} \ldots y_j$ as a multiple of some word from $C^+\setminus P$, and this is impossible because $C^+\setminus P$ is closed under concatenations. Since $P$ is closed under cyclic shifts, we obtain $z\in P$.

Finally, we show that $W_P$ is fairly mixing. Since $W_P$ is prefix-independent, we should care only about the third item of Definition \ref{def:fairly_mixing}. That is, we only have to show the following two claims:

\begin{align}
\label{eq:simple}
&\bigg[( x_1 x_3 x_5\ldots\in W_P) \land (x_2 x_4 x_6\ldots\in W_P) \land(\forall n\ge1\,\, x_n^\omega\in W_P) \bigg] \implies x_1 x_2 x_3\ldots\in W_P,\\
\label{eq:hard}
&\bigg[( x_1 x_3 x_5\ldots\in \overline{W_P}) \land (x_2 x_4 x_6\ldots\in \overline{W_P}) \land(\forall n\ge 1\,\, x_n^\omega\in \overline{W_P}) \bigg] \implies x_1 x_2 x_3\ldots\in \overline{W_P},
\end{align}
for every infinite sequence of words $x_1, x_2, x_3,\ldots \in C^+$. Here, for brevity, by $\overline{W_P}$ we denote $C^\omega\setminus W_P$.

We first show \eqref{eq:simple}. If $x_n^\omega\in W_P$ for every $n$, then $x_n\in \mathsf{per}(W_P) = P$ for every $n$, and hence $x_1x_2x_3\ldots\in W_P$ by definition.

A proof of \eqref{eq:hard} is more elaborate. Assume for contradiction that $x_1 x_2 x_3\ldots \in W_P$. Then we can write $x_1 x_2 x_3\ldots = x y_1 y_2 y_3\ldots $ for some $x\in C^*$ and $y_1, y_2, y_3,\ldots\in P$. One can represent the equality as a sequence  of ``cuts'' inside $x_1 x_2 x_3\ldots$, as on the following picture:

\begin{center}
\begin{tikzpicture}
\draw[draw=black] (0,0) rectangle ++(2,1);

\draw[draw=black] (2,0) rectangle ++(4,1);

\draw[draw=black] (6,0) rectangle ++(3,1);

\draw[draw=black] (9,0) rectangle ++(1,1);

\draw[draw=black] (10,0) rectangle ++(3,1);

\draw[draw=red] (1,0) -- (1,1);

\draw[draw=red] (11.5,0) -- (11.5,1);

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(0,1) -- (2,1) node [black,midway,yshift=0.6cm] {$x_1$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(2,1) -- (6,1) node [black,midway,yshift=0.6cm] {$x_2$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(6,1) -- (9,1) node [black,midway,yshift=0.6cm] {$x_3$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(9,1) -- (10,1) node [black,midway,yshift=0.6cm] {$x_4$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(10,1) -- (13,1) node [black,midway,yshift=0.6cm] {$x_5$};

	\draw [-Latex, color=red](1,-1) -- (1,0);

	\draw [-Latex, color=red](11.5,-1) -- (11.5,0);

\draw [decorate,decoration={brace,amplitude=10pt, mirror},xshift=0pt,yshift=0pt]
(0,0) -- (1,0) node [black,midway,yshift=-0.6cm] {$x$};

\draw [decorate,decoration={brace,amplitude=10pt, mirror},xshift=0pt,yshift=0pt]
(1,0) -- (11.5,0) node [black,midway,yshift=-0.6cm] {$y_1$};

\node[draw=none] (a) at (1.5, 0.5) {$a$};

\node[draw=none] (b) at (10.75, 0.5) {$b$};

\node[draw=none] (a) at (1, -1.2) {\small \textcolor{red}{first cut}};
\node[draw=none] (a) at (11.5, -1.2) {\small \textcolor{red}{second cut}};

\node[draw=none] (a) at (13.7, 0.5) {\huge $\ldots$};
\end{tikzpicture}
\end{center}

Either there are infinitely many cuts inside $x_n$ with odd indices, or there are infinitely many cuts inside $x_n$ with even indices. Without loss of generality, we may assume that we only have cuts inside $x_n$ with odd indices, and at most one for each $n$. Indeed, if necessary, we can join several successive $y_i$'s into one word (this is legal because $P$ is closed under concatenations).

We can now write each $y_i$ as $y_i = ax_{2k} x_{2k + 1}\ldots x_{2m}b$ for some $a,b\in C^*$ and $1\le k\le m$. Now, let $y_i^\prime = a x_{2k + 1} x_{2k + 3} \ldots x_{2m - 1} b$ be a word which can be obtained from $y_i$ by removing $x_n$ with even indices. Additionally, we let $x^\prime\in C^*$ be a word
which can be obtained from $x$ in the same way.
 Since each $x_n$ with an even index lies entirely in some $y_i$ or in $x$, we have that $x_1x_3x_5\ldots = x^\prime y_1^\prime y_2^\prime y_3^\prime\ldots$, as the following picture illustrates:
\begin{center}
\begin{tikzpicture}
\draw[draw=black] (0,0) rectangle ++(2,1);

\draw[draw=black] (2,0) rectangle ++(3,1);

\draw[draw=black] (5,0) rectangle ++(3,1);

\draw[draw=red] (1,0) -- (1,1);

\draw[draw=red] (6.5,0) -- (6.5,1);

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(0,1) -- (2,1) node [black,midway,yshift=0.6cm] {$x_1$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(2,1) -- (5,1) node [black,midway,yshift=0.6cm] {$x_3$};

\draw [decorate,decoration={brace,amplitude=10pt},xshift=0pt,yshift=0pt]
(5,1) -- (8,1) node [black,midway,yshift=0.6cm] {$x_5$};

\draw [decorate,decoration={brace,amplitude=10pt, mirror},xshift=0pt,yshift=0pt]
(0,0) -- (1,0) node [black,midway,yshift=-0.6cm] {$x^\prime$};

\draw [decorate,decoration={brace,amplitude=10pt, mirror},xshift=0pt,yshift=0pt]
(1,0) -- (6.5,0) node [black,midway,yshift=-0.6cm] {$y_1^\prime$};

	\draw [-Latex, color=red](1,-1) -- (1,0);

	\draw [-Latex, color=red](6.5,-1) -- (6.5,0);

\node[draw=none] (a) at (1.5, 0.5) {$a$};

\node[draw=none] (b) at (5.75, 0.5) {$b$};

\node[draw=none] (a) at (1, -1.2) {\small \textcolor{red}{first cut}};
\node[draw=none] (a) at (6.5, -1.2) {\small \textcolor{red}{second cut}};

\node[draw=none] (a) at (8.7, 0.5) {\huge $\ldots$};
\end{tikzpicture}
\end{center}
We will show that $y_i^\prime\in P$ for every $P$. This would contradict a fact that $x_1 x_3 x_5\ldots \in \overline{W_P}$.

First, observe that $x_n\notin P$ for every $n$. Indeed, we are given that $x_n^\omega\in \overline{W_P}$ for every $n$. Hence, $x_n \notin \mathsf{per}(W_P) = P$, as required.

Assume for contradiction that $y_i^\prime = ax_{2k + 1} x_{2k + 3}\ldots x_{2m - 1}b\notin P$.  Using a fact that $C^+\setminus P$ is closed under concatenations and cyclic shifts, we obtain:

\begin{align*}
y_i^\prime = &ax_{2k + 1} x_{2k + 3}\ldots x_{2m - 1} b\notin  P \\
\implies &x_{2k + 1} x_{2k + 3}\ldots x_{2m - 1} ba\notin P \\
\implies &x_{2k}x_{2k + 1} x_{2k + 3}\ldots x_{2m - 1} ba\notin  P && \mbox{because } x_{2k}\notin  P\\
\implies & x_{2k + 3}\ldots x_{2m - 1} ba x_{2k} x_{2k + 1}\notin  P \\
\ \implies &x_{2k + 2}x_{2k + 3}\ldots x_{2m - 1} ba x_{2k} x_{2k + 1}\notin  P && \mbox{because } x_{2k + 2}\notin  P\\
&\vdots\\
\implies &x_{2m} ba x_{2k} x_{2k + 1}\ldots x_{2m - 1} \notin P && \mbox{because } x_{2m}\notin  P\\
\implies &y_i = a x_{2k} x_{2k + 1}\ldots x_{2m} b\notin P,
\end{align*}
contradiction.
\end{proof}

\section{Refuting Conjecture \ref{conj}}
\label{sec:refute}

Consider the free group $F_{\{a, b\}}$ with 2 generators $a, b$. As was proved by Shimbireva~\cite{shimbireva1947theory}, see also~\cite[Page 18]{deroin2014groups}, free groups can be totally ordered. We take an arbitrary total ordering $\le$ of $F_{\{a, b\}}$. We also consider its inverse $\le^{-1}$, which is also a total ordering of $F_{\{a, b\}}$. Define a set of colors $C = \{a, a^{-1}, b, b^{-1}, \varepsilon\}$. Here $a^{-1}, b^{-1}$ are inverses of $a, b$ in $F_{\{a, b\}}$, and $\varepsilon$ is the identity element of $F_{\{a,b\}}$.

Let $W_1\subseteq C^\omega$ be an energy condition over $(F_{\{a, b\}}, \le)$, defined by a (suggestive) valuation of colors which interprets elements of $C$ as corresponding elements of $F_{\{a,b\}}$. Similarly, we let $W_2\subseteq C^\omega$ be an energy condition over $(F_{\{a, b\}}, \le^{-1})$, defined by the same valuation. The only difference between $W_1$ and $W_2$ is that they are defined w.r.t.~different total orderings of $F_{\{a, b\}}$ (one ordering is the inverse of the other one).

We show that the union $W_1\cup W_2$ is not half-positional. It consists of all $\alpha \in C^\omega$ such that $\mathbf{\overline{val}}(\alpha)$  contains either an infinite decreasing subsequence w.r.t.~$\le$ or an infinite decreasing subsequence w.r.t.~$\le^{-1}$. In other words, it consists of  all $\alpha \in C^\omega$ such that $\mathbf{\overline{val}}(\alpha)$  contains either an infinite decreasing subsequence or an infinite \emph{increasing} subsequence w.r.t.~$\le$.

We show that $W_1\cup W_2$ is not half-positional in the following arena.

\begin{center}
\begin{tikzpicture}
\node[draw,  regular polygon, regular polygon sides=4, minimum size=1cm] (1) {};

\node[draw,  circle, right=2cm of 1,minimum size=1cm] (2) {};

\node[draw,  circle, left=2cm of 1, minimum size=1cm] (3) {};

\draw[thick,->] (1) -- (2) node[midway, above] {$\varepsilon$};

\draw[thick,->] (1) -- (3) node[midway, above] {$\varepsilon$};

\path[->]
 (2) edge [thick, in=60,out=120,out distance=1cm,in distance=1cm] node[midway, above] {$b$} (1);

\path[->]
 (2) edge [thick, in=-60,out=-120,out distance=1cm,in distance=1cm] node[midway, below] {$b^{-1}$} (1);

\path[->]
 (3) edge [thick, in=120,out=60,out distance=1cm,in distance=1cm] node[midway, above] {$a$} (1);

\path[->]
 (3) edge [thick, in=-120,out=-60,out distance=1cm,in distance=1cm] node[midway, below] {$a^{-1}$} (1);
\end{tikzpicture}
\end{center}

Here, Alice controls the square and Bob controls the two circles.
Assume that the game starts in the square. We show that Alice has a winning strategy w.r.t.~$W_1\cup W_2$, but not a positional one.

Alice has two positional strategies in this arena: always go to the left and always go to the right. Consider, for example, the first one. Bob has the following counter-strategy which wins against it: alternate the $a$-edge with the $a^{-1}$-edge. We get the following sequence of colors in the play of these two strategies:
\[\varepsilon a \varepsilon a^{-1} \varepsilon a \varepsilon a^{-1}\ldots\]
This sequence does not belong to $W_1\cup W_2$ because
\[\mathbf{\overline{val}}(\varepsilon a \varepsilon a^{-1} \varepsilon a \varepsilon a^{-1}\ldots) = \varepsilon, a, a, \varepsilon, \varepsilon, a, a, \varepsilon,\ldots \]
There are only two distinct elements of $F_{\{a,b\}}$ occurring in $\mathbf{\overline{val}}(\varepsilon a \varepsilon a^{-1} \varepsilon a \varepsilon a^{-1}\ldots)$. Hence, it neither has an infinite decreasing subsequence nor an infinite increasing subsequence. By the same argument, the second positional strategy of Alice (always go to the right) is not winning w.r.t.~$W_1\cup W_2$ either.

On the other hand, Alice has the following winning strategy: alternate the edge to the left circle with the edge to the right circle. Consider any play with this strategy. Its sequence of colors looks as follows:
\[\varepsilon a^{\pm 1} \varepsilon b^{\pm1} \varepsilon a^{\pm 1} \varepsilon b^{\pm1}\ldots\]
We show that this sequence belongs to $W_1\cup W_2$.
A restriction of $\mathbf{\overline{val}}(\varepsilon a^{\pm 1} \varepsilon b^{\pm1} \varepsilon a^{\pm 1} \varepsilon b^{\pm1})$ to elements with even indices looks like this:
\begin{equation}
\label{eq:sec}
a^{\pm 1}, \,\, a^{\pm 1}b^{\pm 1}, \,\,  a^{\pm 1}b^{\pm 1} a^{\pm 1}, \,\, a^{\pm 1}b^{\pm 1} a^{\pm 1}b^{\pm 1} \ldots
\end{equation}

All elements of \eqref{eq:sec} are distinct. Hence, by the Infinite Ramsey Theorem, it either has an infinite decreasing subsequence or an infinite increasing subsequence w.r.t.~$\le$. Indeed, consider an infinite complete graph over $\{1, 2, 3, \ldots\}$, whose edges are colored into green and red as follows. Pick any $i, j\in\{1, 2, 3,\ldots\}$, $i < j$. If the $i$th element of \eqref{eq:sec} is bigger than the $j$th element of \eqref{eq:sec}, then color the edge between $i$ and $j$ into green. Otherwise, color this edge into red (in this case,  the $i$th element of \eqref{eq:sec} is smaller than the $j$th element of \eqref{eq:sec}). Our graph has an infinite induced subgraph in which all edges are of the same color. If they are all green (resp., red), then this subgraph defines an infinite decreasing (resp., increasing) subsequence of \eqref{eq:sec}.

\bigskip

\textbf{Additional remarks.} Energy conditions over free groups are interesting because they are \emph{non-permuting} (if there is more than one generator). A prefix-independent winning condition is permuting if it is closed under permuting periods of periodic sequences. All previously known prefix-independent bi-positional winning condition were permuting. This is because they can be seen as energy conditions over abelian groups (on periodic sequences). In a talk of Colcombet and Niwi\'{n}ski~\cite{slides} it was asked whether there exists a non-permuting bi-positional prefix-independent winning condition. The answer is ``yes''. For example, take $W_1$ as above in this section. Without loss of generality, we may assume that $aba^{-1} b^{-1}$ is negative w.r.t.~$\le$ (otherwise we can consider its inverse). Then $(aba^{-1}b^{-1})^\omega\in W_1$, but $(aa^{-1} b b^{-1})^\omega\notin W_1$.

Additionally, the winning condition $W_1\cup W_2$ is interesting because it refutes  1-to-2-player lifting for half-positionality. Namely, it is easy to see that $W_1\cup W_2$ is positional for Alice in all arenas, where there are no nodes of Bob. This is because she can win in such arenas if and only if there is a reachable non-zero simple cycle. But as we have shown, $W_1\cup W_2$ is not positional for Alice in the presence of Bob. Previously, there were examples that refute 1-to-2-player lifting for half-positionality in stochastic games~\cite{gimbert2014two}.
\section{Discussing Conjecture \ref{our_conj}}
\label{sec:conj}

First, it is useful to understand how prefix-independent bi-positional winning condition are arranged on periodic sequences.
 Luckily, Proposition \ref{prop:technical} gives an answer.

\begin{proposition}
\label{prop:periods}
 Let $C$ be a finite non-empty set. Then for any $P\subseteq C^+$ the following two conditions are equivalent:
\begin{enumerate}[label=\Alph*)]
\item $P = \mathsf{per}(W)$ for some prefix-independent bi-positional winning condition $W\subseteq C^\omega$;
\item $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts;
\end{enumerate}

\end{proposition}
\begin{proof}
The fact that the second item implies the first item follows from Proposition \ref{prop:technical}. Indeed, if $P$ and $C^+\setminus P$  are closed under concatenations and cyclic shifts, then $P = \mathsf{per}(W_P)$ for a prefix-independent fairly mixing winning condition $W_P$, which is bi-positional by Theorem \ref{thm:gztheorem}. We now show that the first item implies the second item. The fact that $P$ and $C^+\setminus P$ are closed under cyclic shifts is a consequence of the prefix-independence of $W$:
\begin{align*}
    c_1 c_2 \ldots c_n \in P &\iff (c_1 c_2 \ldots c_n)^\omega \in W \iff c_n (c_1 c_2 \ldots c_n)^\omega = (c_n c_1\ldots c_{n - 1})^\omega \in W\\
    &\iff c_n c_1\ldots c_{n - 1}\in P.
\end{align*}
We now show that $P$ is closed under concatenations (there is a similar argument for $C^+\setminus P$). Take any $x,y\in P$. Consider the following arena.

\begin{center}
\begin{tikzpicture}

\node[draw,  circle, right=2cm of 1,minimum size=0.5cm] (A) {};

\path[->] (A) edge [anchor=center,loop below, looseness=15, out=60, in=120] node {$x$} (A);

\path[->] (A) edge [anchor=center,loop below, looseness=15, out=-60, in=-120] node {$y$} (A);
\end{tikzpicture}
\end{center}

It has a central circle node which lies on two simple cycles, one of which is colored by $x$ and the other one by $y$. All nodes are controlled by Bob. Since, $x, y\in P$, we have that $x^\omega,y^\omega\in W$. Hence, Bob does not have a positional winning strategy w.r.t.~$W$ from the central circle. Since $W$ is bi-positional, Bob has no winning strategy from the central circle w.r.t.~$W$. Now, assume that Bob alternates the $x$-cycle with the $y$-cycle. He obtains $(xy)^\omega$ as a sequence of colors. Since this strategy is not winning, we have $xy\in P$.
\end{proof}

In turn, periods of periodic sequences of ETOG conditions are arranged as follows.

\begin{proposition}
Let $C$ be a non-empty finite set and $W\subseteq C^\omega$ be an energy condition over a totally ordered group $(G, +, \le)$, defined by a valuation of colors $\mathbf{val}\colon C\to G$. Then $\mathsf{per}(W) = \{x\in C^+\mid \mathbf{val}(x) < 0\}$.
\end{proposition}
\begin{proof}
 Define $P = \{x\in C^+\mid \mathbf{val}(x) < 0\}$. By the argument from the derivation of Theorem \ref{thm:pos}, we have $W = W_P$. Moreover, it was shown there that $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts. Finally, by Proposition \ref{prop:technical}, we have that $P = \mathsf{per}(W_P) = \mathsf{per}(W)$.
\end{proof}

Thus, Conjecture \ref{our_conj} is equivalent to the following conjecture.
\begin{conjecture}
\label{our_conj2}
Let $C$ be any non-empty finite set. Then for any $P\subseteq C^+$ such that $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts there exist a totally ordered group $(G, +, \le)$ and a valuation of colors $\mathbf{val}\colon C\to G$ such that $P = \{x\in C^+\mid \mathbf{val}(x) < 0\}$.
\end{conjecture}

It might be concerning that $P$ and  $C^+\setminus P$ are interchangeable in Conjecture \ref{our_conj2}, while $\mathbf{val}$ treats them asymmetrically. Namely, we require it to be negative on $P$ and \emph{non-negative} on $C^+\setminus P$. However, $\mathbf{val}$ can always be made strictly positive on $C^+\setminus P$. Namely, instead of $G$, consider the direct product $G\times\mathbb{Z}$ with the lexicographic order, and define a new valuation of colors $\mathbf{val}^\prime\colon C\to G\times\mathbb{Z}$,  $\mathbf{val}^\prime(c) = (\mathbf{val}(c), 1)$.

Finally, we notice that our conjecture can be reduced to a reasoning about free groups.
\begin{definition}
A subset $S$ of a group $G$ is called an \textbf{invariant sub-semigroup} of $G$ if the following two conditions hold:
\begin{enumerate}[label=\Alph*)]
\item $xy\in S$ for all $x, y\in S$ (closure under multiplications);
\item $g x g^{-1} \in S$ for all $g\in G, x\in S$ (closure under conjugations with elements of $G$).
\end{enumerate}
\end{definition}
\begin{conjecture}
\label{our_conj3}
Consider an arbitrary non-empty finite set $C$ and  any $P\subseteq C^+$ such that $P$ and $C^+\setminus P$ are closed under concatenations and cyclic shifts. Then there exists an invariant sub-semigroup $S$
of the free group $F_C$ such that, first, $C^+\setminus P$ is a subset of $S$, second, $P$ is disjoint with $S$, and third, for every $g\in F_C$ either $g\in S$ or $g^{-1}\in S$ (in particular, $S$ must have the neutral element).
\end{conjecture}

\begin{proposition}
Conjecture \ref{our_conj2} is equivalent to Conjecture \ref{our_conj3}.
\end{proposition}
\begin{proof}

Consider an arbitrary non-empty finite set $C$.
It is sufficient to show that for any $P\subseteq C^+$ the following two conditions are equivalent:
\begin{enumerate}[label=\Alph*)]
\item  there exist a totally ordered group $(G, +, \le)$ and a valuation of colors $\mathbf{val}\colon C\to G$ such that $P = \{x\in C^+\mid \mathbf{val}(x) < 0\}$.

\item there exists an invariant sub-semigroup $S$
of the free group $F_C$ such that, first, $C^+\setminus P$ is a subset of $S$, second, $P$ is disjoint with $S$, and third, for every $g\in F_C$ either $g\in S$ or $g^{-1}\in S$.
\end{enumerate}
We first establish $\mbox{A)}\implies\mbox{B)}$. Extend $\mathbf{val}$ to a homomorphism from $F_C$ to $G$ by setting $\mathbf{val}(c^{-1}) = -\mathbf{val}(c)$ for $c\in C$.
Set $S = \{g\in F_C\mid \mathbf{val}(g) \ge 0\}$. It is easy to check that all conditions on $S$ are satisfied.

Now we establish $\mbox{B)}\implies\mbox{A)}$. Let $S$ be as in B). Consider a binary relation $\sim$ on $F_C$, defined by $f\sim g\iff fg^{-1}, gf^{-1}\in S$ for $f, g\in F_C$. A fact that $S$ is an invariant sub-semigroup with the neutral element implies that $\sim$ is a congruence on the group $F_C$. Let $G = F_C/\sim$ be the corresponding quotient group. Now, consider a binary relation $\preceq$ on $F_C$, defined by $f\preceq g\iff gf^{-1}\in S$ for $f, g\in F_C$ (observe that $f\sim g \iff f\preceq g, g\preceq f$). It is easy to see that $\preceq$ is correctly defined over $F_C/\sim$, whose elements are equivalence classes of $\sim$. More formally, it holds that if $a\sim b, x\sim y$, then $a\preceq x \iff b\preceq y$ (it can again be derived from the fact that $S$ is an invariant sub-semigroup). It is also routine to check that $\preceq$ defines a total ordering on $G$. We need a condition that either $g\in S$ or $g^{-1}\in S$ for every $g\in F_C$ only to show the totality of our order. Namely, to show that there are no $f, g\in F_C$ with $f\not\preceq g$ and $g\not\preceq f$, we notice that otherwise neither $g f^{-1}$ nor $f g^{-1} = (g f^{-1})^{-1}$ are in $S$. Observe that the equivalence class of $g\in F_C$ w.r.t.~$\sim$ is non-negative in $(G, \preceq)$ if and only if $g\in S$. Now, recall that $C^+\setminus P$ is a subset of $S$ and $P$ is disjoint with $P$.
Hence, if we consider a valuation of colors $\mathbf{val}\colon C\to G$, which maps $c\in C$ to its equivalence class w.r.t.~$\sim$, then $P$ would be the set of words from $C^+$ whose valuation is negative w.r.t.~$\preceq$.
\end{proof}


\section{Introduction}\label{Introduction}
The study of behavior of charged particle beams is important from both theoretical and experimental point of views. In the classical theory of charged particle beams, the purpose is to study the position and motion of the beam along the optical axis while in the quantum theory of charged particles beams, one is interested in time-evolution of the matter-wave or the density matrix describing the beam \cite{1}. The study of the optical properties of quantum charged particle beams has been a significant part of studies and research mainly devoted to the effect of the electromagnetic and gravitational fields on the quantum mechanical phase of the matter-waves \cite{2,3,4,5}.
Quantum beams are directional beams consisting of a stream of charged particles like electrons, protons, ions, etc., moving almost in one direction. From wave-particle duality this stream of particles can be interpreted as a matter-wave propagating in the same direction \cite{6}. The quantum dynamics of the matter-wave is governed by a complicated many body Schrodinger equation containing both the external and internal electromagnetic fields \cite{1}. In the scalar theory of quantum charged beam optics, one deals with a stream of spin-$0$ charged particles \cite{7,8,9}. A natural feature of all beams with the helical phase is the orbital angular momentum that can be generated in the laboratory \cite{10}. Electron microscopy uses electron beams and wave-like properties of the electrons guided through magnetic lenses to achieve a clear and high resolution image of the object. The free electron beams carrying orbital angular momentum (OAM) was investigated in a seminal paper by Bliokh \cite{11}, which triggered extensive experimental works directed toward the generation and applications of structured electron beams \cite{12,13,14}. For a non-interacting beam one can consider the dynamics of a single particle or wave matter. Here for a class of Hamiltonians describing the motion of a neutral or charged particle in the presence of a constant magnetic field and under the influence of a time-dependent external force, we find exact quantum propagators. The scheme is based on the Heisenberg equations of motion and propagator properties.  The initial wave function is taken to be a Gauss-Hermite wave function, though one can also use Laguerre-Gauss wave functions, and for the evolved wave function we will obtain the uncertainties, orbital angular momentum, and the inertia tensor in the center of mass frame of the density function. From the point of view of interference of matter waves, the exact results obtained here give more exact results compared to the approximate methods like the axial approximation. The results obtained here can also be of importance in the non-relativistic quantum electron microscopy.
\section{Matter-wave interacting with a linear time-dependent potential}\label{tdlp}
The Hamiltonian for a matter wave under the influence of a time-dependent force is
\begin{equation}\label{1}
\hat{H}=\frac{\hat{p}_x^2+\hat{p}_y^2+\hat{p}_z^2}{2 m} + \mu(t)\,\hat{x},
\end{equation}
where $\mu(t)\,\hat{x}$ is the external potential exerting a time-dependent force $-\mu (t)$ on the matter wave along the $x$ axis. Hermit-Gaussian beams are a subset of stable laser states with symmetry along the axis of propagation. The initial wave function ($\psi(x,y,z;0)$) is considered to be a Hermite-Gauss wave function. Another important class of initial wave functions are Laguerre-Gauss wave packets having orbital angular momentum and their time-evolution can be obtained similarly. The initial wave function as a wave packet, starts its motion from the origin with the initial momentum $\hbar k_0/m$ in the direction of $z$ axis.
\begin{equation}\label{2}
\psi(x,y,z;0)=(\frac{\alpha}{\pi})^{\frac{1}{4}}\,\sqrt{\frac{2}{\omega_0^2 \pi 2^{m+n}\,n!m!}}\,\exp(-\frac{x^2+y^2}{\omega_0^2})\,\exp(-\frac{\alpha}{2} z^2 + i k_0 z)\,H_n(\frac{\sqrt{2}}{\omega_0} x)\,H_m(\frac{\sqrt{2}}{\omega_0} y),
\end{equation}
where $\omega_0$ is the beam waist and $H_n (x)$ is the Hermite polynomial of order $n$.
Here, we follow a non-relativistic approach to find the density evolution of the matter wave $|\psi(x,y,z:t)|^2$. In the following, In Sec.III, we also consider the effect of a constant magnetic field by modifying the Hamiltonian trough the minimal coupling model and study the density evolution.
\begin{figure}
\centering
\includegraphics[scale=0.4]
{F1.pdf}
\caption{(Color online) The density function $|\psi(x,y,z;0)|^2_{z=0}$ corresponding to Eq. (\ref{2}) on the x-y plane ($z=0$) for $n=m=1$.}
\end{figure}
\subsection{Heisenberg equations of motion}
By making use of the Hamiltonian Eq. (\ref{1}) and the Heisenberg equations for linear momentum operators, we find
\begin{eqnarray}\label{4}
&& i\hbar\, \dot{\hat{p}}_x=[\hat{p}_x,\hat{H}]=-i\hbar\, \mu(t),\nonumber\\
&& i\hbar\, \dot{\hat{p}}_y=[\hat{p}_y,\hat{H}]=0,\nonumber\\
&& i\hbar\, \dot{\hat{p}}_z=[\hat{p}_z,\hat{H}]=0.
\end{eqnarray}
Therefore, the momentums along the $y$ and $z$ axis are constants of motion $\hat{p}_y(t)=\hat{p}_y(0),\,\,\, \hat{p}_z(t)=\hat{p}_z(0)$. The linear momentum along the $x$ direction is
\begin{equation}\label{5}
\hat{p}_x(t)=\hat{p}_x(0)-\nu(t),
\end{equation}
where
\begin{equation}\label{6}
\nu(t)=\int_0^t dt'\,\mu(t').
\end{equation}
The Heisenberg equations for the position operators lead to
\begin{eqnarray}\label{7}
\hat{x}(t) &=& \hat{x}(0)+\frac{t}{m}\hat{p}_x(0)-\frac{\xi(t)}{m},\nonumber\\
\hat{y}(t) &=& \hat{y}(0)+\frac{t}{m}\hat{p}_y(0),\nonumber\\
\hat{z}(t) &=& \hat{z}(0)+\frac{t}{m}\hat{p}_z(0),
\end{eqnarray}
where $\xi(t)=\int_0^t dt'\,\nu(t')$. The classical trajectory of the center of mass of the matter wave can be found from the expectation values of the position operators in the initial state Eq. (\ref{2})
\begin{eqnarray}\label{8}
&& \langle\hat{x}(t)\rangle=\int dx\,dy\,dz\,\psi^*(x,y,z;0)\,\hat{x}(t)\,\psi(x,y,z;0)=\frac{-\xi(t)}{m},\nonumber\\
&& \langle\hat{y}(t)\rangle=0,\nonumber\\
&& \langle\hat{z}(t)\rangle=\frac{\hbar k\,t}{m}.
\end{eqnarray}
In a constant force field, we have $\xi(t)=\mu\,t^2/2$ and the trajectory will be a paraboloid as expected.
\subsection{Orbital Angular Momentum}
The Orbital Angular Momentum (OAM) components corresponding to the matter wave are
\begin{eqnarray}\label{10}
\hat{l}_x(t)=&& \hat{y}(t)\,\hat{p}_z(t)-\hat{z}(t)\,\hat{p}_y(t)=\Big(\hat{y}(0)+
\frac{t}{m}\,\hat{p}_y(0)\Big)\,\hat{p}_z(0)-\Big(\hat{z}(0)+\frac{t}{m}
\,\hat{p}_z(0)\Big)\,\hat{p}_y(0),\nonumber\\
=&& \hat{l}_x(0),\nonumber\\
\hat{l}_y(t)=&& \Big(\hat{z}(0)+\frac{t}{m}\,\hat{p}_z(0)\Big)\,\Big(\hat{p}_x(0)
  -\nu(t)\Big)-\Big(\hat{x}(0)+\frac{t}{m}\,\hat{p}_x(0)-\frac{\xi(t)}{m}
  \Big)\,\hat{p}_z(0),\nonumber\\
\hat{l}_z(t)=&& \Big(\hat{x}(0)+\frac{t}{m}\,\hat{p}_x(0)-\frac{\xi(t)}{m}\Big)\,\hat{p}_x(0)-
\Big(\hat{y}(0)+\frac{t}{m}\,\hat{p}_y(0)\Big)(\hat{p}_x(0)-\nu (t)).
\end{eqnarray}
The OAM along the $x$ axis is a constant of motion $\hat{l}_x(t)=\hat{l}_x(0)$. The expectation values of $\langle\hat{l}_x(t)\rangle$ and $\langle\hat{l}_z(t)\rangle$ in the initial state Eq. (\ref{2}) are zero and for $\langle\hat{l}_y(t)\rangle$ one finds
\begin{equation}\label{11}
  \langle\hat{l}_y(t)\rangle=\frac{\hbar k}{m} \xi(t)+\frac{\hbar k t}{m} \nu(t).
\end{equation}
\subsection{Exact propagator}
Hamiltonian Eq. (\ref{1}) can be decomposed into commuting transverse ($\hat{H}_t$) and longitudinal ($\hat{H}_l$) parts as
\begin{equation}\label{12}
\hat{H}=\underbrace{\frac{\hat{p}_y^2+\hat{p}_z^2}{2 m}}_{\hat{H}_{t}} +\underbrace{\frac{\hat{p}_x^2}{2m}+\mu(t)\,\hat{x}}_{\hat{H}_l}.
\end{equation}
and accordingly, the total quantum propagator can be written as
\begin{equation}\label{13}
k(x,y,z,t|x',y',z',t)=k_t(y,z,t|y',z',0)\,k_l(x;t|x',0),
\end{equation}
where $k_t \,(k_l) $ are transverse (longitudinal) propagators satisfying
\begin{eqnarray}\label{14}
&& \big[i\hbar\,\partial_t+\frac{\hbar^2}{2m}(\partial^2_y+\partial^2_z)\big]\,k_t(y,z,t|y',z',0)
=0,\nonumber\\
&& \big[i\hbar\,\partial_t+\frac{\hbar^2}{2m}\partial^2_x-\mu(t)\,x\big]\,k_l(x,t|x',0)
=0,
\end{eqnarray}
with the explicit solutions \cite{15}
\begin{eqnarray}\label{15}
&& k_t(y,z,t|y',z',0) = \Big(\frac{m}{2\pi i\hbar t}\Big)\exp\Big(\frac{im}{2\hbar t}[(y-y')^2+(z-z')^2]\Big),\nonumber\\
&& k_l(x,t|x',0) = \sqrt{\frac{m}{2\pi i\hbar t}}\,\exp\Big({\frac{-i}{2 m \hbar }\int_0^t dt'\,\Big(\frac{v(t')}{t'}\Big)^2}\Big)\,\exp\Big(\frac{i m}{2 \hbar t}\Big[(x-x')^2-\frac{2\,v(t)}{m}(x-x')-\frac{2\,t\,\nu(t)\,x'}{m}\Big]\Big),
\end{eqnarray}
where $v(t)=\int_{0}^{t} dt'\,t'\,\mu(t')$ and $\nu(t)=\int_{0}^{t} dt'\,\mu(t')$. Having the initial state Eq. (\ref{2}) we can find the evolved state $\psi(x,y,z,t)$ using $\psi(x,y,z,t)=\int dx'dy' dz'$ $k_{t}(y,z,t|y',z',0)k_t (x,t|x',0)\psi(x',y',z',0)$ leading to
\begin{multline}\label{pi0}
\psi(x,y,z;t)=\\
\Big(\frac{m\,\omega_0^2}{2}\Big)\Big(\frac{1}{i\hbar t\,(1+i\lambda)}\Big)\,(\frac{\alpha}{\pi})^\frac{1}{4}\,\frac{1}{\sqrt{2\pi\omega_0^2}}
\sqrt{\frac{m}{m+i\hbar \alpha t}}\, \Big(1-\frac{2}{1+i\lambda}\Big)^\frac{m}{2}\,\Big(1-\frac{2}{1+i\lambda}\Big)^\frac{n}{2}
\,\exp\Big({\frac{-i}{2 m \hbar }\int_0^t dt'\,\Big(\frac{v(t')}{t'}\Big)^2}\Big)\\ .\exp\Big(-\frac{k_0^2}{2\alpha}\Big) \exp\Big(-\frac{y^2}{\omega_0^2}\,(\frac{\lambda^2}{1+i\lambda}+i\lambda)\Big)\, \exp\Big(\frac{(ik_0-(im/\hbar t) z)^2}{2\alpha-2im/\hbar t}\Big)\, \exp\Big(-i\lambda(\frac{x^2}{\omega_0^2})-i\gamma\frac{\sqrt{2}\,x}{\omega_0}\Big)\\
.\exp\Big(-\frac{(\lambda \sqrt{2}\frac{x}{\omega_0} + \gamma + \frac{f}{\sqrt{2}})^2}{2+2i\lambda}\Big)\, H_m\left[\frac{y\,i\,\lambda \sqrt{2}}{\omega_0(1+i\lambda)\sqrt{1-\frac{2}{1+i\lambda}}}\right]\, H_n\left[\frac{i\,\lambda \sqrt{2}\frac{x}{\omega_0} + i\,\gamma + \frac{i\,f}{\sqrt{2}}}{(1+i\lambda)\sqrt{1-\frac{2}{1+i\lambda}}}\right],
\end{multline}
where for simplicity we have defined the time-dependent dimensionless parameters $\lambda=-m\omega_0^2/2\hbar t$, $ \gamma=v(t)\,\omega_0/\sqrt{2}\,t\,\hbar$ and $f=\omega_0 \,\nu(t)/\hbar$.
\begin{figure}
\centering
\includegraphics[scale=0.4]
{F2.pdf}
\caption{(Color online) The density plot $|\psi(x,y,z;t)|^2_{z=0}$ of Eq. (\ref{pi0}), in terms of the scaled dimensionless variables $x/\omega_0$, $y/\omega_0$, $t/\tau$ where $\tau=m\omega_0^2/2\hbar$.}
\end{figure}

\begin{figure}
\centering
\includegraphics[scale=0.4]
{F2-0.pdf}
\caption{(Color online) The density plot $|\psi(x,y,z;t)|^2_{z=0}$ of Eq. (\ref{pi0}) in the presence of a constant force. The scaled dimensionless variables are $x/\omega_0$, $y/\omega_0$, $t/\tau$ and $\tau=m\omega_0^2/2\hbar$.}
\end{figure}

\begin{figure}
\centering
\includegraphics[scale=0.4]
{F3.pdf}
\caption{(Color online) The density plot $|\psi(x,y,z;t)|^2_{z=0}$ of Eq. (\ref{pi0}) in the presence of the time-dependent force $\mu(t)=\mu_0\,\sin(2t/\tau)$ for $n=m=1$ and $\tau=m\omega_0^2/2\hbar$.}
\end{figure}
\subsection{Uncertainties and the Inertia tensor }
By making use of Eq. (\ref{7}), the expectation value of $\hat{x}(t)$ is $-\xi(t)/m$ and for $\hat{x}^2(t)$ we find
\begin{eqnarray}
\Big\langle\hat{x}^2(t)\Big\rangle &=& \int dx\,dy\,dz\,\psi^*(x,y,z;0)\,\hat{x}^2(t)\,\psi(x,y,z;0),\nonumber\\
&=& \frac{3\,\omega_0^2}{4}\,+\,\frac{3\,\hbar^2\,t^2}{m^2\,\omega_0^2}\,+\,\frac{\xi^2(t)}{m^2},
\end{eqnarray}
Therefore, the uncertainty in position is $\triangle x=\sqrt{\langle\hat{x}^2(t)\rangle-\langle\hat{x}(t)\rangle^2}=(3 \omega_0^2/4+3 \hbar^2 t^2/m^2 \omega_0^2)^{1/2}$ which is spreading in time. Similarly, for $\langle\hat{p}_x(t)\rangle$ and $\langle\hat{P}^2_x(t)\rangle$ one easily finds
\begin{eqnarray}
&& \langle\hat{p}_x(t)\rangle\,=\,\int dx\,dy\,dz\,\psi^*(x,y,z;0)\,\Big(\-i\hbar \frac{\partial}{\partial x}-\nu(t)\Big)\,\psi(x,y,z;0)\,=\,-\nu(t),\nonumber\\
&& \langle\hat{p}^2_x(t)\rangle\,=\,\int dx\,dy\,dz\,\psi^*(x,y,z;0)\,\Big(\-i\hbar \frac{\partial^2}{\partial x^2}+\nu^2(t)\Big)\,\psi(x,y,z;0)\,=\,\frac{3\,\hbar^2}{\omega_0^2}+\nu^2(t),
 \end{eqnarray}
 %
Thus
\begin{equation}
\triangle p=\sqrt{\langle\hat{p}^2_x(t)\Big\rangle\,-\,\Big\langle\hat{p}_x(t)\rangle^2}=\frac{\sqrt{3}\,\hbar}{\omega_0},
\end{equation}
and $\triangle x\,\triangle p=(3\,\hbar/2)\sqrt{1+4\hbar^2\,t^2/m^2 \omega_0^4}$. The moment of inertia along the $z$ axis is defined by
\begin{eqnarray}
I_{zz}=m\langle\hat{r}^2\rangle &=& m\langle\psi(0)|\hat{x}^2+\hat{y}^2|\psi(0)\rangle\,=\,m\int dx\,dy\,dz\,\psi^*(x,y,z;0)(\hat{x}^2+\hat{y}^2)\,\psi(x,y,z;0),\nonumber\\
&=& \frac{3m}{2}\,\omega_0^2\,+\,\frac{3\,\hbar^2 t^2}{m\,\omega_0^2}\,+\,\frac{\xi^2(t)}{m},
\end{eqnarray}
and by using the general definition
\begin{equation}
I_{ij}=m\langle\hat{x}^2+\hat{y}^2+\hat{z}^2\rangle_t\,\delta_{ij}-m\langle\hat{x}_i\hat{x}_j\rangle_t,
\end{equation}
one can obtain the inertia tensor corresponding to the density function $|\psi(x,y,z;t)|^2$ as
\begin{align}
I_{ij}=\begin{bmatrix}
\frac{3m}{4}\omega_0^2+\frac{3\hbar^2 t^2}{m\omega_0^2}+\frac{m}{2\alpha}+\frac{\hbar^2 t^2(k^2+\frac{\alpha}{2})}{m} & 0 & \frac{\hbar\,k\,t}{m}\,\xi(t)
\vspace{6mm}
\\0 & \frac{3m}{4}\omega_0^2+\frac{3\hbar^2 t^2}{m\omega_0^2}+\frac{m}{2\alpha}+\frac{\hbar^2 t^2(k^2+\frac{\alpha}{2})}{m}+\frac{\xi^2(t)}{m} & 0
\vspace{6mm}
\\\frac{\hbar\,k\,t}{m}\,\xi(t) & 0 & \frac{3m}{2}\omega_0^2+\frac{6\hbar^2 t^2}{m\omega_0^2}+\frac{\xi^2(t)}{m}
\end{bmatrix}
\end{align}
The inertia tensor with respect to the center of mass is defined by  $I_{ij}^{c}=m[\langle\hat{x}^2+\hat{y}^2+\hat{z}^2\rangle_t\,\delta_{ij}]+ m[R^{2}(t)\,\delta_{ij}-R_i\,R_j]$, where $R_i,\,(i=1,2,3)$ are the center of mass coordinates.
One can easily obtain the inertia tensor in the center of mass frame as
\begin{align}\label{23I}
I_{ij}^{c}=\begin{bmatrix}
\frac{3m}{4}\omega_0^2+\frac{3 \hbar^2 t^2}{m\,\omega_0^2}+\frac{m}{2\alpha}+\frac{\hbar^2 t^2 \alpha}{2m} & 0 & 0
\vspace{6mm}
\\0 & \frac{3m}{4}\omega_0^2+\frac{3 \hbar^2 t^2}{m\,\omega_0^2}+\frac{m}{2\alpha}+\frac{\hbar^2 t^2 \alpha}{2m} & 0
\vspace{6mm}
\\0 & 0 & \frac{3m}{2}\omega_0^2+\frac{6 \hbar^2 t^2}{m\,\omega_0^2}
\end{bmatrix}.
\end{align}
From Eq. (\ref{23I}) we deduce that during the evolution of the matter wave, the shape of the density function remains symmetrical ( $I_{xx} = I_{yy}$) with respect to the $x$ and $y$ axes in the center of mass frame though the volume of the density function is getting larger with respect to the center of mass axes. We also note that the components of $I_{ij}^{c}$ are independent on $\mu(t)$ and the density function just falls only in the direction of the external force.
\section{Propagation kernel in the presence of a constant magnetic field}
Let us consider a charged particle with momentum $\mathbf{P}$ under the influence of a constant magnetic field in the direction of $z$ axis. The Hamiltonian is
\begin{equation}
\hat{H}=\frac{(\mathbf{P}-q\,\mathbf{A})^2}{2 m}\,=\,\frac{\hat{p}_x^2}{2 m}\,+\,\frac{\hat{p}_y^2}{2 m}\,+\,\frac{\hat{p}_z^2}{2 m}\,+\,\frac{1}{2}m \omega^2 \hat{x}^2\,+\,\frac{1}{2}m \omega^2 \hat{y}^2\,+\, \omega \hat{p}_x\hat{y}\,-\,\omega  \hat{p}_y\hat{x},
\end{equation}
where $\mathbf{P}=p_x \hat{i}\,+\,p_y \hat{j}\,+\,p_z \hat{k}$, $\omega=qB/2m$, and the vector potential is chosen as $\mathbf{A}=(-yB\hat{i}+xB\hat{j})/2$. From the Heisenberg equation of motion we find for the position and momentum operators
\begin{eqnarray}\label{posi}
\hat{x}(t) &=& \alpha_R \hat{x}(0)-\alpha_I \hat{y}(0)\,+\,\frac{\beta_R}{m}\hat{p}_x(0)\,-\,\frac{\beta_I}{m}\hat{p}_y(0),\nonumber\\
\hat{y}(t) &=& \alpha_R \hat{y}(0)+\alpha_I \hat{x}(0)\,+\,\frac{\beta_I}{m}\hat{p}_x(0)\,+\,\frac{\beta_R}{m}\hat{p}_y(0),\nonumber\\
\hat{z}(t) &=& \frac{t}{m}\hat{p}_z(0)\,+\,\hat{z}(0),
\end{eqnarray}
\begin{eqnarray}\label{momen}
\hat{p}_x(t) &=& -m\omega^2 \beta_R \hat{x}(0)+m\omega^2 \beta_I \hat{y}(0)+\alpha_R \hat{p}_x(0)-\alpha_I \hat{p}_y(0) ,\nonumber\\
\hat{p}_y(t) &=& -m\omega^2 \beta_R \hat{y}(0)-m\omega^2 \beta_I \hat{x}(0)+\alpha_R \hat{p}_y(0)+\alpha_I \hat{p}_x(0) ,\nonumber\\
\hat{p}_z(t) &=& \hat{p}_z(0),
\end{eqnarray}
where $\alpha_R=\cos^2(\omega t)$, $\alpha_I\,=\,-\sin(\omega t)\cos(\omega t)$, $\beta_R\,=\,\cos(\omega t)\sin(\omega t)/\omega$ and $\beta_I\,=\,-\sin^2(\omega t)/\omega$. The kernel or Feynman propagator in the position space is defined as $\langle\mathbf{r}|\hat{U}(t)|\mathbf{r'}\rangle\,=\, K(\mathbf{r},t|\mathbf{r'},0)$ and can be determined using Eqs. (\ref{posi},\ref{momen}) as
\begin{equation}
K(\mathbf{r},t|\mathbf{r'},0)\,=\,\Big(\frac{m}{2 \pi i \hbar}\Big)^\frac{3}{2} \frac{\omega}{\sin(\omega t)\sqrt{t}}\exp\Big(\frac{i m}{2 \hbar t}(z-z')^2\Big)\exp\Big(\frac{-i m \omega}{\hbar}(x y'-x' y)\Big)\exp\Big(\frac{i m \omega \cot(\omega t)}{2 \hbar}(x-x')^2 +(y-y')^2\Big).
\end{equation}
Having the kernel, the time-evolution of the initial state
\begin{equation}
\psi(x,y,z;0)=(\frac{\alpha}{\pi})^{\frac{1}{4}}\,\sqrt{\frac{2}{\omega_0^2 \pi 2^{m+n}\,n!m!}}\,\exp(-\frac{x^2+y^2}{\omega_0^2})\,\exp(-\frac{\alpha}{2} z^2 + i k_0 z)\,H_\textbf{n}(\frac{\sqrt{2}}{\omega_0} x)\,H_\textbf{m}(\frac{\sqrt{2}}{\omega_0} y),
\end{equation}
is obtained from $\psi(\mathbf{r},t)=\int d^{3}\mathbf{r'}\,k(\mathbf{r},t|\mathbf{r'},0)\, \psi(\mathbf{r'},0)$ as
\begin{multline}\label{70}
\psi(x,y,z;t)=(\frac{ m}{2\pi i\hbar})^\frac{3}{2}(\frac{\alpha}{\pi})^\frac{1}{4}\frac{2\omega\pi \Big(i(m\omega/\hbar)y- i(m\omega/\hbar)\cot(\omega t)x\Big)\Big(-i\beta x-i(m\omega/\hbar)\cot(\omega t)y\Big)}{\omega_0^3\sqrt{2t}\sin(\omega t)\sqrt{\frac{\alpha}{2}-\frac{im}{2\hbar t}}(\frac{1}{\omega_0^2}-\frac{i(m\omega/\hbar)\cot(\omega t)}{2})^3} \\
.\exp\Big(\frac{imz^2}{2\hbar t}+\frac{(ik_0 - \frac{imz}{\hbar t})^2}{2\alpha-\frac{2im}{\hbar t}}\Big) \exp\Big(\frac{(i(m\omega/\hbar)y -i(m\omega/\hbar) \cot(\omega t)x)^2}{(\frac{4}{\omega_0^2} -2i(m\omega/\hbar)\cot(\omega t))}\Big) \exp\Big(\frac{(-i(m\omega/\hbar) x -i(m\omega/\hbar)\cot(\omega t)y)^2}{(\frac{4}{\omega_0^2} -2i(m\omega/\hbar)\cot(\omega t))}\Big)\\
.\exp\Big(\frac{i(m\omega/\hbar)\cot(\omega t)}{2}(x^2+y^2)\Big).
\end{multline}
The probability density $\psi(x,y,z;t)\psi^{*}(x,y,z;t)$ has been depicted in FIG. 5. The applied magnetic field causes the particle to rotate in the direction of the field and gain orbital angular momentum.
\begin{figure}
\centering
\includegraphics[scale=0.4]
{F4.pdf}
\caption{(Color online) The density plot $|\psi(x,y,z;t)|^2_{z=0}$ of Eq. (\ref{70}), in terms of the scaled dimensionless variables $x/\omega_0$, $y/\omega_0$, $t/\tau$ where $\tau=m\omega_0^2/2\hbar$.}
\end{figure}
\section{charged particle in the presence of a constant magnetic field under the influence of external potential}
In this section we generalize the problem investigated in the previous section by adding a linear potential $\mu\hat{x}$ to the Hamiltonian. The linear potential corresponds to a constant force like the gravitational force and we assume that the force is along the $x$ axis. Now the Hamiltonian is
\begin{eqnarray}
\hat{H} &=& \frac{(\mathbf{P}-q\mathbf{A})^2}{2m}+\mu\hat{x},\nonumber\\
&=& \frac{\Big(\hat{p}_x + \frac{q B}{2}\hat{y}\Big)^2}{2m}+\frac{\Big(\hat{p}_y - \frac{q B}{2}\hat{x}\Big)^2}{2m} +\frac{\hat{p}_z^2}{2m}+\mu\hat{x},
\end{eqnarray}
where $\mu>0$ is the absolute value of the force along the $x$ axis.
Following the same process we did in the previous section, we find the position and momentum operators in Heisenberg picture as
\begin{eqnarray}
\hat{x}(t) &=& \alpha_R\hat{x}(0)-\alpha_I\hat{y}(0)+\frac{1}{m}\beta_R\hat{p}_x(0)-\frac{1}{m}\beta_I\hat{p}_y(0)-\frac{\mu}{m}\xi_R,\nonumber\\
\hat{y}(t) &=& \alpha_R\hat{y}(0)+\alpha_I\hat{x}(0)+\frac{1}{m}\beta_I\hat{p}_x(0)+\frac{1}{m}\beta_R\hat{p}_y(0)-\frac{\mu}{m}\xi_I,\nonumber\\
\hat{p}_x(t) &=& \alpha_R\hat{p}_x(0)-\alpha_I\hat{p}_y(0)\,-\,m\omega^2\beta_R\hat{x}(0)\,+\,m\omega^2\beta_I\hat{y}(0)\,-\,\mu\eta_R,\nonumber\\
\hat{p}_y(t) &=& \alpha_R\hat{p}_y(0)+\alpha_I\hat{p}_x(0)\,-\,m\omega^2\beta_I\hat{x}(0)\,-\,m\omega^2\beta_R\hat{y}(0)\,-\,\mu\eta_I.
\end{eqnarray}
where for simplicity we have defined
\begin{eqnarray}
&& \xi(t)\,=\,\int_{0}^{t}\beta(t')dt',\nonumber\\
&& \eta(t)\,=\,\int_{0}^{t} \alpha(t')dt',\nonumber\\
&& \alpha(t)\,=\,\alpha_R\,+\,i\alpha_I\,=\,\cos^2(t)\,-\,i\frac{\sin(2\omega t)}{2},\nonumber\\
&& \beta(t)\,=\,\frac{\sin(2\omega t)}{2\omega}\,-\,i\frac{\sin^2(\omega t)}{\omega}.
\end{eqnarray}
Having the explicit forms for position and momentum operators, the propagation kernel is obtained as
\begin{multline}
K(\mathbf{r},t|\mathbf{r'},0)=\Big(\frac{m}{2 \pi i \hbar}\Big)^\frac{3}{2} \frac{\omega}{\sin(\omega t)\sqrt{t}}\exp\Big(\frac{i m}{2 \hbar t}(z-z')^2\Big)\exp\Big(\frac{-i m \omega}{\hbar}(x y'-x' y)\Big)\exp\Big(\frac{i m \omega\cot(\omega t)}{2 \hbar}(x-x')^2 +(y-y')^2\Big)\\
.\exp\Big(\frac{-i \mu t}{2\hbar}(x+x')\Big)\exp\Big(\frac{-i \mu}{2\hbar}(\frac{1}{\omega}-t\cot(\omega t)(y-y'))\Big),
\end{multline}
and the initial state Eq. (\ref{2}) evolves to
\begin{multline}\label{state}
\psi(x,y,z;t)=(\frac{ m}{2\pi i\hbar})^\frac{3}{2} (\frac{\alpha}{\pi})^{\frac{1}{4}}\exp\Big(\frac{i(m\omega/\hbar)\cot(\omega t)}{2}(x^2+y^2)\Big)\exp\Big(\frac{(i(m\omega/\hbar) y-i(m\omega/\hbar)\cot(\omega t)x-\frac{i\mu t}{2\hbar})^2}{(\frac{4}{\omega_0^2}-2i(m\omega/\hbar)\cot(\omega t))}\Big)\\
.\exp\Big(-\frac{i\mu t}{2\hbar}x\Big)\frac{2 \pi\omega \Big(i(m\omega/\hbar) y-i(m\omega/\hbar)\cot(\omega t)x-\frac{i\mu t}{2\hbar}\Big)\Big(-i\beta x-i(m\omega/\hbar)\cot(\omega t)y-\frac{i\mu}{2\hbar}(\frac{1}{\omega}-t\cot(\omega t))\Big)}{\omega_0^3\sqrt{2t}\sin(\omega t)\sqrt{\frac{\alpha}{2}-\frac{im}{2\hbar t}}(\frac{1}{\omega_0^2}-\frac{i(m\omega/\hbar)\cot(\omega t)}{2})^3}\\
.\exp\Big(\frac{(-i(m\omega/\hbar)x-i(m\omega/\hbar)\cot(\omega t)y-\frac{i\mu}{2\hbar}(\frac{1}{\omega}-t\cot(\omega t)))^2}{(\frac{4}{\omega_0^2}-2i(m\omega/\hbar)\cot(\omega t))}\Big) \exp\Big(\frac{i\mu}{2\hbar}(\frac{1}{\omega}-t\cot(\omega t))y\Big)\exp\Big(\frac{imz^2}{2\hbar t}+\frac{(ik_0 - \frac{imz}{\hbar t})^2}{2\alpha-\frac{2im}{\hbar t}}\Big).
\end{multline}
The probability density of the state Eq. (\ref{state}) is depicted in FIG. 6. The applied magnetic field causes the particle to rotate in the direction of the field and gain an orbital angular momentum and under the influence of external potential, the beam falls down in the $x$ direction.
\begin{figure}
\centering
\includegraphics[scale=0.6]{F5.png}
\caption{(Color online) The density plot $|\psi(x,y,z;t)|^2_{z=0}$ of Eq. (\ref{70}), in the presence of a constant force in terms of the scaled dimensionless variables $x/\omega_0$, $y/\omega_0$, $t/\tau$ where $\tau=m\omega_0^2/2\hbar$.}
\end{figure}
\section{conclusion}
For a class of Hamiltonians describing the motion of a neutral or charged particle in the presence of a constant magnetic field and under the influence of a time-dependent external force, we found exact propagator kernels. The scheme was based on the Heisenberg equations of motion and kernel properties.  The initial wave function was a Gauss-Hermite wave function and for the evolved wave function we studied the uncertainties, orbital angular momentum, and the inertia tensor in the center of mass frame of the density function. From the point of view of interference of matter waves, the exact results obtained here give more exact results compared to the approximate methods like the axial approximation.  Our results can also be of importance in the non-relativistic quantum electron microscopy.\\
\\
\textbf{Data Availability Statement}\\
Data sharing not applicable to this article as no datasets were generated or analysed during the current study.