[Pw_forum] VC-relax questions
Timothy Mason
thmason24 at yahoo.com
Thu Dec 18 17:00:04 CET 2008
Gentlemen and ladies,
I have some rookie questions about VC-relax. I've pasted the input file and greps of Total energy, Total Forces and stresses below. Its the same problem of oscillations that has been posted before. Here are my questions.
1). how can it be that while the forces and stresses are oscillating the total energy continues to decrease albeit slowly? Am I really relaxing fine but just have a long way to go?
2) I've read many posts suggestioning an MD run to determine an effective dt but I don't think I understand. I run molecular dynamics and look for the oscillation period of the total forces???? Then use a dt equal to some fraction of the period?
3) I've seen conflicting suggestions on whether the default cell mass is a good guess. any further thoughts on this?
4) I'm starting off with "! no symmetries" and hoping to relax to a point where I find symmetries. is VC-relax even appropriate for this application?
thanks in advance!
Tim Mason
Graduate student physics
University of Missouri St. Louis
! total energy = -130.99725337 Ry
! total energy = -130.99742389 Ry
! total energy = -130.99765896 Ry
! total energy = -130.99788028 Ry
! total energy = -130.99814044 Ry
! total energy = -130.99838484 Ry
! total energy = -130.99862662 Ry
! total energy = -130.99883535 Ry
! total energy = -130.99906028 Ry
! total energy = -130.99921630 Ry
! total energy = -130.99942817 Ry
! total energy = -130.99962083 Ry
! total energy = -130.99984032 Ry
! total energy = -131.00005894 Ry
! total energy = -131.00027167 Ry
! total energy = -131.00032697 Ry
! total energy = -131.00055195 Ry
! total energy = -131.00075005 Ry
! total energy = -131.00092926 Ry
! total energy = -131.00110579 Ry
! total energy = -131.00128893 Ry
! total energy = -131.00152437 Ry
! total energy = -131.00176666 Ry
! total energy = -131.00201744 Ry
! total energy = -131.00222657 Ry
Total force = 0.025775 Total SCF correction = 0.000281
Total force = 0.020098 Total SCF correction = 0.000361
Total force = 0.015661 Total SCF correction = 0.000081
Total force = 0.016340 Total SCF correction = 0.000253
Total force = 0.012358 Total SCF correction = 0.000853
Total force = 0.012918 Total SCF correction = 0.000221
Total force = 0.011349 Total SCF correction = 0.000287
Total force = 0.013804 Total SCF correction = 0.000249
Total force = 0.012754 Total SCF correction = 0.000124
Total force = 0.016969 Total SCF correction = 0.000230
Total force = 0.011474 Total SCF correction = 0.000262
Total force = 0.013177 Total SCF correction = 0.000210
Total force = 0.010928 Total SCF correction = 0.000265
Total force = 0.009784 Total SCF correction = 0.000097
Total force = 0.010154 Total SCF correction = 0.000534
Total force = 0.023392 Total SCF correction = 0.000539
Total force = 0.018644 Total SCF correction = 0.000303
Total force = 0.013266 Total SCF correction = 0.000284
Total force = 0.012565 Total SCF correction = 0.000137
Total force = 0.011408 Total SCF correction = 0.000406
Total force = 0.013103 Total SCF correction = 0.000212
Total force = 0.011107 Total SCF correction = 0.000195
Total force = 0.009844 Total SCF correction = 0.000492
Total force = 0.010155 Total SCF correction = 0.000692
Total force = 0.014920 Total SCF correction = 0.000309
total stress (Ry/bohr**3) (kbar) P= 0.85
total stress (Ry/bohr**3) (kbar) P= 1.09
total stress (Ry/bohr**3) (kbar) P= 1.44
total stress (Ry/bohr**3) (kbar) P= 2.07
total stress (Ry/bohr**3) (kbar) P= 2.36
total stress (Ry/bohr**3) (kbar) P= 2.74
total stress (Ry/bohr**3) (kbar) P= 2.96
total stress (Ry/bohr**3) (kbar) P= 2.95
total stress (Ry/bohr**3) (kbar) P= 2.92
total stress (Ry/bohr**3) (kbar) P= 3.35
total stress (Ry/bohr**3) (kbar) P= 3.21
total stress (Ry/bohr**3) (kbar) P= 3.13
total stress (Ry/bohr**3) (kbar) P= 2.91
total stress (Ry/bohr**3) (kbar) P= 2.85
total stress (Ry/bohr**3) (kbar) P= 2.61
total stress (Ry/bohr**3) (kbar) P= 1.46
total stress (Ry/bohr**3) (kbar) P= 1.75
total stress (Ry/bohr**3) (kbar) P= 2.29
total stress (Ry/bohr**3) (kbar) P= 2.03
total stress (Ry/bohr**3) (kbar) P= 1.97
total stress (Ry/bohr**3) (kbar) P= 2.01
total stress (Ry/bohr**3) (kbar) P= 2.26
total stress (Ry/bohr**3) (kbar) P= 2.44
total stress (Ry/bohr**3) (kbar) P= 2.66
total stress (Ry/bohr**3) (kbar) P= 2.59
&CONTROL
calculation = "vc-relax",
prefix = 'n1',
tprnfor = .TRUE.
pseudo_dir = '/home/tim/PWSCF/PSPS',
outdir='/home/tim/tmp'
etot_conv_thr = 1.0D-5,
forc_conv_thr = 0.0004,
dt = 20
nstep = 100
/
&SYSTEM
ibrav = 0, nat = 15, ntyp = 5,
ecutwfc = 40
ecutrho = 480
celldm(1) = 1.889725989
/
&ELECTRONS
mixing_mode = 'plain',
mixing_beta = 0.7,
conv_thr = 1.0D-7,
/
&IONS
ion_dynamics = 'damp',
ion_damping = 0.1,
ion_positions = 'from_input',
pot_extrapolation = "second_order",
wfc_extrapolation = "second_order",
upscale = 100,
/
&CELL
press = 0.0,
cell_dynamics = 'damp-w'
/
ATOMIC_SPECIES
Li 1.0 Li.pbe-n-van.UPF
Ca 1.0 Ca.pbe-nsp-van.UPF
N 1.0 N.pbe-rrkjus.UPF
H 1.0 H.pbe-rrkjus.UPF
B 1.0 B.pbe-n-van.UPF
CELL_PARAMETERS
11.03225111 0.49197118 -0.92807585
-0.15038487 5.49373010 -0.38883048
3.04894498 0.25407718 4.19633759
ATOMIC_POSITIONS {crystal}
Li 0.144228066000 0.155479905000 0.105275909000
Ca 0.911201719000 0.659605393000 0.179405853000
N 0.954182418000 0.066667349000 0.203452528000
H 0.879987796000 0.160521652000 0.398773311000
H 0.949093887000 0.150685763000 0.033727675000
H 0.590246767000 0.244749224000 0.483311768000
H 0.563997731000 0.067035740000 0.225636766000
H 0.437455407000 0.306433941000 0.454794956000
H 0.611216941000 0.379448807000 0.154202770000
H 0.231109967000 0.728337975000 0.002905585000
H 0.397236924000 0.644804127000 -0.037067236000
H 0.288752224000 0.877300760000 0.238936122000
H 0.231627436000 0.575300194000 0.301639728000
B 0.550495441000 0.249323405000 0.329678852000
B 0.288025902000 0.707564656000 0.123296526000
K_POINTS automatic
4 4 4 0 0 0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20081218/f7ad6bc0/attachment.html>
More information about the users
mailing list