Virtual-crystal-approximation-

this simple script is to give the combination of two different elements for alloy to calculate the harmonic dynamic properties in almabte

requirements: sencond order force constants for two different elements and put into current directory

for example: FORCE_CONSTANTS_2ND_Si and FORCE_CONSTANTS_2ND_Ge to generate FORCE_CONSTANTS_VCA

then apply phonopy to get what you want.