Geometry optimization and DFT in psi4 is very fast: comparable to ORCA 5, and much faster than NWChem (or at least easier to configure to be fast, i.e., density fitting).
However, calculating the Hessian matrix with DFT using numerical methods is painfully slow. I also cannot find a way to restart a frequency job that stops partway through, so they really waste a lot of computer time.
When are analytic frequency methods going to be implemented in psi4?
And if they aren’t planned, what is the obstacle to their implementation?
Analytic frequencies are desired but not planned. We only have so much time to develop Psi4, and the groups who do the most work on that are focusing on other areas. (To give you a smattering: new SCF algorithms, new SCF convergence techniques, CC properties, DLPNO methods, f12 methods…) None of the major contributors are DFT specialists.
If you or somebody else wants to contribute these to Psi4, hop on our Slack channel and ask for guidance.