<!DOCTYPE html>
<html data-lt-installed="true">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body style="padding-bottom: 1px;">
<p>I would recommend reading the literature instead of asking a
large-language model, but I admit there is a real risk of learning
something useful beyond your query.</p>
<p>Studying the literature to have an idea of the actual magnitude
of the effect, how people have computed it in the past, and if it
may impact your case will save you a lot of time in the long run.<br>
</p>
<p>In any case, whatever Gemini has barfed is not wrong, but very
generic and extremely unpractical for a system with more than a
handful of atoms. If there is a well-separated vibrational mode
associated with the molecule-surface interface, it is possible to
find its energy via the phonon code (option nat_todo), or just
with pw.x using frozen phonons.</p>
hth<br>
<p></p>
<div class="moz-cite-prefix">On 12/08/2025 15:14, Kazume NISHIDATE
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:F35BC4BB-EB30-4041-AF25-56F15B0EB924@iwate-u.ac.jp">
<pre wrap="" class="moz-quote-pre">And I asked it to Gemini, Google’s AI.
Here is her response 😲
=========================
Yes, it's possible to include the zero-point energy (ZPE) correction for both a 2D surface and a 2D surface with adsorbed atoms using Quantum ESPRESSO. This correction is crucial for accurately calculating energies, such as adsorption energies, as it accounts for the vibrational energy of the atoms even at absolute zero temperature (0 K).
To calculate the ZPE correction, you must perform phonon calculations for your systems. Quantum ESPRESSO's phonon code, ph.x, is the primary tool for this. The ZPE is half the sum of all vibrational frequencies, represented by the following formula:
ZPE=21i∑ℏωi
Where:
• ℏ is the reduced Planck constant.
• ωi are the vibrational frequencies (eigenfrequencies) of the system's normal modes.
1. Zero-Point Energy Correction for a 2D Surface
To calculate the ZPE for a 2D surface, you'll follow these general steps:
• Optimize the geometry of your 2D surface slab to find the lowest-energy configuration.
• Run a self-consistent field (SCF) calculation on the optimized structure. This calculation is a prerequisite for the phonon calculation.
• Perform a phonon calculation using ph.x. This involves calculating the vibrational modes across a q-point grid in the Brillouin zone. The output will be the phonon frequencies.
• Process the phonon output to get the total ZPE. This often involves using a post-processing tool like matdyn.x to diagonalize the dynamical matrices and sum up the positive frequencies. The ZPE is then calculated from this sum.
For a 2D surface, it's important to use a supercell that is large enough to minimize the interactions between periodic images and to ensure a sufficient vacuum layer in the z-direction.
2. ZPE Correction for an Adsorbed System
For a system with an adsorbed molecule on a 2D surface, the process is similar but with an important distinction: you need to calculate the ZPE correction for three separate systems to determine the corrected adsorption energy.
The ZPE-corrected adsorption energy (ΔEadsZPE) is calculated as follows:
ΔEadsZPE=(Eslab+adsorbateDFT+ZPEslab+adsorbate)−(EslabDFT+ZPEslab)−(EadsorbateDFT+ZPEadsorbate)
Where:
• EDFT is the ground-state electronic energy from a standard DFT calculation (e.g., using pw.x).
• ZPE is the zero-point energy of each system.
To get the ZPEs, you'll perform a phonon calculation for each of the three components:
• Adsorbed System (Slab + Adsorbate):
• Optimize the geometry of the combined system.
• Run an SCF calculation.
• Run a phonon calculation using ph.x. To save computational time, you can often freeze the atoms of the bulk surface layers and only allow the adsorbate atoms and the top surface layer atoms to move. This is a common approximation.
• Calculate the ZPE from the phonon frequencies.
• Isolated Surface (Slab):
• Optimize the geometry of the clean surface slab.
• Run an SCF calculation.
• Run a phonon calculation. Again, you can freeze the bottom layers to reduce the computational cost.
• Calculate the ZPE.
• Isolated Adsorbate (Molecule):
• Optimize the geometry of the isolated molecule in a large vacuum cell.
• Run an SCF calculation.
• Run a phonon calculation. For a molecule, all atoms are allowed to move.
• Calculate the ZPE from the frequencies.
After obtaining the ZPE values for all three systems, you can substitute them into the equation to get the final ZPE-corrected adsorption energy.
</pre>
<blockquote type="cite">
<pre wrap="" class="moz-quote-pre">2025/08/12 17:28、Chiara Cignarella via users <a class="moz-txt-link-rfc2396E" href="mailto:users@lists.quantum-espresso.org"><users@lists.quantum-espresso.org></a>のメール:
Dear Ashley,
You can use the frequencies obtained by phonon calculations to compute the zero-point energy using the formula (you can find that in textbooks on the quantum harmonic oscillator part - for example Ashcroft pag. 451- ), i.e. doing the sum \sum_i 1/2 * \omega_i * hbar (divided by appropriate normalization factor for your sum)
Note that the computed phonons do not contain themselves the ZPM contributions!
Hope this is useful.
Best regards,
Chiara Cignarella
-------------------------------------
Post-doc researcher
University of Bremen
</pre>
<blockquote type="cite">
<pre wrap="" class="moz-quote-pre">On 12 Aug 2025, at 09:51, Barsha Pal <a class="moz-txt-link-rfc2396E" href="mailto:barsha.pal@aus.ac.in"><barsha.pal@aus.ac.in></a> wrote:
Hello,
I have a small doubt regarding whether or not zero-point energy calculation is possible in Quantum Espresso, I want to perform that for a 2D material and adsorption cases. The property that I aim, includes the zero-point energy correction for the 2D surface and the 2D surface with adsorbed atoms. AI tells me that it is possible to do by phonon calculation. If indeed this is possible, can you guide me towards a short example?
Thank you.
Ashley Cooper
Phd
_______________________________________________________________________________
The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.
--------------------------------------------------------------------------------
Quantum ESPRESSO is supported by MaX (<a class="moz-txt-link-abbreviated" href="http://www.max-centre.eu">www.max-centre.eu</a>)
users mailing list <a class="moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<pre wrap="" class="moz-quote-pre">
_______________________________________________________________________________
The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.
--------------------------------------------------------------------------------
Quantum ESPRESSO is supported by MaX (<a class="moz-txt-link-abbreviated" href="http://www.max-centre.eu">www.max-centre.eu</a>)
users mailing list <a class="moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<pre wrap="" class="moz-quote-pre">
西館数芽
Kazume NISHIDATE Ph.D
Department of Systems Innovation Engineering,
Graduate School of Science and Engineering, Iwate University
4-3-5 Ueda, Morioka, Iwate, 020-8551 JAPAN
Phone:+81-19-621-6391
<a class="moz-txt-link-abbreviated" href="mailto:kazume.nishidate@gmail.com">kazume.nishidate@gmail.com</a>, <a class="moz-txt-link-abbreviated" href="mailto:nisidate@iwate-u.ac.jp">nisidate@iwate-u.ac.jp</a>
<a class="moz-txt-link-freetext" href="https://sites.google.com/site/nisidatelab/">https://sites.google.com/site/nisidatelab/</a>
_______________________________________________________________________________
The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.
--------------------------------------------------------------------------------
Quantum ESPRESSO is supported by MaX (<a class="moz-txt-link-abbreviated" href="http://www.max-centre.eu">www.max-centre.eu</a>)
users mailing list <a class="moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></pre>
</blockquote>
<div class="moz-signature">-- <br>
<small>Dr. Lorenzo Paulatto<br>
IR IMPMC - CNRS UMR 7590 / Sorbonne Université / MNHN <br>
phone: +33 (0)1 442 79822 / telegram: lpaulatto<br>
<a href="http://www.impmc.upmc.fr/~paulatto/"
class="moz-txt-link-freetext">http://www.impmc.upmc.fr/~paulatto/</a>
- <a href="https://anharmonic.github.io/"
class="moz-txt-link-freetext">https://anharmonic.github.io/</a><br>
23-24/423 B115, 4 place Jussieu 75252 Paris CX 05<small></small></small></div>
</body>
<lt-container></lt-container>
</html>