<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Pankaj,<div><br></div><div>&nbsp; &nbsp; you already posted for a similar problem about one month ago. The reply is the same. Have you ever stopped a minute to think how much should be the transmission values that one could expect for this system?</div><div><br></div><div>I see at least one or two additional problems in your input files:</div><div>1. the lead region needs some k-point sampling along z, gamma is not sufficient!</div><div>2. eps_proj looks way too large for me, while nz1 seems too small. Generally you should check convergence of your results w.r.t to these (and other) parameters.</div><div>... but there might be other errors that I could not spot...</div><div>I suggest you to practice with smaller examples in order to understand the method (just running the examples provided is not enough, désolé...)</div><div><br></div><div>HTH</div><div><br></div><div>GS</div><div><br><div><div>Il giorno 21/mag/2012, alle ore 09.30, pankaj sahota ha scritto:</div><blockquote type="cite"><meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"><div><font color="#330099">Hello Everybody,</font></div><div><font color="#330099"><br></font></div><div><font color="#330099">I am trying to reproduce the results of a paper by X. Qian (PHYSICAL REVIEW B 82, 195442 2010). But in the result i am getting ******* in the transmission files. I am pasting the input files for the scattering, lead and the input file for pwcond.x. Can any one have a look on where I am going wrong. I have checked the compilers and ran all the examples in the code. They are all running fine.</font></div>
<font size="4"><b><div><font size="4"><b><br></b></font></div><div><font size="4"><b><br></b></font></div>Scattering region:</b></font><div><font color="#330099"><br></font></div><div><div style="color:rgb(51,0,153)">&nbsp;&amp;CONTROL</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;calculation = 'scf' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; outdir = './tmp' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; pseudo_dir = '/usr/local/src/espresso_potentials/' ,</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; prefix = 'c7-al' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;etot_conv_thr = 1D-7 ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; verbosity = 'high',</div>
<div style="color:rgb(51,0,153)">&nbsp;/</div><div style="color:rgb(51,0,153)">&nbsp;&amp;SYSTEM</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ibrav = 6,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;celldm(1) = 26.46,</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;celldm(3) = 1.57771,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;nat = 39,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ntyp = 2,</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;occupations='smearing' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;smearing = 'mp' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;degauss = 0.01 ,</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ecutwfc = 27 ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ecutrho = 300 ,</div><div style="color:rgb(51,0,153)">&nbsp;/</div><div style="color:rgb(51,0,153)">
&nbsp;&amp;ELECTRONS</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;conv_thr = 1.D-6,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;startingpot = 'file' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;startingwfc = 'atomic' ,</div>
<div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;mixing_mode = 'plain' ,</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;mixing_beta = 0.7 ,</div><div style="color:rgb(51,0,153)">&nbsp;/</div><div style="color:rgb(51,0,153)">
&nbsp;ATOMIC_SPECIES</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; Al &nbsp; 26.981539 &nbsp;Al.pbe-rrkj.UPF</div><div style="color:rgb(51,0,153)">&nbsp; &nbsp; C &nbsp; &nbsp;12.0107 &nbsp; &nbsp;C.pbe-rrkjus.UPF</div><div style="color:rgb(51,0,153)">ATOMIC_POSITIONS crystal</div>
<div style="color:rgb(51,0,153)"><div>Al 0.289286 0.433929 0.0000000</div><div>Al 0.433929 0.289286 0.0000000</div><div>Al 0.578571 0.433929 0.0000000</div><div>Al 0.433929 0.578571 0.0000000</div><div>Al 0.289286 0.289286 0.0916787</div>
<div>Al 0.289286 0.578571 0.0916787</div><div>Al 0.578571 0.289286 0.0916787</div><div>Al 0.578571 0.578571 0.0916787</div><div>Al 0.433929 0.433929 0.0916787</div><div>Al 0.289286 0.433929 0.183357</div><div>Al 0.433929 0.289286 0.183357</div>
<div>Al 0.578571 0.433929 0.183357</div><div>Al 0.433929 0.578571 0.183357</div><div>Al 0.289286 0.289286 0.275036</div><div>Al 0.289286 0.578571 0.275036</div><div>Al 0.578571 0.289286 0.275036</div><div>Al 0.578571 0.578571 0.275036</div>
<div>Al 0.433929 0.433929 0.275036</div><div>C 0.433929 0.433929 0.3203100</div><div>C 0.433929 0.433929 0.380206</div><div>C 0.433929 0.433929 0.440103</div><div>C 0.433929 0.433929 0.500000</div><div>C 0.433929 0.433929 0.559897</div>
<div>C 0.433929 0.433929 0.619794</div><div>C 0.433929 0.433929 0.67969</div><div>Al 0.289286 0.289286 0.724964</div><div>Al 0.289286 0.578571 0.724964</div><div>Al 0.578571 0.289286 0.724964</div><div>Al 0.578571 0.578571 0.724964</div>
<div>Al 0.433929 0.433929 0.724964</div><div>Al 0.289286 0.433929 0.816643</div><div>Al 0.433929 0.289286 0.816643</div><div>Al 0.578571 0.433929 0.816643</div><div>Al 0.433929 0.578571 0.816643</div><div>Al 0.289286 0.289286 0.908321</div>
<div>Al 0.289286 0.578571 0.908321</div><div>Al 0.578571 0.289286 0.908321</div><div>Al 0.578571 0.578571 0.908321</div><div>Al 0.433929 0.433929 0.908321</div><div>K_POINTS gamma</div></div><div style="color:rgb(51,0,153)">
<br></div><div><font size="4"><b>Lead Region:</b></font></div><div style="color:rgb(51,0,153)"><br></div><div style="color:rgb(51,0,153)"><div>&amp;CONTROL</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;calculation = 'scf' ,</div><div>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; outdir = './tmp' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; pseudo_dir = '/home/pankajkumar/potentials/quantum_es/' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; prefix = 'lead2' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;etot_conv_thr = 1D-7 ,</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; verbosity = 'high',</div><div>&nbsp;/</div><div>&nbsp;&amp;SYSTEM</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ibrav = 6,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;celldm(1) = 26.46,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;celldm(3) = 0.289285,</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;nat = 9,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ntyp = 1,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;occupations='smearing' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;smearing = 'mp' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;degauss = 0.01 ,</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ecutwfc = 27 ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ecutrho = 300 ,</div><div>&nbsp;/</div><div>&nbsp;&amp;ELECTRONS</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;conv_thr = 1.D-6,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;startingpot = 'file' ,</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;startingwfc = 'atomic' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;mixing_mode = 'plain' ,</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;mixing_beta = 0.7 ,</div><div>&nbsp;/</div><div>&nbsp;ATOMIC_SPECIES</div><div>&nbsp; &nbsp; Al &nbsp; 26.981539 &nbsp;Al.pbe-rrkj.UPF</div>
<div>ATOMIC_POSITIONS crystal</div><div>Al 0.289285 0.289285 0.0</div><div>Al 0.289285 0.578571 0.0</div><div>Al 0.578571 0.289285 0.0</div><div>Al 0.578571 0.578571 0.0</div><div>Al 0.433928 0.433928 0.0</div><div>Al 0.289285 0.433928 0.5</div>
<div>Al 0.433928 0.289285 0.5</div><div>Al 0.578571 0.433928 0.5</div><div>Al 0.433928 0.578571 0.5</div></div><div style="color:rgb(51,0,153)"><div>K_POINTS gamma</div></div><div style="color:rgb(51,0,153)"><br></div><div style="color:rgb(51,0,153)">
Input file for calculating the conduction:</div><div style="color:rgb(51,0,153)"><br></div><div style="color:rgb(51,0,153)"><div>&amp;inputcond</div><div>outdir='./tmp'</div><div>prefixl='lead2',</div><div>
prefixs='c7-al',</div><div>ikind=1,</div><div>energy0=1.0d0,</div><div>denergy=-0.05d0,</div><div>ewind=1.0d0,</div><div>tran_file='spin_con_up',</div><div>nz1 = 1,</div><div>epsproj = 1.d-3,</div><div>/</div>
<div>&nbsp; &nbsp; 1</div><div>&nbsp; &nbsp; 0.0 &nbsp;0.0 &nbsp;1.0</div><div>41</div></div><div style="color:rgb(51,0,153)"><br></div><div><font size="4"><b style="background-color:rgb(255,255,255)">output file for transmission:</b></font></div><div style="color:rgb(51,0,153)">
<br></div><div style="color:rgb(51,0,153)"><div># E-Ef, T</div><div>&nbsp; 1.000000**********</div><div>&nbsp; 0.950000**********</div><div>&nbsp; 0.900000412.340267</div><div>&nbsp; 0.850000497.711187</div><div>&nbsp; 0.800000802.461207</div><div>
&nbsp; 0.750000440.526889</div><div>&nbsp; 0.700000958.543158</div><div>&nbsp; 0.650000389.683665</div><div>&nbsp; 0.600000908.336839</div><div>&nbsp; 0.550000586.578929</div><div>&nbsp; 0.500000362.790961</div><div>&nbsp; 0.450000268.789968</div><div>&nbsp; 0.400000174.116786</div>
<div>&nbsp; 0.350000265.822178</div><div>&nbsp; 0.300000**********</div><div>&nbsp; 0.250000171.305651</div><div>&nbsp; 0.200000**********</div><div>&nbsp; 0.150000559.105443</div><div>&nbsp; 0.100000**********</div><div>&nbsp; 0.050000**********</div><div>
&nbsp; 0.000000**********</div><div>&nbsp;-0.050000**********</div><div>&nbsp;-0.100000**********</div><div>&nbsp;-0.150000866.786899</div><div>&nbsp;-0.200000**********</div><div>&nbsp;-0.250000**********</div><div>&nbsp;-0.300000378.547020</div><div>&nbsp;-0.350000268.718613</div>
<div>&nbsp;-0.400000162.547897</div><div>&nbsp;-0.450000131.448607</div><div>&nbsp;-0.500000199.233657</div><div>&nbsp;-0.550000532.739318</div><div>&nbsp;-0.600000145.867531</div><div>&nbsp;-0.650000**********</div><div>&nbsp;-0.700000915.016527</div><div>
&nbsp;-0.750000376.695715</div><div>&nbsp;-0.800000**********</div><div>&nbsp;-0.850000**********</div><div>&nbsp;-0.900000**********</div><div>&nbsp;-0.950000**********</div><div>&nbsp;-1.000000**********</div><div><br></div><div>Thanks and regards.&nbsp;</div>
<div>Pankaj.</div></div><div style="color:rgb(51,0,153)"><br></div><div><br></div>-- <br><font color="#000099">Pankaj Kumar<br>DST Project Fellow(Ph.D.)<br>Department of Physics<br>The LNM IIT Jaipur<br>Mob. No. +91 9736694726</font><br>
<br><br>
</div>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero,&nbsp;EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i>&nbsp;&nbsp; PH H2 462, Station 3,&nbsp;CH-1015 Lausanne</i></font></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline"></span>
</div>
<br></div></body></html>