Performance analysis of multi user massive MIMO hybrid beamforming systems at millimeter wave frequency bands

Millimeter-wave (mmWave) and massive multi-input–multi-output (mMIMO) communications are the most key enabling technologies for next generation wireless networks to have large available spectrum and throughput. mMIMO is a promising technique for increasing the spectral efficiency of wireless networks, by deploying large antenna arrays at the base station (BS) and perform coherent transceiver processing. Implementation of mMIMO systems at mmWave frequencies resolve the issue of high path-loss by providing higher antenna gains. The motivation for this research work is that mmWave and mMIMO operations will be much more popular in 5G NR, considering the wide deployment of mMIMO in major frequency bands as per 3rd generation partnership project. In this paper, a downlink multi-user mMIMO (MU-mMIMO) hybrid beamforming communication system is designed with multiple independent data streams per user and accurate channel state information. It emphasizes the hybrid precoding at transmitter and combining at receiver of a mmWave MU-mMIMO hybrid beamforming system. Results of this research work give the tradeoff between multiple data streams per user and required number of BS antennas. It strongly recommends for higher number of parallel data streams per user in a mmWave MU-mMIMO systems to achieve higher order throughputs.


Introduction
The next generation wireless communication systems including 5G NR use mMIMO beamforming techniques to achieve higher SNRs and spatial multiplexing to enhance the data throughput at mmWave frequencies. The major challenge in using mmWave frequencies is propagation loss, therefore 5G NR networks are deployed with mMIMO (large scale antenna arrays) to overcome the path loss. The advantage at mmWave frequency bands is that we can accommodate more antenna elements in the given physical dimensions as wavelengths are very small. To reduce the system cost, antenna elements can be grouped into subarrays and each T-R module is dedicated to an antenna subarray using an emerging technique for mmWave communications called ''Hybrid Beamforming''. Hybrid beamforming designs are capable of transmitting data to multiple users using MU-MIMO with multiplexing gains and these MU-MIMO systems have high potential in mmWave communication networks. Hybrid transceivers consist of lesser RF chains compared to number of transmit antenna elements as they use analog beamformers in the RF domain, and digital beamformers in baseband domain. Hybrid beamforming design balances beamforming gains (to overcome path losses) and power consumption, hardware cost in mmWave mMIMO communication systems.
A common rule-of-thumb in mMIMO systems is that M/K [ 10 (where 'M' is number of antennas and 'K' is number of users) so that user channels are likely to be orthogonal and provides maximum efficiency. It is possible to reduce the computational complexity in hybrid beamforming for mmWave communication systems with lesser number of analog RF chains (compared to number of users) and their performance is close to that of optimal (or fully) digital beamformers [1]. As per 3GPP standard, 5G wireless technologies use frequencies in FR2 band as shown in Table 1 to have short range communication and high data rates [2]. MU-mMIMO wireless communication networks use SDMA to have multiplexing gains by serving multiple UEs with same T-R resources and gives substantial improvements in system throughput. MU-mMIMO system improves spectral efficiency as it allows BS Tx to communicate with many UE Rxs simultaneously using same T-R resources. However, the challenge in MU-mMIMO systems is designing transmit vectors by considering co-channel interference of other users. In UE dense scenarios, MU-MIMO dimensions are increased to fully exploit the spatial multiplexing capabilities. In these scenarios, it is challenging to distinguish UEs in spatial domain as number of pairing users are more. mMIMO increases spatial resolution with more number of narrow beams and gives high degree of freedom for MU paring. It allows number of BS antenna elements of order tens to hundreds, thereby also increasing the number of data streams in a cell to a large value.
MU-mMIMO hybrid beamforming structure shown in Fig. 1 divides beamforming among RF analog domain and baseband digital domain with cost, complexity and flexibility tradeoffs [4]. In RF analog domain, beamforming is achieved by applying phase shift to each antenna element in the antenna subarray. In baseband digital domain, beamforming is achieved using channel matrix to derive precoding and combining weights that helps to transmit and recover multiple data streams independently using a single channel.

Literature review
Hybrid beamforming designs are introduced to reduce the training overhead and hardware cost in mMIMO systems. Hybrid beamforming can be classified based on CSI (average or instantaneous), carrier frequency (mmWave) and complexity (reduced, full or switched complexity). Selection of algorithm to get the best tradeoff between these parameters depends on channel characteristics and application [5]. Acquiring precise and accurate CSI for MU-mMIMO systems is a challenge at mmWave frequencies as the number of BS antennas are high. Accurate channel estimation in MU-mMIMO can be achieved using jointiterative scheme based on step-length optimization [6]. Beamforming neural network (based on deep learning) for mmWave mMIMO systems can optimize the beamforming design and it is robust with higher spectral efficiency (compared to traditional beamforming algorithms) in the  [7]. For a given number of BS antenna elements, optimal number of UEs scheduled simultaneously gives maximum spectral efficiency in mMIMO systems. The spectral efficiency can be same for DL and UL (allows joint network optimization), also it is independent of instantaneous UE locations [8]. Multitask deep learning (MTDL)-based MU hybrid beamforming algorithm for mmWave mMIMO OFDM systems can give better results in terms of sum-rate and lower run-time compared to traditional algorithms [9]. mmWave UL MU-mMIMO system uses lens type antenna array at BS (two-dimensional, both elevation and azimuth angle) and uniform planar array at MS based on arrival/ departure angles of multi-path signals. ''Path delay compensation'' technique at BS transforms MU-MIMO frequency selective channels into parallel smaller MIMO frequency-flat channels at lower hardware cost and increased sum-rates [10]. Effective channel estimation scheme is proposed for time-varying DL channel of mmWave MU-MIMO systems based on angle of arrival/ departure and with minimum number of pilots [11]. Low complex single cell DL MU-mMIMO hybrid beamforming system with perfect CSI gives sum-rate that approaches ideal channel capacity [12]. For a given number of RF chains, performance gap between hybrid and digital beamforming can be reduced by minimizing the number of multiplexed symbols [13]. Clustering and feedback based hybrid beamforming for DL mmWave MU-MIMO NOMA systems give maximum sum-rates compared to OMA systems [14]. A blind MU detection algorithm based on ''markov random field'' to model clustering sparsity and estimate mMIMO channel performs better than the systems that do not exploit clustering sparsity of channel [15]. Manifold optimization, eigenvalue decomposition and OMP algorithms used in designing hybrid beamforming for broadband mmWave MIMO systems offer BER and spectral efficiency closer to fully digital beamforming designs [16]. When SNR or number of RF chains are increasing in mmWave MU-mMIMO system, the optimal hybrid precoding and combining schemes using OMP algorithm gives performance close to fully digital precoding in terms of total sum-rates [17]. Energy efficiency in MU-mMIMO systems is inversely proportional to number of RF chains. Usage of optimal RF and baseband precoding matrices can improve energy and cost efficiency by 76.59% compared to OMP algorithm [18]. Generalized block OMP algorithm for channel estimation in mmWave MU-MIMO system uses different strategies of constructing pilot signals/beamforming weights and this scheme outperforms the existing channel estimation algorithms including OMP [19]. ''Distributed compressive sensing'' method can decrease the feedback overhead, training of CSI estimation at UE and ''Joint OMP algorithm'' performs CSI recovery at BS of a MU-mMIMO systems [20]. Optimal hybrid beamforming scheme for MU-mMIMO relay systems with mixed and fully connected structures based on ''successive interference cancelation'' maximizes the sum-rates [21]. An efficient hybrid beamforming technique for relay assisted mmWave MU-mMIMO system based on ''Geometric Mean Decomposition Tomlinson Harashima Precoding'' algorithm can give performance closer to fully digital beamforming [13]. A channel estimation scheme called ''generalized-block compressed sampling matching pursuit'' for mmWave MU-MIMO systems over frequency selective fading channels can offer better performance than OMP algorithm [22]. Optimal hybrid beamforming design is proposed to minimize Tx power under SINR constraints in a MU-mMIMO system for the cases where number of UEs are less than and greater than the number of RF chains. It gives the optimal Tx powers close that of fully digital beamforming design [23]. Low complex ''hybrid regularized channel diagonalization'' scheme that combines analog beamforming and digital precoding for mmWave MU-mMIMO system performs better than conventional block diagonalization based hybrid beamforming designs even in the presence of lowresolution RF phase shifters [24]. ''Hybrid beamforming with selection'' scheme decreases the computational and hardware cost of small-to-medium bandwidth MU-mMIMO systems with moderate frequency selective channel [25]. Hybrid precoder and combiners for DL frequency selective channels are configured in a mmWave MU-MIMO system based on factorization, iterative hybrid design. BS simultaneously estimates all channels from UEs on each subcarrier using compressed sensing scheme to reduce the number of measurements [26]. Nonconvex hybrid precoding problem in mmWave MU-MIMO systems is addressed using ''penalty dual decomposition'' method under the assumption of perfect CSI. It uses lesser number of RF chains but still, the performance is close to that of fully digital beamforming [27]. ''Atomic Norm  27:1925-1939 1927 Minimization'' method is used for accurate, low-complex channel estimation and spectral efficiency in mmWave MU-MIMO systems which is based on continuous channel representation [13]. Spectral efficiency of a MU-mMIMO hybrid beamforming system can be improved by using low complex manifold optimization algorithm and its efficiency is closer to fully digital beamforming designs [13]. For accurate channel estimation with minimum training overhead and RF chains (compared to beamspace approach), MU-mMIMO hybrid beamforming system is viewed as non-orthogonal angle division multiple access to simultaneously serve multiple users using same frequency band [28]. To improve throughput (close to fully digital beamforming) of mmWave MU-mMIMO OFDM systems, hybrid beamforming with user scheduling algorithm is defined for DL where BS allocates frequency resources for members of OFDM user group (users with identical strongest beams). Analog beamforming vectors are used to find the optimal beam of each user and digital beamforming is used to get best performance gain (by decreasing residual inter-user interference) [29]. Low computational complexity hybrid beamforming scheme for mmWave MU-MIMO UL channel reduces inter-user interference and its performance is close to the corresponding fully digital beamforming design [30]. Optimal decoupling designs for analog precoder and combiner of a DL MU-FDD mMIMO hybrid beamforming systems are defined by selecting the strongest Eigen-beams of receiving covariance matrix with limited instantaneous CSI. Simulation results proved the need of second order channel statistics in designing digital precoder to reduce intergroup interference [31]. Coordinated RF beamforming technique in mmWave mMIMO systems based on ''Generalized Low Rank Approximation of Matrices'' needs only composite CSI instead of complete physical channel matrix. This technique provides competing solution by considering the coordination between BS and UEs to get maximal array gain with no dimensionality constraint in both TDD and FDD systems [13]. Hybrid beamforming is designed for mmWave MU-mMIMO relay system to enhance the sum-rates (decreasing sum MSE between received signal of digital and hybrid beamforming designs) using digital beamforming. The total sum-rates are increased with accuracy of angles of arrival/departure and number of RF chains [32]. Optimal unconstrained precoder and combiner algorithms for mmWave mMIMO systems are designed for the feasibility of low-cost analog RF hardware implementations. Numerical results of the proposed algorithms [33] showed that spectral efficiency of mmWave systems with transceiver hardware constraints approaches the unconstrained performance limits. Low-complexity phased-ZF hybrid precoding is applied in RF domain (to get large power gains) and low dimensional ZF precoding is used in baseband domain (for multi-stream processing) in mmWave MU-mMIMO systems [34]. The tradeoff between energy efficiency (in bits/J) and spectral efficiency (in bits/channel/MS) is quantified for a small-scale fading channel of MU-mMIMO systems to achieve higher spectral and energy efficiencies using ZR or MRC and UL pilot signals at BS [35].

Proposed methodology
In In hybrid beamforming, the number of TR modules (N T RF ) are less than number of antenna elements (N T ) and each antenna element is connected to one or more TR modules for higher flexibility.
The mathematical representation of hybrid beamforming is as follows: Combining weights matrix; W The matrices F RF and W RF define signal phase values. To achieve optimal weights for precoding and combining there exists some constraints during optimization process. For ideal case, resulting combination of matrices F BB *F RF and W RF *W BB are equal to calculation of F and W without any constraints.

Block diagram of MU-mMIMO system
The Fig. 4 represents a block diagram that gives complete steps of data processing in a MU-mMIMO system. At the Tx, multiple user's data is channel encoded using convolutional codes. The channel encoded bits are mapped to equivalent QAM complex symbols and generate mapped symbols from bits/user. The QAM data of each user is divided into multiple transmit data streams. Digital baseband precoding is used to assign weights for subcarriers of transmit data streams. In this paper, precoding weights are computed using ''hybrid beamforming with peak search (HBPS)'' algorithm as it performs better for larger arrays of mMIMO systems and these precoding weights are used to get corresponding combining weights at the Rx. HBPS gives all digital weights and identifies N T RF , N R RF peaks to get corresponding analog beamforming weights instead of searching iteratively for dominant mode (data streams that use most dominant mode of MIMO channel gives higher SNR) in the channel matrix. The resultant digital signal is modulated using OFDM modulation with pilot mapping followed by RF analog beamforming is performed for all Tx antennas. The modulated signal is transmitted through a rich-scattering MU-mMIMO channel and it is demodulated, decoded at Rx side as shown in Fig. 4. Channel sounding and estimation are performed at the Tx and Rx respectively using ''joint spatial division multiplexing (JSDM)'' algorithm as it allows large number of BS antennas with minimum CSI feedback from UEs in a MU-mMIMO downlink channel.

Analysis of mu-mMIMO system design
To perform MU-mMIMO transmission in mmWave cellular communication systems, high-dimensional channels need to be estimated for designing MU precoder. Digital precoding gives high performance at the cost of hardware complexity and power consumption (more number of RF chains and ADCs). On the other hand, analog precoding has less complexity, but with limited performance (it supports only one data stream). Hybrid precoding for MU-mMIMO systems shown in Fig. 5 is a combination of digital precoding and analog precoding. Hybrid precoding design reduces number of RF chains and also maintain spatial multiplexing gain in mmWave MU-mMIMO system. In MU-mMIMO systems, precoding and combining techniques are used to improve signal energy in the direction and channel of interest with help of available channel information at Tx. In a single user MIMO, the benefits of an antenna array are less because of lower channel rank. On the other hand, MU-MIMO system Wireless Networks (2021) 27:1925-1939 1929 creates rich effective channels through spatial separation of users.
The MIMO channel can be characterized asy ¼ Hs þ n ð3Þ where y received vector, s transmitted vector n noise vector The channel matrix, H ¼ where h ij is a complex gaussian random variable that models fading gain between the ith transmitted and jth receiver antenna.

MU-MIMO DL
If CSI is known, diagonalization of channel matrix (H) gives unconstrained optimal precoding weights by taking the first N T RF dominating modes. So, assume that BS has CSI at the Tx (CSIT) then we can perform MU precoding where it is possible to send signals to all users at the same time and frequencies (T-F), still allow users to recover signals with low complexity. The DL link signal or observation vector ð Þ for user where k number of users, x k signal intended for user k, H k channel from BS to user k, n k noise.
2nd term in Eq. (4) represents signals for other users. Instead of sending user data stream x k , we perform precoding as x k = W k s k . Therefore, DL link signal for user W k precoding matrix s k transmitted QAM symbols for the user k.
2nd term in Eq. (5) represents data precoding for other users.
For a case where each user has only one antenna (therefore, one data stream/user), the size of precoding matrix W k is N t 9 1.
The scalar observation of user where h T k DL channel for user k, n k scalar noise. The 1st term in Eq. (6) represents effective channel seen by the user and 2 nd term represents interference due to other users.
After considering all vectors for all users of a MU À MIMO system; y ¼ HWs þ n ð7Þ y ¼ transmitted QAM symbols for all users.
It is essential to design precoding matrix W at the BS using schemes like ZF, MMSE so that the performance of a MU-MIMO system is optimum. The BS sends data to users that change over time and ''user scheduling'' is used for selection of best K users among a large set to have good effective channel H. Precoding is also possible with statistical CSI (like covariance matrix of channel) instead of instantaneous CSI.

MU-MIMO UL
At Rx, all signals are added and the observation vector; The UL observation vector,  Transmitted signals of all users n noise vector.

MU-MIMO precoding
Consider K users with one antenna per user and one BS with M antennas then the DL channel matrix, H of size K x M where each row corresponds to DL channel of a single user.
We consider ZF precoding where precoding matrix; where H ? pseudo inverse of the channel. If the channel matrix is poorly conditioned, W may be very large and it leads to high signal transmission power. Therefore, we need to enforce a constraint that the total power, P Total = W k k 2 One way to achieve this is scale the entire precoding matrix with a low enough value and guarantee the equal SNR per user. On the other hand, scale each column to have fixed power and it gives different SNR per user.
The adjustment of fixed powers can be achieved using . p 1 power assigned to user 1, p k power assigned to user k. The channel capacity (bits/sec) = Available spectrum (Hz) x Spectrum efficiency (dB).
Under no CSIT conditions, the MIMO channel capacity (C) bounds between where q SNR, B channel bandwidth, M t number of Tx antennas, M r number of Rx antennas.
The left hand limit value is channel capacity for LOS channel and the right hand limit is for rich scattering channel.
Case 1: Very large number of Tx antennas and channel matrix, H is almost orthogonal rows where the number of rows represent number of Rx antennas (small) and number of columns represent number of Tx antennas (large).
HH H & M t I Mr scaled identity matrix The capacity of MIMO channel; where B bandwidth, q SNR.

MU-mMIMO
MU-MIMO system alone is not feasible as multiple antennas at UEs is not cost effective and gains are modest due to limited number of antennas (\ 10), also signal processing is highly complex. Alternatively, mMIMO systems can have one antenna per user and 100 0 s of antennas at the BS. Consider one BS with M antennas and K users with one antenna each (K \ M). In order to have favorable propagation conditions, users should be sufficiently separated. Assume that a MU-MIMO at BS with channel reciprocity (use TDD) and nearly orthogonal channels (users are sufficiently separated).
For UL channel, the channel matrix can be decomposed as

MU-mMIMO: UL
The observation model for the UL is given as H has many rows and few columns and x is a vector (size K 9 1) contains signals for each of K users.
BS can apply the following shaper or combiner based on the property The combining matrix; where w~CN (0, N 0 MD), z k = Md k x k ? w k , and w k-* CN (0, N 0 Md k ). The SNR of k th user is given as, SNR k = Md k E x /N 0 where N 0 noise power, N 0 Md k noise variance Therefore, Rate for the kth user; Total sum -rate at BS is sum of all individual rates of 0 K 0 users; R sum ¼ Blog 2 Equation (20) indicates the capacity of the system. This means that, with asymptotically large antennas and simple linear combining at the BS gives optimal results even in the presence of multi-user interference.

MU-mMIMO: DL
The DL is little more complicated because BS needs to do additional processing (precoding) so that users do not see interference from other users.
Observation model for DL with M Â K precoding matrix is given as; where H T downlink channel, W suitable precoding matrix, n noise vector seen by each user, s vector of size K 9 1 and contains data for each of users.
From MU-MIMO, a good precoder is pseudo inverse of channel matrix. Therefore, we choose precoding of the form where H * complex conjugate of the channel HM scaling factor.
D p diagonal matrix that contains square root of powers for each user.
To meet the power constraints, power allocation D p to ensure W k k 2 = trace(W H W) = P Total . We find that W H W = HD p H T H * HD p /M = HD p DDHD p = D p D. D p matrix of powers along the diagonal D matrix of path loss values along the diagonal The observations seen across different users; At user k: the observations can be split into y k The channel in the DL for user i, h i T = Hd i g i T h i T row vector of length M.
We conclude that MU-mMIMO DL channel is a matched filter with asymptotically optimal linear precoder. This means that with simple linear processing at BS, it is possible to remove all MU-interference in DL among all users.

Results analysis
According to 3GPP, future mmWave wireless communication systems recommend 28 GHz frequency band for MU-mMIMO [36]. mmWave MU-mMIMO communication link between the BS and UEs is validated using scattering-based MIMO spatial channel model with ''singlebounce ray tracing'' approximation. This model considers UEs at different T-R spatial locations and randomly placed multiple scatterers. A single channel is used for sounding as well as data transmission, allows path loss modelling with LOS and non-LOS scenarios (these are closer to real scenarios). The channel matrix is updated periodically to mimic variations of MIMO channel over time. The radiation patterns for antenna arrays are isotropic with rectangular or linear geometry. The simulations were performed for a maximum of 256 9 16 MU-mMIMO system for four users and eight users with parameters shown in Tables 2  and 3. At Tx end, 256-element rectangular antenna array is used with 4-RF chains and at Rx, 16-element square array is used with 4-RF chains. Therefore, each antenna element is connected to 4 phase shifters and there are 4-RF chains. To have maximum spectral efficiency of MU-mMIMO, each user is assigned with independent channels and each RF chain is used to send an independent data stream, therefore maximum of 4 streams is supported. The MU-mMIMO Rx at UE of each user is modeled by compensating for path loss, thermal noise. Figures 6,7,8,9, 10 and 11 represent the RMS EVM values in four users and eight users MU-mMIMO systems of 16-QAM, 64-QAM and 256-QAM modulation schemes with multiple number of BS antennas. It is observed that, for users with only one data stream, the RMS EVM is very high compared to users with multiple data streams. For a give modulation scheme, RMS EVM is decreasing as number of BS antennas are increased for users with single data stream. For users with more number of data streams, optimum RMS EVM values are achieved for 128 BS antennas and there is very slight increment in error values for 256 BS antennas. Figures 12, 13, 14, 15, 16 and 7 represent the RMS EVM values in four users and eight users MU-mMIMO systems with 64, 128, and 256 BS antennas respectively for various modulation schemes. The following observations are made from these figures: for the given number of BS antennas, users with single data stream, the RMS EVM is minimizing at higher modulation order and increasing number of BS antennas. When compared with users with more data streams, the reduction rate of RMS EVM is very high in users with single data stream as the number of BS antennas are increased from 64 to 256. For users with more data streams, optimum performance is achieved with 128 BS antennas and lowest RMS EVM values are obtained. RMS EVM values are very slightly increased for 256 BS antennas compared to 64 BS antennas. Interestingly, for a given number of BS antennas there is almost no change in RMS EVM values when the modulation order is increased from 4 to 8. Figures 18, 19 and 20 represent the equalized symbol constellation per data stream in MU-mMIMO system with different combinations of modulation schemes and number of BS antennas. In the constellation diagrams, variance of the recovered streams is more at the users with lower number of independent streams. This is due to the absence of dominant modes in the channel and it causes poor SNR. From the receive constellation diagrams of all combinations, it is very clear that the variance of recovered data streams is very much better for users with multiple data streams as the symbol points are properly located with less dispersion in symbol constellation diagram. The reason is that, data streams use most dominant mode of scattering MIMO channel and these have higher SNR. On the other hand, the symbol points in equalized symbol constellation of users with single data stream are highly dispersed and they have poor SNR values. Figures 21, 22 and 23 represent the signal radiation patterns in MU-mMIMO wireless systems with multiple BS antennas. The stronger lobes of 3D response pattern in mMIMO designs represent distinct data streams of users.
These lobes indicate the spread achieved by hybrid beamforming. From these figures, it is clear that the radiation beams are becoming sharper for increasing number of BS antennas and this increases the reliability of signal there by throughput.

Conclusions and future scope
A mmWave DL MU-mMIMO hybrid beamforming communication system is designed with multiple independent data streams per user. From the overall results, it is         Wireless Networks (2021) 27:1925-1939 1935 observed that for users with lower number of independent data streams, RMS EVM values are higher and for more number of independent data streams, the RMS EVM is less. Therefore, the increasing number of data streams per user leads to decrease in RMS EVM values. For the given modulation scheme, as the number of BS antennas are increasing, the RMS EVM is decreasing. It gives trade-off between the number of Tx/Rx antenna elements and multiple data streams per user. It has been concluded that if the user data is divided into more number of parallel data    streams then it requires less number of active antenna elements to transmit the signals. From the simulation results, it is concluded that 256 9 16 MIMO system is more suitable for eight users and 128 9 16 MIMO system for four users with multiple data streams. It is strongly recommended to have higher number of independent data streams per user in a mmWave MU-mMIMO systems to achieve higher order throughputs. For future research    Wireless Networks (2021) 27:1925-1939 1937 directions, the MU-mMIMO hybrid beamforming designs with an objective of reduction in RMS EVM values for users with single (or less) number of independent data streams and achieve higher throughput are of great investigation interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.
Funding Open access funding provided by Manipal Academy of Higher Education, Manipal.