Driver RPA

(WIP) Here use random-phase approximation to evaluate operators.