“…The IMLS, PIP, NN, and PIP‐NN fitting methods are based on the solution of a weighted least‐squares problem with rigorous implementations of different basis functions, which can be extended to high dimensional PESs easily. The PIP‐NN method was developed very recently by combination of the PIP and NN methods, which leads to a smaller root mean square error and is efficient for fitting the PESs of molecular systems, especially for polyatomic systems . In this review, the PESs of the N + OH, N + CH, and N + C 2 reaction systems are all fitted with the cubic spline interpolation method, which is simple and direct.…”