The adaptive density-guided approach (ADGA) for constructing PESs is an iterative method where a wave function density is used to guide the placement of SPs on a grid of molecular displacements until a convergence is reached. The overall goal of ADGA is to automatically map out an accurate PES with as few SPs as possible. This is done by creating an initial grid for each mode and then iteratively subdivide and extend these grids until all one-dimensional grids have obtained convergence. When the one-dimensional grids are converged higher-dimensional grids for the MCs included are created from the boundaries of the one-dimensional grids. These MCs are then also iteratively subdivided until convergence is achieved.