<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:large">Dear Dr. Lorenzo</div><div class="gmail_default" style="font-size:large">The file size is big (134MB) and it is failing to download because of my poor internet connection. I need to find a good network.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">If you do not mind, I would like to ask you to have a look at my small shell script that can calculate the vibrational entropy using dynmat.out file. </div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">NOTE: It requires only a dynmat.out file and degree of freedom for the molecule.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">I used the Eqn of "S" provided <a href="https://github.com/DavidCdeB/QHA/blob/master/Images_for_README_md/entropy.png">HERE</a> to write this script.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">#!/bin/bash<br>#========Constants=================<br>DOF=5               #Trans. + Rot. DOF for the molecule, Here I use 5 (3+2) for H2, for 3D Molec. I would take it 6<br>THz2Hz='1e+12'                 #Freq in s^(-1)<br>KJPerMol=96                    #KJ/mol<br>KCalPerMol=23                  #Kcal/mol<br>Kb='8.6173303E-5'            #boltzmann's constant in eV/K<br>T=450                                #Temperature in K<br>Pconst='4.135667696e-15'     #planck's constant in eV.S<br>K_to_eV='0.00008617328149741'<br>#=== How many Frequencies<br>Total_freq=$(grep freq dynmat.out | cut -f 2 -d = | cut -f 1 -d [ | wc -l)<br>#==== How many vibrational Frequencies <br>Vibrational_Freq=$(echo "$Total_freq" "$DOF" | awk '{print $1-$2}')<br><br>#==== Summing all vib. Freq., converting into Hz  and storing as a variable<br>FreqHz=$(grep freq dynmat.out | cut -f 2 -d = | cut -f 1 -d [ | tail -n $Vibrational_Freq | awk '{print $1}' | awk -F',' '{sum+=$1;}END{print sum;}' | awk '{print $1*(1e+12)}')<br><br>#===Calculation of vibrational entropy using formula of entropy from [1] in eV/K  and storing as a variable<br>Freq_in_eV_per_K=$(echo "$Kb" "$T" "$Pconst" "$FreqHz" | awk '{print -$1*(log(1-exp(-(($3*$4)/($1*$2))))) + ($3/$2)*($4/(exp(($3*$4)/($1*$2))-1) )}')<br><br>#=== Converting eV/K to eV and storing as a variable<br>eV_per_K_to_eV=$(echo "$Freq_in_eV_per_K" "$K_to_eV" | awk '{print $1*$2}')<br>#converting eV to KCal/mol<br>echo "$eV_per_K_to_eV" "$KCalPerMol" | awk '{print $1*$2}'<br><br>#converting eV to KJ/mol<br>echo "$eV_per_K_to_eV" "$KJPerMol" | awk '{print $1*$2}'<br></div><div class="gmail_default" style="font-size:large">#END</div><font face="georgia, serif" size="4">I look forward to hearing from you or anyone from the  QE  list<span class="gmail_default" style="font-size:large"> whether it is correct or not.</span></font></div><div dir="ltr"><font face="georgia, serif" size="4"><span class="gmail_default" style="font-size:large"></span></font><span style="font-size:large;font-family:georgia,serif"></span></div><div dir="ltr"><div class="gmail_default" style="font-size:large">Regards</div><div class="gmail_default" style="font-size:large">Bhamu</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 24, 2021 at 2:46 PM Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><blockquote><div><div>Let me try D3Q code.</div><div><br></div></div></blockquote><div><div><br><div>If you can compile your own version of QE, you may want to use a recent snapshot of espresso+d3q, as I have added a few features recently. </div><div><a href="https://link.getmailspring.com/link/D84D3DDC-75E3-404B-A829-E30C74FB10F0@getmailspring.com/0?redirect=https%3A%2F%2Fmycore.core-cloud.net%2Findex.php%2Fs%2F1RxbPmcxvHunFVF&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D" title="https://mycore.core-cloud.net/index.php/s/1RxbPmcxvHunFVF" target="_blank">https://mycore.core-cloud.net/index.php/s/1RxbPmcxvHunFVF</a></div><div>cheer</div></div></div><br><div><u></u>--<br><span style="font-size:0.9em;color:gray">Lorenzo Paulatto - Paris</span></div></blockquote></div></div>