<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi,
<div class=""><br class="">
</div>
<div class="">When I was speaking of convergence issues, I was not thinking to the preliminary DFT computation, but mainly to the convergence with the dense grid in BoltzWann - the calculation includes a derivative of a Fermi-Dirac distribution, that is a Dirac-delta
in the zero-temperature limit. Smearing (of the TDF, in BoltzWann) can help but I think will give an incorrect result in practice. </div>
<div class=""><br class="">
</div>
<div class="">Also, the Seebeck coefficient calculation requires to invert a matrix, and at low temperature with the chemical potential in the electronic gap, one risks to ask the code to invert a matrix with basically zero eigenvalues, which is also a problem.</div>
<div class=""><br class="">
</div>
<div class="">Giovanni</div>
<div class=""><br class="">
<div class="">
<div class="">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;">
-- <br class="">
Giovanni Pizzi<br class="">
Theory and Simulation of Materials and MARVEL, EPFL<br class="">
<a href="http://people.epfl.ch/giovanni.pizzi" class="">http://people.epfl.ch/giovanni.pizzi</a><br class="">
http://nccr-marvel.ch/en/people/profile/giovanni-pizzi</div>
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">On 25 Jun 2017, at 05:23, Jun Liu <<a href="mailto:jun.physics@gmail.com" class="">jun.physics@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">Dear Prof. Giovanni,<br class="">
<br class="">
</div>
Thank you very much for your help. I do need to see how conductivity changes across the three cases I mentioned under T=2K. But I should be satisfied with the trend of change if the results look reasonable at higher T. I have submitted jobs for T=50K and see
how they will work out (still not finished yet). <br class="">
<div class=""><br class="">
</div>
<div class="">Could you please elaborate on what you mean by numerical and convergence issues? I think I should have at least reached K-points and ENCUT convergence.<br class="">
<br class="">
</div>
<div class="">Thanks,<br class="">
<br class="">
</div>
<div class="">Sincerely,<br class="">
</div>
Jun</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Thu, Jun 15, 2017 at 12:24 PM, Giovanni Pizzi <span dir="ltr" class="">
<<a href="mailto:giovanni.pizzi@epfl.ch" target="_blank" class="">giovanni.pizzi@epfl.ch</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">Dear Jun,
<div class="">You are working in the range of temperatures 2-10 K. Is this expected? Such low temperatures are very hard to work with in the code, because of numerical and convergence issues.</div>
<div class=""><br class="">
</div>
<div class="">If you work at higher temperature, is the result better?</div>
<div class=""><br class="">
</div>
<div class="">Also, remember that by changing the strain you will change the energy position of the bands (as basically all materials have non-zero deformation potentials) and therefore you cannot really compare values at the same exact chemical potential mu.
E.g. it is possible that when you get close to zero, it is because at that strain there’s no band close to the chemical potential (and since you are at low T, you need bands really close). Changing strain, bands enter in the region.</div>
<div class=""><br class="">
</div>
<div class="">I suggest that you start by T>50K (and without boltz_tdf_smr_fixed_en_width) and you compare the whole plots in a chemical potential range, and once you understand the behaviour you go down in temperature, potentially putting back some smearing
in the TDF, to check also if you are converged enough. Note however that the TDF smearing while probably needed at low T, might introduce some errors in the final results.</div>
<div class=""><br class="">
</div>
<div class="">Hope this helps,</div>
<div class=""><br class="">
</div>
<div class="">Giovanni</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;" class="">
-- <br class="">
Giovanni Pizzi<br class="">
Theory and Simulation of Materials and MARVEL, EPFL<br class="">
<a href="http://people.epfl.ch/giovanni.pizzi" target="_blank" class="">http://people.epfl.ch/<wbr class="">giovanni.pizzi</a><br class="">
<a href="http://nccr-marvel.ch/en/people/profile/giovanni-pizzi" target="_blank" class="">http://nccr-marvel.ch/en/<wbr class="">people/profile/giovanni-pizzi</a></div>
</div>
<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">
<div class="h5">
<div class="">On 31 May 2017, at 18:42, Jun Liu <<a href="mailto:jun.physics@gmail.com" target="_blank" class="">jun.physics@gmail.com</a>> wrote:</div>
<br class="m_-4683542062244942626Apple-interchange-newline">
</div>
</div>
<div class="">
<div class="">
<div class="h5">
<div dir="ltr" class="">
<div class="">
<div class="">
<div class="">Dear wannier users,<br class="">
<br class="">
</div>
<div class="">I would like to ask whether the following input for conductivity tensor calculation looks ok.
<br class="">
<br class="">
num_bands = 128 ! set to NBANDS by VASP<br class="">
num_wann = 88<br class="">
begin projections<br class="">
Mo:d ! 20<br class="">
Te:p ! 24<br class="">
end projections<br class="">
<br class="">
###########BoltzWann ############<br class="">
boltzwann = true<br class="">
kmesh = 400 400 400<br class="">
boltz_relax_time = 0.001<br class="">
boltz_mu_min = 7.95123782 <wbr class=""> #ef-def<br class="">
boltz_mu_max = 8.35123782<br class="">
boltz_mu_step = 0.04<br class="">
boltz_temp_min = 2<br class="">
boltz_temp_max = 10<br class="">
boltz_temp_step = 2<br class="">
boltz_tdf_energy_step=0.01<br class="">
boltz_tdf_smr_fixed_en_width = 0.01<br class="">
boltz_tdf_smr_type = gauss<br class="">
boltz_calc_also_dos = true<br class="">
boltz_dos_energy_min = 5.0<br class="">
boltz_dos_energy_max = 15.0<br class="">
boltz_dos_energy_step = 0.01<br class="">
##############################<wbr class="">###<br class="">
<br class="">
write_hr = .true.<br class="">
# Bandstructure<br class="">
restart = plot<br class="">
bands_plot = true<br class="">
begin kpoint_path<br class="">
....<br class="">
end kpoint_path<br class="">
bands_num_points 100 # bands_plot_format gnuplot<br class="">
<br class="">
<br class="">
spinors = .true.<br class="">
<br class="">
begin unit_cell_cart<br class="">
<div class=""> 3.4716802 0.0000000 0.0000000 /*this and the next line changes*/<br class="">
</div>
<div class=""> 0.0000000 6.3666750 0.0000000 /*accordingly for different cases*/<br class="">
</div>
0.0000000 0.0000000 13.8452386<br class="">
end unit_cell_cart<br class="">
<br class="">
begin atoms_cart<br class="">
...<br class="">
end atoms_cart<br class="">
<br class="">
mp_grid = 12 6 3 /*this line differs from each case*/<br class="">
<br class="">
begin kpoints<br class="">
...<br class="">
end kpoints<br class="">
<br class="">
</div>
<div class="">This input returns reasonable result (in that the conductivity tensor is nearly diagonal for orthorhombic structure). But if I compare results across different cases, they differ too much under as small as 0.5% lattice constant change. More details
are given below.<br class="">
</div>
<div class=""><br class="">
</div>
I tried to calculate the conductivity tensor with postw90.x with version 2.1. I tried it on three cases related with slight elongation along a or b directions defining an orthorhombic structure on the original lattice. The resulting conductivity tensor differs
by nearly 100% for a 0.5% lattice constant change, which cannot be right. To give you some numbers, the following shows different sigma_x,x scanned under different chemical potentials for different lattice structures, (other components of sigma are suppressed
for brevity but the whole sigma matrix does show a diagonal feature expected for an orthorhombic structure.)
<br class="">
<br class="">
</div>
For the base structure (a,b,c)<br class="">
</div>
mu T sigma_x,x<br class="">
<div class=""> 7.963542580 2.000000000 0.4228380355E-05<br class="">
7.963542580 4.000000000 0.6146815602E-01<br class="">
7.963542580 6.000000000 1.261398051 <br class="">
7.963542580 8.000000000 5.250207712 <br class="">
7.963542580 10.00000000 11.71261179 <br class="">
<br class="">
</div>
<div class="">For the elongated structure along a by 0.5%<br class="">
mu T sigma_x,x<br class="">
7.962242580 2.000000000 0.5201871790E-03<br class="">
7.962242580 4.000000000 0.6688530829 <br class="">
7.962242580 6.000000000 6.051666588 <br class="">
7.962242580 8.000000000 16.35207686 <br class="">
7.962242580 10.00000000 27.52636769 <br class="">
<br class="">
</div>
<div class="">For the elongated structure along b by 0.5%<br class="">
mu T sigma_x,x<br class="">
7.951237820 2.000000000 <a href="tel:%28332%29%20903-9665" value="+13329039665" target="_blank" class="">
332.9039665</a> <br class="">
7.951237820 4.000000000 363.9745377 <br class="">
7.951237820 6.000000000 287.1196636 <br class="">
7.951237820 8.000000000 <a href="tel:%28228%29%20918-7157" value="+12289187157" target="_blank" class="">
228.9187157</a> <br class="">
7.951237820 10.00000000 188.4869745 <br class="">
<br class="">
</div>
<div class="">Any idea on how to check the calculation? Thanks all very much for your insightful help!
<br class="">
<br class="">
</div>
<div class="">Sincerely,<br class="">
</div>
<div class="">Jun<br class="">
</div>
<div class=""><br class="">
</div>
</div>
</div>
</div>
______________________________<wbr class="">_________________<br class="">
Wannier mailing list<br class="">
<a href="mailto:Wannier@quantum-espresso.org" target="_blank" class="">Wannier@quantum-espresso.org</a><br class="">
<a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier" target="_blank" class="">http://mailman.qe-forge.org/<wbr class="">cgi-bin/mailman/listinfo/<wbr class="">wannier</a><br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
<br class="">
______________________________<wbr class="">_________________<br class="">
Wannier mailing list<br class="">
<a href="mailto:Wannier@quantum-espresso.org" class="">Wannier@quantum-espresso.org</a><br class="">
<a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier" rel="noreferrer" target="_blank" class="">http://mailman.qe-forge.org/<wbr class="">cgi-bin/mailman/listinfo/<wbr class="">wannier</a><br class="">
<br class="">
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</body>
</html>