You don't really simulate the polaron, it's more an interpretation of the output results. But yes, QuantumATK can do precisely the same calculations as in the paper you mentioned, either with planewave DFT or (probably a lot faster) using LCAO basis sets.