High-order Path Integral Monte Carlo methods for solving quantum dot problems