Abstract
Recently, the single photon avalanche diode optical wireless communication (SPAD OWC) has attracted much attention due to its potential underwater applications. For such system, the channel noise is additive Poisson noise (APN) rather than the commonly encountered additive white Gaussian noise (AWGN) and the corresponding maximum likelihood (ML) detection is hard to provide a useful insight into energy-efficient signal design. By using the previously proposed Hellinger distance design criterion, we design an energy-efficient multi-dimensional constellation within the nonnegative integer set by minimizing the average optical power for a fixed minimum Hellinger distance. Comprehensive simulations indicate that our designed constellation can substantially outperform the currently available pulse amplitude modulation (PAM) and squared PAM for SPAD OWC systems.
© 2017 Optical Society of America
1. Introduction
Optical wireless communication (OWC) has been viewed as an attractive technique for the future wireless communication, due to its unique advantages of high security, freedom from spectral licensing issue and large potential bandwidth [1–7]. For such system, the ubiquitously installed light emitting diodes (LEDs) and photodiodes (PDs) are natural transmitters and receivers using simple intensity modulation and direct detection (IM/DD) [8–10]. Recently, single photon avalanche diode (SPAD), as a new type of detector for OWC [11,12], has drawn much attention due to its higher sensitivity than the commonly used PD [13–15] and its potential applications to underwater communications [16, 17]. For underwater environment, the light signal suffers three main impairments including absorption, scattering and turbulence. The influnece of these impairments can be described by various models such as double-gamma function [18] and log-normal function [19]. After long propagation, the received signal is usually weak, thus high energy-efficiency of transmitted signal is essential to ensure a reliable communication link. In this paper, we consider the energy-efficient constellation design for the SPAD OWC systems.
However, compared with the design criterion and techniques of radio frequency (RF) communication [20–27] and PD OWC, SPAD OWC system has two significant differences. The first difference is that the signal of SPAD OWC or PD OWC is required to be nonnegative which is different from that of the RF communications [28,29]. Therefore, the well-developed radio frequency techniques [28–33] can not be applied to SPAD OWC systems in a straightforward manner. Since the output of SPAD is influenced by photon-counting-induced additive Poisson noise (APN) [13–15,34–38]. For SPAD OWC, the second significant difference is that the signal design criterion and techniques over APN cannot directly follow those over AWGN channels. For example, for PD OWC, pulse amplitude modulation (PAM) is the most energy-efficient scalar constellation based on the Euclidian distance criterion. Unfortunately, according to the results in [39–41], this is no longer true for our considered SPAD OWC, which is under Poisson regime. For SPAD OWC, the received optical signal is usually weak and thus, to ensure energy-efficient transmission, multidimensional signal design is required. However, the error performance analysis of the maximum likelihood (ML) detection in Poisson regime is a challenging task. More recently, the authors in [39] and [42] proposed Anscombe root (AR) receiver and further established a Hellinger distance criterion for signal design. It is also noticed that the designed Hellinger-distance-optimal squared PAM (SPAM) [39,42] is optimal only within one-dimensional integer space. Up to now, the energy-efficient multidimensional constellation design remains unsolved due to the challenge of solving the corresponding max-min optimization problems.
The above factors indeed motivate us to attack this open problem. In this paper, we devise an energy-efficient multidimensional integer constellation by maximizing the Hellinger distance under a normalized average optical power budget. Comprehensive simulations indicate that our designed constellation significantly outperforms the traditional PAM and SPAM schemes for SPAD OWC systems.
2. System model with AR receiver
In this section, we introduce the channel model of SPAD OWC and the previously proposed AR receiver for such system.
2.1. Transmitter
Let us consider an SPAD OWC system over an ideal Poisson channel. Let denote the transmitted signal over L time slots. is the set of all L × 1 vectors whose entries are nonnegative real-valued numbers to satisfy the nonnegativity requirement of intensity modulation. Then, the output vector of the SPADs, r, is given by [43]
where α and β respectively represent the power gain caused by SPAD and the counted photon number caused by the dark current over a symbol interval. Since p is an L-dimensional Poisson noise vector, the entries of r are identically independently Poisson-distributed. Then, the probability density functions (PDFs) of r and p are respectively:where n and n̂ + λℓ are nonnegative integers, and λℓ is the mean of Poisson distribution, which can be expressed as λℓ = αxℓ + β for ℓ = 1, 2, · · · , L [39,43].
2.2. Anscombe root (AR) receiver
A vector-version AR receiver was proposed in [42]. For presentation convenience, now we would like to introduce this receiver in the following.
We first use the existing Anscombe root (AR) operation to transform the entries of r as
Then, when the value of rℓ is sufficiently large [42], the output of AR transform, , is approximately distributed with a Gaussian distribution. At the same time, the proof in [44] tells the fact that where E(.) is the expectation of a random statistic. Therefore, we can finally approximate by where and ξℓ is additive Gaussian noise with zero mean and variance being one. In [42], the scalar-version expression (3) was extended to a vector version. For notation simplicity, we define and x̄ = (x̄1, · · · , x̄L)T. The vector-version expression for the resulting Gaussian channel is expressed as Now, the vector-version AR receiver proposed in [42] for our considered system is given below.AR Receiver: Given the L-dimensional received signal r defined by (1), the estimation of the transmitted x, x̂, at the AR receiver is determined by the following two successive steps:
- AR transform: , where is given by (2).
- Minimum Distance Decision: , where 𝒳 is the constellation of transmitted signal.
In this paper, our main task is to design an energy-efficient L-dimensional constellation 𝒳 for the above-mentioned AR receiver.
3. Energy-efficient multidimensional Hellinger constellation
In this section, we formulate the constellation design problem and then present the main result.
3.1. Problem formulation
In this subsection, our main task is to formulate our design problem for the AR receiver. Here, a natural question comes up: what decides the error performance of the constellation design at the AR receiver? According to the introduction of AR receiver in Subsection 2.2, we know that, when the optical power of the transmitted signal x is sufficiently high, the error performance for our considered channel is dominated by the minimum Hellinger distance between two distinct signal vectors [42]. For convenience, we define HD(𝒳) as the minimum Hellinger distance of the transmitted constellation 𝒳, say
This observation provides us with an explicit Hellinger-distance-based criterion for energy-efficient signal design for AR receiver. In general, for AR receiver, the Hellinger-distance-optimal multidimensional constellation design is as hard as solving a parallel problem of designing the optimal multidimensional constellation over AWGN channels, which is a long-standing and open problem in RF digital communication [20,24,25,45]. Since the optimal design problem for AR receiver is not easy to be transformed into a tractable one, we constrain our discussion into an integer set. Therefore, our design problem is formulated below:
Problem 1 Let . denotes the set of L × 1 vectors with all L entries being nonnegative integers. Then, for arbitrarily given positive integers L and K, devise an L-dimensional size-2K constellation such that the average transmitted optical power is minimized subject to HD(𝒳) = 1, where HD(𝒳) is defined in (5).
3.2. Main result
In this subsection, we present the solution to Problem 1. To make our presentation as clear as possible, let us first introduce some existing mathematical theories related to solving Problem 1. The solution to Problem 1 is related to a Diophantine equation with respect to z. Solving this equation is equivalent to answering the following question: “Under what condition, an integer can be decomposed into the sum of integer squares ?” For such question, there are several classical theorems for specific L. For example, when L = 2, Fermat’s two-square theorem asserts that an odd prime number z can be expressed as z = x2 + y2 for integer x and y, if and only if z has the form 4n + 1, n ∈ ℤ+ [46,47]. Moreover, Legendre’s three-square theorem and Lagrange’s four-square theorem were proved by [48] and [49] respectively. Unfortunately, up to now, the solution to for a general L still remains unsolved [50]. We are interested in the fact that the total number of the solutions to is finite.
For presentation convenience, we denote and define Rn by the cardinality of . Furthermore, we let N̄ be defined by the smallest positive integer satisfying . By our definition of 𝒩, we denote number of all the elements of 𝒩 not larger than N̄ by I. Then, without loss of generality, we arrange all these I numbers in an increasing order as n1 < · · · < nI with n1 = 0 and nI = N̄.
Now, we formally state our main result in this paper.
Theorem 1 An optimal solution to Problem 1 is given by
wherefor 1 ≤ i ≤ I − 1 and where .On Theorem 1, whose proof is provided in the appendix, we would like to make the following remarks:
- In this theorem, we show that our proposed design is closely related to a well-known mathematical problem. By using the fact that the total number of the solutions to is finite, we attain the energy-efficient structure within a nonnegative integer set for AR receiver. Our proposed constellation is given for a general dimension number L. It is observed that, when L = 1, the constellation defined by (6) is reduced to , which is exactly the so-called SPAM, developed in [39].
- To show the significant difference of our proposed constellation from SPAM and commonly used PAM, we illustrate these constellations in Figs. 1 and 2 for L = 2 and L = 3, respectively. As illustrated by these two figures, PAM generated by adding proper direct current to the unipolar version is equally spaced in each dimension, while each dimension of SPAM is the square of modified PAM. By Figs. 1 and 2, our proposed constellation can be viewed as a subset curved from a larger-size SPAM since the energy-efficiency of SPAM is remarkably influenced by the signal points with the largest energy. Therefore, our proposed scheme is the most energy-efficient compared with PAM and SPAM for a fixed Hellinger distance.
4. Simulation results
In this section, we carry out simulations to examine the performance of AR receiver and constellations. Because a proper bit-mapping is still unavailable for the proposed scheme and the relevant study is out of the scope of this paper, average codeword error rate is used as the performance metrics in our simulations, which is defined by the ratio of the number of estimated codewords in error to the total number of transmitted codewords. A codeword error occurs when x ≠ x̂, where x and x̂ are the transmitted codeword and the estimation of x respectively. In addition, we would like to compare the following three schemes.
- Pulse amplitude modulation (PAM): The transmitted signal based on PAM is given by (x1, x2, · · · , xL)T where for 1 ≤ ℓ ≤ L and ;
- Square pulse amplitude modulation (SPAM) [39]: For SPAM, the transmitted signal vector is determined as (x1, x2, · · · , xL)T where for 1 ≤ ℓ ≤ L and .
- Proposed constellation: The proposed constellation 𝒳̃ is defined in Theorem 1.
For comparison fairness, we assure that the above three kinds of constellation have the same bit rate and average transmitted optical power. In addition, the simulation parameters are adopted as α = 4.52 × 1014s/J, β = 7.27Hz [43] and more details are given below.
4.1. Performance of AR and ML receivers
In this subsection, we show the asymptotic statistic behavior of the AR receiver. Robust error performance of AR receiver only requires a reliable estimation of the transmitted symbol x. For AR receiver, the PDF difference between and xぃℓ + ξℓ decides the receiver performance. Thus, we first investigate the PDF of for various optical irradiance and then examine the performance of AR receiver compared with its ML counterpart.
- 1) PDF difference between and x̄ℓ + ξℓ: Without loss of generality, let us consider r̂ℓ which is expressed as . Simulations are carried out to investigate whether the PDF of r̂ℓ approaches that of rℓ which is a Gaussian distributed variable with zero mean and variance being 1. Results are shown in Fig. 3 with various power levels and indicate that the PDF of the AR transformed signals asymptotically approaches that of a Gaussian random variable with increasing optical irradiance.
- 2) Error performance comparisons of ML and AR receivers: To examine the performance of AR receiver by comparing with that of ML receiver for PAM, SPAM and our proposed constellation, we simulate for L = 2 and different modulation order and show the results in Fig. 4, which illustrates that the AR receiver has performance approaching the ML receiver for various modulation order when PAM, SPAM and our proposed multi-dimensional constellation are adopted. Specifically, in Fig. 4(a) the performance loss of AR receiver for PAM constellation are about 0.07dB, 0.05dB and 0.04dB for K = 3, 4, 5 respectively at the error rate of 10−4, while in Fig. 4(b) at the same target error rate, the respective power losses caused by AR receiver are about 0.23dB, 0.20dB and 0.10dB respectively when SPAM is adopted. Fig. 4(c) also indicates the performance gaps between AR receiver and its ML counterpart are respectively 0.20dB, 0.13dB and 0.15dB for K = 3, 4, 5 with the proposed constellation at the error rate of 10−4. These results indicate that our proposed AR receiver can be viewed as an effective alternative to the ML receiver and thus, our proposed Hellinger distance design criterion can be provide a useful insight into the signal design for SPAD optical wireless communications.
4.2. Performance of different constellations
Based on AR receiver, we compare the average codeword error rate of proposed constellation with those of PAM and SPAM for different modulation order with L = 2 and L = 3 respectively in Figs. 5 and 6. Notice that our proposed constellation and SPAM outperform PAM mainly for the reason that PAM is designed based on Euclidian distance and therefore is not energy-efficient any longer for our considered Poisson channels. In addition, the performance comparisons between SPAM and the proposed constellation are carried for L = 4 and L = 5 respectively in Figs. 7 and 8. Combining with the results in Figs. 5 and 6, we can observe that, for L = 2, K = 4, L = 3, K = 6 and L = 4, K = 8, SPAM and our proposed constellations have almost similar performance. It should be noticed that the attained gain by our proposed scheme is dependent on L and K. This phenomenon occurs when K can be exactly divided by L. The reason for this phenomenon is that, under such condition, the signal points of SPAM are evenly distributed on each axis which has been especially illustrated by Fig. 1(b), whose geometrical structure is similar to that of our proposed scheme. We also notice that our proposed constellation can obtain huge power gain over SPAM. Specifically, when K = 3, 5 in Fig. 5, the proposed constellation has respective power advantages of about 1.7 and 1.3 dB over SPAM at the error rate of 10−4. When L = 3, from Fig. 6, the attained power gains of the proposed constellation over SPAM are about 1.75 and 1.2 dB for K = 4, 5, respectively, at the target error rate of 10−4.
For L = 4, as shown in Fig. 7, the power gains by the proposed constellation compared with SPAM at the target error rate of 10−4 are 2.2, 2, 1.4 and 0.4 dB for K = 5, 6, 7, 8, respectively. While, for L = 5, the respective power advantages over SPAM for K = 6, 7, 8, 9 at the target error rate of 10−4 are 1.6, 1.5, 1.9 and 1.1 dB which is shown in Fig. 8. Now, we conclude that our designed constellation outperforms SPAM with the power gain due to its energy-efficient constellation structure.
5. Conclusion
In this paper, we have investigated the energy-efficient multidimensional constellation designs for the SPAD-based OWC systems. For AR receiver, we have designed an energy-efficient multi-dimensional constellation by minimizing the average optical power for a fixed minimum Hellinger distance of two distinct signals within the nonnegative integer set. Comprehensive simulations have shown that our proposed multi-dimensional constellation has significant performance gains over the currently available PAM and SPAM for this system.
This is the first systematic design of multidimensional signals for SPAD OWC systems and will provide a useful insight into the energy-efficient modulation for OWC under Poisson regimes. However, it should be noted that the AR detection complexity of our proposed constellation is exponential with respect to the set size. It is quite interesting to study how to reduce the complexity. In the future, we will attack this problem by utilizing the special structure of the proposed scheme. For example, the regular structure brought by Diophantine equation motivates us to demodulate the signal layer by layer. Moreover, notice that the proposed constellation is curved from a larger-size SPAM while an SPAM is the square of a modified PAM, after AR transform, the constellation is almost equally spaced in transform domain due to the square root operation, which inspires us to design a low-complexity receiver with equally spaced thresholds.
Appendix
Let us consider any size-2K constellation 𝒳 satisfying HD(𝒳) = 1 and for any x ∈ 𝒳. It should be noted that 𝒳 can be always decomposed into I disjoint subsets such that where for 1 ≤ i ≤ I − 1 and where with the following assumptions,
- 1T xi,j ≤ 1T xi,j+1 for j = 1, · · · , Rni − 1,
- 1T xi,Rni ≤ 1T xi+1,1 for i = 1, · · · , I − 1,
- 1T xI,j ≤ 1T xI,j+1 for .
Funding
Henan Major Scientific and Technological Project (161100210200); National Youth Foundation of China (61701536).
References and links
1. J. Li and M. Uysal, “Optical wireless communications: system model, capacity and coding,” in Proceedings of IEEE Vehicular Technology Conference (IEEE, 2003), pp. 168–172.
2. D. K. Borah, A. C. Boucouvalas, C. C. Davis, S. Hranilovic, and K. Yiannopoulos, “A review of communication-oriented optical wireless systems,” EURASIP J. Wirel. Commun. Netw. 91(1), 1–28 (2012).
3. A. Jovicic, J. Li, and T. Richardson, “Visible light communication: Opportunities, challenges and the path to market,” IEEE Commun. Mag. 51(12), 26–32 (2013). [CrossRef]
4. M. A. Khalighi and M. Uysal, “Survey on free space optical communication: A communication theory perspective,” Commun. Surveys Tuts. 16(4), 2231–2258 (2014). [CrossRef]
5. S. Chaudhary and A. Amphawan, “The role and challenges of free-space optical systems,” J. Opt. Commun. 35(4), 558–565 (2014). [CrossRef]
6. A. K. Majumdar, Free-space Optical (FSO) Platforms: Unmanned Aerial Vehicle (UAV) and Mobile (Springer, 2014).
7. S. K. Routray, “The changing trends of optical communication,” IEEE Potentials 33(1), 28–33 (2014). [CrossRef]
8. S. Hranilovic and F. R. Kschischang, “Optical intensity-modulated direct detection channels: signal space and lattice codes,” IEEE Trans. Inf. Theory 49(6), 1385–1399 (2003). [CrossRef]
9. T. Komine and M. Nakagawa, “Fundamental analysis for visible-light communication system using led lights,” IEEE Trans. Consum. Electron. 50(1), 100–107 (2004). [CrossRef]
10. L. Zeng, D. O’Brien, H. Minh, G. Faulkner, K. Lee, D. Jung, Y. Oh, and E. T. Won, “High data rate multiple input multiple output (MIMO) optical wireless communications using white led lighting,” IEEE J. Sel. Areas Commun. 27(9), 1654–1662 (2009). [CrossRef]
11. Y. Li, S. Videv, M. Abdallah, K. Qaraqe, M. Uysal, and H. Haas, “Single photon avalanche diode (SPAD) VLC system and application to downhole monitoring,” in Proceedings of IEEE Global Communication Conference (GLOBECOM) (IEEE, 2014), pp. 2108–2113.
12. T. Mao, Z. Wang, and Q. Wang, “Receiver design for SPAD-based VLC systems under Possion-Gaussian mixed noise model,” Opt. Express 25(2), 799–809 (2017). [CrossRef] [PubMed]
13. H. Elgala, R. Mesleh, and H. Haas, “Indoor optical wireless communication: potential and state-of-the-art,” IEEE Commun. Mag. 49(9), 56–62 (2011). [CrossRef]
14. S. Arnon, J. Barry, G. Karagiannidis, R. Schober, and M. Uysal, Advanced Optical Wireless Communication Systems (Cambridge University, 2012). [CrossRef]
15. Y. J. Zhu, W. F. Liang, J. K. Zhang, and Y. Y. Zhang, “Space-collaborative constellation designs for MIMO indoor visible light communications,” IEEE Photon. Technol. Lett. 27(15), 1667–1670 (2015). [CrossRef]
16. P. A. Hiskett and R. A. Lamb, “Underwater optical communications with a single photon-counting system,” Proc. SPIE 9114, 91140P (2014). [CrossRef]
17. C. Wang, H. Y. Yu, and Y. J. Zhu, “A long distance underwater visible light communication system with single photon avalanche diode,” IEEE Photon. J. 8(5), 1–11 (2016). [CrossRef]
18. S. Tang, Y. Dong, and X. Zhang, “Impulse response modeling for underwater wireless optical communication links,” IEEE Trans. Commun. 62(1), 226–234 (2014). [CrossRef]
19. M. V. Jamali and J. A. Salehi, “On the BER of multiple-input multipleoutput underwater wireless optical communication systems,” in Proceedings of IEEE 4th International Workshop on Optical Wireless Communications (IWOW) (IEEE, 2015), pp. 26–30.
20. R. G. Gallager, Principles of Digital Communication (Cambridge University Press, 2008). [CrossRef]
21. T. S. Rappaport, Wireless Communications: Principles and Practice (Prentice Hall, 1996).
22. V. Tarokh, N. Seshadri, and A. R. Calderbank, “Space-time codes for high data rate wireless communication: performance criterion and code construction,” IEEE Trans. Inf. Theory 44(2), 744–765 (1998). [CrossRef]
23. V. Tarokh, A. Naguib, N. Seshadri, and A. R. Calderbank, “Space-time codes for high data rate wireless communication: performance criteria in the presence of channel estimation errors, mobility, and multiple paths,” IEEE Trans. Commun. 47(2), 199–207 (1999). [CrossRef]
24. G. D. Forney Jr and L.-F. Wei, “Multidimensional constellations–Part I: Introduction, figures of merit, and generalized cross constellations,” IEEE J. Sel. Areas Commun. 7(6), 877–892 (1989). [CrossRef]
25. G. D. Forney Jr, “Multidimensional constellations–Part II: Voronoi constellations,” IEEE J. Sel. Areas Commun. 7(6), 941–958 (1989). [CrossRef]
26. G. Ungerboeck, “Trellis-coded modulation with redundant signal sets Part I: Introduction,” IEEE Commun. Mag. 25(2), 5–11 (1987). [CrossRef]
27. G. Ungerboeck, “Trellis-coded modulation with redundant signal sets Part II: State of the art,” IEEE Commun. Mag. 25(2), 12–21 (1987). [CrossRef]
28. M. Safari and M. Uysal, “Do we really need OSTBCs for free-space optical communication with direct detection?” IEEE Trans. Wireless Commun. 7(11), 4445–4448 (2008). [CrossRef]
29. Y. Y. Zhang, H. Y. Yu, J. K. Zhang, and Y. J. Zhu, “On the optimality of spatial repetition coding for MIMO optical wireless communications,” IEEE Commun. Lett. 20(5), 846–849 (2016). [CrossRef]
30. C. Yuen, Y. L. Guan, and T. T. Tjhung, “Orthogonal space-time block code from amicable complex orthogonal design,” in Proceedings of IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP) (IEEE, 2004), pp. 469–472.
31. M. K. Simon and V. A. Vilnrotter, “Alamouti-type space-time coding for free-space optical communication with direct detection,” IEEE Trans. Wireless Commun. 4(1), 35–39 (2005). [CrossRef]
32. C. Yuen, Y. L. Guan, and T. T. Tjhung, “Power-balanced orthogonal space–time block code,” IEEE Trans. Veh. Technol. 57(5), 3304–3309 (2008). [CrossRef]
33. H. Wang, X. Ke, and L. Zhao, “MIMO free space optical communication based on orthogonal space time block code,” Science in China Series F: Information Sciences 52(8), 1483–1490 (2009). [CrossRef]
34. I. Bar-David, “Communication under the Poisson regime,” IEEE Trans. Inf. Theory 15(1), 31–37 (1969). [CrossRef]
35. J. Cao, S. Hranilovic, and J. Chen, “Capacity-achieving distributions for the discrete-time poisson channel - part I: General properties and numerical techniques,” IEEE Trans. Commun. 62(1), 194–202 (2014). [CrossRef]
36. J. Cao, S. Hranilovic, and J. Chen, “Capacity-achieving distributions for the discrete-time poisson channel - part II: Binary inputs,” IEEE Trans. Commun. 62(1), 203–213 (2014). [CrossRef]
37. E. Fisher, I. Underwood, and R. Henderson, “A reconfigurable single-photon-counting integrating receiver for optical communications,” IEEE J. Solid-State Circuits 48(7), 1638–1650 (2013). [CrossRef]
38. D. Chitnis and S. Collins, “A SPAD-based photon detecting system for optical communications,” J. Lightw. Technol. 32(10), 2028–2034 (2014). [CrossRef]
39. J. Zhang, L.-H. Si-Ma, B.-Q. Wang, J.-K. Zhang, and Y.-Y. Zhang, “Low-complexity receivers and energy-efficient constellations for SPAD VLC systems,” IEEE Photon. Technol. Lett. 28(17), 1799–1802 (2016). [CrossRef]
40. C. Gong, Q. Gao, and Z. Xu, “Analysis and design of amplitude modulation for optical wireless communication with shot noise,” in Proceedings of IEEE International Conference on Communications (ICC) (IEEE, 2016), pp. 1–6.
41. X. Liu, C. Gong, S. Li, and Z. Xu, “Signal characterization and receiver design for visible light communication under weak illuminance,” IEEE Commun. Lett. 20(7), 1349–1352 (2016).
42. L.-H. Si-Ma, J. Zhang, B. Wang, and Y.-Y. Zhang, “Hellinger-distance-optimal space constellations for SPAD underwater MIMO-OWC systems,” IEEE Commun. Lett. 21(4), 765–768 (2017). [CrossRef]
43. Y. Li, M. Safari, R. Henderson, and H. Haas, “Optical OFDM with single-photon avalanche diode,” IEEE Photon. Technol. Lett. 27(9), 943–946 (2015). [CrossRef]
44. F. J. Anscombe, “The transformation of Poisson, binomial and negative-binomial data,” Biometrika 35(4), 246–254 (1948). [CrossRef]
45. C. J. Horton and S. N. J. Alexander, Sphere Packings, Lattices and Groups (Springer-Verlag, 1993).
46. J. H. Silverman, A Friendly Introduction to Number Theory (Pearson Education, Inc., 2006).
47. J. J. Watkins, Number Theory: A Historical Approach (Princeton University, 2014).
48. C. F. Gauss, Disquisitiones Arithmeticae (Yale University, 1966).
49. G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers (Oxford University, 1979).
50. E. Grosswald, Representations of Integers as Sums of Squares (Springer-Verlag, 1985). [CrossRef]