[QE-users] 0.0000 in QHA outputs

Dr. K. C. Bhamu kcbhamu85 at gmail.com
Wed Mar 24 05:10:45 CET 2021


Dear Quantum Espresso/QHA Users
I am trying to calculate the entropy of a molecule (49atoms, C and H) in
the gas phase using QHA provided with QE6.6.
I could run the examples supplied with the QHA code.
But when I run it for H2 molecule, I always get 0.0000 in the
case(H2).QHA.out file.
Similar information (0.000) is written in other files also.

What could be wrong?
I have edited the Edit_Me file properly.
Few header rows from H2.QHA.out file is mentioned below while the entire
directory can be downloaded from here <https://we.tl/t-EfBU7ZFvWE>.

# Zero vibration energy:      0.0000000000  (Ry/cell)
# Phonon DOS norm      :      0.000000      ! 3N for check purpose, N
number of atoms in the unit cell
# T in K, F_vib in Ry/cell, C_v in R (the universal gas constant by 3N
modes), S in k_B
#
#   T         E_internal        F_vibration          Specific heat (C_v)
    Entropy
############################################################################################################
    5.00      0.0000000000      0.0000000000           0.0000000000
 0.0000000000  *>> I used step 0.75 too*
   10.00      0.0000000000      0.0000000000           0.0000000000
 0.0000000000
   15.00      0.0000000000      0.0000000000           0.0000000000
 0.0000000000
   20.00      0.0000000000      0.0000000000           0.0000000000
 0.0000000000

Any help would be appreciated.

Regards

KC Bhamu(Ph.D.)
University of Ulsan
South Korea
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20210324/2ec3dda6/attachment.html>


More information about the users mailing list