[Thermo_pw-forum] nonstoichiometric compounds

[Gunnar Palsson]

Dear Andrea, 

First I would like to tell you that I find your thermo_pw package extremely useful in many ways. I was inspired by your article J. Phys. Condensed Matter 26, 335401 (2014) and have been practising using your PS library v1.0 and the PAW potentials, for pure vanadium and niobium to calculate the thermoelastic properties. I have been thinking for quite a while about calculating the thermophysical properties and changes in free energy of metal hydrides, specifically vanadium hydride, in an attempt to determine the phase boundaries. I am also interested in the changes in the elastic constants with temperature and hydrogen. The question I have is whether one can use thermo_pw to calculate the Helmholtz free energy in the quasi harmonic approximation as a function of temperature for a supercell of vanadium with few hydrogen atoms. The size of the supercell determines what part of the phase diagram one can access in terms of concentration. So for 16 vanadium atoms, one can go from 1/16 to 16/16 in H/V concentration for example. It is my understanding that this requires ibrav=0, which may not be fully implemented? If not it would still be valuable to look into the stoichiometric hydride, V2H, and VH. The V2H belongs to space group 141. Finally, I have some neutron scattering data on hydrogen containing vanadium, and it would be extremely valuable to use the output of thermo_pw to calculate various neutron scattering cross sections and compare with experiment, which I can probably do. 

Best regards,
Dr. Gunnar Palsson,
Researcher,
Uppsala University.