We report the effects of electron lattice coupling on the charge density distribution of two dimensional silicene nanoribbons. By means of an extended tightbinding model with lattice relaxation, we provide theoretical explanation to the unexpected density of states, as well as other electronic properties observed in this kind of material. The methodology follows the Ehrenfest approximation. First, by using the self-consistent field scheme, the Hamiltonian's expectation value leads to an eigenvalue problem in terms of an initial set of coordinates. Under the stationary condition, the Ehrenfest scheme leads to an updated set of coordinates that depends on the eigenstates of the electronic problem. Then, one can use this new set of coordinates as input to the eigenvalue problem. The procedure is repeated until the achieving a given convergence criteria [1,2]. Spin-orbit coupling effects are introduced in order to treat the systems more accurately. The effect of the nanoribbon width is investigated and interesting patterns arise that are expected to greatly influence transport in these materials [3]. Importantly, an extensive analysis on the role played by electron-phonon lattice over the charge distribution is performed. Our results point towards the expected ranges that charge mobility might have in these materials. This contribution is expected to enlighten the behavior of charge distribution in silicone nanoribbon, which are systems of crucial technological appeal [4].