We study a robust and efficient eigensolver for computing the positive dense spectrum of the two-dimensional transmission eigenvalue problem (TEP) which is derived from the Maxwell’s equation with complex media in pseudo-chiral model and the transverse magnetic mode. The discretized governing equations by the N ́ed ́elec edge element result in a large-scale quadratic eigenvalue problem (QEP). We estimate half of the positive eigenvalues of the QEP are on some interval which forms a dense spectrum of the QEP. The quadratic Jacobi-Davidson method with a so-called non-equivalence deflation technique is proposed to compute the dense spectrum of the QEP. Intensive numerical experiments show that our proposed method makes the convergence efficiently and robustly even it needs to compute more than 5000 desired eigenpairs. Numerical results also illustrate that the computed eigenvalue curves can be approximated by the non- linear functions which can be applied to estimate the density of the eigenvalues for the TEP.