One shot force constants calculation: under the implementation