We consider Markovian multiserver retrial queues where a blocked customer has two opportunities for abandonment: at the moment of blocking or at the departure epoch from the orbit. In this queueing system, the number of customers in the system (servers and buffer) and that in the orbit form a level-dependent quasi-birth-and-death (QBD) process whose stationary distribution is expressed in terms of a sequence of rate matrices. Using a simple perturbation technique and a matrix analytic method, we derive Taylor series expansion for nonzero elements of the rate matrices with respect to the number of customers in the orbit. We also obtain explicit expressions for all the coefficients of the expansion. Furthermore, we derive tail asymptotic formulae for the joint stationary distribution of the number of customers in the system and that in the orbit. Numerical examples reveal that the tail probability of the model with two types of nonpersistent customers is greater than that of the corresponding model with one type of nonpersistent customers.We refer to [12,4,8,15,16,17] for effort to find analytical solutions for M/M/c/c retrial queues with more than two servers by the generating function approach. Kim [12] and deal with the case of three servers while Choi and Kim [4] derive analytical solution for a model with feedback. It should be noted that some technical assumptions are imposed in these papers. Using an alternative approach, Phung-Duc et al. [16] show that the joint stationary distribution is expressed in terms of continued fractions for the cases of c = 3 and 4, without any technical assumption as presented in literature. The same authors in [17] further derive analytical solutions for the joint stationary distribution of state-dependent M/M/c/c + r retrial queues with Bernoulli abandonment, where c + r ≤ 4. Pearce [15] presents an expression for the joint stationary distribution in terms of generalized continued fractions for the M/M/c/c retrial queue with any c. Although, the formulae in [15] do not directly yield a numerical algorithm, this is one of the seminal papers providing the most general analytical results for the model. Recently, asymptotic analysis for multiserver retrial queues has been receiving considerable attention. Liu and Zhao [11] use a censoring technique and a level-dependent QBD approach to derive analytical solutions for the M/M/c/c retrial queues for the cases of c = 1 and 2 and investigate the asymptotic behavior for the stationary distribution of the general case with any c. Using the same approach, Liu et al. [10] extend their study to M/M/c/c retrial queues with one type of nonpersistent customers. Kim et al. [13] derive more detailed asymptotic formulae for the joint stationary distribution of M/M/c/c retrial queues in comparison with those obtained by Liu and Zhao [11]. Furthermore, Kim and Kim [14] refine the asymptotic result obtained by Liu et al. [10]. The methodology of [13,14] is based on an investigation of the analyticity of generating functions. However, the asymptoti...