Abstract
In 2020 E. Sakalauskas with coauthors published a paper defining perfectly secure Shannon cipher based on matrix power function, proposing effective parallelization, and ensuring no need for multiple rounds encrypting one data block [1]. In this paper we present computational results with the avalanche effect and bit independence criterion (BIC). These criteria are important when describing the rate of confusion of bits in the ciphertext. It was observed that increasing matrix order and group size enhance BIC and avalanche effect results converging to the desired values. Based on the outputs it is possible to pick appropriate parameters satisfying security needs and available memory in a device where appropriate keys are going to be stored.
1. Introduction
1.1. Avalanche effect and BIC
Cryptography security analysis methods such as avalanche effect and BIC allow us to evaluate block cipher secrecy by computing elements confusion after changing just one bit [2], determine elements confusion and dependance from other elements [3, 4]. The values of these criteria are commonly calculated by considering the avalanche vector ${A}^{{e}_{i}}$, which describes ciphertext bits change after flipping one bit in the plaintext:
where vector has all entries equal to 0 except for the $i$th one which is equal to 1, entry ${a}_{1}^{{e}_{i}}\in \left\{\mathrm{0,1}\right\}$ and function $Enc\left(k,\mu \right)$ is encryption function mapping shared key $k$and plaintext $\mu $ to the ciphertext generally denoted by $c$.
Using expression defined in Eq. (1), we compute the $i$th bit avalanche ${k}_{AVAL}\left(i\right)$ effect as follows:
$W\left({a}_{j}^{{e}_{i}}\right)=\sum _{X\in {\left\{\mathrm{0,1}\right\}}^{n}}{a}_{j}^{{e}_{i}},$
where ${k}_{AVAL}\left(i\right)$ indicates the number of bits changes after flipping $i$th bit. The desired value of the avalanche effect is 0.5 for all the bits, meaning that it is infeasible to distinguish which bit changes occur after flipping a random bit of the original message.
The bit independence of the two entries is being calculated by the maximal absolute correlation coefficient between avalanche vector $j$ and $k$ components. According to [2], BIC can be calculated by the formula:
Furthermore, relying on Eq. (3) we can define the overall BIC for the whole ciphertext block $Enc\left(k,m\right)$ as the maximal correlation by checking all available pairs:
Ideally, the value of BIC should be close to 0 hence ensuring that all the bit changes occur statistically independently.
1.2. Perfectly secure Shannon cipher based on matrix power function
The matrix power function (MPF) was introduced in [5], as the following mapping acting on the Cartesian product of the space of square matrices of order $m$ with itself:
The general notation for this mapping is as follows:
where $W,E\in Ma{t}_{m}\left(\mathbb{S}\right)$ are matrices with entries from semigroup $\mathbb{S}$ and $X,Y\in Ma{t}_{m}\left(\mathbb{R}\right)$ are matrices with entries from a finite ring of integers $\mathbb{R}$. This mapping allows us to raise the base matrix $W$ to the socalled power matrices $X$ and $Y$.
E. Sakalauskas with coauthors used the MPF in 2020 to propose a perfectly secure Shannon cipher defined over ${\mathbb{Z}}_{3}$ [1]. This cipher uses a plaintext matrix $M$, private keys $X$ and $Y$ along with a function $f:{\mathbb{Z}}_{3}\mapsto {\mathbb{G}}_{3}$, which maps elements of ${\mathbb{Z}}_{3}$ to elements of the multiplicative Sylow group ${\mathbb{G}}_{3}=\{\mathrm{1,2},4\}$, which is a subgroup of ${\mathbb{Z}}_{7}^{*}$. Note, that actions in ${\mathbb{G}}_{3}$ are performed modulo 7. A key feature of this mapping is that it does not carry over the addition in ${\mathbb{Z}}_{3}$ to the multiplication in ${\mathbb{G}}_{3}$ and hence is not an isomorphism. The encryption function can be expressed in a following way:
where $F:Ma{t}_{m}\left({\mathbb{Z}}_{3}\right)\mapsto Ma{t}_{m}\left({\mathbb{G}}_{3}\right)$ is an entrywise matrix analogue of the mapping f and ${F}^{1}$ is its inverse. Note that since F is not an isomorphism no cancelations in Eq. (6) are possible. We also use $\u2a00$ to denote Hadamard product of two matrices.
It is worthy noting that the shared key $\left\{X,Y\right\}$ consists of $2{m}^{2}$ entries and hence is at least twice the length of the original plaintext given that extra bits may be added at the end message to make it appropriate length. However, the plaintext and ciphertext are roughly the same size.
To decipher the ciphertext, we denote by ${T}^{H}$ the inverse of matrix $T$ in Hadamard sense i.e., a matrix satisfying the following relation:
where every entry in the matrix is the unit of the group ${\mathbb{G}}_{3}$.
Upon receiving the ciphertext its decryption is performed in the reverse order and can be summarized by the following expression:
Perfect secrecy of the presented block cipher and the statistical independency of the ciphertext $C$ from the plaintext $M$ is proven in [1].
In this paper we investigate the avalanche effect and BIC for the presented block cipher in a more general form i.e., we expand the cardinalities of the algebraic structures considered. In other words, we consider the Sylow group ${\mathbb{G}}_{q}$ of the multiplicative group ${\mathbb{Z}}_{p}^{*}$ and an additive group ${\mathbb{Z}}_{q}$. Hence in Eq. (4) we have $W,E\in Ma{t}_{m}\left({\mathbb{G}}_{q}\right)$ and $X,Y\in Mat\left({\mathbb{Z}}_{q}\right)$. Actions in ${\mathbb{G}}_{q}$ are performed modulo a prime $p=2q+1$.
2. Computational results
The avalanche effect of perfectly secure Shannon cipher defined in Eq. (6) is calculated using Eq. (2). For each fixed pair of parameters $\left\{p,q\right\}$ we investigate the relation between avalanche effect and the matrix order $m$. We executed 1000 experiments and the results averaged for each value of $m$ given the fixed pair $\left\{p,q\right\}$. In Table 1 we present the results of our experiments.
Table 1Avalanche effect with different parameters
$p=\text{7}$, $q=\text{3}$  $p=\text{23}$, $q=\text{11}$  $p=\text{107}$, $q=\text{53}$  $p=\text{4079}$, $q=\text{2039}$  $p=\text{33553799}$, $q=\text{16776899}$  
$m=$ 5  0.4446  0.4615  0.4922  0.5007  0.5001 
$m=$ 8  0.4447  0.4638  0.4919  0.5002  0.5004 
$m=$ 10  0.4451  0.4626  0.4910  0.4997  0.4999 
$m=$ 15  0.4450  0.4627  0.4911  0.4996  0.4998 
$m=$ 16  0.4448  0.4625  0.4921  0.5002  0.4995 
$m=$ 32  0.4442  0.4628  0.4921  0.4999  0.5000 
Analyzing the obtained results, we see that as the parameter $q$ gets larger the avalanche effect increases to 0.5 whereas the matrix order does not have such big of an impact.
We perform the investigation of the BIC in a way similar to the one presented above. As above we performed experiments for each triplet $\left\{p,q,m\right\}$ and using Eq. (4) obtained the BIC values presented in Table 2.
Table 2BIC with the different parameters
$p=\text{7}$, $q=\text{3}$  $p=\text{23}$, $q=\text{11}$  $p=\text{107}$, $q=\text{53}$  $p=\text{4079}$, $q=\text{2039}$  $p=\text{33553799}$, $q=\text{16776899}$  
$m=$ 5  1  0.7391  0.4717  0.2508  0.1731 
$m=$ 8  1  0.5746  0.4074  0.2035  0.1264 
$m=$ 10  1  0.5798  0.3738  0.1555  0.1099 
$m=$ 15  1  0.5678  0.3539  0.1214  0.0785 
$m=$ 16  1  0.5704  0.3436  0.1171  0.0667 
$m=$ 32  1  0.5530  0.3367  0.0782  0.0376 
Note that increasing group size reduces BIC value. However, more importantly we see that small values of the parameter $p$ are clearly not suitable for implementation since the value of BIC approaches the worst possible case. Furthermore, we can see that increasing matrix order has some impact as well and it is more noticeable compared to an analogous result of the analysis of the avalanche effect.
3. Conclusions
In this paper we investigated the previously proposed Shannon block cipher which does not require multiple rounds to encrypt a message. Furthermore, we expanded our research of the initial scheme by introducing a pair of parameters $\{p,q\}$ which makes our cipher more flexible as compared to the original. The obtained results show that even though no information about the plaintext is revealed by the encryption algorithm itself, small values of parameters $q$ cannot be used in practice since the BIC fails even for the largest value of matrix order we considered. However, the avalanche criterion is mostly satisfied and is quite near perfection even for small values of q. Hence, relying on the results presented in Table 1 and Table 2, a good recommendation to choose the system parameters $\{p,q,m\}$ is to find a balance between $q$ and $m$ keeping them reasonably small while also ensuring that BIC is satisfied. Keeping this in mind a triplet {4079, 2039, 15} can be considered a suitable choice for practical implementation.
References

E. Sakalauskas, L. Dindienė, A. Kilčiauskas, and K. Lukšys, “Perfectly secure Shannon cipher construction based on the matrix power function,” Symmetry, Vol. 12, No. 5, p. 860, May 2020, https://doi.org/10.3390/sym12050860

Işıl Vergili, “Avalanche and bit independence properties for the ensembles of randomly chosen n × n SBoxes,” Turkish Journal of Electrical Engineering Computer Sciences, Vol. 9, No. 2, pp. 137–145, 2001.

M. Salman, R. Yugitama, A., and R. F. Sari, “KAMIES: Security optimization of KASUMI algorithm by increasing diffusion level,” International Journal of Security and Its Applications, Vol. 12, No. 3, pp. 29–46, May 2018, https://doi.org/10.14257/ijsia.2018.12.3.04

L. Liu, “Designing a random Sbox with the mixed spatiotemporal chaos,” in Conference series, Vol. 1983, No. 1, 2021, https://doi.org/10.1088/17426596/1983/1/012040

Sakalauskas Eligijus and Lukšys Kęstutis, “The matrix power function and its application to block cipher Sbox construction,” International Journal of Innovative Computing, Information and Control, Vol. 8, No. 4, pp. 2655–2663, 2012.