Optimize the geometry using the L-BFGS-B function minimizer. This is based on the BFGS optimizer, but it does not store the inverse Hessian, instead it is calculated as needed. Because of this, the L-BFGS-B method uses very little storage, and is therefore suitable for optimizing very large systems.

The L-BFGS optimizer is the default if 100 or more variables are to be optimized.  For such systems, L-BFGS is more efficient than EF.  If L-BFGS is not wanted, add keyword EF