[Pw_forum] fqha.x

Eyvaz Isaev eyvaz_isaev at yahoo.com
Wed Apr 4 19:38:59 CEST 2007


Hi again,

I found the fqha.f90 code on the forum and compiled it with mpif77 to get the executable. 

It is not the same I have been talking about, but it does not matter in this particular case. 
Besides, there is no need to compile it using mpif77 or mpif90. Use just a serial version.
 
I am not a physics major, so I don't understand the code completely. 

In this case you will meet lots of troubles. I think it is worth reading some general purpose  textbook  on solid state physics. As example, textbook by Kittel, or Aschroft and Mermin, or  Harrison, or Martin, or Marder  etc.

So, I have a question - please excuse me if it is grossly incorrect. I expect the free energy F to 
equal E at zero temperature. But, the F I obtain from the code is positive at T=0 and becomes negative at larger T. Could you please explain this?

According to quantum-mechanical principles atoms in solid states are not frozen at T=0K, but vibrate due to the uncertinity principle of the quantum mechanics, so there is so called zero vibration energy which is positive by its definition (for a stable structure all frequences are positive). At higher T the phonon contribution to the free energy becomes negative.

Bests,
Eyvaz. 


Thanks,
Shriram

Eyvaz Isaev wrote:
> Hi Shrihram,
>  
>   I'm facing some trouble running this program and hope someone can help 
> me out. I'm using the dos file generated in example 06 to calculate the 
> free energy. But, when I give the input to fqha.x as alas.phdos, for any 
> temperature, the free energy turns out to be zero. I'm not sure what is 
> going wrong here. Can someone please help?
>
> The answer is quite simple: as usual \omega(1) =0, and dos(1)=0, so, you can not 
> calculate  \ln(2sinh(\hbar\omega/2kT)). The starting point has to be 2.
>
> But I think if you find QHA code sent to the forum you will gain much more.
>
> Bests,
> Eyvaz.
>
>
> Thanks,
> Shriram
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
>
>
>
>  
> ____________________________________________________________________________________
> Don't pick lemons.
> See all the new 2007 cars at Yahoo! Autos.
> http://autos.yahoo.com/new_cars.html 
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>   

_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum





 
____________________________________________________________________________________
The fish are biting. 
Get more visitors on your site using Yahoo! Search Marketing.
http://searchmarketing.yahoo.com/arp/sponsoredsearch_v2.php



More information about the Pw_forum mailing list