On the basis of the multi-configuration Dirac–Fock method and the density matrix theory treatment of photoionization (PI) dynamics, a systematic investigation for 4d inner-shell PI of exemplary Ba atom and Ba-like Nd4+, Gd8+ and Yb14+ ions is performed. Starting from the general framework of the relativistic, the relaxation effects are considered by using a set of relaxed one-electron orbitals in solving the coupled Dirac equations. Our results uncover for the first time that the relaxation effects, arising from the overlap integrals between orbitals of the initial state and the relaxed final state, may play a paramount role in determining the magnetic sublevel cross sections, the photoelectron angular-distribution asymmetry parameters, and the fluorescence polarizations of the subsequent x-ray radiation, particularly for neutral atoms close to the threshold regime. Yet the importance of inclusion of the orbital relaxation weakens as the incoming photon energy and/or the degree of atom ionization increases. The nature of these influences is explored, which allows us to get more insight on the physical process and mechanism involved. Our results are consistent with those obtained by the experiments and other theoretical predictions.