Pardiso w/ Quasi Newton method?

I have read through the page describing using PARDISO to solve a nonlinear set of equations after they have been linearized, but I noticed the algorithm is just Newton's method applied with the PARDISO solving routing. 

For faster and more reliable convergence, has anyone created a Quasi Newton solving algorithm using the PARDISO routine during the external iterations?

J