Abstract: In steels, the interaction between screw dislocations and carbon solutes has a great influence on the yield strength. Fe-C potentials used in molecular dynamics (MD) simulations yield a poor description of screw dislocation properties-their core structure and Peierls barrier-compared to ab initio calculations. Here we combine two EAM potentials from the literature, which greatly improves dislocation property accuracy in FeC alloys. Using this hybrid potential, MD simulations of dislocation glide in random solid solutions confirm a powerful solute strengthening, caused by complex interaction processes. We analyze these processes in a model geometry, where a row of carbon atoms is inserted in the dislocation core with varying separations. We use a combination of MD simulations, minimum-energy path calculations, and a statistical model based on the harmonic transition state theory to explain the strengthening induced by carbon. We unveil that carbon disrupts the glide process, as unpinning requires the successive nucleation of two kink pairs. When solute separation is below about 100 Burgers vectors, the activation enthalpy of both kink pairs are markedly increased compared to pure iron, resulting in a strong dependence of the unpinning stress on solute spacing. Our simulations also suggest an effect of carbon spacing on the kink-pair activation entropy. This work provides elementary processes and parameters that will be useful for larger-scale models and, in particular, kinetic Monte Carlo simulations.