Atomistic modeling of light driven electron dynamics are important in studies of photoactive materials. Spin-resolved electronic structure calculations become necessary when dealing with transition metal, magnetic, and even some carbon materials, intermediates, and radicals. An approximate treatment can be pursued in the basis of spincollinear density functional theory. Most transition-metal compounds exhibit open shell nonsinglet configurations, necessitating special treatment of electrons with α/β spin projections. By separate treatment of electronic states with the α/β spin components one is able to describe a broader range of materials, identify new channels of relaxation and charge transfer, and provide knowledge for rational design of new materials in solar energy harvesting and information storage. For this methodology, named spin-resolved electron dynamics, spin-polarized DFT is used as the basis to implement nonadiabatic molecular dynamics. At ambient temperatures, the thermal lattice vibrations results in orbital and energy fluctuations with time. Nonadiabatic couplings are then calculated, which control the dissipative dynamics of the spin resolved density matrix. Different initial excitations are then analyzed and used to calculate relaxation dynamics. Spin-resolved electronic dynamics approach (SREDA) is applied to study vanadium(IV) substitutionally doped bulk anatase in a doublet ground state. The results show that a difference in the electronic structure for α and β spin components determines consequences in optical excitations and electronic dynamics pathways experienced by electrons with α and β spin projections. Specifically, the lone occupied V 3d α-orbital increases the range of absorption and defines the rates and pathways of relaxation for both holes and electrons with α-spin projection. Optical excitations involving occupied V 3d α-orbital are responsible for IR-range absorption, followed by nonradiative relaxation. Certain transitions involving orbitals of α-spin component occur in the visible range and induce localization of a negative charge on the V ion for an extended time period. The slower nonradiative relaxation rate of α-excitations is rationally explained as a consequence of difference of electronic structure for α and β spin projections and specific pattern of energy levels contributed by doping. Specifically, excitations involving orbitals with α-projection of spin experience transitions through larger subgaps in the conduction band compared to the ones experienced by similar excitations involving orbitals with β-projection of spin. It is anticipated that this methodology can be broadly implemented on multiple applications of transition metal based materials, including optoelectronics, information storage, laser crystals, dyes, photovoltaic materials, and metal oxides for photoelectrochemical water splitting.