\section{Introduction}
Capacitance-voltage (CV) measurements are a widely used technique in silicon detector physics. By increasing the reverse bias voltage while scanning the capacitance, it is possible to derive crucial parameters like the full depletion voltage or the doping concentration. However, this only works if the sensor is not too heavily irradiated. After irradiation, the measured CV curves show strong frequency dependencies which is not the case before irradiation and thus, the derived parameters vary over a wide range, indicating that the method is not applicable for such devices.\\
Many different models \cite{BEGUWALA, schibli} have been derived to explain frequency-dependent capacitance values but none of these were able to describe irradiated silicon detectors over the entire frequency range. In this work, an existing model for fitting CV measurements of irradiated sensors is extended. From this, it is investigated where the different frequency behaviour of the current CV analysis technique originates and it is shown how sensor parameters can be determined. Hereby, unirradiated sensors produced by a CMOS foundry \cite{DIEHL} as well as irradiated CMOS, ATLAS R0 and ATLAS R5 sensors have been CV measured over the full, technically available, frequency range and were fitted with the derived model.

\section{Setup and Measurement Procedure}
The setup used in this work is based on a common setup for CV measurements which was extended for the special needs of irradiated devices. The basic setup consists of a Keythley 237 high voltage source used to bias the sensor and simultaneously measure the leakage current. Additionally, a HP 4284 LCR meter can be (IV resp. CV configuration) AC coupled to the high voltage source to measure the capacitance without high voltage at the LCRs inputs. The sensor is placed on a chuck which is directly connected to the HV potential of the HV source and is contacted with micro positioner needles at the frontside to deliver the ground potential to the corresponding pad. The chuck consists of a top layer with embedded holes for vacuum supply and an embedded PT100 element for temperature monitoring. Underneath, peltier elements, allowing to cool the top layer to a constant temperature, are located on a cooling jig which precools the system and removes the heat of the peltier elements wherefore a Julabo chiller is used. The whole chuck stack and the micro positioners are housed in a freezer which is flushed with nitrogen to ensure a dry environment, preventing icing of the sensor while cooling down. The high voltage supply, the LCR, the chiller and the power supply of the peltier elements can be remotely controlled by the PC to acquire data or control and monitor the temperature. The LCR meter was set in parallel mode (CpRp) as this gives the best performance for low capacitances and high resistances. \\
The voltage source offers a bias voltage range up to \SI{1100}{V}, the LCR can measure capacitance, resistance or inductance in a frequency range between \SI{20}{Hz} and \SI{1}{MHz} and the cooling stack reaches temperatures, depending on the size of the sensor, down to \SI{-45}{\celsius} resp. \SI{-50}{\celsius}.\\
The measurements shown in this work were taken by ramping the voltage to a desired value, sweeping the frequency of the LCR meter across the entire offered range before measuring the current  in CV configuration. This induces more noise on the current measurement but, in contrast to ramping the sensor twice in IV and in CV configuration, systematic variations of the current are excluded. These could be caused by two effects: First, the so called hysteresis of the sensor (multiple identical runs show a systematic decrease in current) which would result in a lower current in the second run. Second, when operating the sensor close to breakdown, the current increases for some sensors with every consecutive run. A precise current measurement corresponding to every capacitance is necessary for the later analysis as both are taken into account at the same time.

\section{Theory and Model Construction}
\label{TheoModel}
The depletion width of silicon detectors features a square root dependence on the reverse bias voltage. This depletion width induces a capacitance and, when assuming a parallel plate capacitor, the following relation between the capacitance and the voltage can be deducted \cite{Sze}:

\begin{align}
	\label{C^2}
	C^{-2} = \frac{2 \epsilon_0 \epsilon_r A^2}{q_0 N_\text{eff}}V
\end{align}

From this equation the well used method for deriving the full depletion voltage $V_\text{FD}$ and the effective doping concentration $N_\text{eff}$ can be inferred if the sensor active area $A$ is known: The inverse square of the capacitance is plotted against the voltage, a linear fit is applied to the increasing section, the flat part determines $V_\text{FD}$ and the slope is solely dependent to $N_\text{eff}$.\\
However, the physical explanation and the above shown equation both indicate that the capacitance should only be dependent on the applied voltage and not on the frequency used for the capacitance measurement. This holds true for unirradiated sensors but as the fluence exceeds a threshold of \SI{\approx\,e13}{}-\SI{e14}{\text{n}_\text{eq}}, depending on the base material, a clear dependency can be seen.\\
One reason for frequency dependence is the introduction of defects by radiation. Their description involves trapping and detrapping time constants and can thereby influence the time dependent response to a superimposed AC voltage \cite{BEGUWALA, schibli}. These time constants are in the order of \si{\micro\second} to \si{\milli\second}.  In Ref. \cite{SCHWANDT2019162418} a further explanation which introduces a frequency dependency was discussed. The basic assumption is an increasing resistivity in the non-depleted part of the sensor (which is said to be low enough to be neglected for unirradiated material). Subsequently, the resistivity of the undepleted layer and the capacitance of the depleted layer in series form an RC circuit which shields the bulk material from the oscillating voltage of the LCR. Therefore, for high frequencies, a capacitance measurement always yields the value of the sensor in full depletion. When going lower in frequency, a transition region occurs which leads to a frequency-independent regime following the CV relation given by \cref{C^2}. To further generalize our approach, a finite resistivity of the depleted bulk is introduced.\\
When measuring strip detectors (in contrast to pad diodes) another effect has to be taken into account. Each strip is biased via a bias resistor from the bias ring. This resistance also forms an RC circuit together with the bulk capacitance but, as the resistance is different to the one of the non depleted layer, this results in a further cut off frequency. Combining the beforehand described effects lead to the sensor equivalent circuit shown in \cref{EqCir} with capacitances and resistances being

\begin{align}
	\label{CapsRes}
	C_1 = \dfrac{\epsilon_0 \epsilon_r A}{d-w} \;\; \text{resp.} \;\;	C_2 = \dfrac{\epsilon_0 \epsilon_r A}{w}\\
	R_1 = \dfrac{\rho_\text{nd} (d-w)}{A} \;\; \text{resp.} \;\; R_2 = \dfrac{\rho_\text{d} w}{A}
\end{align}

with $w$ being the depletion width, $d$ the sensor thickness, $A$ the active sensor area, $\rho_\text{nd}$ the resistivity of the non-depleted silicon bulk and $\rho_\text{d}$ the resistivity of the depleted bulk.\\

\begin{figure}
	\centering
	\includegraphics[width=0.4\textwidth]{EqCir.png}
	\caption{Equivalent circuit of an irradiated strip detector.}
	\label{EqCir}
\end{figure}

To construct a fit function from this, one has to compute the impedance of the equivalent circuit of \cref{EqCir}, equate it to the impedance of the simple RC parallel circuit assumed by the LCR meter and solve for $C_2$ and $R_2$ respectively. After this the measured capacitance and resistance can be related to the device parameters $d$, $A$, $\rho_\text{nd}$, $\rho_\text{d}$ and $R_\text{bias}$.

This model has been shown to describe the measurements quite accurately for high frequencies but not for low frequencies (already described in Ref. \cite{Li}). As mentioned before, a frequency-independent regime is expected. However, even for frequencies down to \SI{20}{Hz}, a strong dependency remains. As discovered for unirradiated measurements, it turns out that a current flowing through the sensor results in an additional capacitance contribution being proportional to $e^{-\sqrt{f}}$ which is called low-frequency (LF) increase in the following (\cref{CVincStripUnirrad}). The origin of this contribution is not completely understood but some observations give hints towards diffusion capacitance as illustrated in \cref{UnirradDevs}. It is important to mention that this term cannot be represented by either a single or multiple additional elements in the equivalent circuit as these would only introduce a broken rational function but not an exponential function.

Furthermore, it turns out that the same analytic term is also necessary to be able to fit resistance data. To account for unknown quantities in this term there are two parameters $c_1$ and $c_2$ empirically introduced in the following way

\begin{align}
	C_\text{LF} = c_1 e^{- c_2 \sqrt{f}} \;\; \text{resp.} \;\; R_\text{LF} = c_1 e^{- c_2 \sqrt{f}}
\end{align}

These two terms are added to the capacitance and the resistance deduced from the impedance. The fit functions are analytically computed using Wolfram Mathematica and have a complex analytic structure which makes it not possible to show here. Regarding the capacitance, it has to be noticed that not the whole sensor is biased through the bias resistors but there is a (small) area under the bias ring which is biased directly. This leads to the fact that for very high frequencies, the capacitance does not drop to exactly zero. This remaining capacitance consists of the front to backside capacitance under the bias ring as well as of the capacitance of the bias ring towards the outermost strips and the guard ring. Therefore, it is impossible to find an analytic expression for it and in this work an additional parameter $C_\text{HF}$ was introduced. In the end the fitting functions, $C(f)$ and $R(f)$, are dependent on the global parameters, $d$, $A$, $\rho_\text{nd}$, $\rho_\text{d}$, $R_\text{bias}$ and $C_\text{HF}$, which are expected to be constant for each sensor and independent of the applied voltage, as well as the individual parameters for every voltage, $w$, $c_1$ and $c_2$. The latter two should be ideally voltage-independent but for generality purposes are introduced individually for every voltage and later on their relation to physical quantities is investigated.\\
To illustrate how the frequency dependence looks like and how the above mentioned constituents contribute, a simulation of the capacitance and resistance is shown in \cref{FitFunc}.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{FitFunc.png}
	\includegraphics[width=0.9\textwidth]{FitFuncR.png}
	\caption{Simulation of the frequency dependence of the measured capacitance and resistivity for a partially depleted sensor. The low-frequency part in grey is depletion dependent, i.e. the capacitance decreases and the resistance increases with increasing voltage.}
	\label{FitFunc}
\end{figure}

The red line shows a complete capacitance/resistance sweep for a partially depleted irradiated sensor. In contrast to this, the brown line depicts the same function without the $c_\text{LF}$ term. Hereby, the frequency-independent level, predicted by the equivalent circuit, becomes visible. For the capacitance, this level will move further downwards with increasing depletion, resp. voltage, according to the $\sqrt{V}$ law shown in \cref{C^2}. For the resistance, it will move further up as the more resistive depletion zone grows. Looking at the resistance, the voltage-independent high-frequency limit directly yields the bias resistance divided by the number of strips. This is valid as long as it is high enough to become constant when sweeping the measurement frequency towards higher values. For lower values, the cutoff frequency, and therefore the transition, occurs at higher frequencies which might exceed the experimentally limited frequency range. In orange the function for a vanishing bias resistance (pad diode) is shown which is constant for higher frequencies, for the capacitance this value is equal to the full depletion capacitance. Finally, the line in yellow is for a vanishing resistivity of the non depleted bulk and indicates the capacitance for an unirradiated pad diode with the same area and depletion width which is completely frequency-independent.

\section{Results}

\subsection{Unirradiated Devices}
\label{UnirradDevs}
As an example, the CV measurements of an unirradiated strip detector are shown in \cref{CVincStripUnirrad}. The same behaviour as shown for this detector has been observed and verified for other sensors as well.

\begin{figure}
  \centering
      \includegraphics[width=0.9\textwidth]{CVincStripUnirradLog.png}
  \caption{IV and CV measurement of an unirradiated strip detector with soft breakdown. The line in red is the current and corresponds to the axis on the left, the family of lines are the capacitance at multiple frequencies. The capacitance increase occurs simultaneously with the current increase (\SI{\approx\,60}{V}).}
  \label{CVincStripUnirrad}
  \includegraphics[width=0.9\textwidth]{Cincfit.png}
  \caption{Low-frequency capacitance in the range of the soft breakdown. The lines in green show the fit results.}
  \label{Cincfit}
\end{figure}

The capacitance decreases up to voltages of around \SI{30}{V} which is the full depletion voltage of this sensor. Also, this capacitance decrease, as well as the constant level afterwards, is frequency-independent up to approximately \SI{60}{V}. Above this level, a soft breakthrough occurs which can be seen by the current increase (red line). At the same time, the capacitance starts to increase again albeit with lower magnitude for higher frequencies. This already indicates that it might be caused by the diffusion capacitance which is current driven and therefore commonly said to be neglectable in reverse bias of a diode structure.\\
To analyse the frequency dependence of the capacitance increase, the same data is plotted against the frequency as shown in \cref{Cincfit}.

For higher frequencies and voltages above full depletion (\SI{\approx\,30}{V}), the capacitance is voltage-independent (therefore the capacitance at full depletion $C_\text{FD}$ is introduced), only for frequencies below \SI{100}{Hz} an increase can be seen. It was observed that the latter obeys a $e^{-\sqrt{f}}$ law and, to test the hypothesis of the capacitance increase being influenced by the diffusion capacitance, the fit function was extended by the $\frac{\text{d}I}{\text{d}V}$ value corresponding to each voltage. In total the fit function for every voltage reads:

\begin{align}
	\label{fitfunc}
	C(V_0, f) = C_\text{FD} + c_1 \frac{\text{d}I}{\text{d}V} \bigg\rvert_{V_0} e^{- c_2 \sqrt{f}}
\end{align}

With this formula, the fits of all unirradiated sensors tested are reasonable, the values of the parameters $c_1$ and $c_2$ show no systematics among the voltages and exhibit a standard deviation of 3.8 resp. \SI{0.5}{\%}. The absence of further voltage dependencies indicates that there are no other physical quantities influencing the capacitance increase and, in particular, that the $e^{-\sqrt{f}}$ law is valid and can be taken into account for other current driven capacitance increases.

\subsection{Irradiated Devices}

The derived model in \cref{TheoModel} is applied to full sets of frequency-dependent CV measurements whereby the capacitance and resistance values have been used independently. The result of a capacitance fit is shown, as an example, in \cref{CVFitkomplett}. This fit was applied to different sensors of varying temperatures and the delivered set of parameters as described in \cref{TheoModel} is analysed in the following chapters to prove the reliability of this model.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{CVFitkomplett.png}
	\caption{Full set of frequency-dependent CV measurements with according fit.}
	\label{CVFitkomplett}
\end{figure}

\subsubsection{Depletion Comparison}

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{DepCompCMOS.png}
	\includegraphics[width=0.9\textwidth]{DepCompCMOSR.png}
	\caption{Depletion width of CMOS sensors irradiated to \SI{1e14}{\text{n}_\text{eq}}, \SI{5e14}{\text{n}_\text{eq}} and \SI{1e15}{\text{n}_\text{eq}} at different temperatures from capacitance data (top) and resistance data (bottom).}
	\label{DepCompCMOS}
\end{figure}

First of all the depletion width of three CMOS sensors irradiated to different fluence levels is investigated. Therefore, the depletion width derived from the capacitance data as well as the width derived from the resistance data is shown at the top and bottom of \cref{DepCompCMOS} respectively. Both plots exhibit no systematic differences which indicates a good consistency of the model. It can be seen that the depletion increases with increasing voltage and the slope of the sample irradiated to \SI{1e14}{\text{n}_\text{eq}} is the steepest. At around \SI{90}{V} there is a sharp change to a constant depletion level which is slightly below the sensor thickness of \SI{150}{\micro\meter}. In contrast to expectation, the initial rise of the samples with a fluence of \SI{5e14}{\text{n}_\text{eq}} and \SI{1e15}{\text{n}_\text{eq}} is equal within the measurement uncertainties. However, above \SI{150}{V}, a clear difference in slope becomes visible. From this plot a depletion voltage of around \SI{90}{V}, \SI{300}{V} and \SI{450}{V} can be extracted. For the latter, the sensor could not be biased high enough. Because of that, the increasing depletion was extrapolated to the full depletion level of the \SI{1e14}{\text{n}_\text{eq}} measurement. Additionally, a small dependency on the temperature is visible as a constant offset towards larger depletion for lower temperatures but no significant difference in depletion voltage.\\
To verify the beforehand shown data, the CV measurements are compared to Alibava\cite{Alibava, DIEHL2} measurements of the same sensors in \cref{DepCompCMOSCRAlibava}. A fully depleted sensor of \SI{150}{\micro\meter} width is expected to collect a charge of $\SI{11.6}{\kilo e^-}$. A partially depleted sensor collects proportionally less charge due to the lower active thickness. Irradiated devices have additional charge loss channels like trapping or carrier recombination. As can be seen in \cref{DepCompCMOSCRAlibava}, the measurements for the lower fluences \SI{1e14}{\text{n}_\text{eq}} and \SI{3e14}{\text{n}_\text{eq}} show similar trends within the errors while for the higher fluences \SI{5e14}{\text{n}_\text{eq}} and \SI{1e15}{\text{n}_\text{eq}} the charge collection measurements show significantly less charge than expected from the depletion. This, however, can be explained by the additional charge loss possibilities. In conclusion, it can be stated that the depletion width determined from this CV analysis is in compliance with charge collection measurements.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{DepCompCMOSCRAlibava.png}
	\caption{Depletion width of CMOS sensors irradiated to different fluences determined from capacitance and resistance data of CV measurements compared to charge collection measurements (black) of the according sensors.}
	\label{DepCompCMOSCRAlibava}
\end{figure}

To illustrate the problem that arises from the traditional method of deriving the full depletion capacitance or the effective doping concentration, the capacitance of the sensor irradiated to \SI{1e14}{\text{n}_\text{eq}} is plotted inverse squared in \cref{DepComp1e14}. One can see a linear increase below \SI{100}{V} but the slope is increasing with increasing frequency. Additionally, the lines for higher frequencies show an offset. This gives a lot of freedom to choose the voltage range and the frequency used for the linear fit and therefore a large spread in possible results, a trend even more strongly pronounced for higher fluences.  The reasons for the distorted $C^{-2}$ lines can be seen in \cref{FitFunc}. For low frequencies, there is an additional capacitance term, most likely induced by the diffusion capacitance. This effect reduces the $C^{-2}$ value and therefore the slope is too low. With increasing frequency this influence becomes weaker but unfortunately a further effect kicks in before it is removed completely. At a certain cutoff frequency, the irradiated bulk with its high resistivity starts to screen the depleted layer from the AC signal and therefore the capacitance starts to decrease. By choosing a frequency above the cutoff frequency it would be possible to completely avoid the influence of the low-frequency increase but the cutoff frequency itself depends on the capacitance of the depleted bulk. This is why it shifts towards higher frequencies while depleting the sensor. In contrast to the recommendation given in Ref. \cite{chilingarov}, best results are obtained by choosing the frequency where the steepest slope in the $C^{-2}$ plot is achieved (\SI{531}{Hz} in \cref{DepComp1e14}). Both beforehand mentioned effects result in a decrease of this slope but it is not possible to give a fixed number for this as the cutoff frequency depends on the temperature and sensor fluence (both shift the cutoff frequency to lower values by increasing the bulk resistivity). Nonetheless, at all frequencies an influence of either of the two effects can be observed and, when applying this method, the depletion width as well as the effective doping concentration will always be overestimated.
\begin{figure}
	\includegraphics[width=0.9\textwidth]{CV21e14.png}
	\caption{Inverse squared plot of the capacitance of the sensor irradiated to \SI{1e14}{\text{n}_\text{eq}} at \SI{-20}{\celsius}. The red line indicates the capacitance value of a fully depleted sensor of this geometry.}
	\label{DepComp1e14}
\end{figure}

\subsubsection{Annealing Comparison}
A process used quite often to treat irradiated sensors is annealing. Hereby, the sensor is stored a fixed amount of time at elevated temperature to enable different kinds of defects to rearrange their microscopic structure and therefore also change the electrical properties. To see how large the effect of this treatment is, an irradiated ATLAS R5 sensor with a fluence of \SI{7e14}{\text{n}_\text{eq}} was annealed for \SI{80}{min} at \SI{60}{\celsius} which is expected to bring the sensor to the state of full beneficial annealing. The result can be seen in \cref{AnnComp}.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{AnnComp.png}
	\caption{Depletion of an irradiated ATLAS R5 sensor before and after annealing of \SI{80}{min} at \SI{60}{\celsius}.}
	\label{AnnComp}
\end{figure}

After annealing, the sensor has a higher depletion throughout the whole voltage ramp and a depletion voltage of roughly \SI{1000}{V} (active sensor thickness of \SI{305}{\micro\meter}) can be deducted while before annealing the sensor did not reach full depletion at \SI{1100}{V}. From this it can be concluded that the effective doping concentration was reduced during annealing. However, in contrast to this, in \cref{ResComp} it was found that the bulk resistivity does not change with annealing. This indicates that the mobility based correlation between the effective doping concentration and the bulk resistivity \cite{dopres} is not valid anymore for irradiated silicon.

\subsubsection{Bias Resistance Comparison}
The value of the bias resistor gained from the fits can be compared to either the values given by the sensor manufacturer or direct IV measurements between the bias pad and a single DC pad of a strip as described in Ref. \cite{TechSpec}. Before comparing the value from the fits, it has to be multiplied with the number of strips as effectively all the bias resistors are measured in parallel and therefore a much smaller value is determined. The following table shows example results of fitted bias resistance values and their corresponding comparison values.

\begin{center}
	\begin{tabular}{c||c c}
		           &  CMOS 1e15 &  ATLAS R5 7e14 \\ [0.5ex]
		\hline
		Fit value (\SI{-30}{\celsius}) [\si{\mega\ohm}]  & $8.2\pm0.2$ & $1.42\pm0.12$ \\
		\hline
		Comparison [\si{\mega\ohm}]  & $8.7\pm0.2$ & $1.5\pm0.5$\\
		\label{ResTab}
	\end{tabular}
	\captionof{table}{Bias resistance comparison for CMOS and ATLAS R5 sensors.}
	\label{biasrestable}
\end{center}

For the CMOS sensor the bias resistance is directly measured before irradiation according to Ref. \cite{TechSpec} and the value of the ATLAS R5 sensor was compared to the value given by the manufacturer Hamamatsu. It is observed that sensor annealing as well as using capacitance or resistance data for the fit does not change the result. In contrast to this, a resistance increase with decreasing temperature of $\approx -\frac{1\%}{1 K}$ is found. For irradiated sensors, the method described in Ref. \cite{TechSpec} is not applicable anymore as the irradiation causes a leakage current during depletion of the surface region. This current affects the current between the DC and the bias pad and therefore the measured resistance. To investigate whether irradiation effectively changed the bias resistance the introduced method can be used. As it is explained in \cref{TheoModel}, the most straight-forward way is to apply the highest frequency of the LCR meter as, at high frequencies, the measured resistance is equal to the bias resistance.

\subsubsection{Bulk Resistivity Comparison}
\label{ResComp}
As introduced in \cref{TheoModel}, one general assumption is a change of bulk resistivity for sensors after irradiation, to be more exact, the resistivity of the depleted bulk prior to irradiation is assumed to be infinite and expected to decrease. In contrast the resistivity of the non depleted bulk is determined via the doping level and, for most sensors, given by the manufacturer. It is expected to increase as the radiation damage introduces defects to the silicon. The temperature dependent results of sensors and fluences are shown in \cref{RhoComp}. At first sight, a clear temperature dependence can be observed for the undepleted as well as for the depleted silicon. Plotting the resistivity on log scale vs. temperature, it is well in agreement with a linear trend but also the Arrhenius plot seems to be a valid approach. From an Arrhenius fit, an effective activation energy of \SI{0.60(2)}{\eV} can be determined. This is in good agreement with the position of deep lying defects within the silicon band gap. The resistivity prior to irradiation (\SI{3300}{\ohm\centi\meter}) does not show significant temperature dependence and can be determined from the effective doping concentration which itself can be derived from the full depletion voltage. The same approach applied to the full depletion voltages derived from \cref{DepCompCMOS} yields values between \SI{1.6}{\kilo\ohm\centi\meter} for \SI{1e15}{\text{n}_\text{eq}} and \SI{0.4}{\kilo\ohm\centi\meter} for \SI{1e15}{\text{n}_\text{eq}}. These values are orders of magnitudes lower than the resistivity derived by the fits. From this, it can be concluded that the resistivity increase due to irradiation is caused by deep defects and vastly stronger than the increase caused by shallow defects which only affect the effective doping concentration.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{RhoComp.png}
	\caption{Resistivity of the depleted and undepleted bulk for different sensors, temperatures and fluences. For the ATLAS R5 sensors also the results prior to (red and purple) and post annealing (green and brown) are displayed.}
	\label{RhoComp}
\end{figure}

Additionally, the influence of annealing to the resistivity was investigated. For the undepleted bulk, no change can be observed while the depleted bulk shows a resistivity increase with annealing.\\
The expectation for the fluence dependency cannot fully be confirmed. The resistivity of the non depleted bulk shows a small yet significant increase for the ATLAS sensors but no significant change within the CMOS sensors. In contrast the resistivity of the depleted bulk shows a strong decrease with increasing fluence and can be seen especially for the CMOS sensors.

\subsubsection{Investigation of the LF Parameters}

As previously mentioned in \cref{TheoModel}, the free parameters $c_1$ and $c_2$ are introduced individually for each voltage in the fit. However, only if these values exhibit a clear signature are they likely to have a physical counterpart. Thus the voltage, temperature and fluence dependence of the free parameters of the capacitance and resistance fits are investigated in the following. Ideally, they are constant or directly dependent of another quantity.\\

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{c2_new.png}
	\caption{Voltage and temperature dependence of the free parameter $c_2$ derived from capacitance data of the ATLAS R5 sensor irradiated to \SI{7e14}{\text{n}_\text{eq}}. The inset shows the same data divided by the square root of the bulk resistivity to eliminate the temperature dependency.}
	\label{c2}
\end{figure}

At first, the voltage and temperature dependence of the parameter $c_2$ derived from capacitance data of the ATLAS R5 sensor irradiated to \SI{7e14}{\text{n}_\text{eq}} is shown in \cref{c2}. For very low voltages, values of the order of the ones obtained for the unirradiated sensors are determined but with increasing voltage, a rapid decrease followed by a saturation above \SI{100}{V} occurs. Above that threshold, there is only slight variation which is negligible compared with the errors. These show a similar trend with being small for low voltages and an increase for voltages above \SI{200}{V} and decrease again after roughly \SI{500}{V}. For the temperature dependence, a clear trend can be seen. With decreasing temperature, not only the saturation value increases, but also the errors show a steady increase. The latter one can be explained by the fact that the cutoff frequency moves to lower values with decreasing temperature and therefore the effective frequency range that can be measured decreases as the lower limit of \SI{20}{Hz} is defined by the hardware limit of the LCR. Thus a further reduction of the temperature would result in a complete inapplicability of this method. The increase with decreasing temperature is proportional to the square root of the bulk resistivity  $\rho$ as shown in the inset of \cref{c2}.\\

This general behaviour of high values decreasing rapidly and saturating was seen for all sensors irradiated to more than \SI{1e14}{\text{n}_\text{eq}}. The CMOS sensor with the lowest fluence showed the largest deviation from this behaviour which indicates that for very low fluences this general behaviour is not valid anymore. From the other sensors it can be concluded that with increasing fluence the saturation level of $c_2$ increases linearly (for example the ATLAS R0 sensor irradiated to \SI{1.6e15}{\text{n}_\text{eq}} shows roughly twice the saturation value than the ATLAS R5). A more detailed analysis of the fluence dependence would require further different fluence levels.\\

In the same way, the $c_2$ parameter values derived from the resistance fits were analysed and are shown in \cref{c2res}. These show more of a steady decrease with increasing voltage than a stable plateau. However, the temperature trend remains the same and follows the $\sqrt{\rho}$ trend as well. The decrease in voltage cannot be related uniquely to any physical quantity alone or a combination of those.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{c2res_new2.png}
	\caption{Voltage and temperature dependency of the free parameter $c_2$ derived from the resistance fit for the ATLAS R5 sensor irradiated to \SI{7e14}{\text{n}_\text{eq}}. The inset shows the same data divided by the square root of the bulk resistivity.}
	\label{c2res}
\end{figure}

Also, the dependency of the $c_1$ parameter was investigated. For all sensors, a similar trend was visible and the results of the capacitance fits of the CMOS sensor irradiated to \SI{5e14}{\text{n}_\text{eq}} are shown in \cref{c1} as an example. A strong decrease over two orders of magnitude and a temperature dependence can be seen. However, the latter one has the opposite dependency than for the $c_2$ parameter as it decreases with decreasing temperature. As already observed in \cref{UnirradDevs}, the $c_1$ parameter is dependent on $\frac{\text{d}I}{\text{d}V}$ i.e. the leakage current increase with the bias voltage. For unirradiated detectors, the consideration of this quantity is sufficient to eliminate any voltage dependencies for $c_1$. Therefore, this measure is also included for irradiated detectors now. To compensate the hereby introduced temperature dependence, the bulk resistivity (see \cref{ResComp}) has to be taken into account. Lastly, the increase in depletion width with the bias voltage $\frac{\text{d}w}{\text{d}V}$ was found to completely eliminate any voltage dependence. The resulting plot for the CMOS sensor with a fluence of \SI{5e14}{\text{n}_\text{eq}} is shown in the inset of \cref{c1}.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{c1_new.png}
	\caption{Voltage and temperature dependency of the free parameter $c_1$ derived from capacitance data for the CMOS sensor irradiated to \SI{5e14}{\text{n}_\text{eq}}. A decrease over multiple orders of magnitude can be seen as well as a slight temperature dependency. The inset shows the same data divided by the current increase with voltage, the bulk resistivity and the depletion increase with voltage.}
	\label{c1}
\end{figure}

Furthermore, the dependencies of the $c_1$ parameter derived from resistance fits were under investigation (\cref{c1res}). In contrast to the capacitance fits the voltage dependency has a positive and decreasing slope. Additionally, there is a much stronger pronounced temperature dependence. The consideration of the bulk resistivity $\rho$ nearly compensates this dependency. Taking other quantities like $\frac{\text{d}I}{\text{d}V}$ or $\frac{\text{d}w}{\text{d}V}$ into account can also not fully explain the voltage behaviour.

\begin{figure}
	\centering
	\includegraphics[width=0.9\textwidth]{c1res_new2.png}
	\caption{Voltage and temperature dependency of the free parameter $c_1$ derived from resistance data for the CMOS sensor irradiated to \SI{5e14}{\text{n}_\text{eq}}. An increase over multiple orders of magnitude can be seen and a strong temperature dependence.}
	\label{c1res}
\end{figure}

The hypothesis of the diffusion capacitance being responsible for the LF capacitance increase can be supported by the fact that $c_1$ is dependent on the derivative of the leakage current with respect to the bias voltage. Most sources describe diffusion capacitance for forward bias and therefore carry out the derivative of the diode equation and thus conclude a linear dependence on the current. Therefore, this dependence cannot be expected for our devices in reverse bias. However, the hypothesis is contradicted by the fact that the derived frequency dependency predicts a frequency-independent low-frequency regime \cite{Sze}.

\section{Conclusion}
A novel approach for interpreting irradiated CV measurements was introduced. For this, an equivalent circuit was constructed, the corresponding impedance calculated and the capacitance a LCR instrument would measure from this was derived. In addition with an experimentally observed term, this produced a fit function able to reproduce the whole frequency range of CV measurements. As a proof of concept, multiple sensors have been measured, fitted and the gained parameter values were compared to expectations or results from other measurements and showed no evidence of model deficiencies.\\
The main new discovery is an analytical description for the current induced low-frequency capacitance behaviour. Hints towards a physical explanation are given. The question whether the diffusion capacitance is responsible for the capacitance increase cannot be finally answered.
Furthermore, it has been shown why the present technique of extracting the depletion voltage and the effective doping concentration is affected by irradiation and how to choose the least affected frequency.

\bibliographystyle{elsarticle-num}


\section{Introduction}

Quantum computations is a very perspective and dynamically developing area of science.
Great hopes are placed on quantum computers, primarily due to the expected ability to solve problems that are fundamentally impossible to compute on classical computers due to the asymptotic complexity.
Despite the fact, that the main stress is usually placed  on the problems of quantum cryptography \cite{qc_persp}, it is obvious that since the quantum key distribution and post-quantum cryptography is safe \cite{shohr_preskill}, the main field of applications for the quantum computer tends to be the ``peaceful`` modelling of complex systems, according to the original Feynman's idea \cite{Feynman}. Moreover, the current development of the technologies in the field of quantum computers already enables to solve some problems, which are problematic to be solved on the classical computer. The fields of application are also  wide -- from chemistry \cite{quant_chem} to ``industrial`` problems \cite{quant_industrial}.

In computational problems, it becomes necessary to discretize continuous observables.
The representation considered in the paper can be useful for modeling both classical and quantum systems on both classical and quantum computers.

In this work the representation of continuous quantum observables as a row of powers of observables-digits is considered.
We consider the expansions of the coordinate and momentum into a series in terms of the powers of the base of the corresponding numeral systems, and the digit itself also turns out to be observable.
The momentum operator is introduced as a generator of coordinate shifts. As a result the Fourier expansion in this formalism is equivalent to a series expansion in terms of shift operators.
Explicit expressions are obtained for all systems. Moreover, in the work are considered commutation relations, which also turn out not to be trivial.

Among other things, when constructing this representation, a number of interesting ``physical`` effects naturally arise.
For example, we demonstrate the connection between the choice of a binary digit, gauge transformations and the Aaronov-Bohm effect.

As a part of the research work, various renormalization procedures also naturally arose, allowing one to assign finite values to some divergent series, which may be useful in the context of quantum field theory.

This paper presents the generalization of the digital representation \cite{bin-ivanov} and \cite{ter-bin}, for $q = 2$ and $q = 3$ to arbitrary $q \geq 2$. To make it easier to compare the results of previously published articles to the newest one, the style and contest retained where possible, which sometimes results in blocks of self-citation.
\section{Lattice definition}

To make the notation more compact, we consider here and hereafter the system of units, where the Planck constant is equal to one. In other words, we assume $h = 1$ and hence $1/\hbar = 2\pi$. While speaking about quantum observables we are going to concentrate primary on coordinate and momentum, following the notation of the original articles \cite{bin-ivanov}, \cite{ter-bin}, but it is clear, that we can think about them as about any pair of Fourier-conjugated quantum observables, e.g. time and energy.

\subsection{Coordinate lattice}

We assume that the coordinate is described by $q$ digit-qudit and the coordinate lattice consists of $N = q^n$ nodes, which we assume to be cyclic (after the last one comes the first). If we suppose that the \textit{coordinate lattice constant} is $\Delta x = q^{-n_{-}}$, then the \textit{lattice period} is equal to $\Xi = N\Delta x = q^{n_{+}}$, $n_{+} = n - n_{-}$. We assume that the values $x$ range from 0 to $\Delta x \cdot (N - 1)$.

On the coordinate lattice, a natural addition operation is induced from $\mathds{Z}_{N}$, for which $x = x + \Xi$. It is possible to use other representations of the lattice $\Delta x \cdot \mathds{Z}_{N}$ by real numbers. For example, taking the equivalence of $x$ and $x + \Xi$ into account, we further need a representation in which $x$ ranges from $-\Delta x \cdot (N-1)/2$ to $+\Delta x \cdot (N-1)/2$. The power series for the coordinate on a finite lattice is finite:
\begin{equation}
    x = \sum\limits_{s = -n_{-}}^{n_{+} - 1}x_{s}q^{s} = \sum\limits_{s = -n_{-}}^{n_{+} - 1}\mathbf{d}(s,x)q^{s}.
\end{equation}
Here, $x_{s} = \mathbf{d}(s,x)$ is the $s$-th digit in the {\bf d}igital expansion of $x$. We sometimes specify a range of powers of three that defines a lattice and write $x_{s} = \mathbf{d}_{n_{-}n_{+}}(s,x)$.

We introduce the coordinate basis $\{|x\rangle \}_{x \in \Delta x \cdot \mathds{Z}_{N}}$ for functions defined on the lattice:
\begin{equation}
    \hat{x}|x\rangle = x|x\rangle, \quad \langle x^{\prime}|x^{\prime \prime} \rangle = \delta_{x^{\prime},x^{\prime \prime}}, \quad \psi(x) = \langle x|\psi \rangle, \quad x \in \Delta x \cdot \mathds{Z}_{N}.
\end{equation}
We represent wave functions (ket vectors) in the forms of columns whose rows are ordered in decreasing order of $x$. Thus, if $x$ varies from $0$ to $(N - 1)\Delta x$, then
\begin{equation}\psi (x) =\left(
    \begin{array}{c}
           \psi((N - 1)\Delta x)\\
           \psi((N - 2)\Delta x)\\
           .\\
           .\\
           .\\
            \psi(\Delta x)\\
            \psi(0)\\
    \end{array}\right).
\end{equation}

\subsection{Momentum lattice.}
We define the momentum operator $\hat{p}$ as the generator of the cyclic shifts $\widehat{T}_{A}$ along the coordinate lattice:
\begin{equation}
    \widehat{T}_{A} \psi(x) = \psi(x + A), \qquad \widehat{T}_{A} =\mathrm{e}^{2\pi \mathrm{i} A \hat{p}}, \quad A \in \Delta x \cdot \mathds{Z}.
\end{equation}
Such operators were considered on Weyl's classic book \cite{Weyl} and more detailed later by Schwinger \cite{Schwinger}.

Because the coordinate lattice is periodic, the shift by the period $\Xi$ must be identity transformation, i.e., for eigenvalues of operator $\hat{p}$, we have
$\Xi\cdot p\in\mathds{Z}$.
This gives the \textit{momentum step} $\Delta p$,
\begin{equation}
    \Xi \cdot \Delta p = 1, \quad \Delta p = q^{-n_{+}}, \quad \Delta p \cdot \Delta x = \frac{1}{N} = q^{-n}
\end{equation}
The number of points in the spectrum of momentum is the same as for coordinate,
i.e. for momentum,
we have a periodic lattice with the same number of nodes but a different period $\Pi = \Delta p \cdot N = q^{n_{-}}$, $\Pi \Xi = N$.
The momentum lattice is denoted by $\Delta p \cdot \mathds{Z}_{N}$.
The power series for the momentum is also finite:
\begin{equation}
    p = \sum\limits_{r = -n_{+}}^{n_{-} -1}p_{r}q^{r} = \sum\limits_{r = -n_{+}}^{n_{-} -1} \mathbf{d}(r,p)q^r.
\end{equation}
Here, $p_r=\mathbf{d}(r,p)$ is the $r$-th digit in the \textbf{d}igital expansion of $p$.
We sometimes specify a range of powers of $q$ that defines a lattice and write $p_r = \mathbf{d}_{n_+n_-}(r,p)$.

\subsection{Minimum shift.}
The minimum shift $\widehat{T}_{\Delta x}$ is a shift by the lattice step $\Delta x$; any other shift on a given lattice is a power $\widehat{T}_{A} = (\widehat{T}_{\Delta x})^{A/\Delta x}$, where $A/\Delta x \in \mathds{Z}_{N}$:
\begin{equation}\label{T1}
    \widehat{T}_{A}\psi(x) = \psi(x + A), \, \widehat{T}_A|x\rangle = |x - A\rangle, \,
    \langle x^{\prime}|\widehat{T}_{A}|x^{\prime \prime}\rangle = \delta_{x^{\prime},x^{\prime \prime }- A} = \delta_{x^{\prime} + A,x^{\prime \prime}}.
\end{equation}
Moreover,
\begin{equation}\label{T2}
    \hat T_{\Delta x}\psi(x)=\hat T_{\Delta x}
\left(\begin{array}{c}
\psi((N-1)\Delta x)\\
\psi((N-2)\Delta x)\\
\vdots\\
\psi(2\Delta x)\\
\psi(\Delta x)\\
\psi(0)
\end{array}\right)=
\left(\begin{array}{c}
\psi(0)\\
\psi((N-1)\Delta x)\\
\psi((N-2)\Delta x)\\
\vdots\\
\psi(2\Delta x)\\
\psi(\Delta x)
\end{array}\right)=\psi(x+\Delta x).
\end{equation}
The sum $x + \Delta x$ is taken in the sense $x \in \Delta x \cdot \mathds{Z}_{N}$, i.e., this is a cyclic shift of the function on the lattice down one position.

The eigenvalues of the minimum shift operator are $N$th roots of unity and are related to the eigenvalues of the momentum operator (which has not yet been introduced explicitly):
\begin{equation}
    \lambda^{N} = 1, \qquad \lambda_{p} = \mathrm{e}^{2\pi \mathrm{i} \Delta x p} =  \mathrm{e}^{2\pi \mathrm{i} \Delta x \Delta p p/\Delta p} = (\lambda_{\Delta p})^{p/\Delta p},
\end{equation}
where we take $\Delta x \Delta p = 1/N$ and $p/\Delta p \in \mathds{Z}_{N}$ into account and
\begin{equation}
    \Lambda = \lambda_{\Delta p} = \mathrm{e}^{2\pi \mathrm{i} \Delta p \Delta x} = \mathrm{e}^{2\pi \mathrm{i} /N}.
\end{equation}
The corresponding eigenvalues are obtained from the relation $\psi(x) = \widehat{T}_x\psi(0)$. The normalized eigenvectors have the forms
\begin{eqnarray}
    \psi_{\lambda_{p}}(x^{\prime}) &=& \langle x^{\prime}|\psi_{\lambda_{p}}\rangle = \frac{\lambda_{p}^{x^{\prime}/\Delta x}}{\sqrt{N}} = \frac{\mathrm{e}^{2\pi \mathrm{i} x^{\prime} p}}{\sqrt{N}},\\
    \nonumber
    \langle \psi_{\lambda_{p}}|x^{\prime \prime}\rangle &=& \langle x^{\prime \prime}|\psi_{\lambda_{p}}\rangle^{*} = \frac{\lambda_{p}^{-x^{\prime \prime}/\Delta x}}{\sqrt{N}} = \frac{\mathrm{e}^{-2\pi \mathrm{i} x^{\prime \prime}}}{\sqrt{N}}.
\end{eqnarray}
We can write the projector on the (one-dimensional) eigensubspace of the operator $\widehat{T}_{\Delta x}$ as
\begin{eqnarray}
    \widehat{P}_{\lambda_{p}} &=& |\psi_{\lambda_{p}}\rangle\langle\psi_{\lambda_{p}}|,\\ \nonumber
    \langle x^{\prime}|\widehat{P}_{\lambda_{p}}|x^{\prime \prime} \rangle &=& \frac{\lambda_{p}^{(x^{\prime} - x^{\prime \prime})/\Delta x}}{N} = \frac{\lambda_{p}^{d/\Delta x}}{N} = \frac{\mathrm{e}^{2\pi \mathrm{i} p d}}{N}, \quad d = x^{\prime} - x^{\prime \prime}.
\end{eqnarray}
In this notation $x^{\prime}$ labels rows of matrix, and $x^{\prime\prime}$ labels columns.

The eigenstates of the minimum shift operator are also eigenstates of the momentum operator and can be hence written differently:
\begin{equation}
    |\psi_{\lambda_{p}}\rangle = |\psi_{p}\rangle = |p\rangle, \quad \langle x | p \rangle = \frac{\mathrm{e}^{2\pi \mathrm{i} x p}}{\sqrt{N}}.
\end{equation}

\subsection{Group of shifts.}

We make a trivial remark that might nevertheless be of some interest for an arbitrary positional number system.
We constructed the momentum operator such that it generates a symmetry group with respect to the shifts of the coordinate lattice by an integer number of nodes, i.e. a group isomorphic to the group (with respect to addition) of the residues modulo division by $N:\, \Delta x \cdot \mathds{Z}_{N} \approx \mathds{Z}_{N}$. But we can consider unitary operators of the form $\widehat{T}_{A} = e^{2\pi i A \hat{p}}, \, A \in \mathds{R}.$ Such operators correspond to the cyclic shifts by an arbitrary value (not necessarily a multiple of $\Delta x$). The corresponding group is isomorphic to the group $\mathds{R} / \Xi \approx SO(1) \approx U(1)$ of rotations of a circle by an arbitrary angle. Addition is again understood modulo $\Xi$, ($A = A + \Xi$). In the case $\Xi = \infty$, the group of symmetries coincides with the group $\mathds{R}$ of real numbers with respect to addition.

We see that if the Hamiltonian on the lattice is expressed in terms of the operator $\hat{p}$, then the presence of the lattice does not violate translation invariance under arbitrary translations (not necessarily by an integer number of lattice sites), but the operator $\hat{p}$ (as we see below) turns out to be nonlocal, i.e., matrix elements $\langle x^{\prime} | \hat{p} | x^{\prime\prime} \rangle$ can be nonzero for arbitrary large values $x^{\prime} - x^{\prime \prime}$ (in the lattice).

We can specify a state $|x_{0}\rangle = \widehat{T}_{x_{0}}|0\rangle$ with an arbitrary value $x_{0} \notin \Delta x \cdot \mathds{Z}_{N}$, but such a state is not a state with a certain value of the coordinate, because it decomposes into several basic states $\{ |x\rangle \}_{x \in \Delta x \cdot \mathds{Z}_{N}}$.

\section{Operators of digits and their decomposition by shifts}

We defined the momentum operator such that the Fourier harmonic of the momentum is given by the operator $\widehat{T}_{A} = \mathrm{e}^{2\pi \mathrm{i} A \hat{p}}$ of the coordinate shift. Therefore, if we take Fourier transform for the momentum digits
\begin{equation}
    \mathbf{d}_{n_{+}n_{-}}(r,\hat{p}) = \sum\limits_{A \in \Delta x \cdot \mathds{Z}_{N}}\tilde{\mathbf{d}}_{n_{+}n_{-}}(r, A)\mathrm{e}^{2\pi \mathrm{i} A \hat{p}},
\end{equation}
then we obtain the decomposition of the momentum digit by coordinate shifts%
\begin{equation}
   \mathbf{d}_{n_{+}n_{-}}(r, \hat{p}) = \sum\limits_{A \in \Delta x \cdot \mathds{Z}_{N}} \tilde{\mathbf{d}}_{n_{+}n_{-}}(r, A) \widehat{T}_{A}.
\end{equation}

\subsection{Classical positional numeral systems}

\subsubsection{Base-q numeral system}

Let's consider the numeral system with digits $\{0,1, \hdots, q\}$. After computations we obtain the following expression for decomposition of the momentum digit over coordinate shifts:
\begin{equation}
\hat{p}_{r} = \frac{q - 1}{2} \hat{1} - \frac{\Delta p}{ q^{r}} \sum\limits_{D \in \mathbb{Z}(q^r/\Delta p)} \sum\limits_{\sigma = 1}^{q - 1}\frac{\widehat{T}_{-A}}{1 - \exp(2\pi \mathrm{i} \Delta p A)}, \, A = q^{-r}(D + \sigma/q).
\end{equation}
For decomposition of the coordinate in this case we have:
\begin{equation}
\hat{x}_{s} = \frac{q - 1}{2} \hat{1} - \frac{\Delta x}{ q^{s}}   \sum\limits_{D \in \mathbb{Z}(q^s/\Delta x)} \sum\limits_{\sigma = 1}^{q - 1}\frac{\widehat{S}_{B}}{1 - \exp(2\pi \mathrm{i} \Delta x B)}, \, B = q^{-s}(D + \sigma/q).
\end{equation}
Now let's turn to the examples, considered in previous papers.

\subsubsection{Binary non-symmetric system}

The first system to be considered is the ``classical`` binary system, which was being discussed in the original paper \cite{bin-ivanov}. This system has digits $\{0,1\}$ and will be called \textit{``binary non-symmetric system``}.

\begin{figure}[h]
\label{bns_lattice}
\centering
\begin{tikzpicture}
\draw [->] (-5,0) -- (5,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-2) -- (0,2) node [right] {$\mathbf{b_{ns}}(s,x)$};

\fill[green] (-4,0) circle (2pt);
\fill[green] (-3.5,0) circle (2pt);

\fill[green] (-3,1) circle (2pt);
\fill[green] (-2.5,1) circle (2pt);

\fill[green] (-2,0) circle (2pt);
\fill[green] (-1.5,0) circle (2pt);

\fill[green] (-1,1) circle (2pt);
\fill[green] (-0.5,1) circle (2pt);

\fill[green] (0,0) circle (2pt);
\fill[green] (0.5,0) circle (2pt);

\fill[green] (1,1) circle (2pt);
\fill[green] (1.5,1) circle (2pt);

\fill[green] (2,0) circle (2pt);
\fill[green] (2.5,0) circle (2pt);

\fill[green] (3,1) circle (2pt);
\fill[green] (3.5,1) circle (2pt);

\draw[blue] (0,0) circle (2pt);
\draw[blue] (0.5,0) circle (2pt);
\draw[blue] (1,0) circle (2pt);
\draw[blue] (1.5,0) circle (2pt);

\draw[blue] (2,1) circle (2pt);
\draw[blue] (2.5,1) circle (2pt);
\draw[blue] (3,1) circle (2pt);
\draw[blue] (3.5,1) circle (2pt);

\draw[blue] (-0.5,1) circle (2pt);
\draw[blue] (-1,1) circle (2pt);
\draw[blue] (-1.5,1) circle (2pt);
\draw[blue] (-2,1) circle (2pt);

\draw[blue] (-2.5,0) circle (2pt);
\draw[blue] (-3,0) circle (2pt);
\draw[blue] (-3.5,0) circle (2pt);
\draw[blue] (-4,0) circle (2pt);

\end{tikzpicture}
\caption{Plot of the value of the binary digit number $s$ on a lattice for an ``non-symmetric system'', ($n_{-} = 1$), $s = 0$ -- green filled circles, $s = 1$ -- blue circles.}
\end{figure}

In this case we obtain the following result:

\begin{equation}
\hat{p}_r = \frac{1}{2} \hat{1} - \Delta p \, 2^{-r} \; \sum\limits_{D \in \mathds{Z}_{2^{r}/\Delta p}} \frac{\widehat{T}_{-A}}{1 - \exp(2\pi \mathrm{i} \Delta p A)}, \, A=2^{-r}(D + 1/2).
\end{equation}

\subsubsection{Ternary non-symmetric system}
In paper \cite{ter-bin} was introduced the ternary system with digits $\{0,1,2\}$, called \textit{``ternary non-symmetric system``}.

\begin{figure}[h!]
    \centering
\begin{tikzpicture}
\draw [->] (-6,0) -- (6,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-0.5) -- (0,3) node [right] {$\mathbf{t_{ns}}(s,x)$};

\fill[blue] (-6,0) circle (2pt);
\fill[blue] (-5.666,0) circle (2pt);
\fill[blue] (-5.333,0) circle (2pt);

\fill[blue] (-5,1) circle (2pt);
\fill[blue] (-4.666,1) circle (2pt);
\fill[blue] (-4.333,1) circle (2pt);

\fill[blue] (-4,2) circle (2pt);
\fill[blue] (-3.666,2) circle (2pt);
\fill[blue] (-3.333,2) circle (2pt);

\fill[blue] (-3,0) circle (2pt);
\fill[blue] (-2.666,0) circle (2pt);
\fill[blue] (-2.333,0) circle (2pt);

\fill[blue] (-2,1) circle (2pt);
\fill[blue] (-1.666,1) circle (2pt);
\fill[blue] (-1.333,1) circle (2pt);

\fill[blue] (-1,2) circle (2pt);
\fill[blue] (-0.666,2) circle (2pt);
\fill[blue] (-0.333,2) circle (2pt);

\fill[blue] (0,0) circle (2pt);
\fill[blue] (0.333,0) circle (2pt);
\fill[blue] (0.666,0) circle (2pt);

\fill[blue] (1,1) circle (2pt);
\fill[blue] (1.333,1) circle (2pt);
\fill[blue] (1.666,1) circle (2pt);

\fill[blue] (2,2) circle (2pt);
\fill[blue] (2.333,2) circle (2pt);
\fill[blue] (2.666,2) circle (2pt);

\fill[blue] (3,0) circle (2pt);
\fill[blue] (3.333,0) circle (2pt);
\fill[blue] (3.666,0) circle (2pt);

\fill[blue] (4,1) circle (2pt);
\fill[blue] (4.333,1) circle (2pt);
\fill[blue] (4.666,1) circle (2pt);

\fill[blue] (5,2) circle (2pt);
\fill[blue] (5.333,2) circle (2pt);
\fill[blue] (5.666,2) circle (2pt);

\draw[red] (0,0) circle (2pt);
\draw[red] (0.333,0) circle (2pt);
\draw[red] (0.666,0) circle (2pt);
\draw[red] (1,0) circle (2pt);
\draw[red] (1.333,0) circle (2pt);
\draw[red] (1.666,0) circle (2pt);
\draw[red] (2,0) circle (2pt);
\draw[red] (2.333,0) circle (2pt);
\draw[red] (2.666,0) circle (2pt);

\draw[red] (3,1) circle (2pt);
\draw[red] (3.333,1) circle (2pt);
\draw[red] (3.666,1) circle (2pt);
\draw[red] (4,1) circle (2pt);
\draw[red] (4.333,1) circle (2pt);
\draw[red] (4.666,1) circle (2pt);
\draw[red] (5,1) circle (2pt);
\draw[red] (5.333,1) circle (2pt);
\draw[red] (5.666,1) circle (2pt);

\draw[red] (-0.333,2) circle (2pt);
\draw[red] (-0.666,2) circle (2pt);
\draw[red] (-1,2) circle (2pt);
\draw[red] (-1.333,2) circle (2pt);
\draw[red] (-1.666,2) circle (2pt);
\draw[red] (-2,2) circle (2pt);
\draw[red] (-2.333,2) circle (2pt);
\draw[red] (-2.666,2) circle (2pt);
\draw[red] (-3,2) circle (2pt);

\draw[red] (-6,1) circle (2pt);
\draw[red] (-3.333,1) circle (2pt);
\draw[red] (-3.666,1) circle (2pt);
\draw[red] (-4,1) circle (2pt);
\draw[red] (-4.333,1) circle (2pt);
\draw[red] (-4.666,1) circle (2pt);
\draw[red] (-5,1) circle (2pt);
\draw[red] (-5.333,1) circle (2pt);
\draw[red] (-5.666,1) circle (2pt);

\end{tikzpicture}
    \caption{A plot of the value of the ternary digit number $s$ for a finite lattice ($n_{-} = 1$), $s=0$ -- blue filled circles, $s=1$ -- red circles.}
\end{figure}
For this system we easily get the following expression:
\begin{equation}
    \hat{p}_r =  \Delta p\, 3^{-r}\sum\limits_{D \in \mathds{Z}_{3^r/\Delta p}} \sum\limits_{\sigma = 1}^{2}  \frac{(-1)^{D + \sigma}}{2 \mathrm{i} \sin\left(\pi\,\Delta p \,A\right)}\widehat{T}_{-A}, \,
    A=3^{-r}(D + \sigma/3).
\end{equation}
\subsection{Shifted positional systems}

\subsubsection{Base-q numerical system on the shifted lattice}
Let's now consider the positional system with digits $\{d_1, d_1 + 1, \hdots, d_1 + q - 1\}$, here $d_1$ may be non-integer. The shift of the least digit induces the shift of the lattice. Indeed, the ``least`` number on such lattice $(d_1d_1\hdots d_1, d_1 \hdots d_1)$ can be represented as following expression:
\begin{equation} \label{renorm}
	\sum\limits_{s = -n_{-}}^{n_{+} - 1 } q^{s} d_{1} \,  (\text{mod} \, q^{n_{+}})= -q^{-n_{-}} d_{1} = -d_{1} \Delta x,
\end{equation}
which means that we obtain the shifted on $-\Delta x d_1$ lattice and in general case zero is not obliged to appear the node of the lattice.
For this case we can write the general expression for the momentum digit:
\begin{equation}
\hat{p}_r = \frac{d_1 + q - 1}{2} \; \hat{1} - \frac{\Delta p}{q^{r}} \sum\limits_{D \in \mathds{Z}_{q^r/\Delta p} } \sum\limits_{\sigma = 0}^{q - 1} \frac{\widehat{T}_A }{1 - \exp (2 \pi \mathrm{i} \Delta p  \, A)}  \exp(- 2 \pi \mathrm{i}  \Delta p \, A d_1),
\end{equation}
$\, A = q^{-r}(D + \sigma/q),$ and the following expression for the coordinate digit:
\begin{equation}
\hat{x}_s = \frac{d_1 + q - 1}{2} \; \hat{1} - \frac{\Delta x}{q^{s}} \sum\limits_{D \in \mathds{Z}_{q^s/\Delta x} } \sum\limits_{\sigma = 0}^{q - 1} \frac{\widehat{T}_B }{1 - \exp (2 \pi \mathrm{i} \Delta x  \, B)}  \exp(- 2 \pi \mathrm{i}  \Delta p \, B d_1),
\end{equation}
$\, B = q^{-s}(D + \sigma/q).$
\subsection{Binary system with arbitrary digit}
Hence, for the binary system (q = 2) with arbitrary digits $\{d_{1}, d_{2} = d_{1}+1 \}$ we get the following expression:
\begin{eqnarray} \label{gen_f}
	\hat{p}_r  = \frac{d_{1} + d_{2}}{2} \, \cdot \, \hat{1} - \frac{\Delta p}{2^{r}} \; \sum\limits_{D \in \mathds{Z}_{2^{r}/\Delta p}} \frac{\widehat{T}_{-A}}{1 - \exp(2\pi \mathrm{i} \Delta p A)} \exp(- 2 \pi \mathrm{i} \Delta p A d_{1}) =  \\
=  \frac{d_{1} + d_{2}}{2} \, \cdot \, \hat{1} + \frac{\Delta p}{2^{r}} \; \sum\limits_{D \in \mathds{Z}_{2^{r}/\Delta p}} \frac{\widehat{T}_{-A}}{2\mathrm{i} \sin(\pi \mathrm{i} \Delta p A)} \exp(- \pi \mathrm{i} \Delta p A \left(2d_{1} + 1\right)) \nonumber
\end{eqnarray}
We will later use the case of such binary systems to demonstrate, that the shifted lattices tend to be ``natural`` in some cases.
\begin{figure}[h]
\centering
\begin{tikzpicture}
\draw [->] (-5,0) -- (5,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-1) -- (0,2) node [right] {$\mathbf{b}(0,x)$};

\fill[blue] (-4,0) circle (2pt);
\fill[blue] (-3.5,0) circle (2pt);

\fill[blue] (-3,1) circle (2pt);
\fill[blue] (-2.5,1) circle (2pt);

\fill[blue] (-2,0) circle (2pt);
\fill[blue] (-1.5,0) circle (2pt);

\fill[blue] (-1,1) circle (2pt);
\fill[blue] (-0.5,1) circle (2pt);

\fill[blue] (0,0) circle (2pt);
\fill[blue] (0.5,0) circle (2pt);

\fill[blue] (1,1) circle (2pt);
\fill[blue] (1.5,1) circle (2pt);

\fill[blue] (2,0) circle (2pt);
\fill[blue] (2.5,0) circle (2pt);

\fill[blue] (3,1) circle (2pt);
\fill[blue] (3.5,1) circle (2pt);

\draw[red] (-3.75,-0.5) circle (2pt);
\draw[red] (-3.25,-0.5) circle (2pt);

\draw[red] (-2.75,0.5) circle (2pt);
\draw[red] (-2.25,0.5) circle (2pt);

\draw[red] (-1.75,-0.5) circle (2pt);
\draw[red] (-1.25,-0.5) circle (2pt);

\draw[red] (-0.75,0.5) circle (2pt);
\draw[red] (-0.25,0.5) circle (2pt);

\draw[red] (0.25,-0.5) circle (2pt);
\draw[red] (0.75,-0.5) circle (2pt);

\draw[red] (1.25,0.5) circle (2pt);
\draw[red] (1.75,0.5) circle (2pt);

\draw[red] (2.25,-0.5) circle (2pt);
\draw[red] (2.75,-0.5) circle (2pt);

\fill[green] (-3.75 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (-3.25 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (-2.75 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (-2.25 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (-1.75 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (-1.25 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (-0.75 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (-0.25 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (0.25 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (0.75 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (1.25 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (1.75 - 0.125,0.5 + 0.25) node[cross=3pt,green]{};

\fill[green] (2.25 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};
\fill[green] (2.75 - 0.125,-0.5 + 0.25) node[cross=3pt,green]{};

\end{tikzpicture}
\caption{Plot of the value of the binary digit number $0$, on a lattice ($n_{-} = 1$) for an ``symmetric system'' ($d_1 = -0.5$) -- red circles,  $d_{1} = -0.25$ -- green crosses, and ''non-symmetric system'' $d_1 = 0$ -- blue circles.}
\label{systems}
\end{figure}

\subsection{Ternary symmetric system}
\begin{figure}[h]
    \centering
\begin{tikzpicture}
\draw [->] (-6,0) -- (6,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-1.5) -- (0,1.5) node [right] {$\mathbf{t_o}(s,x)$};

\fill[red] (-4.333,-1) circle (2pt);
\fill[red] (-4,-1) circle (2pt);
\fill[red] (-3.666,-1) circle (2pt);

\fill[red] (-3.333,0) circle (2pt);
\fill[red] (-3,0) circle (2pt);
\fill[red] (-2.666,0) circle (2pt);

\fill[red] (-2.333,1) circle (2pt);
\fill[red] (-2,1) circle (2pt);
\fill[red] (-1.666,1) circle (2pt);

\fill[red] (-0.666,-1) circle (2pt);
\fill[red] (-1,-1) circle (2pt);
\fill[red] (-1.333,-1) circle (2pt);

\fill[red] (-0.333,0) circle (2pt);
\fill[red] (0,0) circle (2pt);
\fill[red] (0.333,0) circle (2pt);

\fill[red] (0.666,1) circle (2pt);
\fill[red] (1,1) circle (2pt);
\fill[red] (1.333,1) circle (2pt);

\fill[red] (1.666,-1) circle (2pt);
\fill[red] (2,-1) circle (2pt);
\fill[red] (2.333,-1) circle (2pt);

\fill[red] (2.666,0) circle (2pt);
\fill[red] (3,0) circle (2pt);
\fill[red] (3.333,0) circle (2pt);

\fill[red] (3.666,1) circle (2pt);
\fill[red] (4,1) circle (2pt);
\fill[red] (4.333,1) circle (2pt);

\draw [color=blue] (-1.666,-1) circle (2pt);
\draw [color=blue] (-2,-1) circle (2pt);
\draw [color=blue] (-2.333,-1) circle (2pt);
\draw [color=blue] (-2.666,-1) circle (2pt);
\draw [color=blue] (-3,-1) circle (2pt);
\draw [color=blue] (-3.333,-1) circle (2pt);
\draw [color=blue] (-3.666,-1) circle (2pt);
\draw [color=blue] (-4,-1) circle (2pt);
\draw [color=blue] (-4.333,-1) circle (2pt);

\draw [color=blue] (-1.333,0) circle (2pt);
\draw [color=blue] (-1,0) circle (2pt);
\draw [color=blue] (-0.666,0) circle (2pt);
\draw [color=blue] (-0.333,0) circle (2pt);
\draw [color=blue] (0,0) circle (2pt);
\draw [color=blue] (0.333,0) circle (2pt);
\draw [color=blue] (0.666,0) circle (2pt);
\draw [color=blue] (1,0) circle (2pt);
\draw [color=blue] (1.333,0) circle (2pt);

\draw [color=blue] (1.666,1) circle (2pt);
\draw [color=blue] (2,1) circle (2pt);
\draw [color=blue] (2.333,1) circle (2pt);
\draw [color=blue] (2.666,1) circle (2pt);
\draw [color=blue] (3,1) circle (2pt);
\draw [color=blue] (3.333,1) circle (2pt);
\draw [color=blue] (3.666,1) circle (2pt);
\draw [color=blue] (4,1) circle (2pt);
\draw [color=blue] (4.333,1) circle (2pt);

\end{tikzpicture}
\caption{A plot of the value of the ternary digit number $s$ for a finite lattice ($n_{-} = 1$), $s=0$ -- red filled circles, $s=1$ -- blue circles.}
\end{figure}
The other case, which was previously discussed in \cite{ter-bin}, where shifted lattice emerged implicitly is the ternary system with digits $\{-1,0,1\}$, also called \textit{ternary symmetric system}. For this system we can obtain the following expression (explicitly or using the general formula):
\begin{equation}
    \hat{p}_r =  \Delta p\, 3^{-r}\sum\limits_{D \in \mathds{Z}_{3^r/\Delta p}} \sum\limits_{\sigma = 1}^{2}  \frac{(-1)^{D + \sigma}}{2 \mathrm{i} \sin\left(\pi\,\Delta p \,A\right)}\widehat{T}_{-A},\quad
    A=3^{-r}(D + \sigma/3).
\end{equation}

\subsection{Renormalizations of infinite and finite sums}

\subsubsection{Motivation for renormalizing infinite sums}
Let's start with the plot of the digit number $0$ for the binary non-symmetric system:
\begin{figure}[h]
\centering
\begin{tikzpicture}
\draw [->] (-5,0) -- (5,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-2) -- (0,2) node [right] {$x_0$};

\draw[|->,thick,blue] (0,0)--(1,0);
\draw[|->,thick,blue] (1,1)--(2,1);
\draw[|->,thick,blue] (2,0)--(3,0);
\draw[|->,thick,blue] (3,1)--(4,1);

\draw[|->,thick,blue] (-4,0)--(-3,0);
\draw[|->,thick,blue] (-3,1)--(-2,1);
\draw[|->,thick,blue] (-2,0)--(-1,0);
\draw[|->,thick,blue] (-1,1)--(0,1);
\end{tikzpicture}
\caption{A plot of the value of the binary digit number $0$ for binary non-symmetric system on the line}
\end{figure}
The digit number $s$ can be obtained by scaling the entire plot over the $x$-axis on multiplier $2^{s}$. Hence, it is easy to see that for any negative number $x$ there is number $n$, for which all the digits with index greater then $n$ will be equal to 1, which means that the sum
\begin{equation}
 x \sim \sum\limits_{s = -\infty}^{\infty} x_{s} 2^{s},
\end{equation}
is diverging. The problem becomes more explicit when we turn to the binary symmetric system:
\begin{figure}[h]
\centering
\begin{tikzpicture}
\draw [->] (-5,0) -- (5,0) node [below] {$x$};
\draw [-|] (0,0) -- (1,0) node [below] {1};
\draw [-|] (0,0) -- (0,1) node [right] {1};
\draw [->] (0,-2) -- (0,2) node [right] {$x_0$};

\draw[|->,thick,red] (0,-0.5)--(1,-0.5);
\draw[|->,thick,red] (1,0.5)--(2,0.5);
\draw[|->,thick,red] (2,-0.5)--(3,-0.5);
\draw[|->,thick,red] (3,0.5)--(4,0.5);

\draw[|->,thick,red] (-4,-0.5)--(-3,-0.5);
\draw[|->,thick,red] (-3,0.5)--(-2,0.5);
\draw[|->,thick,red] (-2,-0.5)--(-1,-0.5);
\draw[|->,thick,red] (-1,0.5)--(0,0.5);
\end{tikzpicture}
\caption{A plot of the value of the binary digit number $0$ for binary symmetric system on the line}
\end{figure}

in this case we obtain the divergence for both negative and positive numbers.

\subsubsection{Renormalization of infinite sums}
To solve the emerging problem, let's consider the arbitrary positional base-$q$ system with digits $x_s$. As we have seen, for such a system the row
\begin{equation}\label{geom}
\sum\limits_{s = -\infty}^{\infty}x_s q^s,
\end{equation}
does not converge in general case. For this case we can introduce the sum ``with prime``, which is determined in the following way:
\begin{equation}
\sum\limits_{s = 0}^{\infty \, \prime}q^s = \frac{1}{1-q} \, ,
\end{equation}
which results in application of the formula for the sum of converging geometrical progression over the field of its application. Such ``renormalization`` results in the expression
\begin{equation}
\sum\limits_{s \in \mathds{Z}}^{\quad \prime}q^s = 0,
\end{equation}
which is, generally speaking, the alternative definition of this renormalization.

It is possible in some situations to consider the row (\ref{geom}) in $q$-adic sence, as, for instance in \cite{p-adic}, it will result in the convergence of the digits with infinite number of non-zero digits after the point, but the divergence of numbers with infinite number of non-zero digits before the point in the other hand.

The alternative way to renormalize the diverging sum, inspired by the examples considered above, can be represented by the following formal calculation:
\begin{equation}
x = \frac{qx - x}{q-1} = \frac{1}{q - 1}\sum\limits_{s \in \mathds{Z}} \, (x_{s - 1} - x_{s})q^s.
\end{equation}
Obviously, both ways to renormalize the sum fix the problems, which occur when we try to write the number as the row of digits in the positional system. Now we are ready to turn back to the lattice.
\subsubsection{Renormalization on the lattice}
As far as the lattice is finite, we can not obtain the divergent sums on it. Therefore, the purpose of the renormalization is no longer to avoid the divergence, but to \textit{change the representation} of $\mathds{Z}_{N}$ from $\{0,1,\hdots,N-1\}$ to, for instance, $\{-k, -k + 1, \hdots , -k + N - 1\}$.

\begin{figure}[h]
\centering
\begin{tikzpicture}
\coordinate (A) at (1.2,-0.3);
\coordinate (B) at (-0.6,-1.7);

\draw[->, red, dashed] (A) to [bend right=45] (B);

\draw[->](-1,0) -- (4,0) node [right] {$x$};
\draw[|-|, blue, thick](-0.8,0) -- (1.2,0) node [below]{\ScissorLeftBrokenTop};

\draw[|-|, red, thick](1.2,0) -- (3.2,0) node [below]{$\Xi$};

\draw[->](-1,-2) --(4,-2) node [right] {$x$};

\draw[|-|, blue, thick](-0.8,-2) -- (1.2,-2);
\draw[|-|, red, thick](-2.8,-2) -- (-0.8,-2);
\end{tikzpicture}
\end{figure}
To demonstrate the concept more precise, let's consider two examples. Let us start with the binary system. For binary non-symmetric system, we obtain the following expressions:
\begin{equation}
x = \sum\limits_{s = -n_{-}}^{n_{+} - 1~\prime}x_s 2^s =
	\sum\limits_{s = -n_{-}}^{n_{+} - 1} (x_{s-1} - x_s)2^{s}, \quad x_{-n_{-} - 1} = 0,
\end{equation}
such renormalization is equal to change of the last digit  $x'_{n_+-1}=-x_{n_+-1}$:
\begin{equation}
\sum\limits_{s = -n_{-}}^{n_{+} - 1~\prime \prime}x_s 2^s
= \sum\limits_{s = -n_{-}}^{n_{+} - 2}x_s 2^s +x'_{n_+-1}2^{n_+-1}.
\end{equation}
and is linear in accordance to the binary digits $x_s$.

It is noteworthy, that in the ternary case the linear renormalization does not work. Indeed, the linear renormalization is described by the expression:
\begin{eqnarray}
    &&x^{\prime} = \sum\limits_{s = -n_{-}}^{n_{+} - 1~\prime}\mathbf{t_{ns}}(s,x) 3^s =
    \frac{1}{2}\sum\limits_{s = -n_{-}}^{n_{+} - 1} (\mathbf{t_{ns}}(s - 1, x) - \mathbf{t_{ns}}(s,x))3^{s}, \nonumber\\
     &&\mathbf{t_{ns}}(-n_{-} - 1, x) = 0.
\end{eqnarray}
Maximal ternary number is renormalized in the right way:
\begin{equation}
    \sum\limits_{s = -n_{-}}^{n_{+} - 1 ~\prime}2\cdot 3^s = \frac{1}{2} (-2\cdot 3^{-n_{-}}) =-\Delta x,
\end{equation}
but the half is not the node of the lattice:
\begin{equation}
\sum\limits_{s = -n_{-}}^{n_{+} - 1 ~\prime} 3^s =-\frac{\Delta x}{2}.
\end{equation}

The alternative (and properly working) way to obtain the lattice
\\ $\{-3^{n-1}\Delta x,\dots,-\Delta x,0,\Delta x,2\Delta x,\dots,(3^n-3^{n-1}-1)\Delta x\}$ from the initial one ($\{0,\Delta x,2\Delta x,3\Delta x,\dots, (3^n-1)\Delta x\}$) is to substract $3^{n_+}$ from the last $3^{n-1}$ nodes. Then we obtain:
\begin{eqnarray}
  \mathbf{t^{\prime\prime}_{ns}}(n_+-1,x)&=&\left\{\begin{array}{cc}
   0,&\mathbf{t_{ns}}(n_+-1,x)=0,\\
   1,&\mathbf{t_{ns}}(n_+-1,x)=1,\\
   -1,&\mathbf{t_{ns}}(n_+-1,x)=2
  \end{array}\right.=\\ \nonumber
  &=&\mathbf{t_{ns}}(n_+-1,x)-\frac32(\mathbf{t_{ns}}(n_+-1,x)-1)\mathbf{t_{ns}}(n_+-1,x).
\end{eqnarray}
\begin{equation}
x^{\prime\prime} = \sum\limits_{s = -n_{-}}^{n_{+} - 1~\prime\prime}\mathbf{t_{ns}}(s,x) 3^s
= \sum\limits_{s = -n_{-}}^{n_{+} - 2}\mathbf{t_{ns}}(s,x)\, 3^s +\mathbf{t^{\prime\prime}_{ns}}(n_+-1,x)\, 3^{n_+-1}.
\end{equation}

Such renormalization works, but is non-linear  to the ternary digits $x_s$.

\subsubsection{Digit on the line}
After defining the infinite sums of digits, we can discuss the limit of $n \rightarrow +\infty$, $\Delta p \rightarrow 0$. In this case we obtain the following expression for the momentum digit on the line:
\begin{equation}
\hat{p}_{r} = \frac{d_1 + q + 1}{2} \; \hat{1} + q^{-r} \sum\limits_{D \in \mathds{Z}_{q^r/\Delta p} } \sum\limits_{\sigma = 0}^{q - 1}  \frac{\widehat{T}_A}{2 \pi \mathrm{i} A}, \; A = q^{-r}(D + \sigma/q).
\end{equation}
This limit corresponds to the infinite number of digits before the point for the coordinate and after the point for the momentum. Hence, we have the infinite coordinate lattice, while the momentum lattice is periodical (and if we additionally consider the limit $\Delta x \rightarrow 0$ -- then we obtain the real lines for both coordinate and momentum). It is noteworthy, that after transition to the limit the sum is taken over both positive and negative integers. On the finite periodical lattice the concept of positive and negative number was not defined.

\section{Commutation relations}
\subsection{Digit-digit commutator}
It is easy to derive the commutation relation between an arbitrary function of the coordinate $f(\hat{x})$ and the shift operator $\widehat{T}_{A}$:
\begin{equation}
    [f(\widehat{x}), \widehat{T}_{A}]\psi(x) = (f(\widehat{x}) - f(\widehat{x} + A))\widehat{T}_{A}\psi(x).
\end{equation}
Hence,
\begin{equation}
    [f(\widehat{x}), \widehat{T}_{A}] = (f(\widehat{x}) - f(\widehat{x} + A))\widehat{T}_{A}.
\end{equation}
Using notation $d(s, \hat{x}) = \hat{x}_s$, we obtain the form of the commutators of the digits of the coordinate and of the momentum on the lattice:
\begin{eqnarray}
       [\hat{x}_s,\hat{p}_r] =
       - \frac{\Delta p}{q^{r}} \,
       \sum\limits_{D \in \mathds{Z}_{q^r/\Delta p}}
        \sum\limits_{\sigma = 1}^{q-1}
        \frac{\{d(s,\hat{x}) - d(s, \hat{x} - q^{-r}(D + 1/q)\} \widehat{T}_{-A}}{1 - \exp(2 \pi \mathrm{i} \Delta p A)}\times
        \nonumber\\ \times\exp(-2\pi \mathrm{i} A \Delta p d_1).
\end{eqnarray}

As in \cite{bin-ivanov} we have that the requirement for commutation is following:
\begin{equation}
	-r - s - 2 \geq 0, \quad s + r \leq -2.
\end{equation}
Hence, the fractional part of the momentum commutes with the fractional part of the coordinate, the lowest digit of momentum does not commute only with the highest digit of the coordinate, and lowest digit of the coordinate does not commute only with the highest digit of the momentum. The fractional parts of the coordinate and momentum can be considered as the full set of observables for one-dimensional  motion.

\subsection{Coordinate - digit commutator}
Because $\hat x = \sum\limits_{s = -n_{-}}^{n_{+} - 1} q^s d(s, \hat x)$, we obtain the following commutator of coordinate and the digit of momentum:
\begin{equation}
	[\hat x, d(r, \hat{p})] = - \frac{1}{q^r} \sum\limits_{D \in \mathbb{Z}_{q^r/\Delta p}}
	 \sum\limits_{\sigma = 1}^{q-1}
	\frac{\Delta p \cdot q^{-r} (D + 1/q)}{1 - \exp(2 \pi \mathrm{i} \Delta p A)} \widehat{T}_{-A} \exp(-2\pi \mathrm{i} A \Delta p d_1).
\end{equation}
\subsection{Coordinate-momentum commutator}
Similarly, we consider $\hat p = \sum\limits_{r = n_{+}}^{n_{-} - 1} q^r d(r, \hat p)$ and hence obtain the following commutator:
\begin{equation}
	[\hat x, \hat p] = -\sum\limits_{r = -n_{+}}^{n_- - 1}
	\sum\limits_{D \in \mathbb{Z}_{q^r/\Delta p}}
	 \sum\limits_{\sigma = 1}^{q-1}
	 \frac{\Delta p \cdot q^{-r} (D + 1/q)}{1 - \exp(2 \pi \mathrm{i} \Delta p A)} \widehat{T}_{-A} \exp(-2\pi \mathrm{i} A \Delta p d_1).
\end{equation}

\subsection{Commutators on the line}
On the line the commutators become independent of $d_1$ (the phase shift becomes infinitely small). The formulas have the following view:
\begin{equation}
	[d(s,\hat x), d(p, \hat r)] = \sum\limits_{D \in \mathbb{Z}}
	 \sum\limits_{\sigma = 1}^{q-1}
	 \frac{d(s, \hat x) - d(s, \hat x - q^{-r}(D + \sigma/q)}{2\pi \mathrm{i} (D + \sigma/q)} \widehat{T}_{-A},
\end{equation}
\begin{equation}
	[\hat x, d(r, \hat p)] = \frac{1}{q^{r}} \frac{1}{2\pi \mathrm{i}} \sum\limits_{D \in \mathbb{Z}}
	 \sum\limits_{\sigma = 1}^{q-1}
	 \widehat{T}_{-A},
\end{equation}
\begin{equation}\label{commutator}
	[\hat x, \hat p] = \frac{1}{2 \pi \mathrm{i}}\sum\limits_{r \in \mathbb{Z}} \sum\limits_{D \in \mathbb{Z}}
	 \sum\limits_{\sigma = 1}^{q-1}
	 \widehat{T}_{-A}.
\end{equation}

\subsection{Renormalization of the commutator on the line}
We obtain a formal decomposition of the commutator in the sum of the shift operators, as $\hbar = 1/2\pi$, we can rewrite (\ref{commutator}) in the following way:
\begin{equation}\label{commutator_2}
[\hat x, \hat p] = - \mathrm{i} \hbar \sum\limits_{r \in \mathbb{Z}} \sum\limits_{D \in \mathbb{Z}}
	 \sum\limits_{\sigma = 1}^{q-1}
	 \widehat{T}_{-A}.
\end{equation}

Let $\mathds{A}$ be a set of numbers whose $q$-nary expansion contains a finite number of nonzero factors with negative powers of $q$ (a finite number of significant $q$-nary digit after $q$-nary point); $\mathbb{A}$ is a group under the summation operation.
Then the set of values of the shifts along which summation occurs has the form $\mathbb{A} / 0$. Given that $\widehat{T}_0 = \hat{1}$, we obtain:
\begin{equation}
    [\hat{x}, \hat{p}] = -\mathrm{i} \hbar \sum\limits_{A \in \mathds{A} \setminus \{0\}}\widehat{T}_{A} = -\mathrm{i} \hbar \left(\sum\limits_{A \in \mathds{A}}\widehat{T}_{A} - \hat{1} \right) = \mathrm{i} \hbar \hat{1} - \mathrm{i} \hbar \sum\limits_{A \in \mathds{A}}\widehat{T}_{A}.
\end{equation}
We know that for the particle coordinate and momentum on the line, there is the canonical commutation relation $[\hat{x}, \hat{p}] = i\hbar \hat{1}$. We thus obtain the renormalization
\begin{equation}
    \sum\limits_{A \in \mathds{A}}\widehat{T}_{A} = \sum\limits_{A \in \mathds{A}} \mathrm{e} ^{2\pi \mathrm{i} A \hat{p}} = 0.
\end{equation}
This renormalization is similar to the formal equality $\int\limits_{\mathds{R}}\mathrm{e} ^{2\pi \mathrm{i} x p}dx = 0$ for all $p \neq 0$ arising in the Fourier transforms.

\section{Physical motivation for shifted lattice}
\subsection{Boundary conditions}
Let us denote the binary digit number $r$ of the momentum $p$ in the system with digits $\{d_1, d_1 + 1\}$ as $b(r,p,d_1)$. From the equation for the digit in the binary system \eqref{gen_f} we can see, that
\begin{equation} \label{phase_0}
b(r,p,d_1) =  d_1+ b(r,p,0) \cdot \mathrm{e}^{-2\pi\mathrm{i}\Delta p A d_1}.
\end{equation}
Therefore, we can mention that the shift of momentum lattice induces the change of boundary conditions for the coordinate. The boundary conditions for the momentum with digits $\{ d_{1}, d_{1} + 1\}$ can be written in the following form:
\begin{equation} \label{phase}
	\Psi(x + 2^{n_{+}}) = \mathrm{e}^{i\phi} \cdot \Psi(x), \, \text{where} ~ \phi = -2\pi d_{1}.
\end{equation}

\begin{figure}[h!]
\begin{center}
\begin{tikzpicture}
\draw[blue,thick] (0,0) circle (2);
\draw[black,fill=green,thick] (0,2) circle (0.1) node [above] {$\mathrm{e}^{\mathrm{i} \phi}  |\psi \rangle \leftarrow 3d_1 + 4=3d_1 \rightarrow |\psi \rangle$};
\draw[black,fill=red,thick] ({2)},{0}) circle (0.1) node [below right] {$3d_1 + 1$};
\draw[black,fill=red,thick] ({0},{-2}) circle (0.1) node [below left] {$3d_1 + 2$};
\draw[black,fill=red,thick] ({-2},{0}) circle (0.1) node [below left] {$3d_1 + 3$};

\end{tikzpicture}
\caption{The rotation over the lattice period gives an additional phase to the wavefunction.}
\end{center}
\end{figure}
\subsection{Aharonov–Bohm effect}
In \cite{shift} was established a relation between flux of the magnetic field passing through the torus, shift of the generalized momentum and boundary conditions of the wavefunction. We analogically can consider the ring, through which there is non-zero flux of the magnetic field $\Phi$.
\\
The stationary state of a charged particle with the charge $e$ and mass $m$ on this ring of length $a$ is the eigenfunction of the operator
\begin{equation}
	\hat{H}_0 = -\frac{\hbar^2}{2m}\left( \frac{\partial}{\partial x} - \mathrm{i}\frac{e}{c\hbar}A_{x}(x)\right)^2,
\end{equation}
where $\mathbf{A}(x)$ is a vector potential.
\begin{figure}[h!]
\begin{minipage}[h]{0.49\linewidth}
\begin{tikzpicture}

\draw[red] (1.5, 1, 1.5) [y={(0,0,1)}] circle (1.5);
\draw [->, thick, blue] (0,0,-0.3) -- (0,2,1);
\draw [->, thick, blue] (1.5,0,-0.5) -- (2,2,0);
\draw [->, thick, blue] (0.7,0,-0.4) -- (1,3,1);
\end{tikzpicture}
\\ a)
\end{minipage}
\begin{minipage}[h]{0.49\linewidth}
\begin{tikzpicture}
\draw[|-|, thick](0,0) --(5,0) node [below] {$a = 0$};
\draw[|-|](5,0) --(0,0) node [left] {$0$};
\draw[->](0,0) --(6,0) node [right] {$x$};
\end{tikzpicture}
\\ b)
\end{minipage}
\caption{a) Ring in the magnetic field, b) Coordinates on the ring}
\end{figure}
\\
Function $\Psi(x)$ satisfies to periodic boundary conditions:
\begin{equation}
\Psi(0) = \Psi(a),
\end{equation}
\begin{equation}
\Psi^{\prime}(0) = \Psi^{\prime}(a),
\end{equation}
After gauge transformation we can nullify $\mathbf{A}$ on the ring. This will lead to the discontinuity of wavefunction at the cuts. It is also clear that $|\Psi(x)|^{2}$ should not change. Then the problem is simplified to the problem of finding eigenfunctions $\psi(x)$ of an operator
\begin{eqnarray} \label{H_1}
	\hat{H}_1 = -\frac{\hbar^2}{2m}\left(\frac{\partial}{\partial x} \right)^2,
\end{eqnarray}
with the phase shift boundary conditions
\begin{eqnarray} \label{bound_1}
	\psi(0) = \mathrm{e}^{\mathrm{i} \phi} \psi(a), \quad \psi^{\prime}(0) =  \mathrm{e}^{\mathrm{i} \phi} \psi^{\prime}(a).
\end{eqnarray}
The gauge transformation, nullifying the vector potential on the ring can be written in the following form:
\begin{eqnarray} \label{gauge}
\mathbf{A}^{\prime}(x) = \mathbf{A} - \nabla f(x) = 0.
\end{eqnarray}
Because we want the probability $|\Psi(x)|^{2}$ not to be changed by gauge transformation \eqref{gauge}, the new wavefunction $\psi(x)$ differs from the old one $\Psi(x)$ only in a phase factor:
\begin{eqnarray} \label{phase_2}
\psi(x) = \Psi(x)\exp\left(
	-\frac{\mathrm{i} e}{\hbar c} f(x)
\right).
\end{eqnarray}
For Hamiltonian \eqref{H_1} we can easily find the eigenfunctions, fitting the boundary conditions \eqref{bound_1}:
\begin{eqnarray}
	\psi(x) = \mathrm{e}^{\mathrm{i} k x}, \, \text{where}~ k = \frac{2\pi n}{a} - \frac{\phi}{a}.
\end{eqnarray}
Now let's consider the relationship between flux of magnetic field $\Phi$ and phase shift $\phi$.
\begin{eqnarray}
	\Phi =  f(a) - f(0) = \underbrace{\frac{\hbar c}{e}}_{\frac{\Phi_0}{2\pi}} \phi.
\end{eqnarray}
From this example we can see, that in a real physical system a shift of (generalized) momentum induces the shift of wavefunction itself (as in \eqref{phase_0}) and the phase factor in boundary conditions (as in \eqref{phase}). In such cases use of the shifted lattice for momentum is quite natural.

\section{Conclusion}
This paper generalizes the results of previously published works \cite{bin-ivanov} and \cite{ter-bin} for the case of the $q$-base numeral systems. The only restriction on the digits of considered positional systems is that the distance between adjacent digits is $1$.

The possible application for the obtained expansions is not only the quantum computations, but also, for instance, the solution of partial differential equations on classical computer. It worth mentioning that since we define the momentum operator as the shift generator, it appears to be non-local (see appendix A), which helps to better describe the symmetries of considered theory.

The natural emergence of renormalizations is particularly interesting. Renormalizing infinite and finite quantities (on the lattice) allows finding the renormalization numerically by passing from the lattice to the limit of a continuous quantity.
Time and energy can also be considered as a coordinate and a momentum \cite{dyn_time}, which allows applying the same renormalization methods to them. Renormalizations in this context are probably related to the quantum theory of measurements (see \cite{book}, \cite{measurement} and the references therein for the quantum theory of measurements).

The digital representation assumes that the coordinate and the momentum are not observed in the experiment. Instead, individual digits of the coordinate and the momentum are directly observed. The measurement of a binary digit of the spatial coordinate corresponds to the passage/nonpassage of a particle through a diffraction grating.
It's quite interestingly, that the system with non-integer digits appears to be in some particular cases (e.g. for non-trivial boundary conditions) more ``natural`` then the common one. It is also noteworthy that the lattices of such systems induce the constructions, which are quite similar to the windings of torus (see appendix B), which may be applicable to ergodicity.
\section{Acknowledgments}
The authors thank the participants of the Conference Phystech-Quant 2020 (Moscow Institute of Physics and Technology, 2020), of the section of theoretical physics of the 62-nd, 63-rd and 64-th Science Conferences of MIPT (2019, 2020 and 2021 respectively),
the seminar of the Department of mathematical physics (Steklov mathematical institute)
and the seminar of the Laboratory of infinite-dimensional analysis and mathematical physics (Moscow State Univercity). Separately, the authors thank I. V. Volovich, Z. V. Khaidukov, V. A. Dudchenko, V. V. Naumov, N. N. Shamarov, V. Zh. Sakbayev and other colleagues for useful discussion.

\section*{Appendix A}
Let us consider the matrices of the momentum digits and operators in the cases, when they are compact enough to be placed on paper.

Everywhere in this section, $\Delta x = 1, \, x \in \mathds{Z}_{N} \in \{0,1, ... ,N - 1  \}$, coordinates and momenta are numbered by ternary numbers, which are marked with a lower index 3, and $\Delta p = 3^{-n} = 1/N$.
\subsection*{The case $n = 1$ and $N = 3^1 = 3$. Symmetric system.} In this case, we have
$$\hat x = \hat x_0 = \left(
\begin{array}{ccc}
     +1&0&0  \\
     0&0&0  \\
     0&0&-1
\end{array}
\right) = \hat{s}_{z}
$$

$$
\hat p_{-1} =\frac{1}{\sqrt{3}} \left(
\begin{array}{ccc}
    0 & \mathrm{i} & -\mathrm{i} \\
    -\mathrm{i} & 0 & \mathrm{i} \\
    \mathrm{i} & -\mathrm{i} & 0
\end{array}
\right) = \frac{1}{\sqrt{3}} \left(\sqrt{2}\hat{s}_{y} + 2\hat{s}_{y}\hat{s}_{x} + \mathrm{i} \hat{s}_{z} \right)
$$

\subsection*{The case $n = 1$ and $N = 3^1 = 3$. Non-symmetric system.}
Here we obtain following results:
$$\hat x = \hat x_0 = \left(
\begin{array}{ccc}
     2&0&0  \\

     0&1&0  \\
     0&0&0
\end{array}
\right) = \hat{1} + \hat{s}_z; \quad
\hat p_{-1} =\frac{1}{6} \left(
\begin{array}{ccc}
     6& 1 - \sqrt{3}\,\mathrm{i}& 1 + \sqrt{3}\,\mathrm{i}  \\
     1 + \sqrt{3}\,\mathrm{i}& 6 & 1 - \sqrt{3}\,\mathrm{i} \\
     1 - \sqrt{3}\,\mathrm{i}& 1 + \sqrt{3}\,\mathrm{i}& 6
\end{array}
\right) =
$$
$$=\hat{1} +  \frac{\sqrt{2}}{6}\hat{s}_{x} + \frac{1}{\sqrt{6}}\hat{s}_{y} - \frac{\sqrt{3}}{6}(2\hat{s}_{y}\hat{s}_{x} + \mathrm{i} \hat{s}_{z})
$$
\subsection*{The case $n = 2$ and $N = 3^2 = 9$. Non - symmetric system.}
$$x = x_0 + 3\cdot x_1=\text{diag}(8;7;6;5;4;3;2;1;0),$$
$$x_0=\text{diag}(2;1;0;2;1;0;2;1;0),$$
$$x_1=\text{diag}(2;2;2;1;1;1;0;0;0).$$

Let's denote:
$$
    E_{n} =  \frac{1}{\exp\left(\frac{-2 \pi \mathrm{i} n}{9}\right) - 1},
    $$
    Then:
$$
\hat{p}_{-1} = \frac{1}{3}\left(
\begin{array}{ccccccccc}
     3&E_{8} &E_{7} & 0& E_{5}& E_{4}&0 &E_{2} & E_{1}\\
     E_{1} & 3 &E_{8} &E_{7} & 0& E_{5}& E_{4}&0 &E_{2}\\
     E_{2} &E_{1} & 3 &E_{8} &E_{7} & 0& E_{5}& E_{4}&0\\
     0&E_{2} &E_{1} & 3 &E_{8} &E_{7} & 0& E_{5}& E_{4}\\
     E_{4} &0&E_{2} &E_{1} & 3 &E_{8} &E_{7} & 0& E_{5}\\
     E_{5} &E_{4} &0&E_{2} &E_{1} & 3 &E_{8} &E_{7} & 0\\
     0 &E_{5} &E_{4} &0&E_{2} &E_{1} & 3 &E_{8} &E_{7} \\
     E_{7} &0 &E_{5} &E_{4} &0&E_{2} &E_{1} & 3 &E_{8} \\
     E_{8} & E_{7} &0 &E_{5} &E_{4} &0&E_{2} &E_{1} & 3
\end{array}
\right).
$$
   $$
   \hat{p}_{-2} = \left(
   \begin{array}{ccccccccc}
        1&0 &0 &3 E_{6} & 0&0 &3E_{3} &0 &0  \\
        0 & 1&0 &0 &3 E_{6} & 0&0 &3E_{3} &0 \\
        0 & 0& 1&0 &0 &3 E_{6} & 0&0 &3E_{3} \\
        3E_{3} &0 & 0& 1&0 &0 &3 E_{6} & 0&0 \\
        0& 3E_{3} &0 & 0& 1&0 &0 &3 E_{6} &0 \\
        0 & 0& 3E_{3} &0 & 0& 1&0 &0 &3 E_{6}\\
        3 E_{6} &0 & 0& 3E_{3} &0 & 0& 1&0 &0\\
        0 & 3 E_{6} &0 & 0& 3E_{3} &0 & 0& 1&0\\
        0 &0 & 3 E_{6} &0 & 0& 3E_{3} &0 & 0& 1\\
   \end{array}
   \right),
   $$
   $\hat{p} = \frac{1}{3} \hat{p}_{-1} + \frac{1}{9} \hat{p}_{-2}$,
$$
\hat{p} = \frac{1}{9}
\left(
\begin{array}{ccccccccc}
     4&E_{8} &E_{7} & 3 E_{6}& E_{5}& E_{4}&3 E_{3} &E_{2} & E_{1}\\
     E_{1} & 4 &E_{8} &E_{7} & 3 E_{6}& E_{5}& E_{4}&3 E_{3} &E_{2}\\
     E_{2} &E_{1} & 4 &E_{8} &E_{7} & 3 E_{6}& E_{5}& E_{4}&3 E_{3}\\
     3 E_{3}&E_{2} &E_{1} & 4 &E_{8} &E_{7} & 3 E_{6}& E_{5}& E_{4}\\
     E_{4} &3 E_{3}&E_{2} &E_{1} & 4 &E_{8} &E_{7} & 3 E_{6}& E_{5}\\
     E_{5} &E_{4} &3 E_{3}&E_{2} &E_{1} & 4 &E_{8} &E_{7} & 3 E_{6}\\
     3 E_{6} &E_{5} &E_{4} &3 E_{3}&E_{2} &E_{1} & 4 &E_{8} &E_{7} \\
     E_{7} &3 E_{6} &E_{5} &E_{4} &3 E_{3}&E_{2} &E_{1} & 4 &E_{8} \\
     E_{8} & E_{7} &3 E_{6} &E_{5} &E_{4} &3 E_{3}&E_{2} &E_{1} & 4
\end{array}
\right).
$$
\subsection*{The case $n = 2$ and $N = 3^2 = 9$. Symmetric system.}
$$x=x_0 + 3 x_1 =
\text{diag}(4;3;2;1;0;-1;-2;-3;-4),$$

$$x_0=\text{diag}(1;0;-1;1;0;-1;1;0;-1),$$

$$x_1=\text{diag}(1;1;1;0;0;0;-1;-1;-1).$$

Let's denote:
$$
 G_{n} =  \frac{(-1)^n}{2\sqrt{3}\sin\left(\frac{9-n}{9}\right)},
$$
 then:
$$
\hat{p}_{-1} = \frac{1}{\sqrt{3}\mathrm{i}}\left(
\begin{array}{ccccccccc}
     0&G_{8} &G_{7} & 0& G_{5}& G_{4}&0 &G_{2} & G_{1}\\
     G_{1} & 0 &G_{8} &G_{7} & 0& G_{5}& G_{4}&0 &G_{2}\\
     G_{2} &G_{1} & 0 &G_{8} &G_{7} & 0& G_{5}& G_{4}&0\\
     0&G_{2} &G_{1} & 0 &G_{8} &G_{7} & 0& G_{5}& G_{4}\\
     G_{4} &0&G_{2} &G_{1} & 0 &G_{8} &G_{7} & 0& G_{5}\\
     G_{5} &G_{4} &0&G_{2} &G_{1} & 0 &G_{8} &G_{7} & 0\\
     0 &G_{5} &G_{4} &0&G_{2} &G_{1} & 0 &G_{8} &G_{7} \\
     G_{7} &0 &G_{5} &G_{4} &0&G_{2} &G_{1} & 0 &G_{8} \\
     G_{8} & G_{7} &0 &G_{5} &G_{4} &0&G_{2} &G_{1} & 0
\end{array}
\right),
$$
$$
   \hat{p}_{-2} = \frac{1}{\sqrt{3}\mathrm{i}} \left(
   \begin{array}{ccccccccc}
        0&0 &0 &1 & 0&0 &-1  &0 &0  \\
        0 & 0&0 &0 &1 & 0&0 &-1  &0 \\
        0 & 0& 0&0 &0 &1 & 0&0 &-1  \\
        -1 &0 & 0& 0&0 &0 &1 & 0&0 \\
        0& -1  &0 & 0& 0&0 &0 &1 &0 \\
        0 & 0& -1  &0 & 0& 0&0 &0 &1\\
        1 &0 & 0& -1  &0 & 0& 0&0 &0\\
        0 & 1 &0 & 0& -1  &0 & 0& 0&0\\
        0 &0 & 1 &0 & 0& -1  &0 & 0& 0\\
   \end{array}
   \right);
   $$
   $\hat{p} = \frac{1}{3} \hat{p}_{-1} + \frac{1}{9} \hat{p}_{-2}$,
$$
\hat{p} = \frac{1}{9\sqrt{3}\mathrm{i}}\left(
\begin{array}{ccccccccc}
     0&3G_{8} &3G_{7} & 1& 3G_{5}& 3G_{4}&-1&3G_{2} & 3G_{1}\\
     3G_{1} & 0 &3G_{8} &3G_{7} & 1& 3G_{5}& 3G_{4}&-1 &3G_{2}\\
     3G_{2} &3G_{1} & 0 &3G_{8} &3G_{7} & 1& 3G_{5}& 3G_{4}-1\\
     -1&3G_{2} &3G_{1} & 0 &3G_{8} &3G_{7} & 1& 3G_{5}& 3G_{4}\\
     3G_{4} &-1&3G_{2} &3G_{1} & 0 &3G_{8} &3G_{7} & 1& 3G_{5}\\
     3G_{5} &3G_{4} &-1&3G_{2} &3G_{1} & 0 &3G_{8} &3G_{7} & 1\\
     1 &3G_{5} &3G_{4} &-1&3G_{2} &3G_{1} & 0 &3G_{8} &3G_{7} \\
     3G_{7} &1 &3G_{5} &3G_{4} &-1&3G_{2} &3G_{1} & 0 &3G_{8} \\
     3G_{8} & 3G_{7} &1 &3G_{5} &3G_{4} &-1&3G_{2} &3G_{1} & 0
\end{array}
\right).
$$
\section{Appendix B}
Let's consider the plot of the nodes of the lattice for the arbitrary digit $d_1$ (Fig. \ref{systems}). We can see that the $d_1$ is the parameter of the translation over the vector $\mathbf{V} = (1/2, -1)$. Thus we obtain something like a torus winding with step $\Delta x$, and $d_1$ defines the section of the torus.
\begin{figure}[h]
\centering
\begin{tikzpicture}
\draw [->] (0,0) -- (6,0) node [below] {$x$};
\draw [->] (0,0) -- (0,3) node [left] {$d_1$};
\draw [->, thick] (0,0) -- (2,0);
\draw [-, thick] (2,0) -- (4,0) node[below] {4};
\draw [->, thick] (0,2) -- (2,2);
\draw [-, thick] (2,2) -- (4,2);
\draw [->, thick] (0,0) -- (0,1);
\draw [->, thick] (0,1) -- (0,2) node[left] {2};
\draw [->, thick] (4,0) -- (4,1);
\draw [-, thick] (4,1) -- (4,2);
\draw[-, thick, red] (0,2) -- (2,0);
\draw[-, thick, red] (1,2) -- (3,0);
\draw[-, thick, red] (2,2) -- (4,0);
\draw[-, thick, red] (3,2) -- (4,1);
\draw[-, thick, red] (0,1) -- (1,0);
\draw[fill = red] (0,0) circle (1pt);
\draw[fill = red] (4,2) circle (1pt);
\draw[-, blue, dashed, thick] (-0.5,2) -- (4.5,2) node[right] {$d_1 = 0$};
\draw[-, blue, dashed, thick] (-0.5,1) -- (4.5,1) node[right] {$d_1 = -1$};
\draw[-, blue, dashed, thick] (-0.5,1.5) -- (4.5,1.5) node[right] {$d_1 = -1/2$};
\end{tikzpicture}
\caption{The torus winded with the red lines represents the view of the lattice for any possible value of $d_1$ for $\Delta x = 1$, $\Xi = 4$}
\end{figure}
\\
Here we can see, that if we change  $d_1$ to $d_1 + 1$, we obtain the same result.

We obtain a less trivial connection with the windings of torus if we plot over the $x$-axis the momentum and over the $y$-axis the momentum with respect of modulo $1$, in other words -- the fractional part of the momentum. Then we obtain the winding of torus with unite slope. Let $\Pi = 2^{n_-} = 4$, $d_1 = 0$, then we have the following plot:
\begin{figure}[h]
\centering
\begin{tikzpicture}
\draw [->] (0,0) -- (5,0) node [below] {$p$};
\draw [->] (0,0) -- (0,1.5) node [left] {$p \, \text{mod} ~ 1$};
\draw [->, thick] (0,0) -- (2,0);
\draw [-, thick] (2,0) -- (4,0) node[below] {4};
\draw [->, thick] (0,1) -- (2,1);
\draw [-, thick] (2,1) -- (4,1);
\draw [->, thick] (0,0) -- (0,0.5);
\draw [-, thick] (0,0.5) -- (0,1) node[left] {1};
\draw [->, thick] (4,0) -- (4,0.5);
\draw [-, thick] (4,0.5) -- (4,1);
\draw[-, thick, red] (0,0) -- (1,1);
\draw[-, thick, red] (1,0) -- (2,1);
\draw[-, thick, red] (2,0) -- (3,1);
\draw[-, thick, red] (3,0) -- (4,1);
\draw[-, thick,  dashed, blue] (0,0.5) -- (0.5,1);
\draw[-, thick,  dashed, blue] (0.5,0) -- (1.5,1);
\draw[-, thick, dashed, blue] (1.5,0) -- (2.5,1);
\draw[-, thick, dashed, blue] (2.5,0) -- (3.5,1);
\draw[-, thick, dashed, blue] (3.5,0) -- (4,0.5);
\end{tikzpicture}
\caption{The plot of fractional part of momentum for $\Pi = 4$, $d_1 = 0$ (red lines), $d_1 = -0.5$ (blue dashed line)}
\end{figure}
\\
Changing the digit $d_1$ we can move this winding without changing it's slope. In particular, the plot is shifted for $-d_1$ over the $p$-axis (or the coordinate axis is shifted for $d_1$).
\addcontentsline{toc}{section}{REFERENCES}