[DFTB-Plus-User] Failure in diagonalisation routine DSYGVR
Jungyoon Kim
jyoonkim0403 at gmail.com
Wed May 29 07:15:06 CEST 2024
Hi all,
I'm currently trying to do SCC-DFTB NVT MD simulations.
Unfortunately, the calculation was failed with this error :
--------------------------------------------------------------------------------
*** Geometry step: 2
iSCC Total electronic Diff electronic SCC error
1 -0.23259089E+03 0.00000000E+00 0.53764967E+00
2 -0.27063312E+03 -0.38042230E+02 0.43113745E+00
3 -0.27064824E+03 -0.15120311E-01 0.12942984E+00
4 -0.27064266E+03 0.55771487E-02 0.11346507E+00
5 -0.27065918E+03 -0.16524719E-01 0.78607192E-01
6 -0.27066028E+03 -0.10932004E-02 0.42746437E-01
7 -0.27066044E+03 -0.16809919E-03 0.25335717E-01
8 -0.27066047E+03 -0.26514960E-04 0.58747631E-02
9 -0.27066040E+03 0.73216544E-04 0.14439062E-02
10 -0.27066034E+03 0.62802498E-04 0.80454294E-03
11 -0.27066033E+03 0.84361088E-05 0.33010731E-03
12 -0.27066034E+03 -0.99345070E-05 0.24382867E-03
13 -0.27066034E+03 -0.22057095E-05 0.51966315E-04
14 -0.27066034E+03 0.17399407E-05 0.27541209E-04
15 -0.27066034E+03 -0.53144510E-06 0.84890399E-05
Total Energy: -266.0198878865 H -7238.7695 eV
Extrapolated to 0K: -266.0203072997 H -7238.7809 eV
Total Mermin free energy: -266.0207267129 H -7238.7923 eV
Force related energy: -266.0207267129 H -7238.7923 eV
Volume: 0.674833E+04 au^3 0.100000E+04 A^3
Electronic Temperature: 0.0031668153 H 1000.0000 K
MD Temperature: 0.3078518889 H 97211.8218 K
MD Kinetic Energy: 32.7862261702 H 892.1586 eV
Total MD Energy: -233.2345005427 H -6346.6337 eV
Pressure: -0.115897E-03 au -0.340981E+10 Pa
Gibbs free energy including KE: -233.2345005427 H -6346.6337 eV
--------------------------------------------------------------------------------
*** Geometry step: 3
iSCC Total electronic Diff electronic SCC error
ERROR!
-> Failure in diagonalisation routine DSYGVR, unable to complete Cholesky
factorization of B 368
I wonder why this error came out and if there is any workaround.
Here is my input file. :
Geometry = GenFormat {
<<< "2tp-test-5.gen"
}
Hamiltonian = DFTB {
SCC = Yes
MaxAngularMomentum = {
H = "s"
C = "p"
N = "p"
S = "d"
}
SlaterKosterFiles = Type2FileNames {
Prefix = "/home/jyoon/DFTBp_input/SK-File/mio-1-1/"
Separator = "-"
Suffix = ".skf"
LowerCaseTypeName = No
}
SpinPolarization = {
collinear = {
UnpairedElectrons = 2
}
}
SpinConstants = {
H {
-0.0717
}
C {
-0.0227
}
N {
-0.0255
}
S {
-3.1592
}
}
KPointsAndWeights = {
0 0 0 1.0
}
# Filling = Fermi {
# Temperature [Kelvin] = 2000
# }
}
Driver = VelocityVerlet{
TimeStep [fs] = 0.5
Thermostat = NoseHoover {
Temperature [Kelvin] = 1000
ChainLength = 10
CouplingStrength [cm^-1] = 3200
AdaptFillingTemp = Yes
}
Steps = 20000
MovedAtoms = 1:-1
MDRestartFrequency = 5
OutputPrefix = "2tp-test-5-out"
}
Options = {
RandomSeed = 0
}
Analysis { CalculateForces = Yes }
ParserOptions { ParserVersion = 10 }
Parallel { UseOmpThreads = Yes }
Thank you.
Best regards,
Jungyoon Kim
답장전달
반응 추가
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20240529/99c871cb/attachment.htm>
More information about the DFTB-Plus-User
mailing list