In this paper, we consider the compressible bipolar Navier-Stokes-Poisson equations with a non-flat doping profile in three-dimensional space. The existence and uniqueness of the non-constant stationary solutions are established when the doping profile is a small perturbation of a positive constant state. Then under the smallness assumption of the initial perturbation, we show the global existence of smooth solutions to the Cauchy problem near the stationary state. Finally, the convergence rates are obtained by combining the energy estimates for the nonlinear system and the L 2 -decay estimates for the linearized equations.x, t/ represent the particle densities, the velocities, and the electrostatic potential, respectively, at time t 0 and position x 2 R 3 . The pressure P i D P i . i / D A i i i with constant A i > 0 and the adiabatic exponent i > 1 for i D 1, 2. The viscosity coefficients and satisfy the usual physical conditions > 0, C 2 3 0. The function b.x/ is the doping profile satisfying b.x/ ! N b, as jxj ! 1, for a positive constant state N b > 0. When only considering the dynamics of one particle in semiconductor devices and plasmas, we have the unipolar Navier-Stokes-Poisson equations (NSP). More details on the (NSP) equations can be found in [1][2][3]. Recently, some important progress has been made for the unipolar (NSP) system. For the pressure law p. / D with the adiabatic exponent > 3=2, the global existence of weak solutions was obtained by [4] when the spatial dimension is three in the framework of Lions and Feireisl for the compressible Navier-Stokes equations [5,6]. When the doping profile is flat, that is, b.x/ D N b, we refer to [7][8][9][10][11] framework. Hao and Li [7] and Zheng [11] framework established the global strong solutions of the initial value problem for the multi-dimensional compressible (NSP) system