\section{Introduction}

Fundamental properties of dynamically evolving plasmas such as plasma channel size, electron and neutral density structure, and resistivity are challenging to measure \cite{Aragon2008}. Quantitative analysis of these basic characteristics is imperative to developing a robust theoretical framework of the dynamics involved in plasma formation and evolution. Common plasma imaging tools such as Schlieren imaging and interferometry are either not quantitative in nature, not directly sensitive to the phase, not spatially (in two dimensions) and temporally resolved, or restricted to weakly scattering plasmas \cite{Weber1997,Traldi2018,Biganzoli2015,Hayasaki2017,Aragon2008}. Furthermore, it is important to note that the hydrodynamics of the plasma and neutral atoms can be quite different. For example the plasma created by a laser or electrical discharge can create shock or sound waves that propagate faster than the plasma expansion.  Therefore, simultaneous imaging of each structure independently is advantageous. By imaging with two wavelengths, the spatial distribution of the electron and neutral background densities may be separated. Dual wavelength interferometry provides a quantitative alternative, but is often not spatially resolved and requires an external reference with interferometric stability \cite{Weber1997}. In this work, we apply recent advances in computational imaging toward improving plasma imaging techniques.

Coherent diffractive imaging (CDI) is a form of computational imaging that simultaneously retrieves the phase and amplitude of diffracting structures \cite{Fienup1978,Fienup1982c,Elser2003b}. In all CDI methods radiation, across the EM spectrum, probes a sample and diffraction is collected with a pixel detector \cite{Seaberg2011a,Porter2017a,Miao1999}. Propagation of the diffracted radiation is known and thus a phase problem is solved in order to reconstruct an image of the sample \cite{Thibault2008b}. While many forms of CDI exist, ptychography is a particularly robust implementation that allows simultaneous phase and amplitude contrast imaging of extended objects \cite{Rodenburg2007,Thibault2009c}.

Ptychography involves probing an object with a coherent illumination function over many partially overlapped regions of the specimen.  The redundancy of information from the partially overlapped regions yields a highly constrained phase problem. This allows ptychographic phase retrieval algorithms to produce reliable, high fidelity, noise-robust reconstructions \cite{Bunk2008}. It has been shown that mixed states can be deconvolved, and partial decoherence can be compensated for in advanced reconstruction algorithms, such as ptychographic information multiplexing (PIM) \cite{Thibault2008b,Batey2014,Thibault2013a}.  PIM extended ptychography to allow for simultaneous imaging with multiple wavelength probe illuminations without significant degradation to reconstruction quality as measured by resolution and signal to noise ratio (SNR) \cite{Wei2019}. While these recent developments in advanced ptychographic phase retrieval algorithms are promising for plasma imaging applications, they require transverse scanning of the object thus preventing the possibility of imaging transient, non-reproducible objects. Single-shot ptychography (SSP) overcomes this scanning requirement by breaking up a coherent illumination into smaller beamlets which simultaneously probe the object \cite{Sidorenko2016e,Chen2018,Wengrowicz2019a,Sidorenko2017}. The diffraction patterns from each beamlet are then collected on a single-detector where they are computationally separated to form a data-set similar to that taken in scanning ptychography. Scanning ptychographic phase-retrieval algorithms are used to reconstruct the object. SSP allows for time-resolved ptychographic imaging, making it an applicable imaging technique to probe plasma dynamics.

In this paper, we propose a novel imaging technique in which SSP is performed simultaneously with probes of multiple wavelength. We call this method wavelength-multiplexed single-shot ptychography (WM-SSP). In the methods section, improvements to SSP and the development of WM-SSP are presented. The results section shows a number of simulated and real experiments that convey the SSP improvements and demonstrate WM-SSP applied to both wavelength dependent and non-wavelength dependent objects. In the WM-SSP Plasma Imaging Simulation section a simulated experiment shows how WM-SSP can be applied to calculate the electron and neutral density structure in a plasma.

\section{Methods}

Before presenting our improvements to single-shot ptychography, we discuss the previously developed method \cite{Sidorenko2016e}. The SSP setup, shown in Figure \ref{SSPDiagram}, consists of a 4f imaging system whereby a diffractive optical element (DOE) is imaged to a detector (camera). The DOE, typically an array of pinholes, breaks up the incident beam into smaller beamlets. These beamlets are collimated and overlap at a point referred to as the crossover point. An object is placed some known distance, $\delta$, away from the crossover point. Beamlets diffract off the object and diffraction pattern intensities are collected around each pinhole image. By computationally segmenting the detector a ptychographic data set is collected in a single-shot. Beamlet locations on the object, analogous to scan positions in scanning ptychography, are calculated using the pinhole image on the detector, the focal length of the second lens, and $\delta$. The diffraction data and calculated probe locations are fed into a conventional scanning ptychographic phase retrieval algorithm, such as the extended  ptychographical  iterative  engine or ePIE algorithm, which simultaneously reconstructs quantitative phase and amplitude images of the object and probe \cite{Gerchberg1971,Maiden2009}.

\begin{figure}[htbp]
\centering
\includegraphics[width=\linewidth]{SSPDiagram.jpg}
\caption{In a) a 3D schematic shows the experimental setup with a red illumination. In b) a detailed diagram of the optical setup is shown.  In this optical setup a 4f imaging system is created using two lenses which images the diffractive optical element (DOE) to the detector. The object is placed some distance, $\delta$, away from the crossover point of the beamlets ($\delta$ can be varied to the desired level of probe overlap and diffraction pattern diversity).}
\label{SSPDiagram}
\end{figure}

\subsection{Improvements to Single-Shot Ptychography}

The development of SSP represented a major advance in single-shot CDI. While SSP creates a highly constrained phase problem for a monochromatic probe, information-multiplexing adds complexity that necessitates more sophisticated methods. Here we present a number of improvements to SSP that provide consistently high-fidelity reconstructions of the object and probe. These improvements facilitate information multiplexing in SSP via probes of different wavelength.

In SSP, diffraction pattern quality has a profound impact on the fidelity and resolution of the reconstructions. Diffraction pattern quality degrades because of factors such as aberrations, back reflections, and SNR. We have designed our SSP microscope with careful consideration of the data collection methods. To mitigate back reflections and spherical aberration we used high quality, achromatic, AR coated lenses as well as a windowless camera. Additional reflections, such as reflections between the object and detector, were removed using tilted ND filters. All of the experimental results in this paper were reconstructed from data taken with a high dynamic range (HDR) algorithm to improve fidelity. The trade off between HDR and single-pulse SSP is a slight reduction in fidelity \cite{Wengrowicz2019a}.

To date many SSP setups use a rectangular DOE, i.e., one that consists of pinholes arranged in a rectangular grid \cite{Sidorenko2016e,Chen2018,Wengrowicz2019a}. In SSP, the DOE pattern is analogous to the scanning pattern in scanning ptychography. It has been shown that a Fermat spiral scanning pattern in ptychography increases the reconstruction fidelity by preventing regularly repeating artifacts and creating more uniformly overlapped probes on the object \cite{Huang2014b}. Furthermore, it is known that the Fermat Spiral offers significantly better performance under conditions with imperfections in the data. Data imperfections are unavoidable in SSP, thus using a DOE with pinholes arranged in a Fermat spiral should be more effective than a rectangular grid DOE \cite{Goldberger2020a}. For our Fermat spiral DOE we considered three design parameters: pinhole diameter, number of pinholes, and pinhole spacing. Here we considered identical pinholes, although theoretically pinholes size and shape variation could be utilized to optimize systematic performance for specific imaging needs. Pinhole diameter was determined from the oversampling condition, given by $\sigma = \frac{\lambda f}{dX D}$ where $dX$ is the pixel size, $D$ is the probe diameter, $\lambda$ is the probe wavelength, and $f$ is the last lens focal length. For our DOE design, we used a pinhole diameter of 55 $\mu$m to give us an oversampling of 8.2. The number of pinholes determines the theoretical lateral resolution given by $Resolution = \frac{\lambda f}{dX N_{px}}$, where $N_{px}$ is the number of pixels in the chopped-out diffraction pattern. Decreasing the number of pinholes increases the chopped-out diffraction pattern size, thus increasing the reconstruction resolution. Therefore, a DOE with only a few pinholes seems favorable, however this must be balanced by the need for redundancy from overlapping probes in a ptychographic data set. Due to these competing effects, we empirically determined that a DOE with 40 pinholes provides a good compromise between resolution and data redundancy for our experimental setup. Finally, pinhole spacing was set such that diffraction signal filled the camera, which minimizes resolution and cross-talk, the coherent superposition of diffraction from adjacent pinholes on the detector. Cross-talk lowers the reconstruction fidelity by introducing mixing of diffraction patterns that is not accounted for in the reconstruction algorithm.

\begin{figure}[htbp]
\centering
\includegraphics[width=\linewidth]{voronoiTessellationImprovementFigure.jpg}
\caption{Images a) and b) show the square and Voronoi tessellations of diffraction data collected from a SSP system with a Fermat Spiral DOE. The amplitude of the reconstructions from the square and Voronoi tessellated Air Force Test Pattern (AFTP) diffraction data are shown in d) and e). Line-outs from the reconstructions are shown in d) and e), showing that the Voronoi tessellation produces higher fidelity images by removing artifacts and flattening the amplitude in the AFTP bars. Images c) and f) show the complex reconstructed probe and AFTP from a reconstruction with all improvements outlined in the methods section of this paper. The brightness of these reconstructions represent the amplitude and the color represents the phase. The linear phase was removed from the AFTP image caused from tilt of the object.}
\label{VTessellation}
\end{figure}

We process the data by tessellating the detector and segmenting (computationally "chopping out") each beamlet's diffraction pattern. With a rectangular DOE, a square tessellation leads to full use of detector space and no double counting (used in multiple chopped out grids) of diffraction data. Tessellating a Fermat Spiral DOE is more complicated as a square tessellation will lead to unused detector space and/or double counted diffraction data, as illustrated in image a) of figure \ref{VTessellation} . To address this complication, we use a Voronoi tessellation. For a given set of points, $p_{j}$, (the centers of each pinhole) the corresponding Voronoi cells are defined as the set of pixels whose distance from a given point in  $p_{j}$ is equal to or less than the distance from any other point in $p_{j}$. The Voronoi tessellation leads to full utilization of the detector and gives preferential treatment to diffraction closest to each pinhole, shown in image b) of figure \ref{VTessellation}.

Finally, we found that SSP is extremely sensitive to accurate \textit{a priori} knowledge of the distance between the object and the crossover points of the beamlets ($\delta$ as shown in figure 1) (more details in the results section). $\delta$  can be measured accurately using translation stages oriented axially, however this relies on an accurate knowledge of the crossover point's axial location in free space. This works for tangible objects, but is not always possible for transient objects. Therefore, we have developed a phase retrieval algorithm that self corrects for inaccurate $\delta$ information. Since $\delta$ globally scales the beamlets' position on the object, our $\delta$ correction algorithm is based on \cite{Guizar-Sicairos2008}. In our $\delta$ correction algorithm, we start by defining the following per iteration error metric:

\begin{equation}
    \epsilon =  \frac{1}{M N J} \sum_{f_x,f_y,j} [|\Tilde{\psi}_{j}(f_x,f_y)| - I_{j}(f_x,f_y)]^2
\label{errMet}
\end{equation}

where $\epsilon$ is the error, $M$ and $N$ are the number of pixels in the $X$ and $Y$ directions respectively, $J$ is the number of beamlets, $f_x$ and $f_y$ are the discrete spatial frequency coordinates in the detector plane, $j$ is the beamlet index, $\Tilde{\psi}_j$ is the Fourier transformed exit surface wave of the $j^{th}$ beamlet, and $I_j$ is the measured diffraction pattern intensities from the $j^{th}$ beamlet. Our algorithm solves for $\delta$ by finding the gradient of this error metric with respect to $\alpha$ where $\alpha = \frac{\delta}{f}$, as:

\begin{equation}
    \frac{d \epsilon}{d \alpha} = \frac{1}{M N J} \sum_{f_x,f_y,j} [|\Tilde{\psi}_{j}| - I_{j}]^2 Re{(\Tilde{\psi}_{j}^* \frac{d \psi}{d \alpha}})
\label{derOfErrMet}
\end{equation}

where $^*$ represents the conjugate. Here $\frac{d \psi}{d \alpha}$ is defined as:

\begin{equation}
\begin{split}
    \frac{d \psi}{d \alpha} = -2 \pi i  \mathcal{F} ( \mathcal{F}^{-1} ( \Tilde{O}(f_x,f_y) (f_x X_j+f_y Y_j) e^{-2\pi i(f_x X_j+f_y Y_j)})P(x,y))
\end{split}
\label{derOfPsi}
\end{equation}

where $\mathcal{F}$ and $\mathcal{F}^{-1}$ are the 2D-discrete Fourier transform and inverse Fourier transform operators, respectively, $\Tilde{P}$ is the Fourier transformed probe, $X_j$ and $Y_j$ are the beamlet positions, x and y are the discrete object space coordinates, and $O$ is the object. The minimization is then implemented using Newton’s method, and after every iteration of the phase retrieval algorithm delta, is updated in the following way:

\begin{equation}

     \delta_{n+1} = (\alpha_{n} + \beta \frac{d \epsilon}{d \alpha}) f
\label{newtMeth}
\end{equation}

 where $\beta$ is a step-size parameter that determines how much $\delta$ changes in each iteration. Use of this algorithm relaxes the required accuracy of the \textit{a priori} $\delta$ information. Through simulated experiments we found that the algorithm retrieved the simulated $\delta$, when the initial $\delta$ was with 5\% of the simulated $\delta$. The fidelity increase from a reconstruction with 5\% error in $\delta$ and 0\% error in $\delta$ is shown in figure \ref{deltaMinimization}.

\begin{figure}[htbp]
\centering
\includegraphics[width=\linewidth]{deltaMinimzationFigure.jpg}
\caption{Graph of Object RMS Error as a function of the Percent Error in $\delta$. This plot was made from doing reconstructions from the same simulated data set over a range of $\delta$'s. Example reconstructions of the cameraman object are shown from a reconstruction at 5\% error in $\delta$ and 0\% error in $\delta$. The upper-right corner shows the corrected $\delta$ as a function of iteration number from a reconstruction from simulated data. The $\delta$ correction algorithm started on iteration 25 (with a 5\% error in $\delta$) and by iteration 80, it had found the correct $\delta$.}
\label{deltaMinimization}
\end{figure}

\subsection{Development of wavelength multiplexed single-shot ptychography (WM-SSP)}

In WM-SSP, SSP is performed simultaneously with probes of multiple wavelengths. Instead of using a single-wavelength phase retrieval algorithm like ePIE, an information multiplexed algorithm, PIM, is used to reconstruct objects for each wavelength \cite{Batey2014}. To reinforce WM-SSP, we developed a novel probe constraint specific to the SSP geometry.

In our SSP geometry, pinholes arranged in a Fermat spiral pattern are used to break up the input illumination into multiple beamlets. The pinholes are small relative to the width of the input beam, such that they effectively act as a spatial filter for the beam. All beamlets for a given wavelength are thus the same and do not need to be treated as different in the algorithm (assuming uniform illumination of the DOE). Probes of different wavelengths, however diffract at different spatial rates. The width of the probe (defined by the zeroth order Jinc) is given by:

\begin{equation}
    W_{probe}= 1.22 \frac{\lambda f}{W_{pinhole}}
\label{wProbe}
\end{equation}

where $w_{probe}$ is the width of the probe, $W_{pinhole}$ is the width of the pinhole. The pixel size in object space, the basic unit of the detector's pixelated grid propagated to the object, is defined as:

\begin{equation}
    dX_{object}= \frac{\lambda f}{W_{detector}}
\label{ObjPx}
\end{equation}

where $dX_{object}$ is the object space pixel size, and $W_{detector}$ is the size of a chopped out diffraction pattern. Both the width of the probes and the object space pixel size scale proportional to $\lambda$. Therefore, in pixel space, probes of different wavelengths will be exactly the same, since the factors of $\lambda$ cancel. In every iteration of the phase retrieval algorithm, we can enforce this constraint by averaging the updated probes for each wavelength. The spectral weights i.e., power contained in each probe wavelength, can still be calculated in each iteration separately, but the spatial profile is constrained to be the same in pixel space for each wavelength. This probe constraint allows the phase retrieval algorithm to give high fidelity reconstructions for each wavelength.

\section{Results}

\subsection{Demonstration of Improvements to Single-Shot Ptychography}

In the previous section, we presented the use of a Fermat spiral DOE, a Voronoi tessellation to segment the detector, and a $\delta$ minimization algorithm to correct for inaccurate \textit{a priori} knowledge of $\delta$. Here we present simulated and experimental demonstrations that illustrates how these additions to the SSP technique improves the method and reconstruction fidelity.

The use of a Fermat Spiral probing pattern in scanning ptychography has been well studied and is known to increase the fidelity of reconstructions. As explained in the Methods section, this implies that using a DOE with a Fermat Spiral pinhole pattern will increase the reconstruction image quality. To show that a Voronoi tessellation, when using a Fermat Spiral DOE, increases reconstruction fidelity, we conducted an experiment using a real world setup. The 4f-imaging system in the setup consisted of two 5 cm air-spaced achromatic lenses, Thorlabs part number ACA254-050-A, with high quality AR coatings at the probing wavelengths. The probing laser was a 532 nm CW-diode laser, which was spatially filtered, expanded and collimated to create relatively uniform illumination on the DOE. The DOE, (a custom made photo-lithography mask: HTA Photomask) was made-up of 40 pinholes each with a 55 $\mu$m diameter arranged in a Fermat spiral. The detector (Thorlabs 8051M-USB) was an 8 mega-pixel camera with the wedge and face-plate removed. The camera sensor array has 3296x2472 square pixels with a pixel size, $dX$ of 5.5 $\mu$m. The DOE under-filled the sensor so that the diffraction filled around 50\% of the total detector space.

To compare square and Voronoi tessellations, an Air Force Test Pattern (AFTP) was imaged at $\delta = 0.9 cm$. Holding everything else constant, the collected diffraction intensities were segmented using both square and Voronoi tessellations. Since the Fermat spiral is not regularly spaced, we chose a square tessellation with cells of 300x300 pixels to balance the effects of double-counting diffraction data and not using enough diffraction data. The Voronoi tessellation will create cells of varying sizes that are chopped out and placed on a square cell of 444x444 pixels (padded with zeros). Reconstructions were performed using both tessellations for 1000 iterations, under the same reconstruction parameters. The amplitude of the retrieved images are shown in d) and e) of figure \ref{VTessellation}. The reconstruction from the square tessellation has artifacts, like the lines between the real lines of the AFTP, and low fidelity as seen in the non-uniform amplitude distribution within the bars of the AFTP. The Voronoi tessellated reconstruction removes these artifacts and has significantly better fidelity of the flat top AFTP bars. This can be seen more easily in the line outs shown in d) and e) of figure \ref{VTessellation}.

The fidelity of a SSP reconstruction is heavily dependent on accurate\textit{a priori} knowledge of $\delta$. To show this, we performed a simulated experiment where reconstructions were done over a range of $\delta$'s. The simulated data was collected with a detector made up of 2048x2048 square pixels of $dX = 10 \mu m$. The 4f-imaging system was made up of two 10 cm focal length lenses and the simulated probe wavelength was 532 nm. The simulated DOE consisted of 30 circular pinholes with 50 $\mu m$ diameters arranged in a Fermat Spiral spaced to fill the detector. The object imaged was the cameraman placed at a $\delta$ of 0.9 cm outside of the crossover point. The field was propagated between lenses using normalized FFTs and between the crossover point and the object plane using spectrum of plane waves. With the simulated diffraction patterns, reconstructions were done from 10\% below to 10\% above the simulated $\delta$. The root mean square (RMS) error was calculated between the reconstructed object at the simulated $\delta$ and the reconstructed object at every sampled $\delta$. A plot of the RMS error vs Percent Error in $\delta$ is shown in figure \ref{deltaMinimization}. The graph shows a drastic increase in object RMS error for small errors in $\delta$, which is reflected in the decrease in fidelity of the cameraman shown in the reconstructions at 5\% error in $\delta$. The $\delta$ correction algorithm, explained previously, was then tested on the simulated data used for this experiment. Starting with a 5\% error in $\delta$, the correction algorithm was able to converge on the correct $\delta$ after 55 iterations. A graph showing $\delta$ change as a function of reconstruction iteration is shown in the inset of figure \ref{deltaMinimization}. Our simulation shows that the $\delta$ correction algorithm was able to retrieve the correct $\delta$ when the guess was within 5\% of the correct $\delta$.

With all these improvements to SSP, an experimental reconstruction of an AFTP was done using the system previously described for the Voronoi Tessellation experiment. The complex probe and AFTP are shown in c) and f) of figure \ref{VTessellation} respectively. In these reconstructions, the brightness represents the amplitude of the object and the color represents the phase. We removed the linear phase on the object which was caused from a tilt of the object with respect to the optical axis. The quadratic phase that is present in the object reconstruction is believed to come from the angled probes present in SSP. These images show the high fidelity reconstructions that are consistently retrieved with our SSP microscope.

\subsection{Experimental Verification of wavelength-multiplexed single-shot ptychography: test pattern}

As an initial test of WM-SSP, we reconstructed an object with no wavelength dependence. The chosen object was an AFTP deposited on glass, which has a negligible amount of spatial dispersion from the beamlets incident on the glass at an angle. The data was collected with 532 nm and 633 nm CW-lasers simultaneously illuminating the experimental SSP system previously described. WM-SSP reconstructions were performed using an algorithm based on PIM with the novel probe constraint. Ghost modes were also used to produce higher fidelity reconstructions. A ghost mode is an additional mode for every wavelength that is summed incoherently at the detector. Every iteration of the phase retrieval algorithm, the ghost mode reconstructed object is flipped across the x-y axis and laterally shifted \cite{rana}. The reconstructed complex AFTP for the 532 nm and 633 nm probes are shown in figure \ref{MultiWavelengthDevFig}. The different field-of-view for the reconstructions is expected due to the difference in probe and object-space pixel sizes shown in equations \ref{wProbe} and \ref{ObjPx}.

\begin{figure}[htbp]
\centering
\includegraphics[width=\linewidth]{MultiWavelengthAFTPReconstructionsFigure.jpg}
\caption{The novel probe constraint is pictorially represented in a). The probe width and the object space pixel size vary linearly with the probe wavelength. The probes in real space units are shown in b), where the red and green probes are not overlapped, due to the wavelength dependence of the probe widths. The probes in pixel space are shown in c), where the red and green probes are spatially overlapped due to the wavelength dependencies of the probe width and object space pixel size cancelling. Complex Experimental WM-SSP AFTP reconstructions are shown in d) and e) (the amplitude is represented by the brightness while the phase is represented by the color). The difference in field of view is due to the wavelength dependence of the object space pixel size.}
\label{MultiWavelengthDevFig}
\end{figure}

\subsection{Experimental verification of wavelength-multiplexed single-shot ptychography: wavelength-dependent object}

We experimentally verified WM-SSP by reconstructing a wavelength-dependent object. The setup is shown in figure \ref{BirefringentMultiWavelength} and is the same as previously described with the addition of an analyzing polarizer. The wavelength-dependent object was made by exploiting the resulting stress-induced birefringence from femtosecond laser exposure of transparent substrates \cite{Champion2013}. The sample imaged is composed of a 4x4 grid of cylinders (imaged along the axis so they appear to look like circles) inscribed in a soda-lime substrate using a femtosecond laser. (Details about the fabrication methods for these structures can be found in \cite{McMillen2014,McMillen2015}). In the case of soda-lime and for the laser exposure conditions used here, the cylinders produce an isotropic stress distribution. Specifically, each cylinder has a diameter of 100 microns and consists itself of a stack of 128 laser-written layers. By simultaneously illuminating the sample with a horizontally polarized 532 nm probe and vertically polarized 633 nm probe and placing an analyzing polarizer after the object, birefringent and non-birefringent objects are imaged separately. The analyzing polarizer was aligned to the vertical polarization such that the 633 nm probe imaged the non-birefringent parts of the sample while the 532 nm probe imaged the birefringent object. Data was taken with both beams simultaneously illuminating the object and each beam individually imaging the object. We reconstructed the single-wavelength data sets, each running 5000 iterations of a standard ePIE algorithm. These reconstructions serve both as a base to compare the WM-SSP reconstructions against but also as object guesses for the WM-SSP retrieval algorithm to start with. Reconstructions of the dual-wavelength data set was done, running 10,000 iterations. The novel probe constraint and ghost modes were used in this reconstruction. The complex single-wavelength reconstructions and multi-wavelength reconstructed objects are shown in figure \ref{BirefringentMultiWavelength}. The 633 nm reconstructions, non-birefringent part of the objects, both show a 4X4 array of circles with dots in the middle. The 532 nm reconstructions, birefringent objects, are more complicated but share the same diagonal features going through the 4X4 array of laser machined circles (cylinders viewed from the top). Since the single-wavelength and multi-wavelength reconstructions compare favorably we are confident that WM-SSP works for wavelength dependant objects. This experiment also provides a novel method for both space and time-resolved quantitative analysis of stress-strain dynamics of laser machining process. This technique has the capability of providing laser machinist valuable information about the stress-induced birefringence they are creating in their machining processes.

\begin{figure*}[htbp]
\centering
\includegraphics[width=\linewidth]{MultiWavelengthBirefringentSample.jpg}
\caption{The WM-SSP setup for imaging the birefringent laser machined sample is pictorially represented. The red probe illumination is vertically polarized and the green probe illumination is horizontally polarized. By aligning an analyzing polarizer after the birefringent object, the red probe images the non-birefringent part of the sample while the green probe images the birefringent part of the sample. The complex object reconstructions from the single-wavelength data are shown in a) and b) for comparison against the complex reconstruction from the multi-wavelength data shown in c) and d) (intensity represents the amplitude while the color represents the phase). A white light image of the laser-written pattern is shown to the left of the reconstructions for reference. Comparing c) to a) and d) to b) shows that WM-SSP reproduces high fidelity reconstructions of wavelength-dependent objects.}
\label{BirefringentMultiWavelength}
\end{figure*}

\section{WM-SSP Plasma Imaging Simulation}

One exciting application of WM-SSP is plasma dynamics analysis. Plasma are created by ionizing gas, creating a mix of neutral unionized atoms, ions, and electrons. The neutral background atoms exhibit different hydrodynamics than the electronic structures and thus in observing the dynamics involved in plasma formation it is favorable to image each species separately. This is difficult due to their spatial overlap, non-reproducible dynamics, and fast time-scales. One way to observing the electron's and neutral's behavior independently is using dual-wavelength interferometry. Dual-wavelength interferometry utilizes the electrons' wavelength dependent index of refraction to separate the contribution to the phase from the electrons and neutrals \cite{Weber1997}. The phase shift from the electrons, relative to vacuum, is given by:

\begin{equation}
    \Delta \phi_e = - \frac{e^2 \lambda_0 \rho_e \Delta Z_e}{2 \pi \epsilon_0 m_e c^2}
\label{ePhase}
\end{equation}

where $e$ is the fundamental charge, $\epsilon_0$ is the vacuum permittivity, $m_e$ is the mass of an electron, $c$ is the speed of light, $\lambda_0$ is the probe wavelength, $\rho_e$is the electronic density, and $\Delta Z_e$ is the plasma thickness (assuming a uniform plasma in the axial direction). The phase shift from the neutral background, relative to vacuum, is given by:

\begin{equation}
    \Delta \phi_N = - \frac{4 \pi (n_0-1) \Delta Z_N}{\lambda_0}  (1 - \frac{\rho_N}{\rho_0})
\label{nPhase}
\end{equation}

where $n_0$ is the index of refraction of the neutral background, $\rho_N$ is the density of the neutral background, $\rho_0$ is the density of the gas the plasma is formed inside, and and $\Delta Z_N$ is the thickness of the uniformly distributed neutrals \cite{Weber1997}. Dual-wavelength interferometry experimentally determines the phase shifts for two wavelengths and uses equations \ref{ePhase} and \ref{nPhase} to setup a system of equations to solve for $\rho_e$ and $\rho_N$. Generally dual-wavelength interferometry (DWI) is not spatially resolved. Spatially (Two-dimensional) and temporally resolved versions of DWI, and generally most plasma imaging methods, are based on pump-probe setups which preclude them from studying non-reproducible plasma dynamics \cite{Hayasaki2017,Aragon2008}. Since WM-SSP numerically solves for the object phase at each wavelength, the same system of equations can be setup and solved to calculate the spatially- and temporally- dependent $\rho_e$ and $\rho_N$ from a single WM-SSP data set. To show the potential of this technique, we performed an experimental simulation.

The simulated WM-SSP setup included two 100 cm-focal length lenses, a detector with 2048x2048 square pixels of $dX = 50 \mu m$, and a DOE with 20 circular pinholes of $130 \mu m$ diameter arranged in a Fermat spiral. The plasma was simulated with $\rho_e = 2.65 * 10^{24}$ electrons/m$^3$, $n_0 = 1.0005$, $\rho_N = 2.386*10^{25}$ atoms/m$^3$, $\rho_0 = 2.65*10^{25}$ atoms/m$^3$, $\Delta Z_e = 10^-4$ m$^3$, and $\Delta Z_N = 10^-3$ m$^3$. These simulation parameters were chosen to represent a plasma being formed in air at STP with an ionization fraction of $10\%$.  The electron's spatial structure was made to resemble an 'e' (everything inside the e was uniformly distributed electrons and everything outside did not contributing to the electron phase shift). The neutral's structure was made to resemble an 'N' (everything inside the N was uniformly distributed neutral atoms and everything outside the N was uniformly distributed background atoms). The plasma was simulated with unit amplitude such that the complex transfer function was written as Object$_{plasma} = e^{i\Delta \phi_e (x,y) + i\Delta \phi_N (x,y)}$. The simulated electron structure, neutral structure, and plasma object are shown in figure \ref{PlasmaDeconFigure}. Diffraction data was simulated by placing the 2D plasma object at $\delta = 4$ cm and probing it simultaneously with 400 nm and 800 nm probe illumination. The plasma objects were reconstructed at each probe wavelength and the electron and neutral densities were calculated. The reconstruction went for 10,000 iterations, the simulated probes were used, the object amplitude for each wavelength was enforced to be one everywhere, and no ghost modes were used. The reconstructed plasma and calculated densities are shown in figure \ref{PlasmaDeconFigure}c and figure \ref{PlasmaDeconFigure}d respectively. The simulated experiment shows the high fidelity electron and neutral density structures that can be calculated from a single WM-SSP data set. The development of WM-SSP represents a major step in computationally studying the dynamics involved in plasma formation.

\begin{figure}[htbp]
\centering
\includegraphics[width=\linewidth]{PlasmaDeconvolutionFigure.jpg}
\caption{The simulated electron and neutral structure is shown in a). This plasma was probed by 400 nm and 800 nm coherent illumination in simulation, shown in b). A simulated WM-SSP diffraction data set was collected and reconstructed. The reconstructions are shown in c). The plasma's electron and neutral structure was calculated, shown in d), from the reconstructed object.}
\label{PlasmaDeconFigure}
\end{figure}

\section{Conclusion}

Through a number of improvements to SSP, we demonstrated wavelength multiplexing in single-shot ptychography. These improvements made SSP a more reliable method of single-shot CDI. Through the development of a novel probe constraint, in addition to the improvements made to SSP, wavelength-multiplexed single-shot ptychography was experimentally realized. By imaging a birefringent object with multiple wavelengths that were orthogonally polarized, we demonstrated our novel computational imaging method by imaging a wavelength dependent object. After reconstructing both wavelengths individually and simultaneously, the reconstructions of each wavelength favorably compared thus verifying that wavelength-multiplex SSP works. WM-SSP was applied to image a simulated plasma composed of both electrons and neutral atoms. Exploiting the electrons' wavelength dependent nature allowed us to calculate the plasma's electron and neutral density structures from a WM-SSP reconstruction. This imaging method is ready to be applied to imaging transient phenomena such as dynamically evolving plasmas. Wavelength-multiplexed SSP allows for quantitative (resolved in both 2D space and time) imaging of both the electron densities and the neutral densities within a single plasma event. Furthermore, this method has the capability of probing quantitatively for fundamental proprieties of plasma, making WM-SSP a powerful novel plasma imaging metrology.

\section*{Funding}
The authors gratefully acknowledge funding from the Air Force through AFOSR FA9550-18-1-0089 and Los Alamos National Laboratory through contract number 501188.

\section*{Disclosures}
The authors declare no conflicts of interest.

\bibliographystyle{ieeetr}


\section{INTRODUCTION}
{\it Introduction.} -- The $\alpha$-particle interaction with nuclei as well as the corresponding optical model potential (OMP) were of special interest from nuclear-physics earliest days.
The widely-used phenomenological OMP parameters have been derived from analysis of either elastic-scattering, which is ruled out below the Coulomb barrier $B$, or $\alpha$-induced reaction data.
They are then used to describe also the $\alpha$--emission from excited nuclei, in terms of the statistical Hauser-Feshbach (HF) \cite{wh52p} and pre-equilibrium emission (PE) \cite{eg92p} models.
However, there are also various assumptions and parameters of these models in addition to the $\alpha$-nucleus potential.
Thus, a definite conclusion on $\alpha$-potential may become possible only using HF+PE consistent parameter sets already validated by analysis of other independent data distinct from either $\alpha$-induced reaction or especially $\alpha$--emission data.

Moreover, there is a so-called $\alpha$-potential mystery of the account at once of both absorption and emission of $\alpha$-particles in nuclear reactions \cite{tr13p}, of equal interest for nuclear astrophysics and fusion technology.
It was referred by Rauscher \cite{tr13p} to competition of the Coulomb excitation (CE) with the compound nucleus (CN) formation in $\alpha$-induced reactions below $B$, while the former does not affect the $\alpha$--emission.
Because the corresponding partial waves and integration radii may provide evidence for distinct account of CE cross section and OM total-reaction cross section \cite{va16ap}, an alternate solution may also concern additional reaction channels but for $\alpha$--emission. Thus, further consideration should be given to eventual role of the pickup direct interaction (DI) leading to increase of the $\alpha$-emission beyond the HF+PE predictions.
Besides, a reaction cross-section enhancement that could be related to the position of a giant quadrupole resonance (GQR) is also worthy of note \cite{ma06p}.

In fact, a former search for new physics in potentials to describe nuclear de-excitation assumed that particle evaporation occurs from a transient nuclear stratosphere of the emitter nucleus (\cite{ma06p} and Refs. therein).
Recent definite conclusion on incident $\alpha$-particle OMP \cite{va14p} made feasible the analysis of a possible difference between the OMPs describing either $\alpha$-particle elastic scattering and reactions or $\alpha$-emission from excited nuclei \cite{va94p}.
On the other hand, OMP \cite{va14p} validation for $\alpha$-emission in low-energy proton-induced reactions on Zn isotopes was related to a surface character \cite{va15p}, at variance with the fast-neutron induced reactions on Zr isotopes \cite{va17p}.
Analysis of $\alpha$--emission from nuclei excited in reactions induced by either neutrons or low-energy protons became thus of interest while there are quite useful recent data of low-lying states feeding \cite{gz10ap,gz10p,yg14p,zw15p,tk18p,hb19p} for $A$$\sim$60 nuclei.

Consequently, it is of interest to find out also under these conditions if the $\alpha$-emission could be described by the same OMP \cite{va14p} which led to a rather good account of $\alpha$-induced reaction data \cite{va14p,ma09p,ma10p,va16p,va19p}.
It provided in the meantime better results also within large-scale nuclear-data evaluation \cite{ajk15p}, being adopted as default option within the world-wide used code TALYS \cite{TALYSp}.
Nevertheless, use of a consistent model parameter set, with no empirical rescaling factors of the $\gamma$ and/or neutron widths, is essential to avoid compensation effects.
Hence the main HF, PE, and DI assumptions and parameters as well as independent data analysis to support them are given in the Supplemental Material  \cite{SM2020p}.

\begin{figure*}
\resizebox{1.8\columnwidth}{!}{\includegraphics{Ni614pa.EPS}}
\caption{\label{Fig:Ni614pa} (Color online) Comparison of (a,b) measured $(p,\alpha)$  excitation functions on $^{61,64}$Ni \cite{exforp}, evaluated (short-dashed curves) \cite{TENDLp}, and PE+CN calculated values using the $\alpha$-particle OMPs of Refs. \cite{va14p} (dashed curves) and \cite{va94p} (short-dash-dotted curves), and (b) total PE $\alpha$-emission (dash-dot-dotted curve), DI pickup (dash-dotted curve), and sum (solid curve); uncertainty bands correspond to error bars of NLD parameter $a$ (with uncertainties in units of the last digit \cite{SM2020p}) and low-lying level cumulative number $N_d$ of residual nuclei (gray band), and (a) for GDH $\alpha$-particle pre-formation probability $\varphi$ (orange band).}
\end{figure*}

{\it DI pickup in addition to PE+CN $\alpha$-emission.} --
The due consideration of DI role for $\alpha$-emission in nucleon--induced reactions is obviously leading to an increase beyond the PE+HF predictions.
Actually, Qaim {\it et al.} \cite{smq95p} extended at incident energies $\leq$15 MeV the conclusion of '90s that the pickup instead of knockout has the main DI contribution to low-lying states in $(p,\alpha)$ and $(n,\alpha)$ reactions (\cite{eg92p,eg89p,eg91p} and Refs. therein).
However, Qaim {\it et al.} normalized the semimicroscopic pickup contribution to their $^{64}$N$(p,\alpha)$$^{61}$Co data at 15 MeV, with results that depend notably on the former PE+HF component \cite{SM2020p}.
Hence, an increased predictive power of the model results demands the primary use of  spectroscopic data.

In the present work, the pickup contribution to $(p,\alpha)$ and $(n,\alpha)$ reactions was determined within the distorted wave Born approximation (DWBA) method using the code FRESCO \cite{FRESCOp} and the same OMPs as within HF+PE analysis.
One--step reaction was considered through the pickup of $^3$H and $^3$He clusters, respectively.
Moreover, the "spectator model" \cite{smits76p,smits79p} was involved, where the two transferred either neutrons or protons in $(p,\alpha)$ and $(n,\alpha)$ reactions, respectively, are coupled to zero angular momentum acting as spectators, while the transferred orbital ($L$) and total ($J$) angular momenta are given by the third, unpaired nucleon of the transferred cluster.
The prior form distorted--wave transition amplitudes, and the finite--range interaction were considered, with more details in \cite{SM2020p}.
Finally, 19 states \cite{smits79p,plj78p,ensdfp}, until the excitation energy of $\sim$5 MeV, were involved for assessment of the pickup excitation function in Fig.~\ref{Fig:Ni614pa}(b).

First, within 5--7 MeV above the effective $(p,\alpha)$ reaction threshold, the DI pickup component is around half the CN one but an order of magnitude larger than PE contribution.
The major role of DI pickup on low-lying residual states, at lower incident energies, is thus proved.
The PE component corresponds to the Geometry-Dependent Hybrid (GDH) model \cite{mb83p}, generalized through inclusion of the angular-momentum and parity conservation \cite{ma90p} and the phenomenological knockout model of $\alpha$--emission based on a pre-formation probability $\varphi$ \cite{eg92p}.
It increases with the incident energy, becoming dominant at proton energies $>$15 MeV.

Second, the DI pickup inclusion is providing a suitable description of the measured data at all energies.
This point is particularly notable below 10 MeV, where there is no effect of nuclear level density (NLD) parameters \cite{SM2020p} as shown by the related uncertainty bands.

Third, there is an obvious data agreement, within several MeV above the effective $(p,\alpha)$ reaction threshold, for the CN component given by the $\alpha$-particle OMP of Ref. \cite{va94p}.
This fact may explain the former results found in this energy range but without DI account \cite{va94p}.

The case of the $(p,\alpha)$ reaction on $^{61}$Ni is also shown in Fig.~\ref{Fig:Ni614pa}(a) making use of the CN isotopic effect of cross--section decreasing with the isotope mass increase \cite{nim77p}.
The lower similar dependence of the DI processes leads to an increased sensitivity of the $(p,\alpha)$ reaction on heavier isotope $^{64}$Ni to the pickup contribution.
Moreover, the larger cross sections for this case underlines the difference between predictions of the $\alpha$-particle OMPs \cite{va14p,va94p}, and thus support the former one.

{\it DI pickup contribution in $\alpha$--emission spectra.} --
An additional check of DI pickup component of the $\alpha$--emission in neutron--induced reactions on $^{59}$Co is worthwhile, as follows.
First, there is a large data basis of total $\alpha$--emission as well as $(n,\alpha)$ reaction cross sections over three decades.
Second, $^{59}$Co may be considered as a benchmark because of its medium  asymmetry parameter $(N-Z)/A$ and corresponding isotopic effect.
Third, there are measured at the incident energy of $\sim$14 MeV valuable $\alpha$--emission angular distributions and spectra \cite{rf86p,smg96p,k99p} that are quite important for DI+PE analysis.

\begin{figure}
\resizebox{0.85\columnwidth}{!}{\includegraphics{Co59naS.EPS}}
\caption{\label{Fig:Co59naS} Comparison of measured $\alpha$-emission spectra from 14--14.5 MeV neutron--induced reactions on $^{59}$Co \cite{rf86p,smg96p,k99p}, and calculated  DI pickup (dash-dotted curve), PE (dashed curves), CN first-- (short-dash curve) and second--emission (dotted curve) components at 14.3 MeV, and their sum (solid curve).}
\end{figure}

Actually, only a tentative attempt was possible in this work, using neutron spectroscopic factors from the $^{55}$Mn$(d,p)^{56}$Mn reaction analysis \cite{NDS56p,comfortp} and the spectator protons pair as picked from 1$f_{7/2}$ subshell  \cite{SM2020p}.
Then, 27 excited states with well-known $J^{\pi}$ and transferred $L$ \cite{ensdfp,ripl3p}, until 2.088 MeV excitation energy, were considered in the pickup cross--section calculation.
The good agreement proved by the analysis of angular distributions for the $\alpha$-energy bins, within center-of-mass system (c.m.), of 6-10 MeV, 10-12 MeV, and 12-14 MeV (Fig. 26 of \cite{SM2020p}), is completed hereafter by $\alpha$-emission spectra and excitation function discussion.

\begin{figure}
\resizebox{0.85\columnwidth}{!}{\includegraphics{Co59na2.EPS}}
\caption{\label{Fig:Co59na} As Fig.~\ref{Fig:Ni614pa} but for $(n,\alpha)$ reaction on $^{59}$Co \cite{exforp} and uncertainty gray band for error bars of $s$-wave neutron-resonance spacing $D_0^{\it exp}$ and $N_d$ of residual nucleus $^{56}$Mn (Table I \cite{SM2020p}), and additionally calculated cross sections of $(n,n'\alpha)$ reaction (thin dashed curve) and total $\alpha$-emission (solid curve).}
\end{figure}

Nevertheless, the results in Fig.~\ref{Fig:Co59naS} support the present attempt of the DI pickup contribution on the lowest-lying states of the odd--odd residual nucleus $^{56}$Mn.
The data agreement could be increased by taking into account that data transformation from laboratory system (lab) to c.m. may provide, as shown for $^{56}$Fe \cite{rf84p}, a shift of up to 1 MeV to higher emission energies.

\begin{figure*}
\resizebox{1.8\columnwidth}{!}{\includegraphics{Fe546nas.EPS}}
\caption{\label{Fig:Fe546nas} As Fig.~\ref{Fig:Co59naS} but for $^{54,56}$Fe target nuclei \cite{exforp,smg79p}, at incident energies of 14.8 and 14.2 MeV, respectively.}
\end{figure*}

\begin{figure*}
\resizebox{1.8\columnwidth}{!}{\includegraphics{Fe546na.EPS}}
\caption{\label{Fig:Fe546na} As Fig.~\ref{Fig:Co59na} but for $^{54,56}$Fe target nuclei and additionally calculated $(n,\alpha)$ cross sections using TALYS-1.9 and its default parameters \cite{TALYSp} (short-dotted curves), and like--GQR components (dash-dot-dotted curves) added to DI+PE+CN (solid curves).}
\end{figure*}

At the same time, this spectrum analysis provides evidence for rather low value of the GDH pre-formation probability $\varphi$=0.04.
An uncertainty band corresponding to half of it is assumed in the excitation function analysis for both total $\alpha$--emission and $(n,\alpha)$ reaction cross section (Fig.~\ref{Fig:Co59na}).
The overlap of the related uncertainty band of calculated cross section and recent data is obvious.

On the other hand, Fig.~\ref{Fig:Co59na} proves also that there are no NLD effects below the incident energy of 10 MeV.
The $\alpha$--particle OMP \cite{va14p} is confirmed, with the additional support of more data since earlier work \cite{va94p}.

{\it Like-GQR in addition to DI+PE+CN $\alpha$-emission.} --
The DI pickup contribution to $(n,\alpha)$ reaction on $^{54,56}$Fe corresponds to 36 excited states up to 5.943 MeV for the residual nucleus $^{51}$Cr, and 18 states up to 5.557 MeV for $^{53}$Cr \cite{SM2020p}.
A first proof of the present approach was given by calculated angular distributions for the $\alpha$-energy bins from 6 to 14 MeV, compared to data of Fischer {\it et al.} \cite{rf84p} for $^{56}$Fe at 14.1 MeV (Fig. 23 of \cite{SM2020p}).
It is of interest in the following to note the best agreement obtained for the medium-energy $\alpha$--emission, i.e. from 10 to 12 MeV.

The most significant feature of the comparison of measured and calculated $\alpha$--emission spectra (Fig.~\ref{Fig:Fe546nas}) is the DI pickup contribution on lowest states.
Then, some spectra overestimation there is at the energies just below the DI maxima, at once with an underestimation of the low-energy side.
However it could be related for $^{54}$Fe to a rough attempt to transform data from lab \cite{smg79p} to c.m. by an average shift of 0.5 MeV to higher energies, at variance to the precise one for $^{56}$Fe \cite{rf84p}.
There are also known question marks at $\alpha$--particle energies below $\sim$6 MeV, pointed out by Fischer {\it et al.} \cite{rf84p}.
At the same time, the corresponding experimental total $\alpha$--emission cross sections are among the lowest data, at their incident energies, within the excitation functions (Fig.~\ref{Fig:Fe546na}).

Nevertheless, the DI pickup contribution is improving the description of these excitation functions except the energy range below $\sim$12 MeV.
The uncertainty bands related to number of low-lying levels, NLD, and PE parameters \cite{SM2020p} have already proved that none of their effects are in order at these energies.
Moreover, the above-mentioned best agreement obtained for the angular distribution for $\alpha$-energy bin of 6-10 MeV, at 14.1 MeV incident energy on $^{56}$Fe, is additionally supporting both the CN approach and the involved  $\alpha$--particle OMP \cite{va14p}.

Following a former suggestion \cite{ma06p}, we have assumed that such reaction cross-section enhancement could be understood as decay from giant resonances populated via neutron capture.
Although generally the decay of the GQR is observed with nucleon emission (\cite{
js81p} and Refs. therein), it has been shown (\cite{mf07p} and Refs. therein) that an appreciable non-statistical decay through $\alpha$--emission can occur.
This assumption has been supported by the position of the extra yield, beyond DI+PE+CN component sum, just at the GQR energy $E_{GQR}$=$65A^{-1/3}$ MeV \cite{js81p}, i.e. 17.092 and 16.89 MeV for $^{55,57}$Fe excited nuclei, respectively.
Therefore, we have obtained a fit of this extra yield by addition of Gaussian distributions at $E_{GQR}$, with FWHM widths of 2.35 and 3.54 MeV, and peak cross sections of 8 and 10 mb, respectively (Fig.~\ref{Fig:Fe546na}).

It is obvious that the widths of these Gaussian distributions are much lower than the systematic 'best' value $\Gamma$=$85A^{-2/3}$ MeV (Fig. 14 of Ref. \cite{js81p}).
They are even lower than the $\Gamma$=$17.5A^{-1/3}$ MeV dependence that seems to describe better the inelastic hadron scattering data which were obtained taking into account also other GRs \cite{js81p}.
Hence we would call these components only like--GQR contributions.
Further conclusions on the physics behind this empirical addition are not yet  evident, while more similar cases to be concerned may help.
Nevertheless, a suitable account of the measured $(n,\alpha)$ reaction and $\alpha$-emission cross sections (Fig.~\ref{Fig:Fe546na}) is obtained by this additional like--GQR contribution which are larger than the DI pickup for incident energies $<$12 MeV.

{\it Conclusions.} --
A consistent set of model parameters, validated by analysis of various independent data, makes possible the assessment of an optical model potential \cite{va14p} also for nucleon-induced $\alpha$-emission for $A$$\sim$60 nuclei.
Particularly,  $\alpha$--emission from $^{55,57,58}$Fe and $^{60}$Co isotopes excited by $(n,\alpha)$ reaction, and $^{62,64,65,66}$Cu and $^{64,65,66,68}$Zn isotopes excited through both $(n,\alpha)$ and $(p,\alpha)$ reactions has been analyzed.
The advantage of rather recent data of low-lying states feeding is essential.
Further consideration of additional reaction channels leading to increase of the $\alpha$-emission cross sections beyond the statistical predictions has concerned the DI pickup and GQR similar features.\\

This work was partly supported by Autoritatea Nationala pentru Cercetare Stiintifica (Project PN-19060102) and carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014-2018 and 2019-2020 under grant agreement No 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission.


\section{Introduction}\label{in}

Schur rings (or $S$-rings) were originally developed by Schur \cite{S33} and Wielandt \cite{W49} in order to study permuation groups as an alternative to character theory and have since been used in many other applications. Schur rings have been widely unsed in the study of association schemes, strongly regular graphs, and different sets.  For the applications of Schur rings to association schemes, see for example Sec. 2.4 and Sec. 4.4 of \cite{CP}.  \medskip

Classification of Shcur rings began with Schur himself, when Schur conjectured that all Schur rings over a group $H$ conincides with a so-called {\it transitivity module,} a submodule of the group ring $\mathbb{Z}[H]$ constructed by a group $G$ acting on its regular subgroup $H$. Such Schur rings are said to be {\it Schurian}. Wielandt \cite{W64} first showed that not all Schur rings are Schurian. Leung and Man in \cite{LM96} and \cite{LM98} showed that all Schur rings over a finite cyclic groups fall into one of four families: partitions induced by orbits of automorphism subgroups ({\it orbit Schur rings}), partitions induced by a group factorization ({\it tensor products}), partitions induced by cosets ({\it wedge products}), and the {\it trivial Schur ring}. We say that a Schur ring is {\it trational } if it is one of these four types.  Thus, the main results of Leung and Man tell us that all Schur rings over finite cyclic groups are traditional. Note that the above mentioned non-Schurian example by Wielandt is non-traditional; see the last paragraph of the fourth page of  \cite{M14}.  \medskip

Rcently, Bastian et al. \cite{B19}  study the Schur rings over infinite groups.  Schur rings over the integers, over a virtually infinite cyclic group, and over torsion-free locally cyclic groups are classified. All these Schur rings are traditional. However, there are non-triditional Schur rings over free groups and free products. \medskip

Let $\mathcal {Z}=\langle z\rangle$ denote the infinite cyclic group, written multiplicative. The finite cyclic group of order $n$ will be denoted by $\mathcal {Z}_n=\langle a\rangle$. In the present paper, all Schur rings over the direct product $\mathcal {Z}$ and $\mathcal {Z}_3$ will be classified. The following is the main result.

\thrml{1451b}All schur rings over the group $\mathcal{Z} \times \mathcal{Z}_{3}$ are of one of the following forms:
\begin{enumerate}
\tm{i}  $F[\mathcal {Z}_3] \wedge F[\mathcal {Z}]$ or $ F[\mathcal {Z}_3] \wedge F[\mathcal {Z}]^{\pm}$,

\tm{ii} $F[H \times \mathcal {Z}_3] \wedge F[\mathcal {Z}]$ or $F[H \times \mathcal {Z}_3]^{\pm} \wedge F[\mathcal {Z}]^{\pm}$,

\tm{iii} $F[K \times \mathcal {Z}_3] \wedge F[\mathcal {Z}]$ or $F[K \times \mathcal {Z}_3]^{\pm} \wedge F[\mathcal {Z}]^{\pm}$,

\tm{iv} $F[\mathcal {Z} \times \mathcal {Z}_3]$ or $F[\mathcal {Z} \times \mathcal {Z}_3]^{\pm}$,

\tm{v} an orbit Schur ring.
\end{enumerate}
where $H \leq  K\leq\mathcal {Z}$.
Hence, all Schur rings over $\mathcal {Z} \times \mathcal {Z}_3$ are either orbit or wedge product Schur
rings, which implies that the group $\mathcal {Z}\times \mathcal {Z}_3$ is traditional.
\end{theorem}

Throughout let $G$ denote an arbitrary group with identity element $1$ and $F$ an arbitrary field of characteristic zero, which contains $\mathbb{Z}$ as a subring. The group algebra over $G$ with coefficients in $F$ is denoted $F[G]$. The torsion subgroup of $G$ is denoted $T(G)$.
\section{Preliminary}\label{Pr}

 For any finite non-empty subset $C\subseteq G$, denote $\overline{C}=\sum_{g\in C}g$ and call $\overline{C}$ a {\it simple quantity}. In addition, we define $C^*=\{g^{-1}: g\in C\}$. A partition ${\mathcal P}$ of $G$ is said to be  {\it finite support} provided any $C\in {\mathcal P}$ is a finite subset of $G$. Our notation for Schur ring is taken from \cite{B19} and \cite{P15}.

 \subsection{Schur Rings}
\dfntnl{1055a} Let ${\mathcal P}$ be a finite support partition of $G$ and $\mathcal {A}$ a subspace of $F[G]$ spanned by $\{\overline{C}: C\in {\mathcal P}\}$. We say that $\mathcal {A}$ is a {\it Schur ring} (or $S$-ring) over $G$ if
 \begin{enumerate}
  \tm{1} $\{1\} \in {\mathcal P}$,
  \tm{2} for any $C\in {\mathcal P}$, $C^*\in {\mathcal P}$,
  \tm{3} for all $C, D\in {\mathcal P}$, $\overline{C}\cdot\overline{D}=\sum_{E\in {\mathcal P}}\lambda_{CDE}\overline{E}$, where all but finitely many $\lambda_{CDE}$ equal $0$.
 \end{enumerate}
\end{definition}

For an $S$-ring $\mathcal {A}$ over $G$, the associated partition ${\mathcal P}$ is denoted ${\mathcal D}(\mathcal {A})$ and each element in ${\mathcal D}(S)$ is called a {\it basic  set}
of $\mathcal {A}$. We say that a subset $C$ of $G$ is an {\it $S$-set} if $C$ is a union of some basic sets of $\mathcal {A}$. An {\it $S$-subgroup} is simutaneously an $S$-set and a subgroup.\medskip

Suppose $\alpha=\sum_{g\in G}\alpha_gg$ and $\beta=\sum_{g\in G}\beta_gg\in F[G]$ (note that here only finitely many nonzero coefficients $\alpha_g$ and $\beta_g$). Then define
$\alpha^*=\sum_{g\in G}\alpha_gg^{-1}$ and $\supp(\alpha)=\{g: \alpha_g\ne 0 \}$. Additionally,  the {\it Hadamard product} is defined as $\alpha\circ \beta=\sum_{g\in G}\alpha_g\beta_gg$. \medskip

The following theorem was proved by Wielandt in \cite{W64} when $G$ is finite. The general case appeared as \cite[Corollary 2.12]{B19}.
\thrml{1505b} Let $G$ be a group and $\mathcal {A}$ a subring of $F[G]$. Then $\mathcal {A}$ is a Schur ring if and only if $\mathcal {A}$ is closed under $\circ$ and $*$, $1\in \mathcal {A}$, and for all $g\in G$ there exists some $\alpha\in \mathcal {A}$ such that $g\in \supp(\alpha)$.
\end{theorem}

Let $f: F\rightarrow F$ be a function such that $f(0)=0$. For any $\alpha=\sum_{g\in G}\alpha_gg$, we set
$$
f(\alpha)=\sum_{g\in G}f(\alpha_g)g.
$$

The following proposition was proved by Wielandt in \cite[Proposition 22.3]{W64}. The general case appeared as \cite[Propostion 2.4]{B19}.
\prpstnl{913a}Let $\mathcal {A}$ be a Schur ring over a group $G$ and the function be as above. Then
$f(\alpha)\in \mathcal {A}$ whenever $\alpha\in \mathcal {A}$.
\end{proposition}

Let $f$ take value $1$ at one non-zero number and $0$ at other numbers.  We get the following statement, which is known as the Schur-Wielandt principle; see \cite[Corollary 1.10]{P15}.

\crllrl{1013c} Let $\mathcal {A}$ be an $S$-ring over $G$. For any $\alpha\in \mathcal {A}$, suppose $g\in \supp({\alpha})$ with $\alpha_g=c$. The set
$$
\{g\in G: \alpha_g=c\}
$$
is an $S$-set.
\end{corollary}

The following two propositions, which were first proved by Wielandt in \cite{W64} and were generalized in \cite{B19},  tell us how to generate $S$-subgroups in any $S$-ring.

\prpstnl{956b} Let $\mathcal {A}$ be a $S$-ring over $G$. Let $\alpha\in \mathcal {A}$ and
$$
{\rm Stab}(\alpha)=\{g\in G: \alpha g=\alpha\}.
$$
Then ${\rm Stab}(\alpha)$ is an $S$-subgroup of $\mathcal {A}$.
\end{proposition}

\prpstnl{956b} Let $\mathcal {A}$ be an $S$-ring over $G$. Let $\alpha\in \mathcal {A}$ and $H=\langle \supp(\alpha)\rangle$. Then~$H$ is an $S$-subgroup of $\mathcal {A}$.

\end{proposition}

Note that under the condition of Proposition \eqref{956b}, $\{C: C\in {\mathcal D}(\mathcal {A}), C\subseteq H\}$
consists of  basic sets of an $S$-ring over $H$. It is denoted $\mathcal {A}_H$.

When we have normal $S$-subgroup, we can construct $S$-ring over the factor group as shown in the following lemma.

 \lmml{1440a}(\cite[Lemma 1.2]{LM90}) Let $\varphi: G\rightarrow H$ be a group homomorphism with ${\rm Ker}(\varphi)=K$ and $\mathcal {A}$ be a Schur ring over $G$. Suppose that $K$ is an $S$-subgroup. Then the image $\varphi(\mathcal {A})$ is a Schur ring over $\varphi(G)$ where ${\mathcal D}(\varphi(\mathcal {A}))=\{\varphi(C): C\in {\mathcal D}(S)\}$.
\end{lemma}

In particular, if $H=G/K$, the factor group of $G$ over $K$, and $\varphi$ is the natural homomorphism,  the corresponding $S$-ring is denoted $\mathcal {A}_{G/K}$.

\medskip

\subsection{Traditional Schur Rings}

The group ring $F[G]$ itself forms a $S$-ring over $G$ with basic set $\{ \{g\}: g\in G\}$. This is called the {\it discrete Schur ring} over $G$. When $G$ is a finite group, the partiotion $\{1, G\setminus 1\}$ produces a Schur ring known as the {\it trivial Schur ring} over $G$. Note that if $G$ is infinite there does not exist trivial $S$-ring over~$G$. \medskip

Assume that $G=H\times K$ and $\mathcal {A}$ and ${\mathcal B}$ are $S$-rings over $H$ and $K$, respectively. Then the set
$$
\{CD: C\in {\mathcal D}(\mathcal {A}), D\in {\mathcal D}({\mathcal B})\}
$$
consists of a basic set  of an $S$-ring over $G$. The corresponding $S$-ring is called the {\it tensor  product} of $\mathcal {A}$ and ${\mathcal B}$. \medskip

If ${\mathcal K}$ is a finite subgroup of $\aut(G)$, the set of elements of $F[G]$ fixed by ${\mathcal K}$ is an $S$-ring over $G$, denoted $F[G]^{{\mathcal K}}$ and called the {\it orbit Schur ring} associated with ${\mathcal K}$. Note that the discrete $S$-ring can be seen as an orbit $S$-ring by taking ${\mathcal K}=1$.  \medskip

An $S$-ring $\mathcal {A}$ over $G$ is a {\it wedge product} if there exist nontrivial proper $S$-subgroups~$H,K$ such that $K\le H$, $K\unlhd G$, and every basic set outside  $H$ is a union of $K$-cosets. In this case, the series
$$
1<K\le H<G
$$
is called a {\it wedge-decomposition} of $\mathcal {A}$. Furthermore, we write $\mathcal {A}=\mathcal {A}_1\wedge \mathcal {A}_2$, where
$\mathcal {A}_1=\mathcal {A}_H$ and $\mathcal {A}_2=\mathcal {A}_{G/K}$. \medskip

An $S$-ring over a group $G$ is called {\it traditional} if it is either a trivial Schur ring (when $G$ is finite), or a tensor product of Schur rings over smaller subgroups, and or a wedge product.  A group $G$ is called {\it traditional} if each $S$-ring over it is traditional.

Several classes of groups are proved to be traditional as shown by the following theorem.
\thrml{1556c} The following classes of groups are traditional:
\begin{enumerate} \tm{1}[\cite{LM96}, \cite{LM98}, \cite{B19}] cyclic groups.
\tm{2}[\cite{B19}]$\mathcal {Z}\times \mathcal {Z}_2$.
\tm{3}[\cite{B19}] Ttorsion-free locally cyclic groups.
\end{enumerate}
\end{theorem}

Actually, it is proved that the only Schur ring over $\mathcal {Z}$ are discrete and symmetric Schur rings; see \cite[Theorem 3.3]{B19}. The symmetric Schur ring over $\mathcal {Z}$ is denoted $F[\mathcal {Z}]^{\pm}$.

\section{Schur rings over $\mathcal {Z}\times \mathcal {Z}_3$}

For $n\in \mathbb{Z}$, define the {\it $n$th Frobenius map} of $F[G]$ as
	$$
	\alpha=\sum_{g\in G}\alpha_gg\, \mapsto \, \alpha^{(n)}=\sum_{g\in G}\alpha_gg^n.
	$$

 \lmml{1951b}(\cite[Theorem 23.9]{W64}, \cite[Theorem 2.20]{B19}) Suppose that $\mathcal {A}$ is a Schur ring over $G$ where $G$ is an abelian group. Let  $m$ be an integer coprime to the orders of all torsion elements of $G$. Then for all $\alpha\in \mathcal {A}$ we have $\alpha^{(m)}\in \mathcal {A}$. In particular, if~$G$ is torsion-free, then all Schur rings over $G$ are closed under all Frobinius maps.
 \end{lemma}

 \lmml{1958c}(\cite[Lemma 2.21]{B19})Let $G$ be an abelian group such that the torsion subgroup $T(G)$ has finite exponent. Let $\mathcal {A}$ be a Schur ring over $G$. Then $T(G)$ is an $S$-subgroup.
 \end{lemma}

The following result is known as the second Schur theorem on multipliers.

\thrml{1535b} Let $\mathcal {A}$ be an S-ring over an abelian group $G$, $p$ a prime divisor of~$|T(G)|$, and $E=\{g \in G: g^{p}=e\}$.  Then for every finite $S$ -set $X$ the set
$$
X^{[p]}=\{x^{p}: x \in X,|X \cap E x| \not\equiv 0 (\bmod p)\}
$$
is an $S$-set.
\end{theorem}
{\bf Proof}.\
For a finite $S$-set $X$, one can see that

\qtnl{1042a}
\overline{X}^{p}=\left(\sum_{x \in X} x\right)^{p} \equiv \sum_{x \in X} x^{p}=\overline{{X}^{(p)}} \, (\bmod p),
\end{equation}

where we write
$$
\sum_{g} a_{g} g \equiv \sum_{g} a_{g}^{\prime} g\, (\bmod p)
$$
 if $a_{g} \equiv a_{g}^{\prime}(\bmod p)$ for all $g\in G$.

 On the other hand, given a coset $ C \in G / E,$ the element $g^p$ does not depend on the choice of $g \in C$.  Denote it by $h_{C}$.  Then the mapping $C \mapsto h_{C}$ is a bijection from $G / E$ onto $G^{p}$.  So
\qtnl{1043b}
\overline{{X}^{(p)}}=\sum_{C \in G / E}|C \cap X| \cdot h_{C} \equiv \sum_{C \in G / E,|C \cap X| \not \equiv 0 \, (\bmod p)}|C \cap X| \cdot h_{C} \, (\bmod p).
\end{equation}
Let $f: F \rightarrow F$ be the function such that $f(x)=1$ if $x$ is an integer not divisible by $p$; $f(x)=0$ for other $x\in F$. The formulas \eqref{1042a} and \eqref{1043b} imply that the set
$$
\left\{h_{C}: C \in G / E,|C \cap X| \not \equiv 0(\bmod p)\right\}
$$
is an $S$-set. Since this set equals $X^{[p]}$,  we are done.
\hfill$\square$

\medskip

The following two propositions are useful in analyzing the structure of $S$-rings over $\mathcal {Z}\times \mathcal {Z}_p$, with $p$ an odd prime.

\prpstnl{1501a}Let $\mathcal {A}$ be an S-ring over $\mathcal {Z} \times \mathcal {Z}_p=\langle z\rangle\times \langle a\rangle$, and let $H$ be the unique maximal $\mathcal {A}$-subgroup contained in $\mathcal {Z}$. If there exists $z^ma^i \in X\in {\mathcal D}(\mathcal {A})$ with~$|X|<p$, then $z^{pm}\in H$.
\end{proposition}
{\bf Proof}.\  Using the notation in Theorem \ref{1535b}, we have $X^{[p]}=X^{(p)}$ by the assumption. Furthermore, it is an $S$-set and hence $\langle X^{(p)}\rangle$ is an $S$-subgroup. According to the maximality of $H$, one can see that
$$
(z^ma^i)^p\in X^{(p)}\subseteq H,
$$
as required.
\hfill$\square$

\medskip

The following basic fact will be used freely in the proof of the main theorem.

{\bf Remark.}\, {\it Let $\mathcal {A}$ be an $S$-ring over $G$ and $K$ a normal  $S$-subgroup. Let $\varphi$ be the natural homomorphism from  $G$ onto $G/K$. For any $X\in {\mathcal D}(\mathcal {A})$, $\varphi^{-1}(X)$ is an~$S$-set since it is equal to $XK$.}

\medskip

\thrml{1521b} Let $\mathcal {A}$ be an S-ring over $\mathcal {Z} \times \mathcal {Z}_p=\langle z\rangle\times \langle a\rangle$, where $p$ is an odd prime. Then for all $X\in {\mathcal D}(\mathcal {A})$ either $X=z^m\mathcal{Z}_{p}$ or $p\neq |X|$.

\end{theorem}
{\bf Proof}.\
Denote $G=\mathcal {Z} \times \mathcal {Z}_p=\langle z\rangle \times\langle a\rangle$. Notice by Lemma \ref{1958c} that $T(G)=\mathcal {Z}_p$ is an $S$-subgroup.

According to Lemma~\ref{1440a}, $\mathcal {A}_{G/T(G)}$ is a Schur ring over $G/T(G)$.
Moreover,~$\mathcal {A}_{G/T(G)}$ is either discrete or symmetric.
By the above remark, it follows that
$$
X\subseteq z^m\mathcal{Z}_{p} \:\:{\rm or}\:\: X\subseteq \{z^m,z^{-m}\}\mathcal{Z}_{p},
$$
for all $X\in {\mathcal D}(\mathcal {A})$, where $z^m\in \mathcal {Z}$ satisfies $X\cap z^m\mathcal {Z}_p\neq \varnothing$.

Towards a contradiction, assume that  there exists  $X\in {\mathcal D}(\mathcal {A})$ with  $|X|=p$ and $X\neq z^m\mathcal {Z}_{p}$ for all $m\in \mathbb{Z}$.  It follows that  $X\subseteq \{z^m,z^{-m}\}\mathcal {Z}_{p}$ for some $m\in \mathbb{Z}$ such that

\qtnl{1039b}
X\cap z^m\mathcal {Z}_{p}\neq \varnothing \quad\text{and}\quad X\cap z^{-m}\mathcal {Z}_{p}\neq \varnothing.
\end{equation}

So, we may assume that
$$
X=\{z^ma^{i_1},z^ma^{i_2}, \ldots,z^ma^{i_t},z^{-m}a^{j_1},z^{-m}a^{j_2},\ldots,z^{-m}a^{j_k}\}.
$$
 Here, ${i_1},\ldots,{i_t},{j_1},\ldots,{j_k}\in \mathbb{N}$ satisfying

\qtnl{1612b}
t+k=p \quad\text{and}\quad  t\neq 0\neq k.
\end{equation}

 It follows that
  $$\overline{X}(\overline{T(G)}-1)= t\overline{z^mT(G)}+k\overline{z^{-m}T(G)}-\overline{X}\in \mathcal {A}.
$$
This implies that  $t\overline{z^mT(G)}+k\overline{z^{-m}T(G)}\in \mathcal {A}$ and hence
$$
t\overline{\{z^m,z^{-m}\}T(G)}+(k-t)\overline{z^{-m}T(G)}\in \mathcal {A}.
$$
By \eqref{1612b}, $k\neq t$ as $p$ is odd. Therefore,  $(k-t)\overline{z^{-m}T(G)}\in \mathcal {A}$ yields that
$\overline{z^{-m}T(G)}\in$~ $\mathcal {A}$.  As a result,  $X\subseteq z^{-m}T(G)$. This is a contradiction  to \eqref{1039b}.
\hfill$\square$

\medskip

{\bf Proof of Theorem \ref{1451b}}\, In the sequel, we fix the following notation:
$$
G=\mathcal {Z} \times \mathcal {Z}_{3}=\langle z\rangle \times\langle a\rangle,
$$
$\mathcal {A}$ is  a Schur ring over $G$, and $H$ is  the unique maximal~$S$-subgroup contained in~$\mathcal {Z}$.

\medskip

Note that $T(G)=\{1, a,a^2\}=\mathcal {Z}_{3}$ is an $S$-subgroup by Lemma \ref{1958c}. Let
 $$
 \pi: \mathcal {Z} \times \mathcal {Z}_{3} \rightarrow \mathcal {Z}
 $$
 be the natural projection map, with  $\ker (\pi)=\mathcal {Z}_{3}$.  Then $\pi(\mathcal {A})$ is a Schur ring over~$\mathcal {Z}$,  which  is either discrete or symmetric by \cite[Theorem 3.3]{B19}.

 \medskip

{\bf Case 1.} \,  Suppose $\pi(\mathcal {A})=F[\mathcal {Z}]$.  Then $\mathcal {A}_H$ is discrete. Furthermore,
 $$
 \left\{z^{m}, a z^{m}, a^2 z^{m}\right\}=\pi^{-1}\left(\left\{z^{m}\right\}\right)
 $$
 is an $S$-set for all integers $m$. Hence
 $$
 \mathcal {A}_{H\times \langle a\rangle}=\mathcal {A}_H\otimes \mathcal{A}_{ \langle a\rangle}.
 $$
  If $H=\mathcal {Z}$,  then $\mathcal {A}=\mathcal {A}_{\mathcal {Z}}\otimes \mathcal {A}_{ \langle a\rangle}$ and hence the $S$-ring $\mathcal {A}$ is traditional.

  \medskip

 Next, we assume that $H<\mathcal {Z}$.
 Choose  $z^{m} \in \mathcal {Z} \setminus H$.   Then
 $$
 \{z^{m}, a z^{m}\}\in {\mathcal D}(\mathcal {A}), \, \, {\rm  or} \, \,  \{z^{m}, a^2 z^{m}\}\in {\mathcal D}(\mathcal {A}), \, \, {\rm  or} \, \, \{z^{m}, a z^{m}, a^2 z^{m} \}\in {\mathcal D}(\mathcal {A}).
 $$
Suppose
$$
\{z^{m}, a z^{m}, a^2 z^{m}\}\in {\mathcal D}(\mathcal{A}).
$$
In this case, we may assume that
$$
\mathcal {A}=\mathcal {A}_{H\times \mathcal {Z}_3} \wedge \mathcal {A}_{G/\mathcal {Z}_3}.
$$
Otherwise, there exists $z^{m'} \in \mathcal {Z} \setminus H$ such that
$$\{z^{m'}, a z^{m'}\}\in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad \{z^{m'}, a^2 z^{m'}\}\in {\mathcal D}(\mathcal {A}).
$$
We may replace $z^{m}$ with $z^{m'}$.

Since both $a$ and $a^2$ are generators of $\mathcal {Z}_3$, without loss of generality we may assume that $\{z^{m}, a z^{m}\}\in {\mathcal D}(\mathcal {A})$.   Therefore,  $a^2z^m\in {\mathcal D}(\mathcal {A})$.  It follows that  $z^{3m}\in \mathcal {A}$ by Proposition  \ref{1501a}.
It is easy to see that $\{a, a^2\}\in {\mathcal D}(\mathcal {A})$.   Let
$$
K=\langle z^{m}, H\rangle,
$$
then $K$ is the unique subgroup of $G$ such that $|K: H|=3$.  Also,
$$
K \times \mathcal {Z}_{3}=\langle z^{m}, a z^{m}, H\rangle
$$
is an $S$-subgroup.

\medskip

{\bf Subcase 1.1.}\,  $K<\mathcal {Z}$.  Then for any  $z^{n} \in \mathcal {Z} \setminus K$,  the basic set containing $z^{n}$ must be $\{z^{n}, z^{n}a, a^2 z^{n}\}$.
 Otherwise, we have $\{z^{3n}\}\in {\mathcal D}(\mathcal {A})$ by Proposition \ref{1501a} and hence $z^{n} \in K$, a contradiction.
 Thus,
  $$
  \mathcal {A}=\mathcal {A}_{K\times \mathcal {Z}_3}\wedge \mathcal {A}_{G/\mathcal {Z}_3}.
  $$
 Consequently,  $\mathcal {A}$ is  traditional whenever $1< K< \mathcal {Z}$.

 \medskip

{\bf Subcase 1.2.}\,  $K=\mathcal {Z}$.  Then
$$
\{z, a z\}\in {\mathcal D}(\mathcal {A}), \, \, {\rm  or} \, \,  \{z, a^2 z\}\in {\mathcal D}(\mathcal {A}), \, \, {\rm  or} \, \, \{z, a z, a^2 z \}\in {\mathcal D}(\mathcal {A}).
$$
Suppose
\qtnl{1619c}
\{z, a z, a^2 z\}\in {\mathcal D}(\mathcal {A}).
\end{equation}

\medskip

{\bf Claim 1.}\, {\it Under the assumption in \eqref{1619c}, we have $\mathcal {A}=\mathcal {A}_{H\times \mathcal {Z}_3}\wedge \mathcal {A}_{G/\mathcal {Z}_3}$. In particular, $\mathcal {A}$ is traditional.}

{\bf Proof}.\  If the claim is false, then there exists $3\nmid k\in \mathbb{Z}$ such that
$$
\{z^k, a z^k\}\in {\mathcal D}(\mathcal {A}),\quad\text{or}\quad \{z^k, a^2 z^k\}\in {\mathcal D}(\mathcal {A}).
$$

If $k\equiv 1\:\:(mod\: 3)$,
then
  $$
\{z,az\}=\{z^{1-k}\}\{z^k, a z^k\}\in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad \{z,a^2z\}=\{z^{1-k}\}\{z^k, a^2 z^k\}\in {\mathcal D}(\mathcal {A}),
 $$
which  is a contradiction to statement \eqref{1619c}.

 If $k\equiv 2\:\:(mod\: 3)$,  then
 $$
 \{z,a^2z\}=\{z^{1+k}\}\{z^{-k}, a^2 z^{-k}\}\in {\mathcal D}(\mathcal {A}) \:\: {\rm or} \:\:\{z,az\}=\{z^{1+k}\}\{z^{-k}, a z^{-k}\}\in {\mathcal D}(\mathcal {A}),
 $$
a contradition.

It follows that $\mathcal {A}=\mathcal {A}_{H\times \mathcal {Z}_3}\wedge \mathcal {A}_{G/\mathcal {Z}_3}$ and consequently $\mathcal {A}$ is traditional.

\hfill$\square$

\medskip

 Now, we may assume that $\{z, a z\}\in {\mathcal D}(\mathcal {A})$. One can see that
 $$
 a^2z\in {\mathcal D}(\mathcal {A}) \quad\text{and}\quad \langle z^3\rangle\leq H\leq \mathcal {Z}.
$$
 Observe that for any $k\in \mathbb{Z}^+$ and $3\nmid k$,
 $$
  \{z^k,a^kz^k\}=\{a^{2k-2}z^k,a^{2k-1}z^k\}=(a^{2}z)^{k-1}\{z, a z\}\in {\mathcal D}(\mathcal {A}).
 $$
Henc, the basic sets in ${\mathcal D}(\mathcal {A})$  are the following:
 $$
 {\mathcal D}(\mathcal {A})=\{1\}; \{z^k,a^kz^k\}, \, k \in \mathbb{Z},3\nmid k; \{a^{2k}z^k\}, \,  k \in \mathbb{Z},3\nmid k; \{az^k,a^2z^k\}, \, k \in \mathbb{Z},3| k.
$$
Let $\psi$ be the automorphism defined as follows:
 \qtnl{1646c}
 \psi: G \rightarrow G,\:\:z\mapsto az,\:\:a\mapsto a^2.
 \end{equation}
 It is easy to see that $\mathcal {A}=F[G]^{\langle \psi\rangle}$. Thus, $\mathcal {A}$ is traditional.

 \medskip

 Similarly, if $\{z, a^2 z\}\in {\mathcal D}(\mathcal {A})$,  then $\mathcal {A}=F[G]^{\langle \delta \rangle}$. Here,  $\delta \in \aut(G)$ defined as:
 \qtnl{1025b}
 \delta: G \rightarrow G,\:\:z\mapsto a^2z,\:\:a\mapsto a^2.
 \end{equation}

 \medskip
 {\bf Case 2.}\, Suppose $\pi(\mathcal {A})=F[\mathcal {Z}]^{\pm}$. Then,
 $$
\{z^{m}, z^{-m}, a z^{m}, a z^{-m},a ^2z^{m}, a ^2 z^{-m}\}=\pi^{-1}\{z^{m}, z^{-m}\}
 $$
  is an $S$-set for all $m\in \mathbb{Z}$.

  Observe that for any nonzero integer $m$, $\{a^kz^m, a^kz^{-m}\}$ can not be a basic set. Otherwise, $\{z^{3m}\}$ would be a basic set by Proposition \ref{1501a}.  This is a contradiction to the assumption.

  \medskip

{\bf Subcase 2.1.}\, $H=\mathcal {Z}$.  Thus, $\{z^m, z^{-m}\}\in {\mathcal D}(\mathcal {A})$ for all $m\in \mathbb{Z}$.
It follows that
$$
\{a z^{m}, a z^{-m},a ^2z^{m}, a ^2 z^{-m}\}
$$
is an $S$-set for all $m\in \mathbb{Z}$.

\medskip

If $\mathcal {A}_{\langle a\rangle}$ is discrete, then $\mathcal {A}=\mathcal {A}_{\mathcal {Z}}\otimes\mathcal {A}_{\langle a\rangle}$ and so $\mathcal {A}$ is traditional.

\medskip

Next, suppose that
$$
\{a,a^2\}\in {\mathcal D}(\mathcal {A}) \quad\text{and}\quad \{az, az^{-1}, a^2z, a^2z^{-1}\}\in {\mathcal D}(\mathcal {A}).
$$
If there exists $k\in \mathbb{Z}^+$ such that $\{az^k,az^{-k}\}\in {\mathcal D}(\mathcal {A})$, then
$$
az^{2k+1}+az^{-1}+az+az^{-2k-1}=(az^k+az^{-k})(z^{k+1}+z^{-k-1})\in \mathcal {A}.
$$
This implies that $\{az^{-1},az\}\in {\mathcal D}(\mathcal {A})$,  a contradiction. Similarly, it is impossible that  $\{az^k,a^2z^{-k}\}\in {\mathcal D}(\mathcal {A})$.   We conclude that
$$
\{az^k,az^{-k},a^2z^k,a^2z^{-k}\}\in {\mathcal D}(\mathcal {A}), \forall  k\in \mathbb{Z}^+.
$$
And therefore,  $\mathcal {A}=\mathcal {A}_{\mathcal {Z}}\otimes\mathcal {A}_{\langle a\rangle}$.  In particular, $\mathcal {A}$ is traditional.

\medskip

Now, suppose that
$$
\{a,a^2\}\in {\mathcal D}(\mathcal {A}) \quad\text{and}\quad \{az,a^2z^{-1}\}
\in {\mathcal D}(\mathcal {A}).
$$
By the observation at the beginning of Case 2, $\{az^{-1}, a^2z\}$ is also a basic set.

Applying Lemma \ref{1951b} to $\{az,a^2z^{-1}\}$ and $\{az^{-1}, a^2z\}$, we obtain that
$$
\{a^kz^{-k}, a^{2k}z^k\}\quad\text{and}\quad \{a^kz^k, a^{2k}z^{-k}\}\in {\mathcal D}(\mathcal {A}), \forall  3\nmid k\in \mathbb{Z}.
$$
Thus, for any $3\nmid k$,
$$
a^{k+1}z^{k-1}+a^{2k+1}z^{-k-1}+a^{k+2}z^{k+1}+a^{2(k+1)}z^{-k-1}=(az^{-1}+a^2z)(a^kz^k+a^{2k}z^{-k})\in \mathcal {A}.
$$
In particular, if $k\equiv 1\:\:(mod\: 3)$, one can see that
$$
\{az^k, a^2z^{-k}\}\in {\mathcal D}(\mathcal {A}), \forall 3|k.
$$
We conclude that $\{az^k, a^{2}z^{-k}\}\in {\mathcal D}(\mathcal {A})$ for all $ k\in \mathbb{Z}$.  It is easy to see that
$$
\mathcal {A}=F[G]^{\langle\xi\rangle},
$$
where $\xi\in \aut(G)$ is defined as
\qtnl{1649b}
\xi:  G \rightarrow G, \:\:z\mapsto z^{-1},\:\: a\mapsto a^2.
\end{equation}

\medskip

Finally, neither  $\{az, az^{-1}\}$ nor  $\{az, a^2z\}$ is a basic set by the observation at the beginning of Case 2.

\medskip

{\bf Subcase 2.2.}\, $H< \mathcal {Z}$.

Choose $z^{m} \in \mathcal {Z} \setminus H$. Then one of the following holds:
\begin{enumerate}
\tm{i} $\{z^{m}, z^{-m}, a z^{m}, a z^{-m},a ^2z^{m}, a ^2 z^{-m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{ii} $\{z^{m}, z^{-m},a z^{m},a ^2 z^{-m}\}\in {\mathcal D}(\mathcal {A})$ and  $\{a z^{-m}, a^2 z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{iii}$\{z^{m}, z^{-m},a z^{-m},a ^2 z^{m}\}\in {\mathcal D}(\mathcal {A})$ and $\{a^2 z^{-m}, a z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{iv}$\{z^{m},  a z^{m}, a z^{-m}\}\in {\mathcal D}(\mathcal {A})$ and $\{ z^{-m}, a^2 z^{-m}, a ^2z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{v} )$\{z^{m},  a^2 z^{-m}, a z^{-m}\}\in {\mathcal D}(\mathcal {A})$ and $\{ z^{-m}, a z^{m}, a ^2z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{vi}$\{z^{m},  a^2 z^{-m}, a^2 z^{m}\}\in {\mathcal D}(\mathcal {A})$ and $\{ z^{-m}, a z^{-m}, a ^2z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{vii}$\{z^{m}, a z^{-m}\}\in {\mathcal D}(\mathcal {A})$, $\{ z^{-m}, a ^2z^{m}\}\in {\mathcal D}(\mathcal {A})$ and $\{ a^2z^{-m}, a z^{m}\}\in {\mathcal D}(\mathcal {A})$,
\tm{viii}$\{z^{m}, a^2 z^{-m}\}\in {\mathcal D}(\mathcal {A})$, $\{ z^{-m}, az^{m}\}\in {\mathcal D}(\mathcal {A})$ and $\{ az^{-m}, a^2 z^{m}\}\in {\mathcal D}(\mathcal {A})$.
\end{enumerate}

If (i)  holds for all $z^m\in \mathcal {Z}\setminus H$ , then
$$
\mathcal{A}=\mathcal{A}_{H\times \mathcal{Z}_3}\wedge\mathcal{A}_{G/\mathcal{Z}_3}.
$$
In particular, $\mathcal {A}$ is traditional.

\medskip

Conditions (iv),  (v),  and (vi) can not happen  by Theorem \ref{1521b}.

\medskip

 Now, assume that statement (ii), (iii), (vii),   or (viii) holds.
In these cases, we conclude that $z^{3m}\in H$ by Proposition \ref{1501a}. Set
 $$
 M=\langle z^{m}, H\rangle=\langle z^{-m}, H\rangle.
 $$
 Then $M$ is the unique subgroup of $G$ such that $|M: H|=3$.
 In addition,
 $$
 M \times \mathcal {Z}_{3}=\langle z^{m},z^{-m}, a z^{m},a^2z^{-m}, H\rangle=\langle z^{m},z^{-m}, a z^{-m},a^2z^{m}, H\rangle
 $$
is an $S$-subgroup.

If $M<\mathcal{Z}$, then for any $z^{l} \in \mathcal{Z} \backslash M$, by the above arguments, one can easily see that  the basic  set containing $z^{l}$ must be
$$
\{z^{l}, z^{-l}, a z^{l}, a z^{-l},a ^2z^{l}, a ^2 z^{-l}\}.
   $$
  Hence,
  $$
  \mathcal {A}=\mathcal {A}_{M\times \mathcal {Z}_3}\wedge \mathcal {A}_{G/\mathcal {Z}_3},
  $$
 and so $\mathcal {A}$ is traditional.

 \medskip

If $M=\mathcal {Z}$,  then $H=\langle z^3\rangle$ by $|M:H|=3$. There are fives cases to consider.
\begin{enumerate}
\tm{1}$\{z, z^{-1}, a z, a z^{-1},a ^2z, a ^2 z^{-1}\}\in {\mathcal D}(\mathcal {A})$,
\tm{2}$\{z, z^{-1},a z,a ^2 z^{-1}\}\in {\mathcal D}(\mathcal {A})$ and $\{a z^{-1}, a^2 z\}\in {\mathcal D}(\mathcal {A})$,
\tm{3}$\left\{z, z^{-1},a z^{-1},a ^2 z\right\}\in {\mathcal D}(\mathcal {A})$,  and $\{a^2 z^{-1}, a z\}\in {\mathcal D}(\mathcal {A})$,
\tm{4}$\{z, a z^{-1}\}\in {\mathcal D}(\mathcal {A})$, $\{ z^{-1}, a ^2z\}\in {\mathcal D}(\mathcal {A})$,  and $\{ a^2z^{-1}, a z\}\in {\mathcal D}(\mathcal {A})$,
\tm{5}$\{z, a^2 z^{-1}\}\in {\mathcal D}(\mathcal {A})$, $\{ z^{-1}, az\}\in {\mathcal D}(\mathcal {A})$,  and $\{ az^{-1}, a^2 z\}\in {\mathcal D}(\mathcal {A})$.
\end{enumerate}

\medskip

Suppose condition (1) holds. Then we have the following claim.

\medskip

{\bf Claim 2.}\,  $\mathcal {A}=\mathcal {A}_{H\times \mathcal {Z}_3}\wedge \mathcal {A}_{\mathcal {Z}_3}$.

{\bf Proof}.\  If the claim is false, then there exists $k\in \mathbb{Z}^+$ not divisible by $3$ such that
$$
\{a^2z^{-k}, az^k\}\in {\mathcal D}(\mathcal {A}) \quad\text{or}\quad  \{az^{-k}, a^2z^k\}\in {\mathcal D}(\mathcal {A}).
$$
First, assume that $\{a^2z^{-k}, az^k\}\in {\mathcal D}(\mathcal {A})$.

\medskip

If $k\equiv 1\:\:(mod\: 3)$,  then $\{z^{k-1},z^{1-k}\}\in {\mathcal D}(\mathcal {A})$ by $3|{k-1}$. It follows that
$$
a^2z^{-1}+ a^2z^{1-2k}+ az^{2k-1}+ az=(a^2z^{-k}+ az^k)\cdot (z^{k-1}+z^{1-k})\in \mathcal {A}.
$$
This yields  that $\{a^2z^{-1},az\}\in {\mathcal D}(\mathcal {A})$, a contradiction.

\medskip

If $k\equiv 2\:\:(mod\: 3)$,  then $\{z^{k+1},z^{-1-k}\}\in {\mathcal D}(\mathcal {A})$ by $3|{k+1}$.  It follows that
$$
az^{-1}+a^2z^{-1-2k}+ az^{2k+1}+ a^2z=(a^2z^{-k}+ az^k)\cdot (z^{k+1},z^{-1-k})\in \mathcal {A}.
$$
Thus,  $\{az^{-1},a^2z\}\in {\mathcal D}(\mathcal {A})$, a contradiction.

\medskip

Therefore,  the claim is valid in this case. Similarly, the claim holds if~$\{az^{-k}, a^2z^k\}$ is a basic set
for some positive integer  $k$ not divisible by $3$.
\hfill$\square$

\medskip

Suppose condition (2) holds. Then, by Lemma \ref{1951b} $\{a^2z^{-2},az^2\}=\{az^{-1},a^2z\}^{(2)}$ is a basic set.  By the statements at the beginning of Subcase 2.2, this implies that
$$
\{z^2, z^{-2},a z^{-2},a ^2 z^2\}\in {\mathcal D}(\mathcal {A}) \quad\text{or}\quad \{z^2, a z^{-2}\}\in {\mathcal D}(\mathcal {A}).
$$
If $\{z^2, a z^{-2}\}\in {\mathcal D}(\mathcal {A})$, then
$$
az+a^2z^{-3}+a^2z^3+z^{-1}=(az^{-1}+a^2z)\cdot(z^2+az^{-2})\in \mathcal {A}.
$$
 It means that $\{az,z^{-1}\}$ is an $S$-set,  a contradiction to the assumption of (2). We conclude that
 \qtnl{1455a}
 \{z^2, z^{-2},a z^{-2},a ^2 z^2\}\in {\mathcal D}(\mathcal {A}).
 \end{equation}
 Note that
 $$
 \{az^3, az^{-3}\}\in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad \{az^3, a^2z^{-3}\} \in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad \{az^3, az^{-3}, a^2z^{-3}, a^2z^3\}\in {\mathcal D}(\mathcal {A}).
 $$

 If $\{az^3, az^{-3}\}\in {\mathcal D}(\mathcal {A})$,  then
 $$
 a^2z^2+z^4+a^2z^{-4}+z^{-2}=(az^3+ az^{-3})(az^{-1}+a^2z)\in \mathcal {A}.
 $$
So, $\{a^2z^2,z^{-2}\}$ is an $S$-set, a contradiction to statement \eqref{1455a}.

 If $\{az^3, a^2z^{-3}\} \in {\mathcal D}(\mathcal {A})$,  then
 $$
 a^2z^2+z^4+z^{-4}+az^{-2}=(az^3+ a^2z^{-3})(az^{-1}+a^2z)\in \mathcal {A}
 $$
 and hence $\{a^2z^2,  az^{-2}\}$ is an $S$-set,  another  contradiction to statement \eqref{1455a}.

As a consequence, we obtain
\qtnl{1509b}
\{az^3, az^{-3}, a^2z^{-3}, a^2z^3\}\in {\mathcal D}(\mathcal {A}).
\end{equation}

Next, we will prove the following claim.

\medskip

{\bf Claim 3.}\, {\it The following statements hold:
\begin{enumerate}
\tm{a} $\{z^n, z^{-n}, a^nz^{n}, a^{2n}z^{-n}\}\in {\mathcal D}(\mathcal {A})$, if  $n\in \mathbb{Z}^+$ and $3\nmid n$,
\tm{b} $\{az^n, az^{-n}, a^2z^{n}, a^2z^{-n}\}\in {\mathcal D}(\mathcal {A})$,  if  $n\in \mathbb{Z}^+$  and  $3| n$.
\end{enumerate}
}
{\bf Proof}.\
We prove the claim by induction on $n$. If $n=1$, $n=2$, or $n=3$, the claims follow by the assumption and statements \eqref{1455a} and \eqref{1509b}.

\medskip

If $n\equiv 0\:\:(mod\: 3)$,  then $n-1\equiv 2\:\:(mod\: 3)$.  Hence,
$$
\{z^{n-1},z^{1-n},az^{1-n},a^2z^{n-1}\}\in {\mathcal D}(\mathcal {A}),
$$
by inductive hypotheses.

Assume that $\{az^n,az^{-n}\}\in {\mathcal D}(\mathcal {A})$. Then
$$
a^2z^{n-1}+z^{n+1}+a^2z^{-n-1}+z^{1-n}=(az^n+az^{-n})(az^{-1}+a^2z)\in \mathcal {A}.
$$
It follows that $\{a^2z^{n-1},z^{1-n}\}$ is an $S$-set,  a contradiction.

Assume that $\{az^n,a^2z^{-n}\}\in {\mathcal D}(\mathcal {A})$. Then
$$
a^2z^{n-1}+z^{n+1}+z^{-n-1}+az^{1-n}=(az^n+a^2z^{-n})\cdot (az^{-1}+a^2z)\in \mathcal {A}.
$$
Thus,  $\{a^2z^{n-1},az^{1-n}\}$ is an $S$-set,  a contradiction.
Hence,
$$
\{az^n, az^{-n}, a^2z^{n}, a^2z^{-n}\}\in {\mathcal D}(\mathcal {A}).
$$

\medskip

If $ n\equiv 1\:\:(mod\: 3)$,  then $n-1\equiv 0\:\:(mod\: 3)$.  So,
$$
\{az^{n-1}, az^{1-n},a^2z^{1-n}, a^2z^{n-1}\}\in {\mathcal D}(\mathcal {A}).
$$
by inductive hypotheses. By Lemma \ref{1951b}, one can see that
$$
\{az^{-n}, a^2z^n\}=\{az^{-1}, a^2z\}^{(n)}
$$
is an $S$-set.  Therefore,  $\{az^{-n}, a^2z^n\}\in {\mathcal D}(\mathcal {A})$.   It  yields that
$$
\{z^{n}, z^{-n},a z^{n},a ^2 z^{-n}\}\in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad  \{z^{n}, a^2 z^{-n}\}\in  {\mathcal D}(\mathcal {A}).
$$
Assume that $\{z^{n}, a^2 z^{-n}\}\in {\mathcal D}(\mathcal {A})$.  Then
$$
az^{n-1}+a^2z^{n+1}+z^{-n-1}+az^{1-n}=(z^{n}+a^2 z^{-n})\cdot(az^{-1}+a^2z)\in \mathcal {A}.
$$
This implies that  $\{az^{n-1},az^{1-n}\}$ is an $S$-set,  a contradiction.

Thus,
$$
\{z^{n}, z^{-n}, a z^{n}, a ^2 z^{-n}\}\in {\mathcal D}(\mathcal {A}).
$$

\medskip

If $n\equiv 2\:\:(mod\: 3)$,  then $n-1\equiv 1\:\:(mod\: 3)$.  Hence,
 $$
 \{z^{n-1},z^{1-n},a^2z^{1-n},az^{n-1}\}\in {\mathcal D}(\mathcal {A})
 $$
 by inductive hypotheses.

 Similarly, we have $\{a^2z^{-n}, az^n\}\in {\mathcal D}(\mathcal {A})$. One can see that
 $$
 \{z^{n}, z^{-n},a z^{-n}, a ^2 z^{n}\}\in {\mathcal D}(\mathcal {A}), \quad\text{or}\quad \{z^{n}, a z^{-n}\}\in {\mathcal D}(\mathcal {A}).
 $$
 Suppose that $\{z^{n}, a z^{-n}\}\in {\mathcal D}(\mathcal {A})$,  then
 $$
 az^{n-1}+a^2z^{n+1}+a^2z^{-n-1}+z^{1-n}=(z^{n}+ a z^{-n})\cdot(az^{-1}+a^2z)\in \mathcal {A}.
 $$
This yields that $\{az^{n-1},z^{1-n}\}\in {\mathcal D}(\mathcal {A})$,  a contradiction.  Hence,
$$
\{z^{n}, z^{-n}, az^{-n}, a ^2 z^{n}\}\in {\mathcal D}(\mathcal {A}).
$$
This completes the proof of Claim 3.
\hfill$\square$

\medskip

Obviously, $\{a,a^2\}\in S(\mathcal{A})$.  As a consequence of Claim 3, one can easily see that
$$
\mathcal {A}=F[G]^{\langle \psi, \xi\rangle},
$$
where $\psi$ and $\xi$ are automorphisms of $G$ defined as in \eqref{1646c} and \eqref{1649b}, respectivelty.

\medskip

Suppose condition (3) holds.  One can similarly get that
$$
\mathcal {A}=F[G]^{\langle \delta, \xi\rangle},
$$
where $\delta$ and $\xi$ are defined as in \eqref{1025b} and \eqref{1649b}, respectively.

\medskip

Suppose condition (4) holds. Applying Lemma \ref{1951b},  for all integers $k$ with $ 3\nmid k$  we obtain that
$$
\{z^k, a^k z^{-k}\}\in {\mathcal D}(\mathcal {A}),  \quad \{ z^{-k}, a ^{2k}z^k\}\in {\mathcal D}(\mathcal {A}),  \quad\text{and}\quad \{ a^{2k}z^{-k}, a^k z^k\}\in {\mathcal D}(\mathcal {A}).
$$
Assume that  $k\equiv 1\:\:(mod\: 3)$.  Observe that
$$
az^{k-2}+a^2z^{k+2}+a^2z^{-k-2}+z^{2-k}=(z^k+a z^{-k})(az^{-2}+a^2z^{2})\in \mathcal {A}.
$$
This yields that  $\{a^2z^{k+2},a^2z^{-k-2}\}\in {\mathcal D}(\mathcal {A})$.

Assume that  $k\equiv 2\:\:(mod\: 3)$.  Then
$$
az^{k+1}+a^2z^{k-1}+z^{1-k}+az^{-1-k}=(z^k+a^2 z^{-k})(az+a^2z^{-1})\in \mathcal {A}.
$$
This implies that $\{az^{k+1},az^{-k-1}\}\in {\mathcal D}(\mathcal {A})$.

We conclude that
$$
\{az^{k},az^{-k}\}\in {\mathcal D}(\mathcal {A}), \quad\text{and}\quad \{a^2z^{k},a^2z^{-k}\}\in {\mathcal D}(\mathcal {A}), \forall  3| k\in \mathbb{Z}^+.
$$
Note that
$$
(z+az^{-1})(z+az^{-1})=z^2+2a+az^{-2}\in \mathcal {A}.
$$
Thus, $\{a\}\in {\mathcal D}(\mathcal {A})$, and hence  $\mathcal{A}_{\langle a\rangle}$ is discrete.

Now let $\rho\in \aut(G)$ be defined as follows:
\qtnl{1105c}
\rho: G \rightarrow G,\:\:z\mapsto az^{-1},\:\: a\mapsto a.
\end{equation}

It is easy to see that
$$
\mathcal {A}=F[G]^{\langle \rho \rangle}.
$$
In particular, $\mathcal {A}$ is traditional.

\medskip

Suppose condition (5) holds.  Let $\sigma\in \aut(G)$ be defined as follows:
\qtnl{1110a}
\sigma: G \rightarrow G, \:\:z\mapsto a^2z^{-1},\:\:a\mapsto a.
\end{equation}
Then one can see that
$$
\mathcal {A}=F[G]^{\langle \sigma \rangle}.
$$