We present a method for the automatic determination of transition states (TSs) that is based on Grimme's RMSD-PP method for the rapid estimation of reaction paths using the GFN-xTB semiempirical tight binding models (Figure ). The RMSD-PP method estimates a reaction path between reactants and products by a geometry optimisation using an energy function augmented by a Gaussian biasing potential that "pushes" and "pull" the structure away from the reactant and towards the product. Our method starts with a series of RMSD-PP calculations with increasingly larger push and pull strengths until reaction completion. The additional structures near the highest point on the reaction path are generated by interpolation and used for DFT single points and the highest energy structure is then used as an initial guess for a TS search. Upon convergence the TS is tested by an IRC calculation and if the TS is found to be incorrect then the initial guess structure is reoptimised with key bond lengths constrained and used as an initial guess for a new TS search. If that fails the entire procedure is repeated but with using an electronic temperature of 6000K for the RMSD-PP calculations.