This protocol calculates and prints Hessian, the matrix of second derivatives of the energy. By default, it also writes it to a file that can be used as an input for the calculation of vibrational frequencies. Translational and rotational degrees of freedom are projected out from the calculated Hessian.