Dear Sydney, You may want to take a look at the Python-based Atomic Simulation Environment (ASE) (https://wiki.fysik.dtu.dk/ase/). You can build simple adsorbate systems using its python interface and directly generate a PWscf input file from a simple Python script.
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 (https://wiki.fysik.dtu.dk/ase/ase/thermochemistry/thermochemistry.html) to generate finite temperature free energies (if this is what you are after). 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. Hope that helps. Steve -- Stephen Weitzner, PhD Postdoctoral Research Scientist Quantum Simulations Group Lawrence Livermore National Laboratory From: users <[email protected]> on behalf of 劉芸瑄 <[email protected]> Reply-To: Quantum ESPRESSO users Forum <[email protected]> Date: Saturday, February 29, 2020 at 3:36 AM To: "[email protected]" <[email protected]> Subject: [QE-users] how to study adsorption energy Dear all I am a beginner in QE. Here are my question about adsorption energy. ( ex H2O adsorb on Cu surface ) I know that I have to optimize the structure of the adsorbent molecule and the surface respectively. And then optimize the adsorbent + surface system. Getting the adsorption energy by the equation Eads= E(Ad+surface)-[(Ead)+(Esurface)] But I am confused about the more specific step of how to use QE to study it. 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 . 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 ? 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 ? I will really appreciate that if you could kindly teach me the more specific instruction of adsorption energy calculation. Thank you so much! Best regards Sydney Liu -- TOK TAIWAN CO.,LTD Sydney Liu E-mail : [email protected]<mailto:[email protected]>
_______________________________________________ Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso) users mailing list [email protected] https://lists.quantum-espresso.org/mailman/listinfo/users
