[Pw_forum] equilibrium volume
Gabriele Sclauzero
sclauzer at sissa.it
Mon Dec 21 10:33:45 CET 2009
Dear Mohaddeseh
mohaddeseh abbasnejad wrote:
>
> Dear all
> When I run an scf calculation for the equilibrium volume obtained by the
> Murnaghan-fit, the stresses are too far from zero.
What does it mean too far? Please specify some numbers, so that we can have an idea of
what's happening. You should also specify some input file of your system.
> Should i expect this volume to be zero pressure volume, or there is a
> problem in the calculation?
I think it will be hardly exactly zero (for obvious numerical reasons) unless you specify
the exact equilibrium lattice constant with many significant digits and use very high
cutoffs (I think that stress is much harder to converge than total energy differences).
You should also keep in mind that the parameters resulting from a fit have an error bar
associated with them (though this might be quite small for the equilibrium lattice
constant, it may be enough to offset the stress a little bit from zero). Try to sample
more points around the minimum and fit with a parabola (such that the chi square of the
fit is very low).
Gabriele
> Thanx for any comment
>
> --
> ---------------------------------------------------------
>
> Mohaddeseh Abbasnejad,
> Room No. 323, Department of Physics,
> University of Tehran, North Karegar Ave.,
> Tehran, P.O. Box: 14395-547- IRAN
> Tel. No.: +98 21 6111 8634 & Fax No.: +98 21 8800 4781
> Cellphone: +989177317514
> E-Mail: m.abbasnejad at gmail.com <mailto:m.abbasnejad at gmail.com>
> Website: http://physics.ut.ac.ir
>
> ---------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
--
o ------------------------------------------------ o
| Gabriele Sclauzero, PhD Student |
| c/o: SISSA & CNR-INFM Democritos, |
| via Beirut 2-4, 34014 Trieste (Italy) |
| email: sclauzer at sissa.it |
| phone: +39 040 3787 511 |
| skype: gurlonotturno |
o ------------------------------------------------ o
More information about the Pw_forum
mailing list