We give the approximate analytic solutions of the Dirac equations for the Rosen-Morse potential including the spin-orbit centrifugal term. In the framework of the spin and pseudospin symmetry concept, we obtain the analytic bound state energy spectra and corresponding two-component upper-and lower-spinors of the two Dirac particles, in closed form, by means of the Nikiforov-Uvarov method. The special cases of the s-wave κ = ±1 (l = l = 0) Rosen-Morse potential, the Eckart-type potential, the PT-symmetric Rosen-Morse potential and non-relativistic limits are briefly studied. Within the framework of the Dirac equation the spin symmetry arises if the magnitude of the attractive scalar potential S(r) and repulsive vector potential are nearly equal, S(r) ∼ V (r) in the nuclei (i.e., when the difference potential ∆(r) = V (r) − S(r) = C s = constant). However, the pseudospin symmetry occurs if S(r) ∼ −V (r) are nearly equal (i.e., when the sum potential Σ(r) = V (r) + S(r) = C ps = constant) [1-3]. The spin symmetry is relevant for mesons [4]. The pseudospin symmetry concept has been applied to many systems in nuclear physics and related areas [2-7] and used to explain features of deformed nuclei [8], the super-deformation [9] and to establish an effective nuclear shell-model scheme [5,6,10].The pseudospin symmetry introduced in nuclear theory refers to a quasi-degeneracy of the single-nucleon doublets and can be characterized with the non-relativistic quantum numbers (n, l, j = l + 1/2) and (n − 1, l + 2, j = l + 3/2), where n, l and j are the single-nucleon radial, orbital and total angular momentum quantum numbers for a single particle, respectively [5,6]. The total angular momentum is given as j = l + s, where l = l + 1 is a pseudoangular momentum and s = 1/2 is a pseudospin angular momentum. In real nuclei, the pseudospin symmetry is only an approximation and the quality of approximation depends on the pseudo-centrifugal potential and pseudospin orbital potential [11]. Alhaidari et al. [12] investigated in detail the physical interpretation on the three-dimensional Dirac equation in the context of spin symmetry limitation ∆(r) = 0 and pseudospin symmetry limitation Σ(r) = 0. Some authors have applied the spin and pseudospin symmetry on several physical potentials, such as the harmonic oscillator [12][13][14][15], the Woods-Saxon potential [16], the Morse potential [17,18], the Hulthén potential [19], the Eckart potential [20][21][22], the molecular diatomic three-parameter potential [23], the Pöschl-Teller potential [24], the Rosen-Morse potential [25] and the generalized Morse potential [26].The exact solutions of the Dirac equation for the exponential-type potentials are possible only for the s-wave (l = 0 case). However, for l-states an approximation scheme has to be used to deal with the centrifugal and pseudo-centrifugal terms. Many authors have used different methods to study the partially exactly solvable and exactly solvable Schrödinger, Klein-Gordon (KG) and Dirac equations in 1D, 3D and/or any D-dimen...