[DFTB-Plus-User] NEGF+XTB
Benjamin Hourahine
benjamin.hourahine at strath.ac.uk
Thu Apr 14 11:02:58 CEST 2022
Hello Luca,
you will need to regenerate the files with contact information before the transport can be run.
The process is described in
https://dftbplus-recipes.readthedocs.io/en/latest/transport/carbon2d-trans.html#contact-calculation
The contact calculations will need something like
KPointsAndWeights = SupercellFolding {
1 0 0
0 1 0
0 0 100
0.0 0.0 0.0
}
and both the electrostatics and Green's function blocks removed.
Regards
Ben
Dr. B. Hourahine, SUPA, Department of Physics,
University of Strathclyde, John Anderson Building,
107 Rottenrow, Glasgow G4 0NG, UK.
+44 141 548 2325, benjamin.hourahine at strath.ac.uk
2013/4 THE Awards Entrepreneurial University of the Year
2012/13 THE Awards UK University of the Year
The University of Strathclyde is a charitable body,
registered in Scotland, number SC015263
________________________________
From: DFTB-Plus-User <dftb-plus-user-bounces at mailman.zfn.uni-bremen.de> on behalf of luca.sementa at pi.ipcf.cnr.it <luca.sementa at pi.ipcf.cnr.it>
Sent: 14 April 2022 07:45
To: dftb-plus-user at mailman.zfn.uni-bremen.de <dftb-plus-user at mailman.zfn.uni-bremen.de>
Subject: [DFTB-Plus-User] NEGF+XTB
CAUTION: This email originated outside the University. Check before clicking links or attachments.
Dear Alessandro
the input file I used is reported here below. As you may check, this is
derived from the test file test/app/dftb+/transport/Au-chain .
The original one runs correctly; the modified one does not.
Best Regards
Luca
----------------------------------------------------------------------
geometry = genformat {
<<< "Au-chain.gen"
}
Transport{
Device{
AtomRange= 1 12
FirstLayerAtoms = 1 5 9
}
Contact{
Id = "source"
AtomRange= 13 20
Potential [eV] = 0.0
FermiLevel [eV] = -6.3051829983324215
}
Contact{
Id = "drain"
AtomRange= 21 28
Potential [eV] = 0.0
FermiLevel [eV] = -6.3051829983324215
}
#Task= ContactHamiltonian{
# contactId = "drain"
#}
}
Hamiltonian = xtb {
method = "GFN1-xTB"
filling = Fermi {
Temperature [K] = 0
}
KPointsAndWeights = SupercellFolding {
1 0 0
0 1 0
0 0 1
0.0 0.0 0.0
}
Electrostatics = Poisson {
#Poissonbox [Angstrom] = 20.0 20.0 20.0
MinimalGrid [Angstrom] = 0.4 0.4 0.4
AtomDensityCutoff [Angstrom] = 7.05
CutoffCheck = No
SavePotential = Yes
PoissonAccuracy = 1e-5
BuildBulkPotential = Yes
#RecomputeAfterDensity = Yes
}
#!Electrostatics = gammafunctional{}
Eigensolver = GreensFunction{
Verbosity = 51
RealAxisPoints = 100
EnclosedPoles = 3
}
#!Eigensolver = DivideAndConquer{}
}
Analysis{
TunnelingAndDOS{
EnergyRange [eV] = -0.0 0.0
EnergyStep [eV] = 0.02
}
}
!Analysis{}
Parallel = {
Blacs{
BlockSize = 4
}
# Allow OMP threads explicitely to test for hybrid parallelisation
with
# MPI-binary. (Check the manual before using this in production runs!)
UseOmpThreads = Yes
}
options = {
writeautotesttag = yes
writeChargesAsText = yes
readChargesAsText = yes
}
parseroptions = {
parserversion = 6
}
_______________________________________________
DFTB-Plus-User mailing list
DFTB-Plus-User at mailman.zfn.uni-bremen.de
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.zfn.uni-bremen.de%2Fcgi-bin%2Fmailman%2Flistinfo%2Fdftb-plus-user&data=04%7C01%7Cbenjamin.hourahine%40strath.ac.uk%7C56f7c397a5414cd9749608da1de27a48%7C631e0763153347eba5cd0457bee5944e%7C0%7C0%7C637855156231668951%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Iu5IjywSec0Wy7OHe8t8WZ%2BxRh3hOnP9LhYhCbG7JpA%3D&reserved=0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20220414/6216631b/attachment.htm>
More information about the DFTB-Plus-User
mailing list