[Pw_forum] same Problem with BFGS

mnunez at ncsu.edu mnunez at ncsu.edu
Tue Jul 11 13:49:03 CEST 2006


I am getting the same bfgs error posted before:
-
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from bfgs : error #         1
     bfgs history already reset at previous step
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

stopping ...
--------------------
the bfgs algorithm seems to get stucked.
the total forces along the minimization  before it stops are:

     Total force =     0.002868     Total SCF correction =     0.034876
     Total force =     0.042610     Total SCF correction =     0.039637
     Total force =     0.008934     Total SCF correction =     0.029486
     Total force =     0.003953     Total SCF correction =     0.016921
     Total force =     0.012592     Total SCF correction =     0.010567
     Total force =     0.004110     Total SCF correction =     0.009995
     Total force =     0.001633     Total SCF correction =     0.012211
     Total force =     0.001412     Total SCF correction =     0.019146
     Total force =     0.001607     Total SCF correction =     0.021538
     Total force =     0.001555     Total SCF correction =     0.012259
     Total force =     0.001361     Total SCF correction =     0.012868
     Total force =     0.001432     Total SCF correction =     0.013882
     Total force =     0.001343     Total SCF correction =     0.011201
     Total force =     0.001422     Total SCF correction =     0.013635
     Total force =     0.001369     Total SCF correction =     0.005947
     Total force =     0.001317     Total SCF correction =     0.011391
     Total force =     0.001485     Total SCF correction =     0.006279
     Total force =     0.001351     Total SCF correction =     0.014118
     Total force =     0.001475     Total SCF correction =     0.011927

the input file
&control
    calculation = "relax",
    title='scf BaOPtAl big slab 1 001'
    restart_mode="from_scratch",
    prefix="BaOPtAlSlabMeth",
    pseudo_dir ="/spin/home/mnunez/pwcheetah/pseudo/",
    outdir="/tmp/work/mnunez/BaOPtAlSlabMethod/",
    forc_conv_thr = 1.0d-4
    etot_conv_thr = 1.0d-5
    tstress = .true.,
    tprnfor = .true.
 /
 &system
    ibrav = 6,
    celldm(1) =7.297,
    celldm(3) = 16.9397,
    nat = 42,
    ntyp = 4,
    nbnd=230,
    ecutwfc =35,
    ecutrho=245,
    occupations ='smearing',
    smearing='gaussian',
    degauss=0.008
 / &electrons
    mixing_mode= 'local-TF'
    diagonalization= 'cg'
    mixing_beta = 0.3,
    mixing_ndim = 11,
    conv_thr =  1.0d-4
 /
 &ions
    upscale=11
 /
ATOMIC_SPECIES
O 15.9994  O.pz-rrkjus.UPF
Ba 137.327 Ba.vdb.UPF
Pt 195.078  Pt.pz-rrkjus.UPF
Al 26.981538  Al.pz-vbc.UPF
ATOMIC_POSITIONS
Pt       0.500000000   0.500000000  -5.843889151
Pt       0.000000000   0.000000000  -5.841443397
Pt       0.500000000   0.000000000  -5.341181221
Pt       0.000000000   0.500000000  -5.341181221
Pt       0.500000000   0.500000000  -4.831963549
Pt       0.000000000   0.000000000  -4.829809502
Pt       0.500000000   0.000000000  -4.316032074
Pt       0.000000000   0.500000000  -4.316032074
Pt       0.500000000   0.500000000  -3.800183390
Pt       0.000000000   0.000000000  -3.799129595
Pt       0.500000000   0.000000000  -3.250938002
Pt       0.000000000   0.500000000  -3.250938002
Al       0.500000000   0.500000000  -2.892741738
Al       0.000000000   0.000000000  -2.774954863
Ba       0.500000000   0.500000000  -2.105588663
O        0.000000000   0.000000000  -2.318966082
Ba       0.000000000   0.000000000  -1.375247423
O        0.500000000   0.500000000  -1.434073949
Ba       0.500000000   0.500000000  -0.681716547
O        0.000000000   0.000000000  -0.728523986
Ba       0.000000000   0.000000000   0.025502608
O        0.500000000   0.500000000  -0.009862356
Ba       0.500000000   0.500000000   0.726855678
O        0.000000000   0.000000000   0.705024595
Ba       0.000000000   0.000000000   1.435097711
O        0.500000000   0.500000000   1.393104543
Ba       0.500000000   0.500000000   2.115982686
O        0.000000000   0.000000000   2.118086625
Pt       0.500000000   0.500000000   2.836917653
Pt       0.000000000   0.000000000   2.696183926
Pt       0.500000000   0.000000000   3.298956976
Pt       0.000000000   0.500000000   3.298956976
Pt       0.500000000   0.500000000   3.810364151
Pt       0.000000000   0.000000000   3.819965483
Pt       0.500000000   0.000000000   4.332138351
Pt       0.000000000   0.500000000   4.332138351
Pt       0.500000000   0.500000000   4.850766432
Pt       0.000000000   0.000000000   4.844433414
Pt       0.500000000   0.000000000   5.360902198
Pt       0.000000000   0.500000000   5.360902198
Pt       0.500000000   0.500000000   5.860968442
Pt       0.000000000   0.000000000   5.865188929
K_POINTS {automatic}
 6 6 1 1 1 0


the pseudopotentials used are the ones in

http://www.pwscf.org/pseudo.htm






More information about the users mailing list