Dear <span style="font-size: 14.6667px;">Matthias,<br /></span><br />Thanks a lot for providing the information.<br /><br />Could you please provide your input file for CP2K for reference? Maybe I made a mistake somewhere.<br /><br />Best regards,<br />Kejiang<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, November 21, 2023 at 11:29:14 PM UTC+8 Krack Matthias wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word">
<div>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt">Hi
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">The EOS curve for FeO (2x2x2) looks fine for me (see attached plot) using PBE and the
<a href="https://pubs.acs.org/doi/10.1021/acs.est.7b01670" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://pubs.acs.org/doi/10.1021/acs.est.7b01670&source=gmail&ust=1700699500156000&usg=AOvVaw1fXBBD3WeZbcuI0ID7nmyy">suggested Hubbard U(eff) value of 1.9 eV for iron</a> with CP2K.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Best<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Matthias<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><img border="0" width="800" height="600" style="width:8.3333in;min-height:6.25in" src="https://groups.google.com/group/cp2k/attach/a252141c4ca/image001.png?part=0.1&view=1" alt="A graph of a function
Description automatically generated"></span><span lang="EN-US" style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<div>
<div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Jürg Hutter <<a href data-email-masked rel="nofollow">hut...@chem.uzh.ch</a>><br>
<b>Date: </b>Tuesday, 21 November 2023 at 09:46<br>
<b>To: </b><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>Re: [CP2K:19537] Re: Significant discontinuity in EOS curver for FeO<u></u><u></u></span></p>
</div>
</div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Hi<br>
<br>
what is the setup for the QE calculation (cell, k-points, magnetization)?<br>
<br>
regards<br>
JH<br>
<br>
________________________________________<br>
From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Kejiang Li <<a href data-email-masked rel="nofollow">lyam_...@126.com</a>><br>
Sent: Monday, November 20, 2023 2:26 AM<br>
To: cp2k<br>
Subject: [CP2K:19531] Re: Significant discontinuity in EOS curver for FeO<br>
<br>
Dear all,<br>
<br>
To add one more point.<br>
As it was said, using a reference cell might solve this problem (<a href="https://groups.google.com/g/cp2k/c/N52jLt2yAIQ/m/JAlV3KTVUtIJ" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/g/cp2k/c/N52jLt2yAIQ/m/JAlV3KTVUtIJ&source=gmail&ust=1700699500156000&usg=AOvVaw0e8FTY2dyhQ4IywiTRyPmV">https://groups.google.com/g/cp2k/c/N52jLt2yAIQ/m/JAlV3KTVUtIJ</a>). You will notice that I used &CELL_REF in my input script,
which did not solve my problem.<br>
<br>
Thanks.<br>
Kejiang<br>
<br>
On Monday, November 20, 2023 at 9:00:11 AM UTC+8 Kejiang Li wrote:<br>
Dear CP2K community,<br>
<br>
I am calculating the EOS for FeO with a 2*2*2 supercell. But there is always a significant discontinuity (sharp decrease) in the range of 4.3-4.4 angstrom lattice constant in the EOS curves, as shown in the figures below.<br>
<br>
I have done the following tests:<br>
- the convergence of the energy with different cutoffs;<br>
- using different functionals, including PBE, HSE, and SCAN;<br>
-including HSE with different HF fractions;<br>
-using different magnetization and +U settings.<br>
However, all the tests produce similar results with an apparent discontinuity (4.3-4.3 A) .<br>
<br>
This discontinuity is not observed in the Quantum Espresso results conducted by myself and in the results of literature that generally used VASP. I believe that the CP2K program might cause this problem.<br>
<br>
Can you give me more suggestions on how to do further tests?<br>
<br>
Here, I also attached one sample code I used during my tests.<br>
<br>
Thanks a lot.<br>
<br>
Best regards,<br>
Kejiang<br>
The University of Science and Technology Beijing<br>
<br>
[cp2k_functionals.jpg]<br>
[cp2k_mag_u.jpg][cp2k_mag_u2.jpg]<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%40googlegroups.com%3chttps:/groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%40googlegroups.com?utm_medium=email&utm_source=footer%3e" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%2540googlegroups.com%253chttps:/groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter%253e&source=gmail&ust=1700699500156000&usg=AOvVaw1EzD65s1UyLGUTpseXOMCX">
https://groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/f6ca6f04-88e1-47cf-aa9f-342a276d7776n%40googlegroups.com?utm_medium=email&utm_source=footer></a>.<br>
<br>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br></span></p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">
To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759D67572D6CC24912369949FBBA%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759D67572D6CC24912369949FBBA%2540ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM&source=gmail&ust=1700699500156000&usg=AOvVaw1g1NO7rhWnKmfGmUwlkXLN">
https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759D67572D6CC24912369949FBBA%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM</a>.<u></u><u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</blockquote></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0182b11e-b2ff-414a-a655-9b4c697e2910n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/0182b11e-b2ff-414a-a655-9b4c697e2910n%40googlegroups.com</a>.<br />