[DFTB-Plus-User] SCC unconvergence problem
Bhatt, Mahesh
M.D.Bhatt at warwick.ac.uk
Thu Dec 6 18:36:00 CET 2012
Dear Vinay,
Thank you for suggestion. Now, the convergence problem has solved, but a new problem occurs when I tried to calculate HOMO and LUMO levels of dicyanovinyl sexithiphene (DCV6T). I got the following message:
WAVEPLOT 0.3
================================================================================
Interpreting input file 'waveplot_in.hsd'
--------------------------------------------------------------------------------
WARNING!
-> The following 3 node(s) had been ignored by the parser:
(1)
Path: waveplot/Basis/O
Line: 86-118 (File: wfc.mio-0-1.hsd)
(2)
Path: waveplot/Basis/P
Line: 172-219 (File: wfc.mio-0-1.hsd)
(3)
Path: waveplot/ParserOptions
Line: 30-32 (File: waveplot_in.hsd)
Processed input written as HSD to 'waveplot_pin.hsd'
Processed input written as XML to 'waveplot_pin.xml'
--------------------------------------------------------------------------------
Doing initialisation
Starting main program
Origin
-10.02280 -1.73955 -1.73695
Box
66.50498 0.00000 0.00000
0.00000 27.16901 0.00000
0.00000 0.00000 23.40492
Spatial resolution [1/Bohr]:
0.75182 1.84033 2.13630
*** glibc detected *** ./waveplot.x86_64-linux: malloc(): memory corruption: 0x00007f0e54021700 ***
======= Backtrace: =========
[0x8ded11]
[0x8e05fd]
[0x83e6eb]
[0x40bfcf]
[0x401e36]
[0x4002bc]
[0x8bb049]
[0x4001b9]
======= Memory map: ========
00400000-009db000 r-xp 00000000 00:31 2694277785 /storage/csc/msslax/dftbplus/mio-0-1/test/waveplot.x86_64-linux
00bda000-00bec000 rwxp 005da000 00:31 2694277785 /storage/csc/msslax/dftbplus/mio-0-1/test/waveplot.x86_64-linux
00bec000-00cbe000 rwxp 00000000 00:00 0 [heap]
4004b000-4004c000 ---p 00000000 00:00 0
4004c000-40444000 rwxp 00000000 00:00 0
406ef000-406f0000 ---p 00000000 00:00 0
406f0000-40ae8000 rwxp 00000000 00:00 0
40e4b000-40e4c000 ---p 00000000 00:00 0
40e4c000-40e5c000 rwxp 00000000 00:00 0
41a6c000-41a6d000 ---p 00000000 00:00 0
41a6d000-41e65000 rwxp 00000000 00:00 0
7f0e48143000-7f0e54000000 rwxp 00000000 00:00 0
7f0e54000000-7f0e54038000 rwxp 00000000 00:00 0
7f0e54038000-7f0e58000000 ---p 00000000 00:00 0
7f0e5a612000-7f0e5a939000 rwxp 00000000 00:00 0
7fffa76b9000-7fffa76ff000 rwxp 00000000 00:00 0 [stack]
7fffa77ff000-7fffa7800000 r-xp 00000000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted
My waveplot_in.hsd and wfc.mio-0-1.hsd files are as follows:
Options = {
TotalChargeDensity = Yes # Total density be plotted?
TotalChargeDifference = Yes # Total density difference plotted?
ChargeDensity = Yes # Charge density for each state?
RealComponent = Yes # Plot real component of the wavefunction
PlottedKPoints = { 1 }
PlottedSpins = { 1 -1 }
PlottedLevels = { 99 } # Levels to plot
PlottedRegion = OptimalCuboid {} # Region to plot
NrOfPoints = { 50 50 50 } # Number of grid points in each direction
NrOfCachedGrids = -1 # Nr of cached grids (speeds up things)
Verbose = Yes # Wanna see a lot of messages?
}
DetailedXML = "detailed.xml" # File containing the detailed xml output
# of DFTB+
EigenvecBin = "eigenvec.bin" # File cointaining the binary eigenvecs
# Definition of the basis
Basis = {
Resolution = 0.01
# Including mio-0-1.hsd. (If you use a set, which depends on other sets,
# the wfc.*.hsd files for each required set must be included in a similar
# way.)
<<+ "wfc.mio-0-1.hsd"
}
ParserOptions {
ParserVersion = 4
}
C = {
AtomicNumber = 6
Orbital = {
AngularMomentum = 0
Occupation = 2.000000
Cutoff = 5.0
Exponents = {
5.000000000000000e-01 1.140000000000000e+00 2.620000000000000e+00
6.000000000000000e+00
}
Coefficients = {
-5.171232639696000e-01 6.773263954720000e-02 -2.225281827092000e-03
1.308444510734000e+01 -5.212739736338000e+00 7.538242674175000e-01
-1.215154761544000e+01 -9.329029568076001e+00 -2.006616061528000e-02
-7.500610238649000e+00 -4.778512145112000e+00 -6.236333225369000e+00
}
}
Orbital = {
AngularMomentum = 1
Occupation = 2.000000
Cutoff = 5.0
Exponents = {
5.000000000000000e-01 1.140000000000000e+00 2.620000000000000e+00
6.000000000000000e+00
}
Coefficients = {
-2.302004373076000e-02 2.865521221155000e-03 -8.868108742828000e-05
3.228406687797000e-01 -1.994592260910000e-01 3.517324557778000e-02
1.328563289838000e+01 -7.908233500176000e+00 6.945422441225000e+00
-5.876689745586000e+00 -1.246833563825000e+01 -2.019487289358000e+01
}
}
}
Please let me know what is the problem here.
Mahesh
More information about the DFTB-Plus-User
mailing list