<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"MS Gothic";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:"\@MS Gothic";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
span.EmailStyle19
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Dear Sydney,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">You may want to take a look at the Python-based Atomic Simulation Environment (ASE) (<a href="https://wiki.fysik.dtu.dk/ase/">https://wiki.fysik.dtu.dk/ase/</a>). You can build simple adsorbate systems using its python interface and directly
 generate a PWscf input file from a simple Python script.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Your expression for the binding energy is correct. You should perform vibrational calculations for the adsorbate bound to the surface and for the isolated molecule in vacuum to include at the very least zero point energy corrections. You
 can also use ASE to compute enthalpic and entropic corrections using its thermochemistry tools (<a href="https://wiki.fysik.dtu.dk/ase/ase/thermochemistry/thermochemistry.html">https://wiki.fysik.dtu.dk/ase/ase/thermochemistry/thermochemistry.html</a>) to
 generate finite temperature free energies (if this is what you are after).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">You should perform convergence tests to see what properties are important for your application, because it depends on what it is you are trying to calculate. You should try these calculations with and without van der Waals corrections,
 change the slab cross-sectional area (only if you are trying to eliminate lateral adsorbate-adsorbate interactions), and try several vacuum heights in addition to the usual convergence tests you should perform. Be aware that you may have a dipole across your
 supercell if your slab is not symmetric. You may need to include dipole corrections.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hope that helps.<o:p></o:p></p>
<p class="MsoNormal">Steve<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<p class="MsoNormal">Stephen Weitzner, PhD<o:p></o:p></p>
<p class="MsoNormal">Postdoctoral Research Scientist<o:p></o:p></p>
<p class="MsoNormal">Quantum Simulations Group<o:p></o:p></p>
<p class="MsoNormal">Lawrence Livermore National Laboratory<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">users <users-bounces@lists.quantum-espresso.org> on behalf of
</span><span style="font-size:12.0pt;font-family:"MS Gothic";color:black">劉芸瑄</span><span style="font-size:12.0pt;color:black"> <danosunny@gmail.com><br>
<b>Reply-To: </b>Quantum ESPRESSO users Forum <users@lists.quantum-espresso.org><br>
<b>Date: </b>Saturday, February 29, 2020 at 3:36 AM<br>
<b>To: </b>"users@lists.quantum-espresso.org" <users@lists.quantum-espresso.org><br>
<b>Subject: </b>[QE-users] how to study adsorption energy<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal" style="font-stretch:normal">Dear all <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:13.0pt">  I am a beginner in QE. Here are my question about adsorption energy. ( ex H2O adsorb on Cu surface )</span><o:p></o:p></p>
</div>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">  I know that I have to optimize the structure of the adsorbent molecule and the surface respectively. And then optimize the adsorbent + surface system.<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">Getting the adsorption energy by the equation<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">Eads= E(Ad+surface)-[(Ead)+(Esurface)] <o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt"><o:p> </o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">  But I am confused about the more specific step of how to use QE to study it.<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">First , I think I have to build the surface with the molecule above it. I try to use the BURAI to build the slab but I fail to build the molecule above it. Because
 I can not built the molecule with accurate bond length or angle which just like the Gaussview can do . <o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">  Alternatively, I hope to use VESTA to build the slab and then output the xyz file to Avogadro. And then build the more accurate molecule structure above the slab
 by Avogadro. After that, output the xyz file to BURAI to create the QE input file. Is that alright ? Or is there any software or method that is more easily to build the molecule + slab surface system ?  <o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">  Also, which part or parameter that I should be aware of ? I think there are still a lot to be concerned. Like the size of the slab, the distance between adsorbent
 and surface , the QE input file parameter, all of these could affect the convergence of the SCF and the accuracy of the energy. As a beginner, how to start studying the adsorption energy ?<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">  I will really appreciate that if you could kindly teach me the more specific instruction of adsorption energy calculation.<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">Thank you so much!<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt"><o:p> </o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">Best regards<o:p></o:p></span></p>
<p style="margin:0in;margin-bottom:.0001pt;font-stretch:normal"><span style="font-size:13.0pt">Sydney Liu<o:p></o:p></span></p>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<p class="MsoNormal">TOK TAIWAN CO.,LTD<br>
Sydney Liu<br>
E-mail : <a href="mailto:danosunny@gmail.com">danosunny@gmail.com</a><o:p></o:p></p>
</div>
</div>
</body>
</html>