A Closed Form Selected Mapping Algorithm for PAPR Reduction in OFDM Multicarrier Transmission †

: Nowadays, the demand for communication multi-carriers’ channels, where the sub-channels are made mutually independent by using orthogonal frequency division multiplexing (OFDM), is widespread both for wireless and wired communication systems. Even if OFDM is a spectrally efﬁcient modulation scheme, due to the allowed number of subcarriers, high data rate, and good coverage, the transmitted signal can present high peak values in the time domain, due to inverse fast Fourier transform operations. This gives rise to high peak-to-average power ratio (PAPR) with respect to single carrier systems. These peaks can saturate the transmitting ampliﬁers, modifying the shape of the OFDM symbol and affecting its information content, and they give rise to electromagnetic compatibility issues for the surrounding electric devices. In this paper, a closed form PAPR reduction algorithm is proposed, which belongs to selected mapping (SLM) methods. These methods consist in shifting the phases of the components to minimize the amplitude of the peaks. The determination of the optimal set of phase shifts is a very complex problem; therefore, the SLM approaches proposed in literature all resort to iterative algorithms. Moreover, as this calculation must be performed online, both the computational cost and the effect on the bit rate (BR) cannot be established a priori. The proposed analytic algorithm ﬁnds the optimal phase shifts of an approximated formulation of the PAPR. Simulation results outperform unprocessed conventional OFDM transmission by several dBs. Moreover, the complementary cumulative distribution function (CCDF) shows that, in most of the packets, the proposed algorithm reduces the PAPR if compared with randomly selected phase shifts. For example, with a number of shifted phases U = 8, the CCFD curves corresponding to the analytical and random methods intersect at a probability value equal to 10 − 2 , which means that in 99% of cases the former method reduces the PAPR more than the latter one. This is also conﬁrmed by the value of the gain, which, at the same number of shifted phases and at the probability value equal to 10 − 1 , changes from 2.09 dB for the analytical to 1.68 dB for the random SLM.


Introduction
Orthogonal frequency division multiplexing (OFDM) is an efficient multi-carrier modulation technique both for wireless communication and power line communication (PLC).
For its high data rate, OFDM has been indeed widely used in many wireless communication standards, such as digital video broadcasting (DVB) and based mobile worldwide interoperability for microwave access [1]. It offers a considerable high spectral efficiency, multipath delay spread tolerance, immunity to the frequency selective fading channels and power efficiency [2,3].
Moreover, the rapid development of smart grids and electric transportation has made PLC technology a promising complementary scheme to the wireless communications [4][5][6][7][8].
Signal distortion techniques (SD), such as clipping and filtering, which means the non-linear treatment of amplitude near the peak power [21][22][23]; 3.
Multiple signaling and probabilistic techniques (MSP), such as selected mapping methods (SLM) [24][25][26], partial transmission sequence (PTS) [27,28], phase optimization [29], tone reservation and tone injection [30], and constellation shaping [31,32], which generate a permutation of the multi-carrier signal and choose the signal with the minimum PAPR for transmission, or modify different parameters in the OFDM signal, and optimize them to minimize the PAPR.
Over the last few years, several studies have been presented, which combine two or more PAPR reduction methods to use the characteristics of each method in a hybrid scheme to improve performance and to improve the cost efficiency in terms of better PAPR reduction [33]. Among these techniques, in literature there are some that combine the approach of selective mapping (SLM) with other approaches, often demonstrating advantages in terms of complexity and/or information integrity. In [34] two parallel forms of multiple signal representation (MSR) approaches, the SLM methodology and the modified-partial transmit sequence (PTS) technique called the cyclic shift sequences, PTS are mixed to enhance the PAPR efficiency and computational complexity level better than the conventional PTS methodology. In [35,36], SLM methods are used with compounding methods. The latter are applied to speech signal having high peaks and are also applied to the OFDM signal. In particular in [36], a new hybrid method has been proposed considering the cascade of two stages. In the first stage, the SLM is combined with PTS reduction technique and in the second stage the PTS technique is combined with compounding methods. This hybrid method is able to reduce the system complexity and to increase the bit error rate performance. In order to reduce the complexity of SLM and achieve a significant PAPR reduction, in [37] the authors propose a hybrid approach called low complexity hybrid selective mapping for PAPR mitigation for asymmetrically clipped optical OFDM. The new approach utilizes the properties of IFFT and Hermitian symmetry to reduce the complexity of the conventional SLM. Despite the great number of methods, until now, no one has been able to solve the problem at low computational cost and, at the same time, preserving the information content of the signal.
In [38], a detailed analysis of advantages/limitations of several approaches that combine SLM and PTS is reported. This analysis shows that there is a trade-off between the efficiency in the PAPR reduction and the computational complexity level of these hybrid methods.
Our proposal focuses on improving the performance of the SLM method, which is one of the most used in hybrid approaches. The SLM [25] approach consists in applying a phase shift to the components of the OFDM symbol, changing the shape of the OFDM symbol in such a way to limit the peaks. The determination of the optimal set of phase shifts is a very complex problem; thus, the SLM approaches proposed in literature all resort to iterative algorithms to solve it. In this paper, the problem of determining the optimal set of phase shift is solved in closed form on an approximated formulation of the PAPR. Such approximation uses the n-norm instead of infinity norm, with n = 4, in the definition of the PAPR. The gradient expression of the objective function has been analytically defined, and its solution consists simply in the calculation of the roots of a sixth degree polynomial, which represent the optimal solutions of the optimization problem. To better approximate the PAPR value, the analytical evaluation can be easily extended to a greater value of the n-norm. We have limited ourselves to n = 4 in order not to burden the reading of the analytical steps of the algorithm reported in the Appendix A of the paper. However, the proposed approach can be conceptually extended to a higher value of n at the expense of a higher number of analytical steps. Simulation results demonstrate the validity of the n-norm approximation and the advantage with respect to SLM iterative methods.
The paper is organized as follows. Section 2 introduces the OFDM system model. In Section 3, the PAPR problem is depicted together with the selected mapping method. In Section 4, the closed form algorithm is proposed, whereas the simulation results are reported in Section 5. The last Section 6 provides discussion and conclusions. Finally, in the Appendix A, the analytical steps are reported. Figure 1 shows the basic scheme of an OFDM system. An OFDM modulation scheme allows digital data to be efficiently and reliably transmitted over a transmitter channel and it performs well even in multipath environments with reduced receiver complexity. Using ODFM, it is possible to exploit the time domain, the space domain, the frequency domain and even the code domain to optimize transmitter channel usage. As shown in Figure 1, the bit data must be firstly converted from serial stream into parallel streams, each of which are coded and modulated on to a subcarrier. In fact, OFDM allows one to subdivide the available band into mutually orthogonal subcarriers so that each of them can be considered independently. These subcarriers are regularly spaced in frequency, forming a block of spectrum. The frequency spacing and time synchronization of the subcarrier is chosen in such a way that the subcarriers are orthogonal, which means they do not cause interference with one another. Moreover, using narrow subcarriers, the frequency response can be considered constant and then the signal does not undergo distortion effects. The signal sent through the single sub channel is the sum of two iso frequential in quadrature sinusoidal signals. The two degrees of freedom, given by the amplitude of the two sinusoids, allow one to send a two-dimensional message, corresponding to points in the plane ( Figure 2). The set of points constitutes a constellation, and each point is associated with a binary string; therefore, the more points, the more bits are associated to each point, and then the higher the bit rate. In Figure 2, two different kinds of constellation are shown: the amplitude phase shift keying (APSK), and the quadratic amplitude modulation (QAM). Each constellation is associated with a maximum power level, corresponding to the farthest points from the origin. To increase the number of points and then the length of the string associated to each of them, one can augment the power available for the sub channel or put the points closer each other. Note that, the first solution can give rise to saturation of the transmission devices and interferences, whereas the second solution can cause receiving errors, due to the unavoidable presence of noise. The QAM scheme is more diffused than the APSK because it is easier to implement, and the average signal power is lower. However, in this work, to cope with possible heavy noise conditions, we preferred to adopt the APSK constellation scheme.

PAPR in OFDM Systems
The OFDM typically suffers a high peak to average power ratio (PAPR), defined for transmitted signal x(t), as where T is the interval of OFDM signal (OFDM symbol period). PAPR has negative consequences, both on the interference caused, for example, by PLC transmission and on the allowed peak of power the system can support. In particular, if the instantaneous signal overcomes the linear range of the amplifier, this clips the peaks of the signal, which will result in distorted [39].

Selected Mapping Approach for PAPR Reduction
Selected mapping approach basically consists of shifting the phase of the sinusoidal components in order to reduce the value of the peaks. As previously highlighted, although SLM is one of the most used PAPR reduction techniques, it has some drawbacks. In fact, typically this approach is applied online, because a different set of phase shifts must be determined for any stream of bit and it is not reasonable calculating a priori the shifting for all the possible streams, but this implies that the set of phases before the shifting has to be sent as side information, by using to this purpose a part of the available channel. Note that, the reduction of the PAPR, whichever approach is used, is paid in terms of bit rate.
Most of the SLM approaches presented in literature consist in iteratively generating a large set of phase shift vectors, and then selecting the OFDM symbol having the lowest PAPR. The quality of the solution in terms of low PAPR depends on the number of phase shift vectors. However, the required calculation time must not exceed the rate the streams are generated, so as not to affect the transmission rate.
In [40], it is shown that the iterative procedure usually converges towards a limited set of phases to shift, because the incremental benefit that can be obtained by shifting the most parts of the components is very low. This fact suggests that a possible way to reduce the computational cost is to select a reduced number of phases to shift. Limiting the number of phases to shift has a benefit in terms of bit rate, because it reduces the side information of the transmission. Indeed, if the point of the constellation undergoes a phase shifting, this shift must be attached to the transmission, to make it possible for the receiver to reconstruct the original message. Therefore, the higher the number of shifted components, the higher is the part of the channel used to transmit side information.
In this work, the number of phases to be shifted is established a priori. The components whose phase is shifted are those having the highest values of power, which are potentially those mainly affecting the PAPR. In this way, the receiver can deduce the set of shifted components without any side information. To overcome the computational complexity of the algorithm, in this paper, the problem of finding the optimal phase shift vector is solved in closed form on an approximated formulation of the PAPR.
To compare the results of a PAPR reduction technique it is common to refer to the complementary cumulative distribution function (CCDF) of the PAPR, which provides a statistical description of the power levels in the OFDM signal. In particular, it gives the probability that PAPR of the OFDM signal x(t) will be above a given threshold PAPR 0 [10].
A PAPR reduction algorithm may satisfy the multiple objectives of reducing the PAPR without affecting the bit rate (BR) and the bit error rate (BER), hence limiting the computational complexity and the side information to be transmitted. Several performance indexes can be introduced. As SLM algorithms do not affect BER, in this paper the following performance index has been proposed [13]: 10 PAPR a f ter PAPR be f ore (2) which measures the relative PAPR reduction after the application of the algorithm at a fixed CCDF. A large value for the "gain" Γ implies better PAPR reduction.

Analytical Solution of SLM Problem
The calculation of the optimal phase shift vector is a hard task, which could require a high computational cost even in the case of a small number of phases to shift. The complex envelope of the transmitted OFDM signal can be written as: where X = {X k , k = −M, · · · , M} is the block of M symbols, M is the number of subcarriers, ω is the angular frequency of the fundamental harmonic, and j = √ −1. The equivalent real envelope can be obtained by summing M sinusoids corresponding to as many subcarriers: The aim of the optimization problem consists in minimizing the PAPR value given by (1) of the OFDM symbol (4) acting on a prefixed number U of phases ϕ k . The initial values of both amplitudes A k and phases ϕ k depend on the current frame of bits transmitted with the OFDM symbol, so that a new optimization problem must be solved online for each symbol.
The first objective is to remove the time variable in the formulation of the optimization problem. In fact, it is a difficult task to foresee the position of the peaks in the OFDM symbol, and this would be a further problem to address inside that of optimizing the PAPR. To this end, we resort to the norm of the signal, which is defined for the generic order n as: As it is well known, as n → ∞ , E n tends to the maximum value of x(t), but, in general, a small order n > 2 is sufficient to put in evidence the presence of an instantaneous peak. In particular, we can observe that the integral inside the root in (5), for n = 2 is equal to the energy of x(t), which is invariant with respect to the phases ϕ k of its components.
Therefore, even for n = 3, the norm could be used to estimate the presence of a peak in the signal. For the sake of simplicity, the norm n = 4 has been adopted, to avoid introducing the absolute value of the objective function. Furthermore, the root operation in (5) does not affect the position of the minimum, so that the optimization problem is formalized as follows: where Φ is the set of phases to shift and P is a function whose stationary points are assumed in place of those of the PAPR[x(t)]. Figure 3 reports PAPR[x(t)] (continuous red line) and P (continuous blue line) when changing only one phase angle from 0 • to 360 • . As expected, the two curves present an offset, due to their different definitions, but their minimums (highlighted with vertical dashed lines) are quite close and, assuming the stationary point evaluated on the blue curves, non-significant variation is introduced in the PAPR. Note that, a further improvement of the PAPR approximation can be achieved if a norm n, with n greater than 4, is considered, at the cost of a greater effort in the analytic determination of the optimal phase shifts. In the same Figure 3, the continuous green line reports the trend of P evaluated using the norm n = 16 instead of n = 4. As can be noted, the offset greatly reduces and the two minima points almost coincide. Considering (4) and (6), the general expression of the norm P can be calculated, where the fundamental frequency ω has been normalized to 1: In this paper, a line search algorithm [41] is used to seek the minimum of the objective function. To this end, the gradient expression of the objective function must be defined, in order to calculate the stationary points. It is worth noticing that this does not affect the generality of the results. By deriving (7) with respect to ϕ r ∈ Φ, the first-order condition can be determined: Let us re-write Equation (8) highlighting the r-th term in the sum: Let us define s k = sin(kt + ϕ k ) and c k = cos(kt + ϕ k ). Then, the first order condition becomes: The four terms in Equation (10) can be evaluated separately in closed form by simple analytical steps, reported in Appendix A, leading to: Second Term where where where Determination of ϕ r In the following, the phases ϕ r , corresponding to the stationary points, are deduced. Once the coefficients, A, B, C, D, E and F, whose derivation is reported in the Appendix A, are evaluated in closed form, the first-order condition (10) becomes: ∂P ∂ϕ r = A cos(3ϕ r ) + B sin(3ϕ r ) + C cos(2ϕ r ) + D sin(2ϕ r ) + E cos(ϕ r ) +F sin(ϕ r ) = 0 Applying the Weierstrass substitution (or tangent half-angle substitution) [42], which introduces the variable z = tan ϕ r 2 , the following parametric trigonometric relations can be stated: Then, using the duplication and triplication formulas, it follows that: Hence, the first-order condition becomes: Solving this equation with respect to z, the values of ϕ r are obtained as: It is worth to point out that the phase ϕ r , corresponding to the stationary point, is not unique, because it is obtained as the argument of a trigonometric function. This is in accordance with the fact that the stationary point could be either a minimum or a maximum point.

Simulation Results
Following the analytical procedure described in previous Sections, a first simulation has been conducted by considering a random bit stream composed of 10 7 bits and OFDM symbols with 64-APSK modulation to present the statistical analysis of the PAPR reduction by varying the number of shifted phases U. In the simulation, an OFDM base-band signal with N = 64 subcarriers and 390,625 frames has been considered.
In Figure 4, the comparison between the PAPR CCDF curves, obtained with the conventional OFDM and the SLM method with a number U of shifted phases ranging from 2 to 32, is reported. It can be observed that a small number of shift vectors allows the PAPR to be strongly reduced, but at the same time at the improvement tends to saturate very soon. Figure 5 reports the marginal improvement of the gain Γ of the random SLM method with the increase of the number U of shifts from 2 to 32. The gain has been calculated for CCDF equal to 10 −3 . As can be noted, the PAPR reduction tends to saturate as the number of shift vectors increases, and no further measurable advantage can be achieved with a number of shift vectors greater than eight.
The proposed analytic SLM has been tested on the same bit stream used for the above analysis. In Figure 6, the CCDF diagrams are shown for different numbers of shifted carriers for the analytical SLM and the random SLM. As can be observed, for U = 8, the analytical SLM outperforms the random one in the overwhelming majority of the cases, while the latter works better in the rare cases where the PAPR is very high.    Table 1, which presents the comparison between the analytic and the random SLM in terms of the gain Γ. Three cases are reported corresponding to three different values of CCDF. The larger values of Γ for analytic SLM for higher values of CCDF implies its better PAPR reduction in the large majority of cases. However, as already highlighted, with the random SLM algorithm, the maximum value of the PAPR reduction can be obtained with a limited number of phase shift vectors and, due to the low computational cost, it would be possible to apply in parallel the two algorithms and, for each frame, the one corresponding to the maximum reduction can be assumed to calculate the OFDM symbol. In this case, the effective CCDF curve is represented by the leftmost envelope of the random and analytical diagrams. Although the example referred to in Figure 6 is very basic, the results allow some important considerations to be made. The first is that the improvements achieved with the random approach tend to saturate for low values of U. In the real cases, in which the number of subcarriers is much higher than in the example, the saturation phenomenon will be even more evident, because the relative weight of the proposed random solutions will be much smaller. Another aspect that appears evident is that, with the same number of shifted phases, the proposed method allows one to obtain a greater reduction in PAPR in most cases. For example, with U = 8, the curves corresponding to the random and analytical methods intersect at a probability value equal to 10 −2 ; this means that in 99% of cases the analytical method reduces the PAPR more. In 1% of cases, the random method allows for a better result. This fact is attributable to the fact that a low-order norm was adopted. As shown in Figure 3, as the order of the norm increases, the minimum point of the norm tends to coincide with the global minimum, so that the growth of the norm necessarily involves a shift to the left of the curve in Figure 6. As explained previously, the adoption of a higher order norm does not involve a greater computational burden during transmission, but only a greater analytical development of the model.

This is confirmed by the values reported in
For the sake of comparison, in the same Table 1, the values of the gain provided in recent literature [43], by using the dispersive selection mapping (DSLM) scheme proposed in [44], have been added for the same number N of subcarriers and the same number U of phases shifted. Even if the DSLM presents limited advantages with respect to analytic and random SLM at a probability value equal to 10 −2 , the analytic SLM outperforms both random and DSLM at a probability value equal to 10 −1 . Note that DSLM has a higher computational complexity with respect to analytic and random SLM.

Discussion and Conclusions
The OFDM is a very appealing technique both for wireless and wired communications because of its spectrum efficiency and robust channel. One of its main disadvantages is that the composite transmitted signal may present a very high peak to average power ratio (PAPR) when the input sequences are strongly correlated. Hence, PAPR reduction is mandatory. A PAPR reduction algorithm may satisfy the multiple objectives of reducing the PAPR without affecting the bit rate (BR) and the bit error rate (BER), hence limiting the computational complexity and the side information to be transmitted.
In this paper, a closed form algorithm for the PAPR reduction in OFDM modulation scheme is proposed, which can be used both in wired and wireless communications. It represents an analytical solution of the optimization problem associated with the selective mapping (SLM) method, where the peaks occurring in the multi carrier signals are reduced by applying a phase shifting to the sinusoidal components of the signal. The analytical solution is obtained by an approximation of the PAPR definition, which uses the n-norm instead of infinity-norm, with n = 4, in the definition of the PAPR, and it presents several advantages with respect to the random SLM iterative approaches presented in literature.
Firstly, the analytically obtained solution represents the optimal solution of the minimization problem if the approximation introduced by replacing 4-norm with the infinitynorm in the definition of the PAPR is disregarded. Note that, the closed-form algorithm can be formulated using n-norm with increased value of n, achieving a better approximation of PAPR at the cost of increased effort in off-line analytical steps. In this paper, the order n = 4 has been adopted just so as not to burden the analytical steps reported in the Appendix A.
Moreover, the computational time of the random SLM method to define the phase shift is not fixed. In fact, a certain number of random attempts are performed, until the value of the PAPR is acceptable or no improvement is obtained after a prefixed maximum number of runs. Conversely, the computational complexity of the proposed closed form method is limited to the solution of a 6th degree polynomial equation, which guarantees that the best solution is found in a fixed time.
The potentiality of the proposed approach has been tested referring to a case study and comparison has been made with the random SLM PAPR reduction gain with respect to the conventional OFDM.
As shown in the results, with the same number of candidate phase vectors U with U > 4, the random SLM method exceeds the analytical SLM in the rare cases where the PAPR is very high, while the former works better in the vast majority of cases. Even with a number of shifted phases U = 8, the CCFD curves corresponding to the analytical and random methods intersect at a probability value equal to 10 −2 . This means that in 99% of cases the first method reduces the PAPR more than the second one. This is also confirmed by the value of the gain Γ, which measures the relative PAPR reduction after the application of the algorithm at a fixed CCDF. The larger the value of the gain, the better the PAPR reduction. With a number of shifted phases U = 8 and at the probability value equal to 10 −2 the gain assumes the same value Γ = 2.43 dB for both the analytic and random methods, whereas at the probability value equal to 10 −1 , Γ = 2.09 dB for the analytical and Γ = 1.93 dB for the random SLM. From the comparison with literature [43], analogue conclusions can be drawn.
The results could be improved if a better approximation of the PAPR is considered by increasing the norm order, at the cost of further offline effort in the analytic determination of the stationary points.
Note that, while with random SLM the computational complexity increases with the increase of U, and, at the same time, the relative improvement decreases, with the analytical approach, the computational complexity does not change with the increase of U while maintaining significant margins for improvement.
Recently, several hybrid approaches have been proposed to minimize the PAPR value, at the expense of increasing the computational complexity. Literature clearly demonstrates that the most promising approach is to combine two or more PAPR reduction methods to use the characteristics of each method in a hybrid scheme. However, there is a trade-off between the efficiency in the PAPR reduction and the computational complexity level of these hybrid methods. Improving the performance of each of the methods in the hybrid approach will improve the efficiency of the overall hybrid method. Many of them combine the approach of selective mapping (SLM) with other approaches, such as the partial transmit sequence (PTS) methods [34,36]. Our proposal focuses on improving the performance of the SLM method, which is in fact one of the most used in hybrid approaches. The proposed closed form SLM algorithm, which per se is not superior to other methods, is not intended as an alternative to these recent hybrid methods, but, if ever, as an improvement of the SLM approach present in the hybrid methods with consequent impact on both the computational complexity of hybrid methods and the computational times.
A comparison with the most performing techniques will be done in future works by implementing the proposed analytical SLM approach within hybrid techniques. Moreover, in future works, the analytic approach will be extended to a higher value of norm order. Funding: This work has been supported by Fondazione di Sardegna under project "SISCO-ITC methodologies for the security of complex systems", CUP: F74I19001060007.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A
In this appendix, the four terms in Equation (10) The first term of (10) is null ∀r because it is equal to the integral of two sinusoids within a multiple of their period.
Second Term The first term of (A3) leads to the sum of integrals whose argument is the product of two sinusoids with different frequency (because k = r); therefore, it is null ∀r. Concerning the second term, applying the Werner formulas to c 2r c r , the following formula is obtained: The first addend in (A4) leads to the sum of integrals whose argument is the product of two sinusoids with different frequency (because k = r); therefore, it is null ∀r. For what concerns the second addend: Then, the integral (A5) is not null for k = 3r. In fact, writing the expression (A5) in the extended form, and using the Werner formulas, and considering k = 3r: where: Let us now analyze the two integrals in (A8) separately. For what concerns the first integral: The first integral in (A9) is null ∀r because it is the integral of a sinusoid within a multiple of its period. The second integral leads to the sum of integrals whose argument is the product of two sinusoids with different frequency (because k = r); therefore, it is null ∀r.
Let us now analyze each term one by one.

1.
Calculation of the term: The first addend in (A14) is null ∀r because it leads to the sum of integrals whose argument is the product of two sinusoids with different frequency (because k = r).
Concerning the second addend: The second integral in (A15) is null ∀r because it leads to the sum of integrals whose argument is the product of two sinusoids with different frequency (because k = r). Then, only the first addend survives leading to: Then, the integral in (A16) is not null only for k = r 3 : Calculation of the term: The integral (A19) is not null for two cases, i.e., k ± 2l = r: where: 3. Calculation of the term: The integral (A22) is not null in four cases, i.e., k + l ∓ m = r ; k − l ± m = r. In particular, for: k + l ∓ m = r [sin(2rt + ϕ k + ϕ l ∓ ϕ m + ϕ r ) The same exact procedure holds for the second case: i.e., for k − l ± m = r [sin(2rt + ϕ k − ϕ l ± ϕ m + ϕ r ) Summarizing the previous results, for k ± ( l ∓ m) = r 2π 0 6 M ∑ k, l, m = r k > l > m A k A l A m s k s l s m dt = E cos(ϕ r ) + F sin(ϕ r ) where