<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="">
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="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 31 May 2017, at 18:42, 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="">
<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>
_______________________________________________<br class="">
Wannier mailing list<br class="">
<a href="mailto:Wannier@quantum-espresso.org" class="">Wannier@quantum-espresso.org</a><br class="">
http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier<br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>