<style class="ke-style">
[list-style-type] {padding-left:20px;list-style-position:inside}
[list-style-type] li {margin:0}
[list-style-type] li:before, span.ke-list-item-matter {font-family:"sans serif",tahoma,verdana,helvetica}
[list-style-type] li p,[list-style-type] li h1,[list-style-type] li h2,[list-style-type] li h3,[list-style-type] li h4,[list-style-type] li h5,[list-style-type] li div,[list-style-type] li blockquote{display:inline;word-break:break-all}
[list-style-type] li table {display:inline-block;vertical-align:top}
p{margin:0}
td {word-break: break-word}
.default-font-1717765926100{
font-size:16px;
}
</style><div class="default-font-1717765926100" dir="ltr"><p>Dear Developers:</p><p><br></p><p>Thanks for your reply to my previous inquiry.</p><p><br></p><p>After encountering this issue, I manually keep tracking on the WF center to restore the continuity of the indices. However, I am still curious if there is a more permanent solution to address this matter. Because when I review historical data generated by VASP v5.4.1 combining with Wannier90 v2.1.0, I observed that the WF center could remain close to the atom center even when the fraction coordinates exceed 0.5. For instance, the atomic position section in the .win file is:</p><p><br></p><p>Begin Projections</p><p>Cd:sp3:z=-1,-1,-1,x=1,0,-1</p><p>Se:p:z=-1,-1,-1,x=1,0,-1</p><p>End Projections</p><p>...</p><p>begin unit_cell_cart</p><p> 6.0500002 0.0000000 0.0000000</p><p> 0.0000000 6.0500002 0.0000000</p><p> 0.0000000 0.0000000 6.0500002</p><p>end unit_cell_cart</p><p><br></p><p>begin atoms_cart</p><p><br></p><p>Cd 1.5125000 1.5125000 1.5125000</p><p>Cd 4.5375001 4.5375001 1.5125000</p><p>Cd 4.5375001 1.5125000 4.5375001</p><p>Cd 1.5125000 4.5375001 4.5375001</p><p>Se 0.0000000 0.0000000 0.0000000</p><p>Se 0.0000000 3.0250001 3.0250001</p><p>Se 3.0250001 0.0000000 3.0250001</p><p>Se 3.0250001 3.0250001 0.0000000</p><p>end atoms_cart</p><div>...</div><p><br></p><p>And the .wout file indicates: </p><p><br></p><p>...</p><p>Final State</p><p> WF centre and spread 1 ( 2.073500, 0.951423, 0.951533 ) 2.87731075</p><p> WF centre and spread 2 ( 2.074079, 2.074121, 2.074076 ) 2.87774052</p><p> WF centre and spread 3 ( 0.951538, 0.951427, 2.073501 ) 2.87730840</p><p> WF centre and spread 4 ( 0.950970, 2.073917, 0.950972 ) 2.87765217</p><p> WF centre and spread 5 ( 5.098500, 3.976423, 0.951533 ) 2.87731073</p><p> WF centre and spread 6 ( 5.099079, 5.099121, 2.074076 ) 2.87774049</p><p> WF centre and spread 7 ( 3.976538, 3.976427, 2.073501 ) 2.87730839</p><p> WF centre and spread 8 ( 3.975970, 5.098917, 0.950972 ) 2.87765219</p><p> WF centre and spread 9 ( 5.098500, 0.951423, 3.976534 ) 2.87731074</p><p> WF centre and spread 10 ( 5.099079, 2.074121, 5.099076 ) 2.87774051</p><p> WF centre and spread 11 ( 3.976538, 0.951427, 5.098501 ) 2.87730836</p><p> WF centre and spread 12 ( 3.975970, 2.073917, 3.975972 ) 2.87765218</p><p> WF centre and spread 13 ( 2.073500, 3.976423, 3.976534 ) 2.87731076</p><p> WF centre and spread 14 ( 2.074079, 5.099121, 5.099076 ) 2.87774049</p><p> WF centre and spread 15 ( 0.951538, 3.976427, 5.098501 ) 2.87730838</p><p> WF centre and spread 16 ( 0.950970, 5.098917, 3.975972 ) 2.87765217</p><p> WF centre and spread 17 ( -0.142515, -0.115089, -0.142653 ) 1.86289944</p><p> WF centre and spread 18 ( -0.000300, 0.202594, 0.000294 ) 1.87615749</p><p> WF centre and spread 19 ( 0.142692, -0.088809, 0.142242 ) 1.86859312</p><p> WF centre and spread 20 ( -0.142515, 2.909912, 2.882347 ) 1.86289944</p><p> WF centre and spread 21 ( -0.000300, 3.227594, 3.025294 ) 1.87615749</p><p> WF centre and spread 22 ( 0.142692, 2.936191, 3.167242 ) 1.86859313</p><p> WF centre and spread 23 ( 2.882485, -0.115089, 2.882347 ) 1.86289944</p><p> WF centre and spread 24 ( 3.024700, 0.202594, 3.025294 ) 1.87615749</p><p> WF centre and spread 25 ( 3.167693, -0.088809, 3.167242 ) 1.86859313</p><p> WF centre and spread 26 ( 2.882485, 2.909912, -0.142653 ) 1.86289944</p><p> WF centre and spread 27 ( 3.024700, 3.227594, 0.000294 ) 1.87615749</p><p> WF centre and spread 28 ( 3.167693, 2.936191, 0.142242 ) 1.86859312</p><p> Sum of centres and spreads ( 66.549857, 66.548337, 66.549865 ) 68.47064743</p><p><br></p><p>Currently, the center of WF5 is positioned at (5.10, 3.98, 0.95) within the home cell. Given that the previous version maintained the WF center within the home cell and ensured continuity, I am hopeful that the latest version can achieve the same outcome. I am eager to understand if there are any ignored things on my end or if manual intervention is always required in Wannier v3.1.0.</p><p><br></p><p>I appreciate your insights on this matter and look forward to your guidance.</p><p><br></p><p>Thank you for your attention and assistance.</p><p>Best regards,</p><p>Tenghui Li </p><p>Department of Chemistry, Zhejiang University, China</p><p><br></p><p><br></p><blockquote name="replyContent" class="ReferenceQuote" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0;">-----原始邮件-----<br>
<b>发件人:</b> <span id="rc_from">"Giovanni Pizzi" <giovanni.pizzi@epfl.ch></span><br>
<b>发送时间:</b> <span id="rc_senttime">2024-06-07 20:21:18 (星期五)</span><br>
<b>收件人:</b> 李滕辉 <lith@zju.edu.cn><br>
<b>抄送:</b> "wannier@lists.quantum-espresso.org" <wannier@lists.quantum-espresso.org><br>
<b>主题:</b> Re: [Wannier] Centre of Wannier functions<br><br>
<div style="overflow-wrap:break-word;-webkit-nbsp-mode:space;line-break:after-white-space;">Hi, I don't think there is a very simple solution by just changing some code flags.<div><br>
<div>My suggestion would be, assuming (as it should be) that you have as small tilmestep in your MD, that you post-process the data and shift at every step all WFs by integer lattice vectors, so as to keep their position as "continuous" as possible (essentially,
given a WF center at step i, you look for all possible replicas by shifting by lattice vectors e.g. in a small supercell, and find the closest one to the corresponding WF in MD step (i-1), and consider that as the WF position).</div><div><br>
</div><div>Best,</div><div>Giovanni</div><div><br>
</div><div><div><br>
<blockquote type="cite"><div>On 7 Jun 2024, at 10:58, 李滕辉 <lith@zju.edu.cn> wrote:</div><br class="Apple-interchange-newline">
<div><div class="default-font-1717743508832" dir="ltr" style="font-size:16px;caret-color:#000000;font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;">Dear Wannier90 developers:<div>I meet an issue while using Wannier90 v3.1.0, together with VASP v6.4.1. I've done some AIMD calculations to generate a continuous trajectory. Then I proceeded with Wannier analysis and subsequently examined the Hamiltonian. Within the input file, I have
specified the atom centers, aiming to obtain WFs centered on these atom locations, or at least in close proximity to them.</div><div><br>
</div><div><div style="margin:0px;">...</div><div style="margin:0px;">Begin Projections</div><div style="margin:0px;"> Cd:sp3:z=-1,-1,-1,x=1,0,-1</div><div style="margin:0px;"> Se:p:z=-1,-1,-1,x=1,0,-1</div><div style="margin:0px;">End Projections</div><div style="margin:0px;">...</div><div style="margin:0px;">begin unit_cell_cart</div><div style="margin:0px;"> 6.0500000 0.0000000 0.0000000</div><div style="margin:0px;"> 0.0000000 6.0500000 0.0000000</div><div style="margin:0px;"> 0.0000000 0.0000000 6.0500000</div><div style="margin:0px;">end unit_cell_cart</div><div style="margin:0px;">begin atoms_cart</div><div style="margin:0px;">Cd 1.5178271 1.7957993 1.4240287</div><div style="margin:0px;">Cd 4.4307325 4.4894561 1.7053815</div><div style="margin:0px;">Cd 4.5160066 1.2418392 4.5483831</div><div style="margin:0px;">Cd 1.5803075 4.4594387 4.5691762</div><div style="margin:0px;">Se 0.0050033 0.0332671 5.9680866</div><div style="margin:0px;">Se 6.0433290 3.0636897 3.0562890</div><div style="margin:0px;">Se 2.8807786 0.0724863 3.0164791</div><div style="margin:0px;">Se 3.2493695 3.0420933 5.8999131</div><div style="margin:0px;">end atoms_cart</div></div><div>...</div><div><br>
</div><div>However, I have observer that the WFs are not aligning with the atom centers of the home cell. Instead, they are centered on atoms within a unit cell whose center is located at(0.0, 0.0, 0.0). For instance, since I use<span class="Apple-converted-space"> </span><em>sp<sup>3</sup></em><sup><span class="Apple-converted-space"> </span></sup>orbitals
for Cd, the center of WF5 corresponds to the second Cd atoms, while the atom center positioned at (4.43. 4.49. 1.71) and the WF center at (-1.01, -2.12, 1.14) in the .wout file.</div><div><br>
</div><div>...</div><div><div style="margin:0px;"> Final State</div><div style="margin:0px;"> WF centre and spread 1 ( 2.085684, 1.260598, 0.798043 ) 3.47876053</div><div style="margin:0px;"> WF centre and spread 2 ( 2.049201, 2.344944, 1.921311 ) 3.44168714</div><div style="margin:0px;"> WF centre and spread 3 ( 0.904277, 1.241577, 1.954697 ) 3.43610592</div><div style="margin:0px;"> WF centre and spread 4 ( 1.042490, 2.326357, 0.905988 ) 3.58232243</div><div style="margin:0px;"> WF centre and spread 5 ( -1.010143, -2.123145, 1.141211 ) 3.49051492</div><div style="margin:0px;"> WF centre and spread 6 ( -1.128292, -0.972610, 2.258101 ) 3.52243681</div><div style="margin:0px;"> WF centre and spread 7 ( -2.121837, -2.038340, 2.220324 ) 3.46813098</div><div style="margin:0px;"> WF centre and spread 8 ( -2.218949, -1.012121, 1.103897 ) 3.51827962</div><div style="margin:0px;"> WF centre and spread 9 ( -1.015362, 0.681962, -2.070276 ) 3.60904653</div><div style="margin:0px;"> WF centre and spread 10 ( -0.890789, 1.882241, -0.973578 ) 3.63289332</div><div style="margin:0px;"> WF centre and spread 11 ( -2.108611, 0.719542, -0.996879 ) 3.46637567</div><div style="margin:0px;"> WF centre and spread 12 ( -2.115134, 1.897953, -2.162619 ) 3.66683836</div><div style="margin:0px;"> WF centre and spread 13 ( 2.150342, -2.128650, -2.068309 ) 3.51187251</div><div style="margin:0px;"> WF centre and spread 14 ( 2.156870, -1.008390, -0.926044 ) 3.53284787</div><div style="margin:0px;"> WF centre and spread 15 ( 1.050800, -2.157773, -0.945039 ) 3.46859217</div><div style="margin:0px;"> WF centre and spread 16 ( 0.958142, -1.003749, -2.060564 ) 3.59315863</div><div style="margin:0px;"> WF centre and spread 17 ( 0.037244, 0.001632, -0.049092 ) 2.50422447</div><div style="margin:0px;"> WF centre and spread 18 ( 0.015322, 0.080868, -0.134619 ) 2.21017329</div><div style="margin:0px;"> WF centre and spread 19 ( -0.059586, 0.030045, -0.095360 ) 2.28365022</div><div style="margin:0px;"> WF centre and spread 20 ( 0.037435, -3.029794, -2.977580 ) 2.46391151</div><div style="margin:0px;"> WF centre and spread 21 ( 0.000015, -2.932250, -3.012771 ) 2.31121349</div><div style="margin:0px;"> WF centre and spread 22 ( -0.056325, -2.986181, -2.999162 ) 2.38773168</div><div style="margin:0px;"> WF centre and spread 23 ( 2.893882, 0.030034, 3.037144 ) 2.38000522</div><div style="margin:0px;"> WF centre and spread 24 ( 2.928631, 0.123782, 2.984530 ) 2.38126096</div><div style="margin:0px;"> WF centre and spread 25 ( 2.821335, 0.075610, 3.032140 ) 2.43241519</div><div style="margin:0px;"> WF centre and spread 26 ( -2.791587, -3.058092, -0.116545 ) 2.45592149</div><div style="margin:0px;"> WF centre and spread 27 ( -2.744926, -2.945582, -0.212516 ) 2.32087924</div><div style="margin:0px;"> WF centre and spread 28 ( -2.851460, -3.004112, -0.153946 ) 2.45045905</div><div style="margin:0px;"> Sum of centres and spreads ( 0.018668,-17.703645, -0.597512 ) 85.00170923</div><div>...</div></div><div><br>
</div><div>While this discrepancy does not present an issue in static calculations, it becomes problematic in MD trajectories. For instance, When an atom moves from (0, 0, 0.49) to (0, 0, 0.51), the center shifts from (0, 0, 0.49) to (0, 0, -0.49), resulting in a discontinuity
in the Hamiltonian non diagonal element indices.</div><div><br>
</div><div>I've made attempts to address this issue by adjusting parameters in the .win file, such as setting translate_home_cell = true, or explicitly setting the center of each orbital. But none of them works. It seems that the WF centres are predetermined in ab
initial calculations. Therefore, I am reaching out to inquire if there exist specific parameters or solutions to solve this problem.</div><div><br>
</div><div>Thank you for your attention and assistance. </div><div><br>
</div><div>Best regards,</div><div>Tenghui Li</div><div>Department of Chemistry, Zhejiang University, China</div></div><span style="caret-color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;float:none;display:inline !important;">_______________________________________________</span><br style="caret-color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;">
<span style="caret-color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;float:none;display:inline !important;">Wannier
mailing list</span><br style="caret-color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;">
<a href="mailto:Wannier@lists.quantum-espresso.org" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;orphans:auto;text-align:start;text-indent:0px;text-transform:none;white-space:normal;widows:auto;word-spacing:0px;-webkit-text-stroke-width:0px;">Wannier@lists.quantum-espresso.org</a><br style="caret-color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;-webkit-text-stroke-width:0px;text-decoration:none;">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/wannier" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;orphans:auto;text-align:start;text-indent:0px;text-transform:none;white-space:normal;widows:auto;word-spacing:0px;-webkit-text-stroke-width:0px;">https://lists.quantum-espresso.org/mailman/listinfo/wannier</a></div></blockquote></div><br>
</div></div></div></blockquote></div>