[QE-users] epsilon.x with nspin=2 and "bad bands number"
mkondrin at hppi.troitsk.ru
mkondrin at hppi.troitsk.ru
Tue Mar 31 12:51:43 CEST 2020
Dear QE Developers and Users!
When launching epsilon.x after scf run with nspin=2 I have got an error
"bad bands number" although the same calculations for spin unpolarized
calculations completes OK. My QE version is 6.4. The example input files
are below.
I have found that increasing nbnd parameter to twice as large number (that
is nbnd=240 in the example below) aloows one to avoid this error. So I
simply commented out the check of bands number in the epsilon.x code. This
also solves my problem.
Sincerely yours,
M. Kondrin
#scf.in
&control
title='GeV3H2 TDDFT props',
calculation='scf',
prefix='gev3h2',
disk_io='low',
pseudo_dir = '../../../../QE/pseudo',
outdir='./tmp'
/
&SYSTEM
ibrav = 0
A = 7.56616
nat = 53
ntyp = 3
tot_charge=-2,
ecutwfc=70,
nspin=2,
tot_magnetization=0,
occupations='fixed',
nbnd=120
/
&electrons
mixing_beta = 0.7,
conv_thr = 1.0d-6
/
&ions
/
&cell
/
ATOMIC_SPECIES
Ge 73. Ge.pbe-mt_fhi.UPF
C 12. C.pbe-mt_fhi.UPF
H 1.01 H.pbe-mt_fhi.UPF
K_POINTS {Gamma}
CELL_PARAMETERS (alat= 14.29797028)
1.017268700 0.000966298 -0.001276172
0.509471188 0.880497388 -0.001276172
0.507868198 0.293217841 0.838703908
ATOMIC_POSITIONS (crystal)
Ge 0.030176730 0.053156603 0.041666667
C 0.339908341 -0.004633174 -0.011768993
C 0.660599947 -0.006910063 -0.006729167
C -0.000457701 0.001992166 0.668790467
C 0.663957422 -0.007183033 0.353505398
C 0.002429970 0.661897060 0.336749457
C 0.333251116 0.667546805 0.332866839
C 0.666870377 0.667773394 0.332963326
C 0.335418462 0.331461365 0.000834511
C 0.663274875 0.354933130 -0.007069496
C 0.333102298 -0.001079778 0.330162200
C 0.002942686 0.661030547 0.670114003
C 0.004630479 0.338493709 0.657033452
C 0.331879118 0.332848285 0.666332211
C 0.665920813 0.333088498 0.667953085
C -0.013745535 -0.003011939 0.348514247
C 0.332436215 0.002337467 0.663352392
C 0.334816038 0.662930722 0.002745506
C 0.675328893 0.659398752 0.001209447
C 0.666876571 -0.000038910 0.667481514
C 0.002724484 0.670973412 0.004393570
C 0.333135665 0.666691617 0.667510819
C 0.667302097 0.666041411 0.667309037
C 0.001571177 0.334667054 0.333270037
C 0.336541956 0.330626528 0.331882294
C 0.665918954 0.332907549 0.332005117
C 0.412359921 0.080608849 0.078939763
C 0.086345273 0.097078869 0.734819086
C 0.744839624 0.078702855 0.426299882
C 0.080995867 0.750897118 0.419980941
C 0.416641716 0.750197668 0.415822514
C 0.750485049 0.751454216 0.417001122
C 0.090243396 0.422733386 0.090062500
C 0.423934581 0.408004440 0.082123887
C 0.728400204 0.420058458 0.090402830
C 0.422302786 0.080390647 0.413219331
C 0.084413111 0.750231035 0.753171133
C 0.090516367 0.419375912 0.729827935
C 0.415559939 0.416462956 0.750370007
C 0.750425784 0.417414379 0.751328216
C 0.081341168 0.083791035 0.414542867
C 0.421436273 0.080903363 0.746583877
C 0.420402612 0.748517295 0.080587828
C 0.751871968 0.747914872 0.082498823
C 0.748666279 0.081762157 0.750063296
C 0.087966508 0.743424992 0.095102326
C 0.415786528 0.750082218 0.750466495
C 0.752706805 0.746791377 0.751451039
C 0.083372243 0.416456762 0.415851819
C 0.417291922 0.416031236 0.416024296
C 0.750244835 0.417412520 0.415380248
H 0.031525882 0.371888089 0.028431477
H 0.711445245 0.051807451 0.054901857
#epsilon.in
&inputpp
calculation = 'eps'
prefix='gev3h2',
outdir='./tmp'
/
&energy_grid
smeartype='gaussian'
intersmear=0.1
wmin=0.0
wmax=15.0
nw=1000
/
More information about the users
mailing list