<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="color:rgb(7,55,99)">Dear Wannier developers and users,</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)"> I am new to the Wannier90 code and trying to calculate the polarization of bulk BaTiO3.</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)"> I think the bulk polarization should not depend on the choice of initial projections once the convergence has reached but, strangely enough, the results seem to depend on the initial projections schemes. </div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)"> For example, I tested two initial projections for BaTiO3.</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Case #1: atomic orbital projections</div><div class="gmail_default"><div class="gmail_default"><font color="#073763">dis_num_iter = 80000</font></div><div class="gmail_default"><font color="#073763">dis_mix_ratio = 0.9d0</font></div><div class="gmail_default"><font color="#073763">dis_conv_tol = 1.0d-9</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">conv_window  = 100</font></div><div class="gmail_default"><font color="#073763">conv_tol = 1.0d-9</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">search_shells = 50</font></div><div class="gmail_default"><br></div><div class="gmail_default"><font color="#073763"><b>Begin Projections</b></font></div><div class="gmail_default"><font color="#073763"><b>Ba: s;px;py;pz</b></font></div><div class="gmail_default"><font color="#073763"><b>O: sp3</b></font></div><div class="gmail_default"><font color="#073763"><b>End Projections</b></font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">guiding_centres = true</font></div><div class="gmail_default"><font color="#073763">#num_no_guide_iter = 20</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">num_iter = 80000</font></div><div class="gmail_default"><font color="#073763">num_print_cycles = 100</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">use_bloch_phases = F</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">write_r2mn = true</font></div><div class="gmail_default"><font color="#073763">write_xyz = true</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">num_bands = 32</font></div><div class="gmail_default"><font color="#073763">num_wann = 16</font></div><div style="color:rgb(7,55,99)">===========================================<br></div><div style="color:rgb(7,55,99)"><br></div></div><div class="gmail_default" style="color:rgb(7,55,99)">Case #2: random projections</div><div class="gmail_default"><div class="gmail_default"><font color="#073763">dis_num_iter = 80000</font></div><div class="gmail_default"><font color="#073763">dis_mix_ratio = 0.9d0</font></div><div class="gmail_default"><font color="#073763">dis_conv_tol = 1.0d-9</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">conv_window  = 100</font></div><div class="gmail_default"><font color="#073763">conv_tol = 1.0d-9</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">search_shells = 50</font></div><div class="gmail_default"><br></div><div class="gmail_default"><font color="#073763"><b>Begin Projections</b></font></div><div class="gmail_default"><font color="#073763"><b>Random</b></font></div><div class="gmail_default"><font color="#073763"><b>End Projections</b></font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">guiding_centres = true</font></div><div class="gmail_default"><font color="#073763">#num_no_guide_iter = 20</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">num_iter = 80000</font></div><div class="gmail_default"><font color="#073763">num_print_cycles = 100</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">use_bloch_phases = F</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">write_r2mn = true</font></div><div class="gmail_default"><font color="#073763">write_xyz = true</font></div><div class="gmail_default"><font color="#073763"><br></font></div><div class="gmail_default"><font color="#073763">num_bands = 32</font></div><div class="gmail_default"><font color="#073763">num_wann = 16</font></div><div style="color:rgb(7,55,99)"><br></div></div><div class="gmail_default" style="color:rgb(7,55,99)">===================================================</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">The only difference in the two input files is the Projections block.</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Although both cases are converged, the sum of wannier centers for the occupied bands are different for case 1 and case 2.</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Can anyone help me with this?</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Any comment and/or suggenstions would be very much appreciated.</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Thank you. </div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div class="gmail_default" style="color:rgb(7,55,99)">Best,</div><div class="gmail_default" style="color:rgb(7,55,99)"><br></div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="font-size:small"><b>|</b></span><b style="font-size:small"> Duk-Hyun Choe</b><span style="font-size:small">, Post-doctoral fellow, </span><br></div><div>  Department of Physics, Applied Physics, and Astronomy <br><div>  Rensselaer Polytechnic Institute, Troy, NY, USA<br></div>  E-mail: <a href="mailto:hoalad@gmail.com" target="_blank">hoalad@gmail.com</a>,  Tel: +1-518-687-4389<div>  Home: <a href="https://sites.google.com/site/dhchoeprofile/" target="_blank">https://sites.google.com/site/dhchoeprofile/</a></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>