Hi all:
Recently days, I am learning how to consider vibrational contribution in cluster expansion. When I get the static calculation result and fit successfully, I run ‘foreachfile str_relax.out pwd ; svsl -d’. However the program inform me a lot of unstable configurations. I know from manual that we need to find out this unstable mode and generate perturbations along the unstable directions. However what I concern is, does those modes matters in vibrational free energy of a structure using the Stiffness VS Length method? Because the configurations I selected in fitting stiffness are all ground states, they are tought to be stable by program. Those configurations who are unstable are acutally above ground states. In principle, they should not influence the fitting process and vibrational energy calculation, right? So, in this case, can I ignore those unstable modes and force svsl to continue on all structures including those unstable modes by svsl -fn?