Novel Spectrum Allocation Among Multiple Transmission Windows for Terahertz Communication Systems
Abstract
This paper presents a novel spectrum allocation strategy for multiuser terahertz (THz) band communication systems when the to-be-allocated spectrum is composed of multiple transmission windows (TWs). This strategy explores the benefits of (i) allowing users to occupy sub-bands with unequal bandwidths and (ii) optimally avoiding using some spectra that exist at the edges of TWs where molecular absorption loss is high. To maximize the aggregated multiuser data rate, we formulate an optimization problem, with the primary focus on spectrum allocation. We then apply transformations and modifications to make the problem computationally tractable, and develop an iterative algorithm based on successive convex approximation to determine the optimal sub-band bandwidth and the unused spectra at the edges of TWs. Using numerical results, we show that a significantly higher data rate can be achieved by changing the sub-band bandwidth, as compared to equal sub-band bandwidth. We also show that a further data rate gain can be obtained by optimally determining the unused spectra at the edges of TWs, as compared to avoiding using pre-defined spectra at the edges of TWs.
Index Terms:
Terahertz communication, spectrum allocation, adaptive bandwidth, transmission windows.I Introduction
The terahertz (THz) band communication (THzCom), operating in the electromagnetic spectrum ranging from 0.1 to 10 THz, is envisioned as an enabling technology to support emerging wireless applications that require an explosive amount of data in the sixth-generation (6G) and beyond systems. Specifically, the enormous available bandwidths in the order of a few tens to a hundred gigahertz (GHz) and the extremely short wavelengths offer enormous potentials for wireless transmission [1]. Despite the promise, THzCom encounters unique challenges that have never been addressed at lower frequencies, e.g., higher channel sparsity, severe spreading loss, and frequency-dependent molecular absorption loss [2]. These challenges need to be wisely tackled for developing ready-to-use THzCom systems.
To fully harness the potentials of the THz band, the huge available bandwidth that spans across multiple ultra-wide transmission windows (TWs) needs to be optimally useds [3, 4]. This necessitates the exploration of novel and efficient spectrum allocation schemes for THzCom system when the to-be-allocated spectrum of interest is composed of multiple TWs. In this work, we consider the entire spectrum regions that exist between neighboring molecular absorption coefficient peaks at the THz band as TWs [4]. Multi-band-based spectrum allocation scheme has recently been studied for micro- and macro-scale THzCom systems with multiple TWs in the to-be-allocated spectrum of interest. In this scheme, the spectrum of interest is divided into a set of non-overlapping sub-bands which have a relatively small bandwidth, and then the sub-bands are utilized to serve the users in the system. This scheme has the ability to efficiently allocate spectral resources while ensuring intra-band interference-free transmission. Triggered by this ability, multi-band-based spectrum allocation has been widely explored for micro- and macro-scale multiuser THzCom systems [4, 5, 6, 7, 8].
A novel sub-band assignment strategy for multi-band-based spectrum allocation was proposed in [4] to improve the throughput fairness among users in multiuser THzCom systems. Moreover, spectrum allocation problems in the THz band backhaul network and non-orthogonal multiple access (NOMA) assisted THzCom systems were addressed in [5], and [6], respectively. Furthermore, efficient transmit power allocation and sub-band assignment algorithms were developed in [7] and [8], respectively, for multiuser THzCom systems. It is noted that the studies [4, 5, 6, 7, 8] focused on multi-band-based spectrum allocation with equal sub-band bandwidth (ESB), where the spectrum of interest is divided into sub-bands with equal bandwidth. However, since the molecular absorption loss varies considerably within the THz band, the variation among the sub-bands would be very high when ESB is considered. This variation can be reduced by exploring spectrum allocation with adaptive sub-band bandwidth (ASB), which allows to change the sub-band bandwidth. This can eventually lead to an overall improvement in the data rate performance.
Motivated by the potentials of ASB, the multi-band-based spectrum allocation with ASB was first introduced in our previous study [9]. However, the design in [9] is only applicable when the spectrum of interest exists entirely within a specific region of a TW. When the to-be-allocated spectrum of interest spans across multiple TWs, on one hand, the number of sub-bands that exist within each TW needs to be optimally decided. On the other hand, the set of users that are served by each TW and the sub-band assignment policy need to be identified. These challenges have not been touched by [9], but need to be carefully addressed to exploit the potentials of the huge available bandwidth that spans across multiple TWs at the THz band. This motivates this work.
We note that since molecular absorption loss is high at the edges of TWs, it may be beneficial if some spectra at the edges of TWs are not used during spectrum allocation. Despite so, previous studies including ours [9, 4, 5, 6, 7, 8] considered that the entire spectrum existing within the to-be-allocated spectrum of interest is fully utilized. Several recent studies [10, 11] considered to avoid some spectra at the edges of TWs. However, the designs in [10, 11] considered the unused spectra at the edges of TWs to be pre-defined, i.e., the unused spectra are arbitrarily chosen and fixed, which may not guarantee the optimal performance.
In this work, we design a novel multi-band-based spectrum allocation strategy when the to-be-allocated spectrum of interest is composed of multiple TWs. Specifically, we explore the benefits of ASB and optimally avoiding some spectra that exist at the edges of TWs during spectrum allocation. With these considerations in mind and focusing on an indoor THzCom scenario where a single access point (AP) supports the downlink of multiple users, we formulate a spectrum allocation problem. To analytically solve this problem, we propose transformations, as well as modifications, and arrive at an approximate convex problem. Subsequently, we develop an iterative algorithm based on the successive convex approximation technique to solve the approximate convex problem. Using numerical results, we show that when the to-be-allocated spectrum of interest is composed of multiple TWs, enabling ASB achieves a significantly higher data rate (between and ) as compared to adopting ESB. We also show that an additional gain of up to can be obtained by optimally determining the unused spectra at the edges of TWs, as compared to avoiding using pre-defined spectra that exist at the edges of TWs, especially when the power budget constraint is stringent. We finally illustrate that the feasibility region of the spectrum allocation problem can be improved by considering ASB and optimally avoiding some spectra that exist at the edges of TWs.
II System Model and Problem Formulation
We consider a three-dimensional (3D) indoor THzCom system where a single AP, which is located at the center of the ceiling of the indoor environment, supports the downlink of users which demand high data rates. We consider that the users are stationary and distributed uniformly on the floor in this environment. We denote as the set for the users and further denote and as the horizontal and 3D distances of the link between the AP and the th user, respectively, where is the fixed difference between the heights of AP and the users [12].
II-A THz Spectrum
The intermittent molecular absorption loss peaks that are observed throughout the THz band at different frequencies divides the entire THz band into multiple ultra-wideband transmission windows (TWs) [9]. In this work, the entire spectrum regions that exist between neighboring molecular absorption coefficient peaks are denoted by TWs [4]. We focus in this work on the scenario where the to-be-allocated spectrum of interest is composed of multiple TWs, as shown in Fig. 1. For analytical purposes, we consider that each TW is composed of multiple positive absorption coefficient slope regions (PACSRs) and negative absorption coefficient slope regions (NACSRs). Here, PACSRs and NACSRs are defined as the regions within TWs with increasing and decreasing molecular absorption coefficients, respectively. Considering this, we find that the to-be-allocated spectrum of interest is made of multiple PACSRs and NACSRs, as shown in Fig. 1. We denote and as the sets of the PACSRs and the NACSRs in the spectrum of interest, respectively. We further denote as the set of all regions in the spectrum of interest.
We focus on multi-band-based spectrum allocation. Thus, each region within the spectrum of interest is divided into sub-bands that are separated by guard bands. An illustration of the arrangement of sub-bands within a PACSR and an NACSR is depicted in Fig. 2. We denote as the set of the sub-bands that exist in the region and further denote and as the bandwidth and the center frequency of the th sub-band that exist in the region , respectively. For notational convenience, we consider that sub-bands are labeled such that when and when . Thus, can be expressed as
(1) |
where , is the end-frequency of the region when and the start-frequency of the region when , as shown in Fig. 2. In addition, when when , is the bandwidth of the unused spectrum that exists in the region , and is the bandwidth of the guard band that exists between the th and the th the sub-band in the region . Moreover, , , and , . In this work, we optimally determine the value of by solving the optimization problem that will be formulated in the subsequent section. Also, we consider that the guard band bandwidth is fixed and is equal to , i.e., 111 To further improve the performance of the spectrum allocation strategy, it may be beneficial to optimally design the values of , which will be considered in our future work..
Remark 1
As shown in Figs. 1 and 2, the values of the molecular absorption coefficient are very high at the edges of TWs. Due to this, it may be beneficial if some spectrum that exists at the edges of TWs are not utilized during spectrum allocation. To optimally determine the unused spectra at the edges of TW, we introduce a design variable in our spectrum allocation strategy, where denotes the unused spectra at the edge of each region where the molecular absorption coefficient is very high. It is reflect in (II-A) and depicted in Fig. 2.
With all the aforementioned considerations, we express the total bandwidth available within the region as
(2) |
We assume that the users in the system are served by separate sub-bands [4, 9]. This assumption ensures intra-band interference-free data transmission and eliminate the hardware complexity and signal processing overhead caused by frequency reuse in the system222A higher spectral efficiency can be achieved by allowing (i) users to occupy multiple sub-bands and/or (ii) sub-band reuse during spectrum allocation, which will be considered in our future work.. Under this assumption, we set the total number of sub-bands in the spectrum of interest to be equal to the number of users in the system, i.e.,
(3) |
Although the total number of sub-bands in the spectrum of interest, i.e., , is defined by (3), the number of sub-bands that exist within each region, i.e., , is unknown. Thus, it should be carefully decided such that the spectra available within the regions are utilized optimally, while ensuring that the constraint (3) is satisfied.
Let us denote as the sub-band assignment indicator variable such that if the th sub-band in the region is assigned to the th user and , otherwise, . Under the assumption that the users in the system are served by separate sub-bands, we have
(4a) | ||||
(4b) |
where (4a) ensures that each user in the system is assigned with one sub-band and (4b) ensures that each sub-band in the spectrum of interest is assigned to one user.
II-B Achievable Data Rate
Considering the spreading loss and the frequency-selective molecular absorption loss at the THz band, the data rate achieved by the th user in the th sub-band in region is obtained as an integral of capacity on frequency, given by [4]
(5) |
where , is the capacity when the frequency is , is the transmit power allocated to the th user, and are the antenna gains at the AP and users, respectively, is the speed of light, and is the molecular absorption coefficient at 333The impact of inter-band interference (IBI) is not considered in this work due the fact that there exist IBI suppression schemes that can suppress IBI with minimal throughput degradation and waveform designs that can minimize power leakages to adjacent bands [9, 10].. Following [13, 9], the noise power is calculated as a function of the sub-band bandwidth. Specifically, noise power of the th sub-band in region , , is obtained as , where is the noise power density [13, 9]. We assume that users in the system are equipped with tunable front-ends such that they can dynamically adjust their operating frequency and bandwidth to match those of their assigned sub-bands.
We note that although the values of can be found by using the HITRAN database [14], line by line radiative transfer model (LBLRTM) [15], or the International Telecommunication Union (ITU) model [16], there does not exist a tractable expression that maps to within the entire THz band and can be used to analytically model spectrum allocation problems with ASB [9]. With this in mind, through curve fitting, we model the values of corresponding to each region using separate exponential functions of . Specifically, we approximate in the region as
(6) |
where , and are the model parameters obtained for the region . It is noted that and . Thereafter, by substituting (6) into (5), is obtained as
(7) |
where and is given in (II-A). When behaves according to (6), within the PACSR (or the NACSR) decreases (or increases) monotonically with frequency. We approximate this monotonically decreasing (or increasing) behaviour of using a piecewise linear function. This approximation allows to modify in (7) into a simplified expression, given by
(8) |
Finally, the achievable data rate of the th user is obtained as
(9) |
III Optimal Spectrum Allocation
In this section, we present the multi-band-based spectrum allocation problem with ASB when the to-be-allocated spectrum of interest is composed of multiple TWs. We aim to maximize the sum data rate among all users under given sub-band bandwidth, power, and rate constraints. Mathematically, this problem is formulated as
(10a) | |||
(10b) | |||
(10c) | |||
(10d) | |||
(10e) | |||
(10f) | |||
(10g) |
The problem consists of sub-band assignment, sub-band bandwidth allocation, determination of the unused spectra at the edges of TW, and power control. The constraint (10b) reflects the power budget at the AP, (10c) ensures that the power allocated to each user is upper bounded by , (10d) ensures that the data rate achieved by each user is greater than the threshold , and (10e) ensures that the bandwidth of each sub-band is upper bounded by .
III-A Solution Approach
We now present the solution to the optimization problem . We note that it is extremely difficult to use traditional optimization theory techniques to analytically solve in its current form. This is because the number of sub-bands that exist within each region, i.e., , is unknown. This leads to the fact that the number of design variables and are unknown. Therefore, we need to render the problem computationally tractable. Considering this, we introduce modifications to the problem , as follows:
We first consider that there exist sub-bands within each region, i.e., . This implies that there would be sub-bands in the entire spectrum of interest. Out of these sub-bands, we ensure that only some of the sub-bands are used by at most one user, while the remaining sub-bands are not assigned to any user. To reflect this modification, we replace (3) and (4b) in with
(11a) | |||
(11b) |
After this modification, although the intractability of caused by the unknown number of design variables corresponding to each region can be solved, a new challenge arises due to the existence of guard bands in each region. Specifically, our spectrum allocation strategy considers that guard bands with bandwidth of exist between sub-bands. Thus, guard bands appear in each region as each region consists sub-bands after the modification. This would result in guard bands existing even between sub-bands which have zero bandwidth, which must to be avoided. To overcome this challenge, we first consider that the bandwidths of guard bands, , are variables such that . Then, we ensure that there exists exactly one guard band with the bandwidth of between sub-bands that have non-zero bandwidth, and the remaining guard bands have zero bandwidth. This is achieved by ensuring that the guard band between the th and the th sub-band in a given region has non-zero bandwidth only if the th sub-band of this region has non-zero bandwidth. Mathematically, this is expressed as
(12) |
We then note that when . Considering this, (12) can be equivalently expressed as
(13) |
With these aforementioned modifications, we transform into the following tractable optimization problem, given by
(14) | ||||
subject to |
Next, we note that although is tractable, it is non-convex in its current form. In particular, the non-convexity in arises due to (i) in (9) is non-concave w.r.t. the design variables , ,, , and (ii) the design variable is binary. Considering this, first, to handle the non-concavity of , we introduce the following substitutions for , , , given by
(15a) | ||||
(15b) | ||||
(15c) |
where are real constants with and , are the new design variables that replace , , in the optimization problem, with and . Considering these replacements, we transform into the following equivalent problem, given by
(16a) | |||
(16b) | |||
(16c) |
where , is a very small positive number, , , and . In , is given by
(17) |
with . We also note that, in , is selected such that , , , to ensure the convexity of the objective function in and the constraint function (10d) w.r.t. , since in (17) is concave w.r.t. only when [9]. Here, , , and is the maximum of the link distances, i.e., . We clarify that for any selected combination of real values for , the feasible region of the term is , as long as . Thus, the constraint introduced on does not add any restriction on the domain of the problem .
We now handle the non-convexity of caused by the binary constraint function (10d). To this end, we first relax the binary variables in into real variables [17, 18, 9]. Let us denote as the real variables that would replace in the optimization problem. We then rewrite the binary constraint function (10g) equivalently as
(18a) | |||
(18b) |
Considering that the constraint (18b) is non-convex w.r.t , we relax the constraint (18b) and include it as a penalty function in the objective function [17, 18, 9]. In doing so, we transform as
(19) | ||||
where is the objective function of , is the constant penalty factor, , and is an approximate of . Moreover, , , , and are the constraints obtained by replacing with in , , , and (16b), respectively. We note that is obtained as a standard convex optimization problem. Thus, by using the successive convex approximation technique, we propose an iterative algorithm to solve , which is summarized in Algorithm 1.
We clarify that in Step 3 of Algorithm 1, is optimized for the given . This yields [19, 20]. Thus, Algorithm 1 produces a monotone sequence of improved feasible solutions, , to the problem once is initialized from a feasible point, e.g., . Moreover, is bounded by constraints (10d) and (16c). Based on these, it can be concluded that the convergence of Algorithm 1 is guaranteed, i.e., as [19, 20].
By checking the operations of Algorithm 1, we find that the same steps are executed in each iteration; thus, we focus only on the complexity of one iteration. The dominating computational complexity of the algorithm arises from Step 3, where is solved to obtain and . Thus, solving each iteration of Algorithm 1 requires a complexity of [18, 9]. We find from [9, 4] that solving the algorithm for spectrum allocation with ESB requires a complexity of approximately . Based on this and the fact that when , we conclude that the proposed algorithm for spectrum allocation with ASB is times more complex than that for spectrum allocation with ESB.
IV Numerical Results
In this section, we present numerical results to examine the performance of the proposed multi-band-based spectrum allocation scheme. Specifically, we consider an indoor environment of size , where an AP at the center of the indoor environment serves 30 users. We also consider that the spectrum that exists between and is used to serve the users. We clarify that the considered spectrum spans across two TWs and each of the TW is composed of a PACSR and an NACSR. We denote the two TWs and the four regions by and , respectively, as depicted in Fig. 1. Unless otherwise stated, the values of the other parameters used for numerical results are as follows: , , , , , , , , , , and .
In Fig. 3, we plot the performance gain achieved by (i) considering ASB when the to-be-allocated spectrum of interest is composed of multiple TWs, i.e., , and (ii) optimally determining the unused spectra that exists at the edge of TWs, i.e., . In Fig. 3(a), we first observe the aggregated multiuser data rate, , achieved by the proposed scheme with ASB while considering the state-of-the-art ESB-based spectrum allocation scheme in [9] as the benchmark. We clarify that in the benchmark scheme, (i) ESB is considered, (ii) the sub-bands are optimally assigned to users while satisfying the constraints (4a), (4b), and (10g), and (iii) the transmit power is allowed to vary while satisfying the constraints (10b) and (10c). We clarify that for the two schemes, is fixed at , where denotes the bandwidth of the spectrum at the edges of each region where the molecular absorption coefficient is higher than 0.3 [10, 9]. We observe that for all power budget values, a significantly higher (between and ) can be achieved by adopting ASB by varying the sub-band bandwidths. Next, we compare achieved by the proposed spectrum allocation scheme with ASB when and , and then plot the data rate gain of considering relative to considering , . We observe that although is marginal when the power budget is high, a noticeable data rate gain () exists when the power budget constraint is stringent. These two observations show the significance of considering and . We further observe in Fig. 3(b) that the average value of the unused bandwidth within regions, , increases when the power budget decreases, which explains the increasing when the power budget decreases.
In Fig. 4, we show how the spectral resources are allocated according to the proposed multi-band-based allocation scheme. We plot the allocation of spectral resources within the transmission window, . We note that the gray bars represent the sub-bands and the width and the height of the bars represent the bandwidth of the sub-bands and the transmission distance of the users assigned to the sub-bands, respectively. We first observe that, according to the optimal sub-band assignment configuration, the users with longer distances are assigned to center sub-bands where the absorption loss is low while the users with shorter distances are assigned to edge sub-bands where the absorption loss is high. This sub-band assignment configuration is similar to that considered in [4, 5] for the scenario with ESB. This shows that, since the absorption loss increases exponentially with distance, the data rate gain achieved by the users with longer distances, through occupying the center sub-bands as compared to occupying the edge sub-bands, is higher than the data rate loss incurred by the users with shorter distances, through occupying the edge sub-bands as compared to occupying the center sub-bands. We also observe that although the edge sub-bands experience higher absorption losses, a wider bandwidth is allocated to the edge sub-bands compared to the center sub-bands, since the edge sub-bands are occupied by the users with shorter distances.
In order to examine the impact of different spectrum of interest on the proposed spectrum allocation scheme, we plot versus the number of regions, , in Fig. 5. Here, the regions corresponding to of 1, 2, 3, and 4 are , , , and , respectively. For the sake of fairness, we made the number of users in the system to be proportional to the total available bandwidth within the spectrum of interest. More precisely, the number of users corresponding to of 1, 2, 3, and 4 are 10, 17, 27, and 30, respectively. We first observe that achieved by our proposed scheme with ASB is equal to that achieved by the scheme in [9], which is only applicable when the spectrum of interest is entirely within a region, i.e., . This shows that the scheme investigated in [9] is a special case of our proposed multi-band based spectrum allocation scheme in this work. We also observe that the proposed spectrum allocation scheme with ASB outperforms the state-of-the-art spectrum allocation scheme with ESB by for all values of . Furthermore, we observe that an additional gain can be achieved by considering as compared to . Finally, we observe that for the proposed spectrum allocation scheme improves when the upper bound on sub-band bandwidth, , increases. This is expected since a larger gives more room to exploit the ASB capability, which improves .

Finally, we illustrate the improvement in the feasibility region of the problem that occurs as a result of considering ASB and . To this end, Fig. 6 plots the percentage of simulation trails for which the problem is feasible, , for the benchmark scheme with ESB and the spectrum allocation scheme with ASB when and . We first observe that increases when power budget, , and rate threshold, , increases and decreases, respectively. This is expected because the increase in and the decrease in will lead to the constraints (10b) and (10d) in to be less stringent, thereby leading to more opportunities to attain a feasible solution to . Second, we observe that our consideration of ASB increases as compared to considering ESB. Third, we observe that our consideration of increases as compared to considering . The second and third observations are due to the fact that considering ASB and allows sub-band bandwidth to vary and some spectra that exist at the edges of TWs to be avoided from utilization during spectrum allocation, respectively, if that can yield a feasible solution point. This again shows the importance of our considerations of ASB and during spectrum allocation.
V Conclusion
This paper investigated the benefits of ASB and optimally determining the unused spectra at the edges of TWs during multi-band-based spectrum allocation in THzCom systems where the to-be-allocated spectrum of interest is composed of multiple TWs. We first formulated an optimization problem, with the main focus on spectrum allocation. Thereafter, to analytically solve the formulated problem, we proposed transformations, as well as modifications, and developed iterative algorithms based on the SCA technique. Using numerical results, we showed that the proposed spectrum allocation scheme with ASB achieves a significantly higher data rate compared to the multi-band based spectrum allocation with ESB. We also showed that an additional gain can be obtained by optimally determining the unused spectra at the edges of TWs as compared to avoiding using pre-defined spectra that exist at the edges of TWs, especially when the power budget constraint is stringent.
We clarify that the solution proposed in this work is only applicable when the molecular absorption coefficient within the TWs can be modeled using exponential functions of frequency with minimal approximation errors. Nevertheless, to harness the potentials of the huge available bandwidths of THz band, it would be beneficial to develop a solution for the spectrum allocation with ASB when the molecular absorption coefficient varies in a highly non-linear manner within the TWs and cannot be modeled using exponential functions of frequencies, which is the case in many spectrum regions within the THz band.
References
- [1] Z. Chen, C. Han, Y. Wu, L. Li, C. Huang, Z. Zhang, G. Wang, and W. Tong, “Terahertz wireless communications for 2030 and beyond: A cutting-edge frontier,” IEEE Comm. Mag., vol. 59, no. 11, pp. 66–72, Nov. 2021.
- [2] M. Polese, J. M. Jornet, T. Melodia, and M. Zorzi, “Toward end-to-end, full-stack 6G terahertz networks,” IEEE Commun. Mag., vol. 58, no. 11, pp. 48–54, Nov. 2020.
- [3] I. F. Akyildiz, A. Kak, and S. Nie, “6G and beyond: The future of wireless communications systems,” IEEE Access, vol. 8, pp. 133 995–134 030, July 2020.
- [4] C. Han and I. F. Akyildiz, “Distance-aware bandwidth-adaptive resource allocation for wireless systems in the terahertz band,” IEEE Trans. THz Sci. Technol., vol. 6, no. 4, pp. 541–553, July 2016.
- [5] M. Yu, A. Tang, X. Wang, and C. Han, “Joint scheduling and power allocation for 6G terahertz mesh networks,” in Proc. Int. Conf. Comput. Netw. Commun. (ICNC), Big Island, HI, USA, Feb. 2020, pp. 631–635.
- [6] X. Zhang, C. Han, and X. Wang, “Joint beamforming-power-bandwidth allocation in terahertz NOMA networks,” in Proc. Int. Conf. Sensing, Commun., Netw. (SECON), Boston, MA, USA, Sept. 2019, pp. 1–9.
- [7] H. Zhang, H. Zhang, W. Liu, K. Long, J. Dong, and V. C. M. Leung, “Energy efficient user clustering and hybrid precoding for terahertz MIMO-NOMA systems,” in Proc. IEEE Int. Conf. Commun. (ICC), Dublin, Ireland, June 2020, pp. 1–5.
- [8] H. Zhang, Y. Duan, K. Long, and V. C. M. Leung, “Energy efficient resource allocation in terahertz downlink NOMA systems,” IEEE Trans. Commun., vol. 69, no. 2, pp. 1375–1384, Feb. 2021.
- [9] A. Shafie, N. Yang, S. A. Alvi, C. Han, S. Durrani, and J. M. Jornet, “Spectrum allocation with adaptive sub-band bandwidth for terahertz communication systems,” IEEE Trans. Commun., vol. 70, no. 2, pp. 1407–1422, Jan. 2022.
- [10] H. Yuan, N. Yang, K. Yang, C. Han, and J. An, “Enabling massive connections using hybrid beamforming in terahertz micro-scale networks,” in Proc. IEEE Wireless Commun. Netw. Conf. (WCNC), Seoul, South Korea, June 2020, pp. 1–7.
- [11] Z. Hossain and J. M. Jornet, “Hierarchical bandwidth modulation for ultra-broadband terahertz communications,” in Proc. IEEE Int. Conf. Commun. (ICC), Shanghai, China, May 2019, pp. 1–7.
- [12] A. Shafie, N. Yang, S. Durrani, X. Zhou, C. Han, and M. Juntti, “Coverage analysis for 3D terahertz communication systems,” IEEE J. Sel. Areas Commun., vol. 39, no. 6, pp. 1817–1832, June 2021.
- [13] A. Saeed, O. Gurbuz, A. O. Bicen, and M. A. Akkas, “Variable-bandwidth model and capacity analysis for aerial communications in the terahertz band,” IEEE J. Sel. Areas Commun., vol. 39, no. 6, pp. 1768–1784, June 2021.
- [14] L. Rothman, I. Gordon, A. Barbe, D. Benner, P. Bernath, M. Birk, V. Boudon, and L. Brown, “The HITRAN 2008 molecular spectroscopic database,” J. Quantitative Spectroscopy Radiative Transfer, vol. 110, no. 9, pp. 533–572, July 2009.
- [15] A. Saeed, O. Gurbuz, and M. A. Akkas, “Terahertz communications at various atmospheric altitudes,” Physical Commun., vol. 41, pp. 101–113, Aug. 2020.
- [16] ITU-R, “Attenuation by atmospheric gases and related effects,” Recommendation ITU-R P.676-12, Aug. 2019.
- [17] D. W. K. Ng and R. Schober, “Secure and green SWIPT in distributed antenna networks with limited backhaul capacity,” IEEE Trans. Wireless Commun., vol. 14, no. 9, pp. 5082–5097, Sept. 2015.
- [18] E. Che, H. D. Tuan, and H. H. Nguyen, “Joint optimization of cooperative beamforming and relay assignment in multi-user wireless relay networks,” IEEE Trans. Wireless Commun., vol. 13, no. 10, pp. 5481–5495, Oct. 2014.
- [19] T. T. Vu, D. T. Ngo, M. N. Dao, S. Durrani, and R. H. Middleton, “Spectral and energy efficiency maximization for content-centric C-RANs with edge caching,” IEEE Trans. Commun., vol. 66, no. 12, pp. 6628–6642, Dec. 2018.
- [20] B. R. Marks and G. P. Wright, “A general inner approximation algorithm for nonconvex mathematical programs,” Oper. Res., vol. 26, no. 4, pp. 681–683, Aug. 1978.