Not sure of the reason. Have you looked at the log from PATH? Is the time being spent in JuMP? Or is PATH taking many more iterations to converge? You might consider making the lower bound of the denominator larger.
If you’re going to manually code the gradient, you should almost just use PATHSolver.solve_mcp directly: