The radial dependent potential and neutron production rate in spherical inertial electrostatic confinement fusion (IECF) devices is investigated. The electrostatic potential is determined by solving the Poisson equation for various deuteron and electron distribution functions. The fusion reaction rates are determined using energy distribution function. Also, dependence of potential structure and neutron production rate on some important parameters as the ion and electron convergence, working pressure, kinetic energy of the secondary electrons emitted from the cathode and the fraction of secondary electrons drawn inside the cathode are studied. Total produced neutrons as a function of input power at different working conditions are also obtained.