An integral and expensive part of coupled climate model simulations is the gas-phase chemistry which gives rise to hundreds of coupled differential equations. We propose a method which improves the convergence and robustness properties of commonly used Newton–Raphson solvers. The method is flexible and can be appended to most algorithms. The approach can be useful for a broader community of computational scientists whose interests lie in solving systems with intensive interactive chemistry.