<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Thank you Giuseppe<br></span></div><div>&nbsp;</div><div>DEBBICHI Mourad<br>Unité de Recherche Physique des Solides,99/UR/13-19,<br>Département de Physique, Faculté des Science de Monastir,<br>Avenue de l'Environnement 5019, Monastir Tunisie.<br>tél:+21697487042<br>mourad_fsm@yahoo.fr<br></div>  <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1">  <b><span style="font-weight: bold;">De&nbsp;:</span></b> Giuseppe Mattioli &lt;giuseppe.mattioli@ism.cnr.it&gt;<br> <b><span style="font-weight: bold;">À&nbsp;:</span></b> debbichi mourad &lt;mourad_fsm@yahoo.fr&gt;; PWSCF Forum &lt;pw_forum@pwscf.org&gt; <br> <b><span style="font-weight: bold;">Envoyé
 le :</span></b> Mardi 19 juin 2012 12h27<br> <b><span style="font-weight: bold;">Objet&nbsp;:</span></b> Re: [Pw_forum] (no subject)<br> </font> </div> <br><br>Dear Mourad<br>No special keywords, only a bit of physical chemistry...<br>Let us suppose we are calculating The formation heat (aka formation enthalpy, DeltaH_f) of wurtzite <br>ZnO. First of all, what is it the formation enthalpy? Suppose that we are not very rigorous, we may <br>say "the energy gained to form ZnO from Zn and O in their standard state". We may consider therefore <br>this chemical process:<br><br>Zn(metal)+O2(gas)--&gt;ZnO(wurtzite)<br><br>How can we calculate such an energy? Well, all chemical "Delta" quantities are defined as a difference <br>between some final state quantity (free energy, enthalpy, ...) and the related initial state <br>quantity. If we consider negligible the PV contribution to the above chemical process (H = U+PV), we <br>can estimate formation heats as
 differences of total energies calculated in supercells:<br><br>DeltaH_f (ZnO) = E(ZnO) - E(Zn) - E(O2)<br><br>You must only be sure that the number of atoms in all the supercells is fully consistent with your <br>calculation:<br><br>ZnO(wurtzite): 4 atoms (2 ZnO units)<br>Zn(metal, hcp): 1 atom<br>O2 molecule: 2 atoms<br><br>Your final equation is therefore<br><br>DeltaH_f (ZnO) = E(ZnO) - 2*E(Zn) - E(O2)<br><br>Modify it according to your purpose...<br><br>HTH<br><br>Giuseppe<br><br><br>On Tuesday 19 June 2012 11:21:54 debbichi mourad wrote:<br>&gt; Dear pwscf users.<br>&gt; How to calculate the formation enthalpy, which command&nbsp; I can added to the<br>&gt; input? best regards<br>&gt;&nbsp; <br>&gt; DEBBICHI Mourad<br>&gt; Unité de Recherche Physique des Solides,99/UR/13-19,<br>&gt; Département de Physique, Faculté des Science de Monastir,<br>&gt; Avenue de l'Environnement 5019, Monastir Tunisie.<br>&gt; tél:+21697487042<br>&gt; <a
 ymailto="mailto:mourad_fsm@yahoo.fr" href="mailto:mourad_fsm@yahoo.fr">mourad_fsm@yahoo.fr</a><br><br>-- <br>********************************************************<br>- Article premier - Les hommes naissent et demeurent<br>libres et ègaux en droits. Les distinctions sociales<br>ne peuvent être fondèes que sur l'utilitè commune<br>- Article 2 - Le but de toute association politique<br>est la conservation des droits naturels et <br>imprescriptibles de l'homme. Ces droits sont la libertè,<br>la propriètè, la sùretè et la rèsistance à l'oppression.<br>********************************************************<br><br>&nbsp;  Giuseppe Mattioli&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <br>&nbsp;  CNR - ISTITUTO DI STRUTTURA DELLA MATERIA&nbsp;  <br>&nbsp;  v. Salaria Km 29,300 - C.P. 10&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <br>&nbsp;  I 00015 - Monterotondo Stazione (RM)&nbsp;
 &nbsp; &nbsp; &nbsp; &nbsp; <br>&nbsp;  Tel + 39 06 90672836 - Fax +39 06 90672316&nbsp; &nbsp; <br>&nbsp;  E-mail: &lt;<a ymailto="mailto:giuseppe.mattioli@ism.cnr.it" href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>&gt;<br><br><br> </div> </div>  </div></body></html>