<div dir="ltr"><div>Dear Dr. Gabriele Sclauzero<br><br>Very thanks for your reply and suggestions.<br></div><br>Sincerely yours<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jan 28, 2014 at 2:48 PM, Sclauzero Gabriele <span dir="ltr"><<a href="mailto:gabriele.sclauzero@mat.ethz.ch" target="_blank">gabriele.sclauzero@mat.ethz.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><br>
<br>
<br>
Thank you very much for your reply, it is really appreciated. I will changed all parameters of my input file. I have a planer molecule in X-Y plane and I want to use a planar lead in Y-Z. I realized I have to change the transmission direction to Z.<br>
<br>
</div>Good!<br>
<div class="im"><br>
Could I have a lead in X-Y plane and calculate transmission for all structure (as an extended molecule or scattering region)?<br>
Reference article for leads: Physica B: Condensed Matter, Volume 406, Issue 20, 15 October 2011, Pages 3969-3974<br>
<br>
</div>You can also calculate transport for a system geometry with leads made of "thick" nanowires or extend 2D surfaces as in the paper you cite. The latter case requires to sample the k_perpendicular-dependent transmission in the 2D Brillouin zone. Moreover, the calculation for each k_perp. is much more demanding than in the monatomic chain case.<br>
I would suggest you to start with the simple case of the monatomic chain and then move to the massive leads.<br>
The two PRB papers I have mentioned in my first reply can give you a good idea for the two cases. Those calculations have been performed with PWscf and PWCOND.<br>
<br>
GS<br>
<br>
<br>
Wishes<br>
<div><div class="h5"><br>
<br>
On Tue, Jan 28, 2014 at 1:41 PM, Gabriele Sclauzero <<a href="mailto:gabriele.sclauzero@epfl.ch">gabriele.sclauzero@epfl.ch</a><mailto:<a href="mailto:gabriele.sclauzero@epfl.ch">gabriele.sclauzero@epfl.ch</a>>> wrote:<br>
Dear Khadije,<br>
<br>
On top of the possible memory problem pointed out by Paolo, I see several issues in your input.<br>
<br>
1. Unless things have changed recently, you must use fully-relativistic pseudopotentials in combination with lspinorb=.true. (I wonder why the code does not stop earlier with an error message). Anyway, I suggest you to start with scalar-relativistic collinear calculations. They are much faster and the relaxed structure should not change much. Once you got familiar with collinear calculations in PWscf and PWCOND you might consider investigating spin-orbit effects.<br>
<br>
2. You are using PAW datasets, which is a relatively recent feature. I don't think that it has been extended to the transmission calculations with PWCOND yet (please check).<br>
<br>
3. The new atomic geometry is still far from what you need for the scattering region. You are still simulating an isolated (0D) object, not a 1D object. I would suggest you to use XCrysDen (<a href="http://www.xcrysden.org/" target="_blank">http://www.xcrysden.org/</a>) to visualize your input geometry. If you repeat your unit cell along z (using shift-n key), you will probably realize it immediately. If your left/right leads are monatomic Au chains, then you have to put a short linear chain on each side of the molecule (not a square as you do now) and then appropriately choose the length of the unit cell along z.<br>
<br>
4. Beware: in the transmission calculations the transport direction is always along z, not along y as you instead seem to assume now. Also, the scatterer has to be at the center of the unit cell, and the continuations of the leads (atomic chains) at the edges. Please have a VERY CAREFUL look at the examples in PWCOND/examples/, also with the aid of xcrysden.<br>
<br>
5. You will need to select the k-points in a more suitable way: one k point in the xy plane should be sufficient, while you need to check convergence along z. But before doing this, you need to solve issue #3.<br>
<br>
<br>
HTH<br>
<br>
GS<br>
ETHZ, Materials Theory<br>
<br>
Dear Gabriele Sclauzero<br>
<br>
Thanks alot for your answer. I changed my molecule and added 4 Au atoms instead of 2. I am getting an error after running it. could you help me please?<br>
<br>
input file:<br>
&control<br>
calculation = 'relax'<br>
restart_mode='from_scratch',<br>
prefix='adn',<br>
tprnfor = .true.<br>
pseudo_dir = '/home/khalili/espresso-5.0.2/pseudo',<br>
outdir='./'<br>
/<br>
&system<br>
ibrav= 0, celldm(1)=6.0,<br>
nat= 23, ntyp= 5,<br>
noncolin=.true.,<br>
lspinorb=.true.,<br>
starting_magnetization(1)=0.0,<br>
ecutwfc = 35.0,<br>
ecutrho = 350.0,<br>
occupations='smearing',<br>
smearing='methfessel-paxton',<br>
degauss=0.01<br>
/<br>
&electrons<br>
diagonalization='david'<br>
electron_maxstep = 500,<br>
mixing_mode = 'plain'<br>
mixing_beta = 0.3<br>
conv_thr = 1.0d-8<br>
/<br>
&ions<br>
ion_dynamics = 'bfgs'<br>
/<br>
ATOMIC_SPECIES<br>
C 12.0107 C.pbe-n-kjpaw.UPF<br>
H 1.0 H.pbe-kjpaw.UPF<br>
N 14.0067 N.pbe-kjpaw.UPF<br>
Au 196.96655 Au.pbe-dn-kjpaw.UPF<br>
S 32.065 S.pbe-n-kjpaw.UPF<br>
ATOMIC_POSITIONS<br>
N 0.170976000 0.000000000 0.000000000<br>
C -0.170976000 0.000000000 0.000000000<br>
N -0.382005000 0.268863000 0.000000000<br>
C -0.210257000 0.563407000 0.000016000<br>
C 0.142155000 0.600825000 -0.000027000<br>
C 0.334612000 0.302284000 0.000800000<br>
N 0.676126000 0.302006000 0.001939000<br>
N 0.233362000 0.942208000 -0.004019000<br>
C -0.064351000 1.105626000 -0.011570000<br>
N -0.330002000 0.881899000 0.000315000<br>
S -0.270976000 -0.450000000 0.000000000<br>
H 0.798025000 0.039749000 0.001986000<br>
H 0.805314000 0.502022000 -0.000734000<br>
S -0.064351000 1.568242000 -0.011570000<br>
H -0.553279000 0.914575000 0.131427000<br>
Au -0.464351000 1.914575000 -0.011570000<br>
Au 0.335649000 1.914575000 -0.011570000<br>
Au -0.064351000 1.914575000 -0.411570000<br>
Au -0.064351000 1.914575000 0.388430000<br>
Au -0.670976000 -0.950000000 0.000000000<br>
Au 0.129024000 -0.950000000 0.000000000<br>
Au -0.270976000 -0.950000000 0.400000000<br>
Au -0.270976000 -0.950000000 -0.400000000<br>
CELL_PARAMETERS cubic<br>
6.0 0.0 0.0<br>
0.0 6.0 0.0<br>
0.0 0.0 6.0<br>
K_POINTS {Automatic}<br>
2 2 1 1 1 1<br>
<br>
</div></div>mpirun -np 20 /home/khalili/espresso-5.0.2/bin/pw.x -in <a href="http://adn.rx.in" target="_blank">adn.rx.in</a><<a href="http://adn.rx.in/" target="_blank">http://adn.rx.in/</a>> | tee adn.rx.out<br>
<div class="im"><br>
<br>
output:<br>
.<br>
.<br>
.<br>
per-process dynamical memory: 627.4 Mb<br>
<br>
Self-consistent Calculation<br>
<br>
iteration # 1 ecut= 35.00 Ry beta=0.70<br>
Davidson diagonalization with overlap<br>
<br>
===================================================================================<br>
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES<br>
= EXIT CODE: 9<br>
= CLEANING UP REMAINING PROCESSES<br>
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES<br>
===================================================================================<br>
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)<br>
This typically refers to a problem with your application.<br>
Please see the FAQ page for debugging suggestions.<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
</div><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><mailto:<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>><br>
<div class="im"><a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
<br>
--<br>
Khadije Khalili<br>
Ph.D Student of Solid-State Physics<br>
Department of Physics<br>
University of Mazandaran<br>
Babolsar, Iran<br>
</div><a href="mailto:kh.khalili@stu.umz.ac.ir">kh.khalili@stu.umz.ac.ir</a><mailto:<a href="mailto:kh.khalili@stu.umz.ac.ir">kh.khalili@stu.umz.ac.ir</a>><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><mailto:<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
Dr. Gabriele Sclauzero<br>
Materials Theory - ETHZ<br>
ETH Zurich, HIT G 43.2<br>
Wolfgang-Pauli-Str. 27<br>
8093 Zürich, Switzerland<br>
<br>
Phone +41 44 633 94 10<br>
Fax +41 44 633 14 59<br>
<a href="mailto:gabriele.sclauzero@mat.ethz.ch">gabriele.sclauzero@mat.ethz.ch</a><mailto:<a href="mailto:gabriele.sclauzero@mat.ethz.ch">gabriele.sclauzero@mat.ethz.ch</a>><br>
<a href="http://www.theory.mat.ethz.ch" target="_blank">www.theory.mat.ethz.ch</a><br>
<div class="HOEnZb"><div class="h5"><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div>Khadije Khalili</div><div>Ph.D Student of Solid-State Physics</div><div>Department of Physics</div><div>University of Mazandaran</div><div>
Babolsar, Iran</div><div><a href="mailto:kh.khalili@stu.umz.ac.ir" target="_blank">kh.khalili@stu.umz.ac.ir</a></div><div> </div></div>
</div>