Phase field fracture modelling using quasi-Newton methods and a new adaptive step scheme