Modeling and Analysis of MmWave Communications in Cache-enabled HetNets

01/26/2018 ∙ by Wenqiang Yi, et al. ∙ 0

In this paper, we consider a novel cache-enabled heterogeneous network (HetNet), where macro base stations (BSs) with traditional sub-6 GHz are overlaid by dense millimeter wave (mmWave) pico BSs. These two-tier BSs, which are modeled as two independent homogeneous Poisson Point Processes, cache multimedia contents following the popularity rank. High-capacity backhauls are utilized between macro BSs and the core server. A maximum received power strategy is introduced for deducing novel algorithms of the success probability and area spectral efficiency (ASE). Moreover, Monte Carlo simulations are presented to verify the analytical conclusions and numerical results demonstrate that: 1) the proposed HetNet is an interference-limited system and it outperforms the traditional HetNets; 2) there exists an optimal pre-decided rate threshold that contributes to the maximum ASE; and 3) 73 GHz is the best mmWave carrier frequency regarding ASE due to the large antenna scale.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

To support the explosive data traffic of future fifth-generation (5G) cellular networks, numerous researches [1, 2, 3, 4] have paid attention to an innovative framework that densifies the traditional networks with massive small base stations (BSs) or exploiting power domain. However, the improvement of these heterogeneous networks (HetNets) is mainly restricted to the capacity of the backhauls [1]. A recent study [5] has relaxed the limitation by equipping caches at all BSs to store the most popular files as only the minority of multimedia contents is frequently requested by the majority of customers in the real world. Accordingly, the aforementioned cache-enabled HetNets have been studied in various papers. The primary work [6] analyzed the energy efficiency and throughput of cache-enabled cellular networks with a regular hexagonal grid. Since stochastic geometry is a useful tool to acquire the networks’ randomness, modeling a tier of BSs in small cell networks or HetNets with a Poisson Point Process (PPP) is more accurate than the traditional hexagonal scenario [7]. Under this condition, the throughput of multi-tier cache-enabled HetNets was discussed in [1], where macro BSs connected to the core networks via backhauls and small cell devices cached the content through wireless broadcasting. Unfortunately, the further analysis on the impact of backhaul capacity was omitted, which is the key parameter when comparing with the conventional HetNets.

In addition to the network densification, another key capacity-increasing technology is exploiting new spectrum bands, such as millimeter wave (mmWave) [8]. Two distinctive characteristics of mmWave are small wavelength and the sensitivity to blockages [9, 10]. Thanks to short wavelength, steerable antennas can be deployed at devices to enhance the directional array gain [9]. On the other side, the sensitivity gives rise to severe penetration loss for mmWave signals when passing through building exteriors [10], so it is unrealistic to expect the outdoor-to-indoor coverage from macro mmWave BSs. An ingenious hybrid network is created to solve this issue, where mmWave transmitters contribute to the ultra-fast data rate in short-range small cells, and sub-6 GHz BSs provide the universal coverage [11]. The same with cache-enabled HetNets, stochastic geometry has also been wildly utilized in mmWave networks, where the locations of transceivers were modeled following PPPs [12]. With the aid of such structure, the primary article [12] introduced a stochastic blockage model to represent the actual mmWave communication environment, but the antenna pattern was over-simplified as a flat-top model. Then the authors in [13] proposed an actual antenna pattern for increasing the accuracy. Considering the hybrid HetNets, a tractable structure combining mmWave with sub-6 GHz was analyzed in [11], which performed close to the reality.

As discussed above, although HetNets with caches have been fully analyzed under traditional sub-6 GHz networks, there is still lack of articles on a hybrid system with mmWave small cells. Since mmWave is able to provide fast data rate in short-distance networks [14], adopting mmWave into a dense pico tier of cache-enabled HetNets is an promising way to increase the throughput of 5G cellular networks. The other benefit of such hybrid HetNets is no mutual interferences due to applying distinctive carrier frequencies between tiers. These advantages motivate us to create this paper. In contrast to [1], we introduce fiber-connections between macro BSs and the multimedia server to evaluate the impact of backhaul capacity. Then, due to the employment of mmWave, the propagation environment and antenna beamforming pattern in the small cells are replaced by Nakagami fading and actual antenna arrays, respectively. Moreover, we have compared the performance of various mmWave frequencies in this paper. On the other hand, unlike noise-limited assumption in [11], we demonstrate that in the dense mmWave networks, the system becomes an interference-limited scenario. With the aid of the content placement, we conclude that the throughput of our system is conditionally decided by the storage capacity. The main contributions are summarized as follows: 1) we discuss the Laplace transform of interference for traditional sub-6 GHz macro cells and mmWave small cells with the actual antenna pattern; 2) the novel algorithms for the success probability and area spectral efficiency (ASE) under Maximum Received Power (Max-RP) scheme are derived; 3) our model is an interference-limited system and there is an optimum value of rate requirement for obtaining the maximum ASE; and 4) 73 GHz is the best mmWave carrier frequency under the considered association strategy.

Ii System Model

Ii-a Network Architecture

Fig. 1: Layouts of the proposed cache-enabled two-tier hybrid HetNets with traditional macro cells and mmWave small cells.

In this paper, we present a two-tier cache-enabled hybrid HetNet consisting of numerous macro and pico BSs as shown in Fig. 1. The locations of BSs in -th tier are modeled following an independent homogeneous PPP with a density , where and represent the macro and pico tier, respectively. In our model, a typical user is located at the origin of a plane so that the distribution of distance from the typical user to its nearest BS in -th tier will follow . Apparently, the number of pico BSs in the real HetNets is much more than that of macro BSs, so we consider . In order to compare the performance of our proposed networks with the traditional HetNets, we provide a server to supply the less-popular contents. This server connects to each macro BS through a high-capacity wired backhaul.

Moreover, different carrier frequencies are employed in our system. The macro BSs connect to receivers with traditional sub-6 GHz, while mmWave is utilized between the pico BS and its corresponding user. As a result, there is no interference between these two tiers. However, in the same tier, other active BSs apart from the corresponding BS that serves the typical user will be the source of interferences. Additionally, the quantity of users is assumed to be large enough to ensure that all BSs are active when the corresponding BS is transmitting messages to the typical user.

Ii-B Blockage Model

In tier 1, since macro BSs employ sub-6 GHz as its carrier frequency, the path loss law with a distance is same with that in traditional cellular networks as shown below

(1)

where is the path loss exponent and is the intercept for the macro tier.

In tier 2, the communication environment is changed due to the employment of mmWave. More specifically, the transmitting routes of pico BSs are divided into line-of-sight (LOS) group and non-line-of-sight (NLOS) group and each of them has its unique path loss law. Note that the density of the pico tier is huge, which indicates that tier 2 can be visualized as dense mmWave networks. One blockage may obstruct all BSs behind it under this case. Therefore, we consider an LOS ball to model the blockage process. A recent study [15] has advocated that such blockage pattern fits the real-world scenarios better than other models. In terms of the NLOS links, the authors in [12] have demonstrated that the impact of NLOS signals is so weak that it can be ignored in mmWave networks, only LOS signals will be considered in this paper.

As a consequence, we define the LOS ball with a radius , which represents the departure from nearby obstacles. The probability of LOS links will be one inside the ball and zero outside the area. The path loss law in tier 2 with a distance is shown as follows

(2)

where and are the path loss exponent and the intercept of LOS links in tier 2, respectively. is the unit step function.

Ii-C Directional Beamforming

In -th tier, we employ antenna arrays composed of elements at all cache-enabled BSs and the transmit power is assumed to be a constant with . The uniform linear array formed by antenna elements is utilized at all pico BSs. However, we only consider an omnidirectional antenna pattern at macro BSs () and the typical user for tractability of the analysis [11].

Directional antenna arrays deployed at the pico BSs will supply substantial beamforming gains to compensate the path loss, so the received signal for -th tier at the typical user can be expressed as follows

(3)

where interfering BSs are located at when the typical user is receiving the message from the serving BS at

. The channel vector from the BS to the typical user and the beamforming vector of the BS in

-th tier are denoted by and , respectively. represents the thermal noise with power .

Combining with the aforementioned assumptions, the product of the fading gain and beamforming gain of the BS located at in -th tier is shown as below [16]

(4)

where is the small fading term for -th tier. is the spatial angle of departure (AoD) from the interfering BS to the typical user, and is the spatial AoD between the BS at location and its corresponding receiver. is the array gain function. More specifically, an actual array pattern is employed at pico BSs so that , where

is a uniformly distributed random variable over

. and are the antenna spacing and wavelength, respectively [16]. On the other hand, the array gain function for macro BSs is due to the omnidirectional antenna pattern.

Ii-D Channel Model

Assuming that all active BSs have the full knowledge of the AoD from itself to the typical user, the serving BS will align the antenna beam towards the typical user for achieving maximum directivity gain . In our model, since the corresponding BS is only interfered by the active BSs located in the same tier, the interference for -th tier can be expressed as follows

(5)

Therefore, the signal-to-interference-plus-noise-ratio (SINR) of -th tier at the typical user is given by

(6)

where follows independent Nakagami fading due to utilizing mmWave and the parameter of Nakagami fading is considered to be a positive integer for simplifying the analysis [12]. Therefore, is a normalized Gamma random variable. On the other side, we assume a Rayleigh fading model for the macro tier so that the fading parameter .

Ii-E Cache-enabled Content Access Protocol

In this paper, we assume that a static multimedia content catalog containing files is stored at the server. Every macro BSs has a restricted storage with files, and each of pico BSs has a weaker ability to cache files, namely, . High-speed backhauls are employed for connecting the core server to macro BSs like traditional HetNets, and the backhaul capacity is denoted by . When the data traffic load of our proposed networks becomes low, the content will be cached in a sequence of its popularity rank at all BSs via broadcasting until the storage is fully occupied. The content with -th popularity rank can be represented by the Zipf distribution as shown below

(7)

where

is the skew parameter of the probability distribution.

With the aid of Max-RP scheme, the access protocol in our paper is defined as follows.

Access Protocol: When the typical user demands a multimedia file, it communicates with the proximate macro or pico BS with maximum received power. However, if the demanded content is absent from the corresponding BS due to limited storage capacity, the typical user will request that file from the server via the nearest macro BS.

Iii Laplace Transform of Interference

The expected value of interference can be derived by Laplace Transform of Interference, using which the success probability and ASE for the considered HetNets will be deduced.

Iii-a Laplace Transform of Interference Analysis in Tier 1

In the macro tier, we utilize sub-6 GHz as the carrier frequency, and the fading channel is assumed to be Rayleigh fading. Therefore, the exact expression for Laplace transform of interference can be expressed as below.

Lemma 1

The Laplace transform of interference in the macro tier is given by

(8)

where is the Gauss hypergeometric function.

Proof:

The proof procedure is similar as Appendix B in [1], but fading parameter equals to 1 in our case due to Rayleigh fading assumption.

Iii-B Laplace Transform of Interference Analysis in Tier 2

Since the path loss exponent of LOS links is no less than 2 for the practical mmWave networks, we will divide the analysis into two conditions ( and ) in order to achieve several closed-form equations.

Lemma 2

Under the condition , the -th Laplace transform of interference in the pico tier is as follows

(9)

where

(10)

are Gauss-Chebyshev nodes over , and is a tradeoff parameter between the accuracy and complexity [17, 18]. When , the equality is established.

Numerous actual channel measures [19, 20] have indicated that the path loss exponent of LOS link is under various carrier frequencies, e.g. GHz, GHz and GHz, so we are more interested in the Laplace transform of interference under the condition of . The equation (2) will be changed into

(11)

where

(12)
Proof:

See Appendix A.

Iv Success Probability and Area Spectral Efficiency Analysis

Before discussing the success probability, we define our Content Placement as follows: All BSs in different tiers will choose the most popular contents to store. As a result, the probability of the condition that -th tier BSs cache the same copy of -th ranked file can be expressed as below

(13)

Therefore, a set of -th tier BSs containing the -th ranked file will form an independent non-homogeneous PPP with the density

From the costumer’s perspective, the success probability is an important parameter to appraise the quality of service. In our cache-enabled HetNets, the instantaneous data rate at the typical user exceeding the pre-decided rate threshold will contribute to the success probability [1]. As discussed in the previous sections, we conclude that the proposed system has two different processes in sending the multimedia contents: 1) User Association Mode; and 2) Server Mode. We will discuss these two modes in details as below.

Iv-a User Association Mode

The typical user will choose the macro or pico BS with maximum received power as its serving BS. We define the Max-RP association probability when the typical user connects with -th tier as follows

(14)

where and

Therefore, the probability density function (PDF) of the distance

between the typical user and its serving BS with -th ranked file in -th tier is changed into [21]

(15)

where and .

Iv-B Server Mode

In server mode, the nearest Macro BS will act as the relay to retransmit the multimedia file from the server to the typical user. Under this condition, the backhaul capacity will restrict the performance of our system. To simplify the notation, we first derive the coverage probability of tier 1 without considering the cache capacity.

Lemma 3

Since in various articles [11, 1], the noise can be ignored in the traditional cellular networks with sub-6 GHz, we only consider the signal-to-interference-ratio (SIR) instead of SINR for analyzing the performance of tier 1, namely, . Under this condition, the coverage probability for no-caching tier 1 is given by

(16)
Proof:

Note that the SIR coverage probability where due to Rayleigh fading assumption. Thus such coverage probability can be expressed as . Note that , the coverage probability will be simplified as above.

Corollary 1

The success probability in server mode is given by

(17)
Proof:

In server mode, if is smaller than , the will be zero as the system rate is not large enough for transmitting the content. Moreover, if is larger than , the success probability for -th ranked file under this mode will be . We sum up the -th to -th elements to calculate as shown above.

Iv-C Success Probability

Based on the Laplace transform of the interference, we are able to calculate the content-related coverage probabilities of two tiers, which is the basement for analyzing the success probability.

Lemma 4

For -th ranked content, the coverage probability for -th tier under Max-RP scheme is shown below

(18)
(19)

where

(20)
Proof:

Note that the probability of distance under this user association is shown in equation (15), with the similar calculation process with Appendix A from [12] and Lemma 3, coverage probability for two tier can be expressed as above.

As mentioned in the beginning of this section, the universal success probability can be defined as below

(21)
Theorem 1

With the aid of Corollary 1 and Lemma 4, the success probability of Max-RP will be as follows

(22)
Proof:

When the required rate is , the pre-decided SINR threshold is , so the success probability for -th ranked content in the first popularity rank is . Considering the whole multimedia contents, the success probability will be calculated as above.

Iv-D Area Spectral Efficiency

The ASE is the average instantaneous data rate transmitted in unit bandwidth and unit area. Assuming that Gaussian Codebooks are utilized for all transmissions, we are able to define ASE with the aid of Shannon’s Capacity Formula. It is expressed as follows , where is the SINR coverage probability of the considered networks and denotes the active BSs’ density [22].

Proposition 1

The ASE for Max-RP user association strategy is given by

(23)

V Numerical Results

The general network settings for our system are shown in Table. I [11, 1, 12] and the reference distance for the intercept is one meter.

LOS ball range m
Density of tier1 m
Density of tier2 m
Bandwidth MHz; GHz
Path loss law ; ,
Number of antennas ;
Carrier frequency for two tiers GHz; GHz
Transmit Power at BSs dBm; dBm
Transmit Power at the typical user dBm
Backhaul capacity Mbps
Number of content , ,
Skew of the popularity distribution
TABLE I: General Settings of the Network

Fig. 2 illustrates the impact of and on success probability. Comparing the analytical results of the success probability with the simulation results, we note that they match each other ideally, thereby certifying the analysis. It is obvious that our cache-enabled HetNets outperform the traditional HetNets where the macro BSs have no caching capacity, especially in the area . Moreover, the SIR scenario fits the analytical results perfectly so that our model can be regarded as an interference-limited system. Lastly, with the rise of antenna scales , the success probability will correspondingly increase.


Fig. 2: Success probability versus pre-decided rate threshold with  MHz.

Fig. 3: Success probability versus cache capacity in tier 2 with different , and bps.

Fig. 4: Area spectral efficiency versus pre-decided rate threshold with different , carrier frequencies and MHz.

Fig. 3 shows the success probability versus cache capacity with different . The success probability has a negative correlation with the transmit power . Massive contributes to large received power, when increases to 100 dBm, the success probability becomes independent on . In this case, the typical use will associate with tier 1 all the time. In terms of the cache capacity, the success probability is a monotonic increasing function with pico BSs’ storage capacity . For the cache capacity of macro BSs , when , the success probability has no relationship with since the less-popular content that only contained in the server can be transmitted freely through the backhaul. Under this condition, the proposed HetNet is same with the traditional one. On the other side, when , the server will be blocked so that the success probability will be benefited by the large , which represents that more multimedia files can be stored at macro BSs.

Fig 4 plots the ASE versus with different and carrier frequencies. The optimum pre-decided rate threshold for achieving the maximum ASE can be easily figured out due to the convex property. When the density of pico tier increases from m to m, the optimal number decreases. Then, we compare the carrier frequencies at 28 GHz, 38 GHz, 60 GHz and 73GHz with antenna scales 10, 20, 40 and 80 [14], respectively. All of them has the same LOS path loss exponents (), except that the path loss exponent for 60 GHz is 2.25 [19, 20]. Fig. 4 illustrates that 73 GHz is the best choice thanks to the largest antenna scales. The simulation of 60 GHz matches the analytical results with an insignificant difference, thereby verifying all analytical expressions under the condition .

Vi Conclusion

In this treatise, we utilize the stochastic geometry to analyze the performance of our cache-enabled hybrid HetNet. More specifically, the proposed network, which performs better than the traditional HetNet, can be regarded as an interference-limited system due to the high density of mmWave tier and the nature of sub-6 GHz tier. As discussed in the previous sections, our system has a positive correlation with pico BSs’ antenna scales and the cache capacity of both tiers. Additionally, there exists an optimum value of pre-decided rate threshold contributing to the maximum ASE. Lastly, for Max-RP association strategy, 73 GHz is the best carrier frequency of mmWave tier.

Appendix A: Proof of Lemma 1

The Laplace transform of interference in tier 2 is given by

(A.1)

where (a) follows the Gamma random variable’s moment generating function 

[12]; (b) is computing the expectation of tier 2 antenna gain . When , (A.1) can be simplified as follows

(A.2)

(c) follows Gauss hypergeometric function. When , (A.1) can be simplified as below

(A.3)

(d) follows (2.117-1), (2.117-3) and (2.118-1) in [23]. With the aid of Gauss-Chebyshev Quadrature, we obtain Lemma 1.

References

  • [1] D. Liu and C. Yang, “Caching policy toward maximal success probability and area spectral efficiency of cache-enabled HetNets,” IEEE Trans. Commun., vol. 65, no. 6, pp. 2699–2714, Jun. 2017.
  • [2] J. Zhao, Y. Liu, K. K. Chai, A. Nallanathan, Y. Chen, and Z. Han, “Spectrum allocation and power control for non-orthogonal multiple access in HetNets,” IEEE Trans. Wireless Commun., vol. 16, no. 9, pp. 5825–5837, Sep. 2017.
  • [3] Y. Liu, Z. Qin, M. Elkashlan, A. Nallanathan, and J. A. McCann, “Non-orthogonal multiple access in large-scale heterogeneous networks,” IEEE J. Sel. Areas Commun., vol. 35, no. 12, pp. 2667–2680, Dec. 2017.
  • [4] Z. Ding, Y. Liu, J. Choi, Q. Sun, M. Elkashlan, C. L. I, and H. V. Poor, “Application of non-orthogonal multiple access in LTE and 5G networks,” IEEE Commun. Mag., vol. 55, no. 2, pp. 185–191, Feb. 2017.
  • [5] X. Wang, M. Chen, T. Taleb, A. Ksentini, and V. C. M. Leung, “Cache in the air: Exploiting content caching and delivery techniques for 5G systems,” IEEE Commun. Mag., vol. 52, no. 2, pp. 131–139, Feb. 2014.
  • [6] D. Liu and C. Yang, “Energy efficiency of downlink networks with caching at base stations,” IEEE J. Sel. Areas Commun., vol. 34, no. 4, pp. 907–922, Apr. 2016.
  • [7] R. W. Heath, M. Kountouris, and T. Bai, “Modeling heterogeneous network interference using Poisson point processes,” IEEE Trans. Signal Process., vol. 61, no. 16, pp. 4114–4126, Aug. 2013.
  • [8] Y. Liu, Z. Qin, M. Elkashlan, Z. Ding, A. Nallanathan, and L. Hanzo, “Nonorthogonal multiple access for 5G and beyond,” Proc. IEEE, vol. 105, no. 12, pp. 2347–2381, Dec. 2017.
  • [9] T. S. Rappaport, F. Gutierrez, E. Ben-Dor, J. N. Murdock, Y. Qiao, and J. I. Tamir, “Broadband millimeter-wave propagation measurements and models using adaptive-beam antennas for outdoor urban cellular communications,” IEEE Trans. Antennas Propag., vol. 61, no. 4, pp. 1850–1859, Apr. 2013.
  • [10] A. V. Alejos, M. G. Sanchez, and I. Cuinas, “Measurement and analysis of propagation mechanisms at 40 GHz: Viability of site shielding forced by obstacles,” IEEE Trans. Veh. Technol., vol. 57, no. 6, pp. 3369–3380, Nov. 2008.
  • [11] H. Elshaer, M. N. Kulkarni, F. Boccardi, J. G. Andrews, and M. Dohler, “Downlink and uplink cell association with traditional macrocells and millimeter wave small cells,” IEEE Trans. Wireless Commun., vol. 15, no. 9, pp. 6244–6258, Sep. 2016.
  • [12] T. Bai and R. W. Heath, “Coverage and rate analysis for millimeter-wave cellular networks,” IEEE Trans. Wireless Commun., vol. 14, no. 2, pp. 1100–1114, Feb. 2015.
  • [13] D. Maamari, N. Devroye, and D. Tuninetti, “Coverage in mmWave cellular networks with base station co-operation,” IEEE Trans. Wireless Commun., vol. 15, no. 4, pp. 2981–2994, Apr. 2016.
  • [14] W. Yi, Y. Liu, and A. Nallanathan, “Modeling and analysis of D2D millimeter-wave networks with poisson cluster processes,” IEEE Trans. Commun., vol. 65, no. 12, pp. 5574–5588, Dec. 2017.
  • [15] J. G. Andrews, T. Bai, M. N. Kulkarni, A. Alkhateeb, A. K. Gupta, and R. W. Heath, “Modeling and analyzing millimeter wave cellular systems,” IEEE Trans. Commun., vol. 65, no. 1, pp. 403–430, Jan. 2017.
  • [16] X. Yu, J. Zhang, M. Haenggi, and K. B. Letaief, “Coverage analysis for millimeter wave networks: The impact of directional antenna arrays,” IEEE J. Sel. Areas Commun., vol. 35, no. 7, pp. 1498–1512, Jul. 2017.
  • [17] Y. Liu, Z. Qin, M. Elkashlan, Y. Gao, and L. Hanzo, “Enhancing the physical layer security of non-orthogonal multiple access in large-scale networks,” IEEE Trans. Wireless Commun., vol. 16, no. 3, pp. 1656–1672, Mar. 2017.
  • [18] Y. Liu, Z. Ding, M. Elkashlan, and H. V. Poor, “Cooperative non-orthogonal multiple access with simultaneous wireless information and power transfer,” IEEE J. Sel. Areas Commun., vol. 34, no. 4, pp. 938–953, Apr. 2016.
  • [19] S. Deng, M. K. Samimi, and T. S. Rappaport, “28 GHz and 73 GHz millimeter-wave indoor propagation measurements and path loss models,” in Prof. IEEE Int. Conf.on Commun. Workshop (ICCW), Jun. 2015, pp. 1244–1250.
  • [20] T. S. Rappaport, E. Ben-Dor, J. N. Murdock, and Y. Qiao, “38 GHz and 60 GHz angle-dependent propagation for cellular & peer-to-peer wireless communications,” in IEEE Proc. of International Commun. Conf. (ICC), Jun. 2012, pp. 4568–4573.
  • [21] H. S. Jo, Y. J. Sang, P. Xia, and J. G. Andrews, “Heterogeneous cellular networks with flexible cell association: A comprehensive downlink SINR analysis,” IEEE Trans. Wireless Commun., vol. 11, no. 10, pp. 3484–3495, Oct. 2012.
  • [22] M. Afshang, H. S. Dhillon, and P. H. J. Chong, “Modeling and performance analysis of clustered device-to-device networks,” IEEE Trans. Wireless Commun., vol. 15, no. 7, pp. 4957–4972, Jul. 2016.
  • [23] A. Jeffrey and D. Zwillinger, Table of integrals, series, and products.   Academic press, 2007.