<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-2022-jp">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Iurii, </p>
<p><br>
</p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;">Thank you for the reference links. I followed them and ran some tests. It works with U prediction. However, with CoPS3, I wonder whether even Hubbard's inclusion
is working. </span><br>
</p>
<p><br>
</p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><b>So there is some change. Is it not what you expect?</b></span><br>
</p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><b>Problem:
</b>Without Hubbard U inclusion lattice parameter comes </span><b style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;">{a ≠ sqrt(3)*b}, </b><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;">but
with the inclusion of U, it should be corrected to more like <b>{a = sqrt(3)*b}.
</b>According to the literature. </span><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont;"><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont;">The output structure does not hold the monoclinic symmetry; Hubbard U should help it achieve that. However, </span><span style="font-size: 12pt; font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont;">I
do not get it despite using any Hubbard U. </span><br>
</p>
<p><br>
</p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;">I also ran calculations with higher convergence criteria but ended up with similar lattice parameters. I am now running a calculation with ibrav= -12, but I am
wondering if it is the right approach in this system. </span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><br>
</span></p>
<p><span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;">Regards<span style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"><b style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px;"></b></span></span></p>
<p><font face="Calibri, Helvetica, sans-serif, serif, EmojiFont">Vivek</font></p>
<p><font face="Calibri, Helvetica, sans-serif, serif, EmojiFont">TU Dresden</font></p>
<p><font face="Calibri, Helvetica, sans-serif, serif, EmojiFont"><br>
</font></p>
</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> users <users-bounces@lists.quantum-espresso.org> on behalf of Iurii TIMROV via users <users@lists.quantum-espresso.org><br>
<b>Sent:</b> Friday, November 25, 2022 5:11:30 PM<br>
<b>To:</b> users@lists.quantum-espresso.org<br>
<b>Subject:</b> Re: [QE-users] Problem in vc-relax calculation with Hubbard U in magnetic system</font>
<div> </div>
</div>
<div><style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Vivek,</p>
<p><br>
</p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">> I want to predict the self-consistent Hubbard parameter using HP code.</span></font><br>
</p>
<p><br>
</p>
<p>Please have a look here:</p>
<p><a href="https://www.materialscloud.org/learn/sections/tctKdH/advanced-quantum-espresso-tutorial-2022" class="OWAAutoLink" id="LPlnk199888" previewremoved="true">https://www.materialscloud.org/learn/sections/tctKdH/advanced-quantum-espresso-tutorial-2022</a></p>
<p><br>
</p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">> ...but the structure remains the same.
<br>
</span></font></p>
<p><br>
</p>
<p>From your input and output files I see this:</p>
<p><br>
</p>
<p>OLD</p>
<p></p>
<div>CELL_PARAMETERS (angstrom)<br>
5.8969998360 0.0000000000 0.0000000000<br>
0.0000000000 10.2159996033 0.0000000000<br>
-1.9672647655 0.0000000000 6.3670060311</div>
<p></p>
<p><br>
</p>
<p>NEW</p>
<p></p>
<div>CELL_PARAMETERS (angstrom)<br>
5.771620672 0.000000000 0.036554788<br>
0.000000000 10.142807293 0.000000000<br>
-1.885622677 0.000000000 6.187226171</div>
<p></p>
<p><br>
</p>
<p>So there is some change. It is not what you expect?</p>
<p><br>
</p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">>
<span>conv_thr = 1.d-6</span></span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px"><br>
</span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">Try 1.d-10 or even smaller<br>
</span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px"><br>
</span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">>
<span>forc_conv_thr = 1.d-3</span><br>
</span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px"><br>
</span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px">Try
<span>forc_conv_thr = 1.d-4 or 1.d-5</span></span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px"><span>and</span></span></font></p>
<p><font size="2" face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols"><span style="font-size:16px"><span><span>etot_conv_thr</span> = 1.d-5 or 1.d-6<br>
</span></span></font></p>
<p><br>
</p>
<p>Have a look also here:</p>
<p><a href="https://www.materialscloud.org/work/tools/qeinputgenerator" class="OWAAutoLink" id="LPlnk903279" previewremoved="true">https://www.materialscloud.org/work/tools/qeinputgenerator</a></p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</p>
<p><span style="font-family: Constantia, "Hoefler Text", serif;">Iurii</span><br>
</p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, "EmojiFont", "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
<span style="font-family:Cambria,Georgia,serif">Dr. Iurii TIMROV</span><br>
<span style="font-family:Cambria,Georgia,serif">Senior Research Scientist</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">Theory and Simulation of Materials (THEOS)</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-size:12pt; color:rgb(128,128,128); font-family:Cambria,Georgia,serif">Swiss Federal Institute of Technology Lausanne (EPFL</span><font style="font-size:12pt; color:rgb(128,128,128); font-family:"Times New Roman",Times,serif" color="808080"><font face="'Times New Roman', Times, serif"><span style="font-family:Cambria,Georgia,serif">)</span></font></font><br>
</div>
<span style="font-family:Cambria,Georgia,serif"></span><font color="808080"></font><span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">CH-1015 Lausanne, Switzerland</span><br>
<span style="font-family:Cambria,Georgia,serif">+41 21 69 34 881</span></font></div>
<span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><a href="http://people.epfl.ch/265334" id="LPNoLP" tabindex="0"><span style="font-family:Cambria,Georgia,serif">http://people.epfl.ch/265334</span></a><br>
</div>
</div>
</div>
</div>
</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> users <users-bounces@lists.quantum-espresso.org> on behalf of Vivekanand Shukla <vivekanand.shukla@chalmers.se><br>
<b>Sent:</b> Friday, November 25, 2022 4:34:20 PM<br>
<b>To:</b> users@lists.quantum-espresso.org<br>
<b>Subject:</b> [QE-users] Problem in vc-relax calculation with Hubbard U in magnetic system</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Dear all, </span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;"><br>
</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">I am trying to run a <b>vc-relax</b> calculation for CoPS3 in bulk with the
Hubbard U parameter. CoPS3 results in off-lattice parameters without Hubbard U. I checked it with VASP, and it seems realistic. However, in QE, the Hubbard U parameter has no effect on the lattice parameter. It remains off with or without U. I can see the
Hubbard effect in forces in the output file, but it does not reflect on the lattice parameter in QE output.</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;"><br>
</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">I want to predict the self-consistent Hubbard parameter using HP code. It
does not make sense if I do not get the inclusion of the U parameter right in the first place. I used <b>QE versions 6.4.1, 6.7, and 7.0</b> and various convergence strategies, but the structure remains the same. </span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;"><br>
</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, serif, EmojiFont; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Could anyone please suggest anything wrong with the calculation? Please find
one of the input and output files in the attached link. </span></font></div>
<p></p>
<p><br>
</p>
<p><a href="https://cloudstore.zih.tu-dresden.de/index.php/s/8cCf3N8qd3swW5d" class="OWAAutoLink" id="LPlnk326444" previewremoved="true">https://cloudstore.zih.tu-dresden.de/index.php/s/8cCf3N8qd3swW5d</a></p>
<div id="LPBorder_GT_16693904047850.23621807749044854" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_16693904047820.789985927596291" role="presentation" cellspacing="0" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);">
<tbody>
<tr valign="top" style="border-spacing: 0px;">
<td id="ImageCell_16693904047830.7133422280175312" colspan="1" style="width: 250px; position: relative; display: table-cell; padding-right: 20px;">
<div id="LPImageContainer_16693904047830.06500177644709404" style="background-color: rgb(255, 255, 255); height: 200px; position: relative; margin: auto; display: table; width: 200px;">
<a id="LPImageAnchor_16693904047830.4778745008730201" href="https://cloudstore.zih.tu-dresden.de/index.php/s/8cCf3N8qd3swW5d" target="_blank" style="display: table-cell; text-align: center;"><img aria-label="Preview image with link selected. Double-tap to open the link." id="LPThumbnailImageID_16693904047830.3217748667410407" width="200" height="200" style="display: inline-block; max-width: 250px; max-height: 250px; height: 200px; width: 200px; border-width: 0px; vertical-align: bottom;" src="https://cloudstore.zih.tu-dresden.de/core/img/favicon-fb.png"></a></div>
</td>
<td id="TextCell_16693904047830.9653261746673563" colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;">
<div id="LPRemovePreviewContainer_16693904047830.5335458730522691"></div>
<div id="LPTitle_16693904047830.8349895077318461" style="top: 0px; color: rgb(0, 120, 215); font-weight: 400; font-size: 21px; font-family: wf_segoe-ui_light, "Segoe UI Light", "Segoe WP Light", "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; line-height: 21px;">
<a id="LPUrlAnchor_16693904047840.7668073479167283" href="https://cloudstore.zih.tu-dresden.de/index.php/s/8cCf3N8qd3swW5d" target="_blank" style="text-decoration: none;">Cloudstore</a></div>
<div id="LPMetadata_16693904047840.699146346872868" style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: 400; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 14px;">
cloudstore.zih.tu-dresden.de</div>
<div id="LPDescription_16693904047840.3291956348988161" style="display: block; color: rgb(102, 102, 102); font-weight: 400; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;">
Cloudstore - based on <a href="https://nextcloud.com/\">Nextcloud</a></div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<div style="font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Sincerely</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Vivek</span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Postdoctoral Researcher </span></font></div>
<div style="font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px; margin-top: 0px; margin-bottom: 0px;">
<font face="Calibri,Helvetica,sans-serif,EmojiFont,Apple Color Emoji,Segoe UI Emoji,NotoColorEmoji,Segoe UI Symbol,Android Emoji,EmojiSymbols" size="2"><span style="font-size: 16px;">Technical University Dresden </span></font></div>
<p></p>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p><br>
</p>
</div>
</div>
<style type="text/css" style="display:none"> <!-- p {margin-top:0; margin-bottom:0} --> </style></div>
</div>
</div>
</body>
</html>