<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.2963" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Dear users</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I have a problem with decharging during my
transport calculations. I have read the suggestions in the previous mails and
acted accordingly.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I set my EnergyLowestBound to -100 Ry and I even
start the transport calculation from a converged Siesta
calculation.</FONT></DIV>
<DIV><FONT face=Arial size=2>I use the same pseudo´s, basis and kpoints in x-
and y-direction for lead and scattering region.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Still i get the following output in .CHR
:</FONT></DIV>
<DIV><FONT face=Arial size=2> negf: iter
Charge Bias [V]
Temperatu<BR>negf: 1
554.00000 0.00000 0.00190
Ry<BR>negf: 2 454.87244
0.00000 0.00190 Ry<BR>negf:
3 556.36230 0.00000 0.00190
Ry</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>What could be wrong? See my input file
below.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>/Anders</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>d14n36$ cat Au_bridge.fdf
<BR>SystemName
Au_bridge<BR>SystemLabel
Au_bridge<BR>NumberOfAtoms
64<BR>NumberOfSpecies
4<BR>WriteKpoints
.false.<BR>#ProcessorY
1<BR>DiagMemory
3<BR>DiagScale
3<BR>TryMemoryIncrease .true. <BR>MD.TypeOfRun
CG<BR>MD.UseSaveXV
.false.<BR>MD.NumCGsteps
0<BR>MaxSCFIterations 100</FONT></DIV>
<DIV> </DIV><FONT face=Arial size=2>
<DIV><BR>DM.MixingWeight 0.05 #Default value=0.25<BR>DM.NumberPulay 4
#0,1=linear
mixing<BR>XC.functional
GGA
#default=LDA<BR>XC.authors
PBE<BR>#SpinPolarized
.false
#default<BR>MeshCutoff
300 Ry
#default<BR>ElectronicTemperature 300.00 K
#default<BR>NeglNonOverlapInt
F<BR>DM.UseSaveDM
T<BR>OccupationFunction
MP<BR>OccupationMPOrder 5</DIV>
<DIV> </DIV>
<DIV>%block ChemicalSpeciesLabel<BR> 1 1
H<BR> 2 6 C<BR> 3 16 S <BR> 4 79 Au
#Species index, atomic number, species label<BR>%endblock
ChemicalSpeciesLabel</DIV>
<DIV> </DIV>
<DIV>LatticeConstant 1 Ang
</DIV>
<DIV> </DIV>
<DIV>%block LatticeVectors<BR>6.285 6.285 0.0<BR>6.285 -6.285
0.0<BR>0.0 0.0 26.065<BR>%endblock
LatticeVectors</DIV>
<DIV> </DIV>
<DIV>AtomicCoordinatesFormat ScaledCartesian<BR>#AtomCoorFormatOut Ang</DIV>
<DIV> </DIV>
<DIV>%block AtomicCoordinatesAndAtomicSpecies<BR>
2.0950
0 0
4<BR> 4.1900
2.0950 0
4<BR> 4.1900
-2.0950 0
4<BR> 6.2850
4.1900 0 4
<BR> 6.2850
0 0
4<BR> 6.2850
-4.1900 0
4<BR> 8.3800
2.0950 0
4<BR> 8.3800
-2.0950 0
4<BR> 10.4750
0 0
4<BR>
0 0
2.0950 4<BR> 2.0950
2.0950 2.0950 4<BR>
2.0950 -2.0950 2.0950
4<BR> 4.1900 4.1900
2.0950 4<BR>
4.1900 0
2.0950 4<BR> 4.1900
-4.1900 2.0950 4<BR>
6.2850 2.0950 2.0950
4<BR> 6.2850 -2.0950
2.0950 4<BR>
8.3800 0
2.0950 4<BR>
2.0950 0
4.1900 4<BR> 4.1900
2.0950 4.1900 4<BR>
4.1900 -2.0950 4.1900
4<BR> 6.2850
0 4.1900 4<BR>
4.9396 0.1407 8.1433
2.0000<BR> 5.0536 1.2630
8.9487 2.0000<BR> 5.6764
1.0098 10.2181 2.0000<BR>
6.0217 -0.3372 10.3530
2.0000<BR> 5.9735 2.0971
11.1574 2.0000<BR> 4.6868
2.2425 8.6435 1.0000<BR>
6.1706 3.0544 10.6651
1.0000<BR> 6.5267 -0.8115
11.1905 1.0000<BR> 5.6130
-1.2688 8.9509 3.0000<BR>
4.1900 0
6.5400 3.0000<BR> 6.0715
2.1131 12.5198 2.0000<BR>
5.7612 1.1048 13.5374
2.0000<BR> 4.8817 0.0224
13.4308 2.0000<BR> 6.2842
1.2154 14.8716 2.0000<BR>
5.8148 0.2491 15.7471
2.0000<BR> 6.4105 3.0628
12.9459 1.0000<BR> 6.9876
1.9875 15.1801 1.0000<BR>
4.3034 -0.2781 12.5600
1.0000<BR> 4.7096 -0.8433
14.9205 3.0000<BR>
6.2850 0
17.4370 3.0000<BR>
4.1900 0
19.7800 4.0000<BR>
8.3800 0
19.7800 4.0000<BR> 6.2850
2.0950 19.7800 4.0000<BR>
6.2850 -2.0950 19.7800
4.0000<BR>
2.0950 0
21.8750 4.0000<BR> 4.1900
2.0950 21.8750 4.0000<BR>
4.1900 -2.0950 21.8750
4.0000<BR> 6.2850 4.1900
21.8750 4.0000<BR>
6.2850 0
21.8750 4.0000<BR> 6.2850
-4.1900 21.8750 4.0000<BR>
8.3800 2.0950 21.8750
4.0000<BR> 8.3800 -2.0950
21.8750 4.0000<BR>
10.4750 0
21.8750
4.0000<BR>
0 0
23.9700 4.0000<BR> 2.0950
2.0950 23.9700 4.0000<BR>
2.0950 -2.0950 23.9700
4.0000<BR> 4.1900 4.1900
23.9700 4.0000<BR>
4.1900 0
23.9700 4.0000<BR> 4.1900
-4.1900 23.9700 4.0000<BR>
6.2850 2.0950 23.9700
4.0000<BR> 6.2850 -2.0950
23.9700 4.0000<BR>
8.3800 0
23.9700 4.0000<BR>%endblock
AtomicCoordinatesAndAtomicSpecies</DIV>
<DIV> </DIV>
<DIV>PAO.BasisSize
SZ<BR>#PAO.BasisType
split # Type of
basis ('nones', 'nonodes',
'split')<BR>#PAO.SplitNorm
0.15 # Amount
of norm carried by the second zeta</DIV>
<DIV> </DIV>
<DIV>%block PAO.Basis<BR>Au 3
0.27141<BR> n=6 0 1 #E
166.43341 3.66795<BR>
4.53756 #2.16868<BR>
1.00000 #1.00000<BR> n=6 1
1 #E 12.34073
4.02728<BR> 4.54026<BR> 1.00000<BR> n=5
2 1 #E 151.22986
4.34561<BR> 4.95606 #3.19253<BR>
1.00000 #1.00000<BR>%endblock PAO.Basis</DIV>
<DIV> </DIV>
<DIV>%block Ps.lmax<BR> H
1 <BR> C
2<BR> S
3<BR> Au 2<BR>%endblock
Ps.lmax<BR> <BR>%block kgrid_Monkhorst_Pack<BR>10 0 0 0<BR>0 10 0 0<BR>0 0
1 0<BR>%endblock kgrid_Monkhorst_Pack</DIV>
<DIV> </DIV>
<DIV>#Smeagol flags:<BR>#BulkTransport
.true.<BR>#BulkLeads
LR </DIV>
<DIV> </DIV>
<DIV>EMTransport
T<BR>NEnergReal
1000<BR>NEnergImCircle
50
#50<BR>NEnergImLine
40
#20<BR>NPoles
20
#10<BR>VInitial
0.0d0
eV<BR>VFinal
1.0d0
eV<BR>NIVPoints
100
<BR>#Delta
0.d-6<BR>EnergLowestBound -100.0d0
Ry<BR>SpinConfLeads
0<BR>NSlices
1<BR>TrCoefficients
F<BR>#AtomLeftVCte
3
<BR>#AtomRightVCte
10
<BR>NTransmPoints
100<BR>InitTransmRange -5.0 eV
<BR>FinalTransmRange 5.0
eV<BR>HartreeLeadsBottom -7.87933588
eV<BR>HartreeLeadsLeft
00.0
Ang<BR>HartreeLeadsRight 26.065
Ang</DIV>
<DIV> </DIV>
<DIV>#%block SaveBiasSteps<BR># 0 <BR>#%endblock
SaveBiasSteps<BR>SaveMemTranspK
T<BR>PeriodicTransp
T<BR>ParallelOverK
T<BR>BlockSize
8<BR>TransmissionOverk
F<BR>InitTransport
T</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV></BODY></HTML>