<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0">
Hi Aaron,</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2 ContentPasted3">
It seems that SCAN0 in the current `develop` branch is not recognized early enough as a hybrid functional to trigger the needed the exx initialization. A simple fix can be found in the attached diff file.
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2">
However, I am not sure if it is worth checking into the official repository since Fabrizio might have a much better solution to be released soon.
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2">
<br>
</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2">
Best,</div>
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted2">
Hsin-Yu<br>
</div>
<div class="elementToProof">
<div style="font-family: "Times New Roman", Times, serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<font color="#888888">--
<div dir="ltr">
<div dir="ltr">
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Hsin-Yu Ko</span></span></div>
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Postdoctoral Research Fellow<br>
</span></span></div>
<span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Department of Chemistry and Chemical Biology<br>
</span></span><span><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Cornell University</span></span></span></div>
</div>
</font><br>
</div>
</div>
</div>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Aaron Kaplan <kaplan@temple.edu><br>
<b>Sent:</b> Thursday, February 2, 2023 3:51 PM<br>
<b>To:</b> Hsin-Yu Ko <hsin-yu.ko@cornell.edu>; General discussion list for Quantum ESPRESSO developers <developers@lists.quantum-espresso.org><br>
<b>Cc:</b> Paolo GIANNOZZI <paolo.giannozzi@uniud.it><br>
<b>Subject:</b> Re: [External] Re: [QE-developers] [External] Re: Hybrid meta-GGAs in QE</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Dear Dr. Ko,</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Running SCAN0 with PWscf also results in a segfault - the EXX calculation is initialized, but the calculation immediately stops afterwards.</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
For example, using ACE for EXX and setting input_dft = SCA0 (same as XC-000L-000L-000L-000L-264L-267L) prints the following lines in pw.x ' s output:</div>
<blockquote itemscope="" itemtype="https://schemas.microsoft.com/QuotedText" style="border-color:rgb(200,200,200); border-left:3px solid rgb(200,200,200); padding-left:1ex; margin-left:0.8ex">
<div class="x_elementToProof x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
convergence has been achieved in  42 iterations
<div class="x_ContentPasted0">     Writing meta-gga kinetic term</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">     Using ACE for calculation of exact exchange</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">     EXX grid:    31103 G-vectors     FFT dimensions: (  40,  40,  40)</div>
     ACE projected onto    58 (nbndproj) and applied to    58 (nbnd) bands<br>
</div>
</blockquote>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
After which the segfault is thrown (guessing that at least one variable needed for the EXX calculation isn't allocated for meta-GGAs).<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Best,</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Aaron<br>
</div>
<div id="x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Hsin-Yu Ko <hsin-yu.ko@cornell.edu><br>
<b>Sent:</b> Thursday, February 2, 2023 9:23 AM<br>
<b>To:</b> General discussion list for Quantum ESPRESSO developers <developers@lists.quantum-espresso.org><br>
<b>Cc:</b> Aaron Kaplan <kaplan@temple.edu>; Paolo GIANNOZZI <paolo.giannozzi@uniud.it><br>
<b>Subject:</b> [External] Re: [QE-developers] [External] Re: Hybrid meta-GGAs in QE</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_x_elementToProof x_x_ContentPasted0" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Hi Aaron,</div>
<div class="x_x_elementToProof x_x_ContentPasted0" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_x_elementToProof x_x_ContentPasted0 x_x_ContentPasted1 x_x_ContentPasted2" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
While there is not a systematic implementation of hybrid meta-GGA from libxc, as Paolo pointed out, I think there is an instance of SCAN0 (25% EXX mixing with SCAN meta-GGA) implemented (<a href="https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgitlab.com%2FQEF%2Fq-e%2F-%2Fblob%2Fdevelop%2FXClib%2Fdft_setting_routines.f90%23L166&data=05%7C01%7Ckaplan%40temple.edu%7C5082fea1d5d1417dfdfd08db05290c54%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109446082237004%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=B2LtR6Jaqpi8lyQJnr95a8cvwNI6Oq4JlUdDAkk5hHE%3D&reserved=0" originalsrc="https://gitlab.com/QEF/q-e/-/blob/develop/XClib/dft_setting_routines.f90#L166" shash="X2u+EJcj6R/sDYR1EvnrfmnhBDG+puhOWBdL90PqY6icQJlx/QA2gojCG0ZVTkQmQzpHG1LX5BwIqiAiL8Vj17q6uFB1uP5O/MqgqJj7S9RiSdeigBTQdjlw4LFQh52fGKQJ5yyNg53Z2UTgg+6ut/52+jxMKCjhFONeAeWujg4=" id="LPNoLPOWALinkPreview_1">https://gitlab.com/QEF/q-e/-/blob/develop/XClib/dft_setting_routines.f90#L166</a>).
 You can probably follow a similar procedure to introduce the needed meta-GGA part from libxc.<br>
</div>
<div class="x_x__Entity x_x__EType_OWALinkPreview x_x__EId_OWALinkPreview x_x__EReadonly_1">
</div>
<br>
<div class="x_x_elementToProof">
<div class="x_x_elementToProof" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
Best,</div>
<div class="x_x_elementToProof" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
Hsin-Yu<br>
</div>
<div class="x_x_elementToProof" style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div id="x_x_Signature">
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<font color="#888888">--
<div dir="ltr">
<div dir="ltr">
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Hsin-Yu Ko</span></span></div>
<div><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Postdoctoral Research Fellow<br>
</span></span></div>
<span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Department of Chemistry and Chemical Biology<br>
</span></span><span><span style="font-family:georgia,serif"><span style="color:rgb(102,102,102)">Cornell University</span></span></span></div>
</div>
</font><br>
</div>
</div>
</div>
</div>
<div id="x_x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> developers <developers-bounces@lists.quantum-espresso.org> on behalf of Aaron Kaplan <kaplan@temple.edu><br>
<b>Sent:</b> Thursday, February 2, 2023 8:55 AM<br>
<b>To:</b> Paolo Giannozzi <paolo.giannozzi@uniud.it>; General discussion list for Quantum ESPRESSO developers <developers@lists.quantum-espresso.org><br>
<b>Subject:</b> Re: [QE-developers] [External] Re: Hybrid meta-GGAs in QE</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Thank you, Dr. Giannozzi.</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Agreed that meta-GGAs + unscreened-EXX may be unusable in a PW code, but I'm tasked with determining that.
<br>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
I will try to follow the implementation of hybrid-GGAs and see if it's feasible to add the functionality for hybrid meta-GGAs.
<br>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Best,</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Aaron<br>
</div>
<div id="x_x_x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_x_x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Paolo Giannozzi <paolo.giannozzi@uniud.it><br>
<b>Sent:</b> Thursday, February 2, 2023 4:17 AM<br>
<b>To:</b> General discussion list for Quantum ESPRESSO developers <developers@lists.quantum-espresso.org><br>
<b>Cc:</b> Aaron Kaplan <kaplan@temple.edu><br>
<b>Subject:</b> [External] Re: [QE-developers] Hybrid meta-GGAs in QE</font>
<div> </div>
</div>
<div class="x_x_x_BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="x_x_x_PlainText">Hybrid meta-GGA functionals are not currently implemented. I don't think
<br>
it is difficult, but: hybrid functionals are computationally heavy, <br>
meta-GGAs are numerically unstable, the combination of the two may turn <br>
out to be unusable<br>
<br>
Paolo<br>
<br>
On 2/1/23 22:17, Aaron Kaplan wrote:<br>
> <br>
> Hello,<br>
> <br>
> I'm a former student in John Perdew's group, currently working with Mike <br>
> Klein, and one of the r2 SCAN developers.<br>
> <br>
> I'm curious if it's possible to use hybrid meta-GGAs with PW in QE 7.1 - <br>
> I've tried many variations of this calculation but a segfault is always <br>
> thrown when the ACE EXX calculation starts.<br>
> <br>
> I saw in the mail archives that PW in QE 6.* could not perform hybrid <br>
> meta-GGA calculations <br>
> (<a href=""></a>https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fpipermail%2Fusers%2F2020-February%2F044020.html&data=05%7C01%7Ckaplan%40temple.edu%7C15aba46157ac483cf47308db04fe55cf%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109262658111891%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=5TiSFsCKjOl5xoQo%2BZBRqmqiAk1Skbg85nqWXZxb%2Bzo%3D&reserved=0
 <<a href="https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fpipermail%2Fusers%2F2020-February%2F044020.html&data=05%7C01%7Ckaplan%40temple.edu%7C5082fea1d5d1417dfdfd08db05290c54%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109446082237004%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=Z8HkFisRFHA6Ei53Y%2BpxkP6001mZyolS0MfYW0jmWj4%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/pipermail/users/2020-February/044020.html" shash="kj2u3OPLkXhtXOa3VtuhSFnVbDNYJ3I08YHKKQL4DOlHb5K2/QPWLVxG/5YXLzbpviOR9CcWb345s1rLAVSys4G26Kj7u56s618eobwNWYHVSaZNUzjkTD4LUj7Y++fPH1md5iwB+Aw1jgRT6aSixn8I9RN/dWD/+CFL7uN4pUw=">https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fpipermail%2Fusers%2F2020-February%2F044020.html&data=05%7C01%7Ckaplan%40temple.edu%7C15aba46157ac483cf47308db04fe55cf%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109262658111891%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=5TiSFsCKjOl5xoQo%2BZBRqmqiAk1Skbg85nqWXZxb%2Bzo%3D&reserved=0</a>>),
 but am not sure if that still applies to QE 7.<br>
> <br>
> I'll note a few other issues: it seems like LibXC's hybrid GGAs don't <br>
> work either (e.g., the global hybrid of PBEsol, ID = 274), and newer <br>
> versions of LibXC cannot be used because they removed identifiers like <br>
> HYB_GGA_* and HYB_MGGA_* (e.g., in LibXC v. 6 and on, a hybrid GGA is <br>
> identified as a GGA).<br>
> <br>
> Best,<br>
> Aaron Kaplan<br>
> <br>
> ________________________________________________<br>
> The Quantum ESPRESSO community stands by the Ukrainian people<br>
>   and expresses its concerns about the devastating effects that<br>
> the Russian military offensive has on their country and on the<br>
> free and peaceful scientific, cultural, and economic cooperation<br>
> amongst peoples.<br>
> _______________________________________________<br>
> developers mailing list<br>
> developers@lists.quantum-espresso.org<br>
> <a href="https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fdevelopers&data=05%7C01%7Ckaplan%40temple.edu%7C5082fea1d5d1417dfdfd08db05290c54%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109446082237004%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=q0ZWprC6r2np6wT255sXTbLqbw%2BC9EYyIo8285Istm0%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/mailman/listinfo/developers" shash="LMkWjlT2torTYGoWHdyyeGI5r6GE4USn//l2ncbV6KIktgoZLjspG3fROdPdwBuICMFm5SC9AYWJqdD+7vSxVp1j58lF+OOFd0hvu0HFt9S/LSlVliRemyhS6I24+gebTFuYpO0U3ellku3dGIC0yynhSGi/tBTZGm2aZrtJWjk=">
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fdevelopers&data=05%7C01%7Ckaplan%40temple.edu%7C15aba46157ac483cf47308db04fe55cf%7C716e81efb52244738e3110bd02ccf6e5%7C0%7C0%7C638109262658111891%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=%2B7BSyzGDTsUHs9kykpGCctwmfRhDljGpsjp6e9d%2FYLo%3D&reserved=0</a><br>
<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
</div>
</span></font></div>
</div>
</div>
</div>
</body>
</html>