Research Article | Open Access

Eunchul Yoon, Duckdong Hwang, Chanki Jang, Janghyun Kim, Unil Yun, "Blind Selected Mapping with Side Information Estimation Based on the Received Pilot Signal", *Wireless Communications and Mobile Computing*, vol. 2018, Article ID 8523680, 9 pages, 2018. https://doi.org/10.1155/2018/8523680

# Blind Selected Mapping with Side Information Estimation Based on the Received Pilot Signal

**Academic Editor:**Ernestina Cianca

#### Abstract

Blind selected mapping (BSLM) is a technique to reduce the peak to average power ratio (PAPR) of OFDM by using a set of scrambling sequences. Since the side information (SI) on the scrambling sequence selected by the transmitter is not directly sent to the receiver in BSLM, the SI needs to be estimated by the receiver. In this paper, a BSLM scheme that applies a new SI estimation approach is proposed. The proposed BSLM scheme estimates the SI by selecting the minimum of the candidate metrics that are computed with the received pilot signal, the pilot symbol sequence, and the scrambling sequences. It is shown by simulation that the proposed BSLM scheme yields better SI error rate (SIER) and bit error rate (BER) performance and lower implementation complexity than four traditional BSLM schemes.

#### 1. Introduction

Orthogonal frequency division multiplexing (OFDM) is an attractive transmission technique for broadband communication. OFDM can overcome the detrimental impact of intersymbol interference efficiently without requiring complex equalization [1]. Since OFDM is considered as appropriate for high spectral efficiency and scalable to low complexity receivers, it is well suited for the 5G cellular network [2]. However, OFDM has a critical problem of high peak to average power ratio (PAPR). When the signal with high PAPR passes nonlinear devices such as digital-to-analog converter and high power amplifier, the problems of quantization noise, in-band signal distortion, and out-of-band radiation can be severely aggravated [3, 4]. Therefore, various PAPR reduction schemes have been proposed (e.g., [5–17]). Selected mapping (SLM) introduced in [18] has been known as an efficient PAPR reduction technique with low implementation complexity. SLM generates multiple candidate transmission sequences by multiplying the data symbol sequence by multiple scrambling sequences and selects a candidate transmission sequence that can lead to the lowest PAPR. To recover the data symbol sequence, the receiver needs the side information (SI) on the scrambling sequence selected by the transmitter. Since correct delivery of the SI is particularly crucial for reliable recovery of the data symbol sequence, the SI should be protected with strong channel coding, which decreases the effective data transfer rate. In order to apply the SLM technique without directly transmitting the SI from the transmitter to the receiver, several blind SLM (BSLM) schemes have been studied. In [19], the maximum likelihood (ML) detection was introduced for BSLM and a simplified ML detection approach was suggested to reduce the complexity of the ML detection. Although the simplified ML detection substantially reduced the complexity, it still had high complexity. In [20], the SI was embedded in the scrambling sequence by partitioning each scrambling sequence into multiple subsequences and multiplying those subsequences by a distinct set of phase-offsets. While the BSLM scheme in [20] had low complexity, it showed slightly worse bit error rate (BER) performance than other BSLM schemes at a low signal to noise ratio (SNR). Reference [21] proposed a BSLM scheme that embedded the SI in the scrambling sequence by allocating power-extended scrambling symbols. While it considerably reduced complexity, its BER performance was deteriorated with a large symbol constellation size. Reference [22] proposed a BSLM scheme that estimated the SI by assuming high autocorrelation between the channel frequency response (CFR) coefficients over the adjacent pilot subcarriers. The BSLM scheme in [22] showed low complexity and performed comparatively well even with a large symbol constellation size. However, our simulation results showed that, given a multipath fading channel consisting of a large number of effective channel taps with nontrivial tap powers, the BER performance degraded severely. It is because a greater number of channel taps can cause more fluctuations in the CFR amplitudes leading to more weakened autocorrelation between the CFR coefficients over the adjacent pilot subcarriers. Reference [23] proposed a BSLM scheme to decrease the number of the subcarriers that were used to indirectly help to estimate the SI. However, at the expense of reducing the number of subcarriers, the BER performance of the BSLM scheme became deteriorated.

In this paper, a BSLM scheme that applies a new SI estimation approach is proposed. The proposed BSLM scheme estimates the SI by selecting the minimum of the candidate metrics that are computed with the received pilot signal, the pilot symbol sequence, and the scrambling sequences. It is shown by simulation that the proposed BSLM scheme yields better SI error rate (SIER) and BER performance and lower implementation complexity than four traditional BSLM schemes.

The remainder of this paper is organized as follows. In Section 2, the system model is explained. In Section 3, the previous BSLM schemes are introduced. In Section 4, the proposed BSLM scheme is explained and its complexity is analyzed. In Section 5, the SIER and BER performances of the proposed BSLM scheme and its implementation complexity are compared with those of the conventional BSLM schemes. Finally, in Section 6, this paper is concluded.

*Notations. * and denote the transpose operator and the complex conjugate transpose operator, respectively. denotes the uniform distribution with its random variable located between and . denotes the Frobenius norm operator. denotes a diagonal matrix whose diagonal components are given by the components of a vector . denotes the component of a matrix in the -th row and the -th column. Finally, denotes an identity matrix.

#### 2. System Model

We consider an OFDM system with a block size , where and denote the number of data subcarriers and the number of pilot subcarriers, respectively. The channel impulse response is modeled by a finite impulse response filter with maximum length aswhere are uncorrelated complex Gaussian random variables with zero mean. The average power of for is given by the exponential power delay profile (PDP) [24, 25]:where denotes a power normalization factor that guarantees and denotes the sampling period with being the system bandwidth, and denotes the root mean square channel delay spread. With a large valued , the average channel tap power decays slowly with increasing , which means there are a larger number of effective channel taps with nontrivial tap powers. The CFR coefficient vector,is given by the -point discrete Fourier transform of the channel impulse response. If a matrix is defined byfor and , the CFR coefficient vector is given byThe data symbols are chosen from an -ary QAM constellation, , which has unit average symbol power, and the pilot symbols are chosen from the BPSK constellation, . Data and pilot symbols are mapped to an OFDM block to form an OFDM symbol sequence:With the pilot subcarrier indices given by , and , the pilot symbol sequence can be written asand the CFR coefficient vector for the pilot subcarriers can be written asWith the data subcarrier indices given by , and , the data symbol sequence and the CFR coefficient vector for the data subcarriers are defined analogously. Note that can be written in terms of aswhere is an matrix given byfor and . The time-domain sample sequence,is given by the inverse discrete Fourier transform (IDFT) of . The PAPR of OFDM is defined byIn SLM, the transmitter generates scrambled sequences in frequency domain by multiplying by scrambling sequences:for . The scrambling sequences are assumed to be known to both the transmitter and the receiver. The transmitter generates candidate time-domain sample sequences by taking the IDFT of each scrambled sequence and selects a candidate time-domain sample sequence that can lead to the minimum PAPR. In BSLM, the receiver estimates the SI, i.e., the index of the scrambling sequence selected by the transmitter. If the SI is , the received signal over the -th subcarrier for can be written aswhere denotes the additive Gaussian noise with mean zero and variance .

#### 3. Previous BSLM Schemes

In this section, traditional BSLM schemes are described for comparison purposes.

##### 3.1. The Conventional ML-Based BSLM Scheme

The conventional ML-based BSLM scheme can be written aswhere denotes the CFR coefficient over the -th data subcarrier estimated at the receiver. In the above, was defined by , which satisfied the condition of . This scheme requires computations of the ML metric .

##### 3.2. The Simplified ML-Based BSLM Scheme

Reference [19] suggested a simplified ML-based BSLM scheme. The scheme can be summarized aswhereIn [19], was defined by , which satisfied the condition of . In order to reduce implementation complexity, all the values of can be computed and stored in memory to be used in search of . Once the optimal is determined in (16), can be estimated by taking the optimal arguments of the minimization problems in (17) in search of . This scheme requires computations of the ML metric .

##### 3.3. The BSLM Scheme Using Phase Offset Factors

Reference [20] proposed an BSLM scheme that partitioned the scrambling sequence into subsequences and multiplied each subsequence by a phase offset factor to enable the receiver to estimate effectively. In [20], was defined by , which satisfied the condition of . The phase offset sequences were given byfor , which were designed to have the largest minimum Hamming distance. Since the -th scrambling subsequence is multiplied by the phase offset factor , the received signal over the -th data subcarrier for , can be written aswhere . With the definition of , the BSLM scheme in [20] can be summarized aswherefor . In order to reduce implementation complexity, all the values of can be computed and stored in memory to be used in search of . Once the optimal is determined in (20), can be estimated by taking the optimal arguments of the minimization problems in (21) in search of with , and then by converting into based on . This scheme requires computations of the ML metric .

##### 3.4. The BSLM Scheme Using CFR Autocorrelation

Reference [22] proposed an BSLM scheme that exploited CFR autocorrelation to estimate under the assumption of high correlation between two CFR coefficients over two neighboring pilot subcarriers. In [22], the scrambling sequences were defined byfor , and used to scramble instead of . In [22], was defined by , which satisfied . The received signal carrying the -th pilot symbol for , can be written asThe BSLM scheme in [22] can be summarized aswhereNote that corresponds to the CFR coefficient estimated over the -th pilot subcarrier, whereas with is irrespective of the CFR coefficient over the -th pilot subcarrier. If the autocorrelation between two CFR coefficients over two neighboring pilot subcarriers is very high, it follows that . Therefore, it is highly probable that in (25) can be maximized with . In order to reduce implementation complexity, all the values of can be computed and stored in memory to be used in search of . Once the optimal is determined in (24), can be estimated byfor . This scheme performs well when correlation between two CFR coefficients over two neighboring pilot subcarriers is high. However, its performance can be seriously degraded with a large number of channel multipaths because the autocorrelation between two CFR coefficients over two neighboring pilot subcarriers decreases with an increasing number of channel multipaths. This phenomenon will be certified by simulation in Section 5.

#### 4. Proposed BSLM Scheme

In the proposed BSLM scheme, is chosen as under a constraint of . The scrambling sequence for the pilot symbols is defined asThe received signal on the pilot subcarriers can be written aswhereIn the proposed scheme, the received pilot signal, , is multiplied by and then multiplied by every matrix of for to yield If , the expression in (31) reduces towhereWe assume that the components of are independent and identically distributed Gaussian random variables with mean zero and variance . If , based on the least square method [26], the channel impulse response can be approximated in terms of asFrom (9) and (34), the estimated CFR coefficient vector for the pilot subcarriers can be written asDisregarding the noise term and regarding as the right SI , the received pilot signal, , can be approximated in terms of asBy substituting in (31) into (35) and substituting in (35) into (36) sequentially, it follows thatwhereThe error vector between in (37) and in (29) is given byThe total error between and for is defined bySince , in the case of , the first term in the right-hand side of (39) reduces to zero. It is easy to show that . Therefore, in the case of , is given byThe averaged total error with respect to reduces to In the case of , is given by (37) and the error vector between and is given by (39). Then, the averaged total error with respect to reduces toNote that because can be written in a positive semidefinite matrix structure. Since , , and , it is highly probable that the value of which minimizes corresponds to . From (37), it follows thatFrom (38), (40), and (44), it is easy to know that the SI can be estimated at the receiver by minimizing the candidate metric aswhere To reduce computational complexity, the matrices, , can be computed off-line and saved in memory for later use. Then, the SI estimation based on (45) requires complex multiplications (CMs), complex additions (CAs), operations of , and real additions (RAs). Since one CM is translated into 4 real multiplications (RMs) and 2 RAs, one CA is translated into 2 RAs, and one operation of is translated into 2 RMs and 1 RA [22], the computational complexity of the SI estimation based on (45) corresponds to RMs and RAs. Once is estimated, the -ary QAM symbols can be estimated by the symbol-based ML detection, which requires additional computation of RMs and RAs.

#### 5. Numerical Results

For simulation, the system parameters are chosen as , , , and . In order to accommodate the nonlinear amplification characteristics, Rapp’s solid-state power amplifier (SSPA) model [27] is adopted with smoothness factor 3 and input back-off factor 6dB. In this section, the proposed BSLM scheme and four traditional BSLM schemes are compared by simulation:(i)PROPOSED: the BSLM scheme presented in the previous section.(ii)BSLM-I: the BSLM scheme in [19], which applies the simplified ML decoding.(iii)BSLM-II: the BSLM scheme in [20], which applies partition and phase offset techniques with its partition number given by .(iv)BSLM-III: the BSLM scheme in [21], which applies the power-extended scrambling.(v)BSLM-IV: the BSLM scheme in [22], which uses the autocorrelation between the CFR coefficients over the adjacent pilot subcarriers.

Figure 1 shows the SIER results of the five BSLM schemes in two cases of and when . It can be seen that PROPOSED yields lower SIERs than other BSLM schemes in both the cases of and . In the case of , BSLM-IV yields SIERs comparable to those of PROPOSED. However, in the case of , BSLM-IV degrades severely yielding the worst SIER performance. It is because a large valued incurs a greater number of effective channel taps leading to more severely fluctuating CFT amplitudes and more weakened CFR autocorrelation. Since the SI estimation in BSLM-IV is based on the high CFR autocorrelation assumption, the SIER performance of BSLM-IV degrades as increases. Contrary to BSLM-IV, PROPOSED shows good SIER performance regardless of .

**(a)**

**(b)**

In Figure 2, the complementary cumulative distribution functions (CCDFs) of the five BSLM schemes are plotted for various values of when the transmitted symbols are chosen from the 16QAM constellation. The CCDF is defined as which means the probability that the PAPR is greater than a given quantity, . For referential purpose, NONSLM is included in the figure, which denotes the performance of the original OFDM that does not apply any SLM scheme. From the figure, it can be seen that all the BSLM schemes reduce the PAPR substantially when compared with NONSLM. Note that the five BSLM schemes show almost identical PAPR reduction. The amounts of PAPR reduction by the five BSLM schemes tend to increase as the number of scrambling sequences increases.

Figure 3 shows the BER results of PROPOSED for various values of when and the transmitted symbols are chosen from the 16QAM constellation. For referential purpose, NONSLM is included in the figure, which denotes the performance of the original OFDM that does not apply any SLM scheme. It can be seen that PROPOSED improves the BER performance substantially for all the values of when compared to NONSLM. The BER improvement by PROPOSED with increasing appears to reach saturation when . It implies that increasing the number of the scrambling sequences over brings only trivial BER improvement while it further complicates SI estimation. Therefore, for the remaining part of this section, is chosen in simulation.

Figure 4 shows the BER results of the five BSLM schemes in two cases of and when the transmitted symbols are chosen from the 16QAM constellation. In the case of , it can be seen that PROPOSED and BSLM-IV yield lower BERs than other BSLM schemes. In the case of , BSLM-IV yields severely degraded BERs even worse than the BERs of NONSLM. On the contrary, in the case of , PROPOSED outperforms all the other BSLM schemes. In both the cases of and , BSLM-I shows BERs comparable to those of PROPOSED at a high SNR. Since PROPOSED has lower implementation complexity than BSLM-I, PROPOSED can be regarded as more efficient than BSLM-I. BSLM-II shows slightly worse BER performance than other BSLM schemes at a low SNR.

**(a)**

**(b)**

Figure 5 shows the BER results of the five BSLM schemes under the same simulation conditions as in Figure 4 except that the transmitted symbols are chosen from the 64QAM constellation. By comparing Figure 4 and Figure 5, it can be seen that BSLM-I, BSLM-II, BSLM-IV, and PROPOSED show the same BER trends in both figures. However, BSLM-III degrades severely as the symbol constellation size increases. In addition, it can be seen that BSLM-III degrades further as increases. It is because a large symbol constellation size or a large valued incurs a large variation of the received signal power, which makes it difficult for BSLM-III to estimate the SI correctly by using the locations of the power-extended scrambling symbols.

**(a)**

**(b)**

In Figure 6, the implementation complexities of the BSLM schemes are compared. As in [20], the detection complexity reduction ratio (DCRR) is defined asNote that a positive DCRR value indicates that PROPOSED has lower implementation complexity than the conventional BSLM scheme being considered. In computing the complexity, both the complexity of SI estimation and that of symbol detection are considered together. The RM-based and RA-based DCRRs of BSLM-I, BSLM-II, BSLM-III, and BSLM-IV are computed for various QAM constellation sizes by referring to the complexity analysis in [20, 22] and plotted in Figure 6. It can be seen that the RM-based and RA-based DCRRs of BSLM-I, BSLM-II, and BSLM-III are positive for all the QAM constellation sizes. In addition, the RM-based and RA-based DCRRs of BSLM-IV are positive when . This implies that PROPOSED has lower implementation complexity than the four traditional BSLM schemes. The RM-based and RA-based DCRRs of BSLM-I, BSLM-II, BSLM-III, and BSLM-IV tend to increase as the QAM constellation size increases. This implies that the amount of complexity reduction by PROPOSED increases as the symbol constellation size increases.

**(a)**

**(b)**

#### 6. Conclusion

A BSLM scheme that applied a new SI estimation approach based on the pilot signal was proposed. The proposed BSLM scheme estimated the SI by selecting the minimum of the candidate metrics that were computed with the received pilot signal, the pilot symbol sequence, and the scrambling sequences. It was shown by simulation that the proposed BSLM scheme yielded better SIER and BER performance and lower implementation complexity than four traditional BSLM schemes.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This paper was supported by Konkuk University in 2016.

#### References

- Y. Wu and W. Y. Zou, “Orthogonal frequency division multiplexing: A multi-carrier modulation scheme,”
*IEEE Transactions on Consumer Electronics*, vol. 41, no. 3, pp. 392–399, 1995. View at: Publisher Site | Google Scholar - N. Bhushan, T. Ji, O. Koymen et al., “5G Air Interface System Design Principles,”
*IEEE Wireless Communications Magazine*, vol. 24, no. 5, pp. 6–8, 2017. View at: Google Scholar - J. Tellado,
*Multicarrier Modulation with Low PAR: Applications to DSL and Wireless*, Kluwer, Norwell, Mass, USA, 2000. - T. Jiang and Y. Wu, “An overview: peak-to-average power ratio reduction techniques for OFDM signals,”
*IEEE Transactions on Broadcasting*, vol. 54, no. 2, pp. 257–268, 2008. View at: Publisher Site | Google Scholar - A. E. Jones, T. A. Wilkinson, and S. K. Barton, “Block coding scheme for reduction of peak to mean envelope power ratio of multicarrier transmission schemes,”
*IEEE Electronics Letters*, vol. 30, no. 25, pp. 2098-2099, 1994. View at: Publisher Site | Google Scholar - J. Zakaria and M. F. M. Salleh, “PAPR reduction scheme: Wavelet packetbased PTS with embedded side information data scheme,”
*IET Communications*, vol. 11, no. 1, pp. 127–135, 2017. View at: Publisher Site | Google Scholar - X. Li and L. J. Cimini Jr., “Effects of clipping and filtering on the performance of OFDM,”
*IEEE Communications Letters*, vol. 2, no. 5, pp. 131–133, 1998. View at: Publisher Site | Google Scholar - J. Tellado,
*Peak to Average Power Ratio Reduction for Multicarrier Moulation [Ph.D. thesis]*, University of Stanford, Stanford, 1999. - X. Wang, T. T. Tjhung, and C. S. Ng, “Reduction of peak-to-average power ratio of OFDM system using a companding technique,”
*IEEE Transactions on Broadcasting*, vol. 45, no. 3, pp. 303–307, 1999. View at: Publisher Site | Google Scholar - V. Tarokh and H. Jafarkhani, “On the computation and reduction of the peak-to-average power ratio in multicarrier communications,”
*IEEE Transactions on Communications*, vol. 48, no. 1, pp. 37–44, 2000. View at: Publisher Site | Google Scholar - Y. Kou, W.-S. Lu, and A. Antoniou, “New peak-to-average power-ratio reduction algorithms for multicarrier communications,”
*IEEE Transactions on Circuits and Systems I: Regular Papers*, vol. 51, no. 9, pp. 1790–1800, 2004. View at: Publisher Site | Google Scholar - S. H. Han and J. H. Lee, “PAPR reduction of OFDM signals using a reduced complexity PTS technique,”
*IEEE Signal Processing Letters*, vol. 11, no. 11, pp. 887–890, 2004. View at: Publisher Site | Google Scholar - H. Kim, E. Hong, C. Ahn, and D. Har, “A pilot symbol pattern enabling data recovery without side information in PTS-based OFDM systems,”
*IEEE Transactions on Broadcasting*, vol. 57, no. 2, pp. 307–312, 2011. View at: Publisher Site | Google Scholar - Y. Rahmatallah and S. Mohan, “Peak-to-average power ratio reduction in ofdm systems: a survey and taxonomy,”
*IEEE Communications Surveys & Tutorials*, vol. 15, no. 4, pp. 1567–1592, 2013. View at: Publisher Site | Google Scholar - X. Huang, G. Tan, Q. Xu, N. Xu, and S. Wang, “A kind of PAPR reduction method based on pruning WPM and PTS technology,”
*Journal of Electronics (China)*, vol. 30, no. 3, pp. 261–267, 2013. View at: Google Scholar - L. Li, D. Qu, and T. Jiang, “Partition optimization in LDPC-coded OFDM systems with PTS PAPR reduction,”
*IEEE Transactions on Vehicular Technology*, vol. 63, no. 8, pp. 4108–4113, 2014. View at: Publisher Site | Google Scholar - H.-S. Joo, K.-H. Kim, J.-S. No, and D.-J. Shin, “New PTS Schemes for PAPR Reduction of OFDM Signals Without Side Information,”
*IEEE Transactions on Broadcasting*, vol. 63, no. 3, pp. 562–570, 2017. View at: Publisher Site | Google Scholar - R. W. Bäuml, R. F. H. Fischer, and J. B. Huber, “Reducing the peak-to-average power ratio of multicarrier modulation by selected mapping,”
*IEEE Electronics Letters*, vol. 32, no. 22, pp. 2056-2057, 1996. View at: Publisher Site | Google Scholar - A. D. S. Jayalath and C. Tellambura, “SLM and PTS peak-power reduction of OFDM signals without side information,”
*IEEE Transactions on Wireless Communications*, vol. 4, no. 5, pp. 2006–2013, 2005. View at: Publisher Site | Google Scholar - H.-S. Joo, S.-J. Heo, H.-B. Jeon, J.-S. No, and D.-J. Shin, “A new blind SLM scheme with low decoding complexity for OFDM systems,”
*IEEE Transactions on Broadcasting*, vol. 58, no. 4, pp. 669–676, 2012. View at: Publisher Site | Google Scholar - S. Y. Le Goff, S. S. Al-Samahi, B. K. Khoo, C. C. Tsimenidis, and B. S. Sharif, “Selected mapping without side information for PAPR reduction in OFDM,”
*IEEE Transactions on Wireless Communications*, vol. 8, no. 7, pp. 3320–3325, 2009. View at: Publisher Site | Google Scholar - E. Hong, H. Kim, K. Yang, and D. S. Har, “Pilot-aided side information detection in SLM-based OFDM systems,”
*IEEE Transactions on Wireless Communications*, vol. 12, no. 7, pp. 3140–3147, 2013. View at: Publisher Site | Google Scholar - T. Kitazawa and K. Ohuchi, “On blind estimation with reduced complexity in CSS-OFDM systems,” in
*Proceedings of the 2016 IEEE International Symposium on Signal Processing and Information Technology, ISSPIT 2016*, pp. 120–125, cyp, December 2016. View at: Google Scholar - B. Wang, Y. Wu, F. Han, Y.-H. Yang, and K. J. R. Liu, “Green wireless communications: A time-reversal paradigm,”
*IEEE Journal on Selected Areas in Communications*, vol. 29, no. 8, pp. 1698–1710, 2011. View at: Publisher Site | Google Scholar - M. Emami, M. Vu, J. Hansen, A. J. Paulraj, and G. Papanicolaou, “Matched filtering with rate back-off for low complexity communications in very large delay spread channels,” in
*Proceedings of the Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers*, pp. 218–222, usa, November 2004. View at: Google Scholar - M. Morelli and U. Mengali, “A comparison of pilot-aided channel estimation methods for OFDM systems,”
*IEEE Transactions on Signal Processing*, vol. 49, no. 12, pp. 3065–3073, 2001. View at: Publisher Site | Google Scholar - C. Rapp, “Effects of HPA-nonlinearity on 4-DPSK OFDM Signal for a Digital Sound Broadcasting System,” in
*Proceedings of 2nd European Conference on Satellite Commununications*, pp. 179–184, Liege, Belgium, 1991. View at: Google Scholar

#### Copyright

Copyright © 2018 Eunchul Yoon et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.