[Pw_forum] convergence of phonon

Nicola Marzari marzari at mit.edu
Sat Jan 15 08:48:50 CET 2005


Dear Eduardo,

in this order, I would consider

1) wfc cutoff. For a reasonable ultrasoft psp should be between
25 and 30 Ry, and unlikely to be more than 40 Ry.

2) dual (i.e. ecutrho). 6 to 12 times the wfc cutoff.

3) k-points - you might want to use m-p meshes containing gamma,
since the phonon q-mesh will also need to contain gamma (last
I checked, the q-to-r and r-to-q interpolators require that
you have gamma in your q mesh). Usually shifted k-point meshes
are slightly better for convergence, and I would use those if
only a total energy and force calculation was needed.

4) degauss. I tend to use ~0.5eV for the cold smearing (note that
in the code degauss is given in Ry). Something slightly larger
(up to 0.7eV) can be used for certain systems (the larger the
smearing, the fewer the k-points needed to achieve convergence).
Occasionally you might need to go down to ~0.3eV. It's system
dependent (how complex the fermi energy surface is, how parabolic
the dos is...). Note that this smearings are good for integrated
quantities (like the forces); if you needed to study a propery dependent
only on the fermi surface, you would need to be more carefull.

5) convergence in scf cycle. As accurately ocnverged as the code
would go....

6) convergence in the nscf/phonon. You might want to try a few tests here.

In addition, you need to make sure that all your atomic positions
are accuarately relaxed. Volume is often the theoretical volume,
but there are cases in which it might be wise to use the
experimental parameters (e.g. the c/a in graphite - we have
a recent preprint on arxiv/cond-mat on the matter).

At the end, you want a convergence in the frequencies with not more
than 1% error due to all the parameters above, so you can roughly
estimate (from sqrt(k/m), where k is the force constant, and
you know the experimental frequencies) the error on the forces that
you can tolerate.

Anyone on the list with additional suggestions is welcome to add them...

All the best,

			nicola

Eduardo Ariel Menendez P wrote:

> Dear colleagues,
> 
> I am doing a phonon calculation for a rather large system, with 12 atoms
> in the unit cell and is metallic (Ti_3Si_C). The phonon calculations are
> long. I started one 15 days ago and  for most q-points it is still running.
> I need some advice to perform the convergence check.
> 
> 1) Which parameters are relevant for a phonon calculation?
> I guess: number of k-points,degauss,ecut,ecutrho
> anyone else?
> 
> 2) What output should be checked for convergence in previous SCF
> calculations: the forces, or the energy is enough?
> ?What convergence of forces or energy can be considered safe for a phonon
> calculation? Considering the following data, which ecut can I use?
> 
> ecut(Ry)  total_E(Ry) 	P(kbar)  Force(Ry/au)	Force(eV/A)
> 
> 30	 -761.90543963 	-98.68 	-0.00675419 	-0.173648
> 40 	 -761.96493584 	-92.26 	-0.00574380 	-0.147672
> 50 	 -761.97745068 	-80.15 	-0.00621080 	-0.159678
> 60 	 -761.97911012 	-80.46 	-0.00611748 	-0.157279
> 70 	 -761.98131155 	-80.07 	-0.00611850 	-0.157305
> 80 	 -761.98190159 	-79.20 	-0.00615275 	-0.158186
> 90 	 -761.98206846 	-79.33 	-0.00613964 	-0.157849
> 100 	 -761.98234169 	-79.42 	-0.00613564 	-0.157746
> 
> Above, Force is taken as the z force component on a specific atom wich is
> not at the theoretical equilibrium position. P(ressure)
> is negative as I am using the experimental lattice constant
> (should I relax the cell volume?).
> 
> Thank you
> 
> 
> Eduardo A. Menendez Proupin
> Department of Physics
> Faculty of Science
> University of Chile
> Las Palmeras 3425
> Ñuñoa, Santiago
> Chile
> Phone: 56+2+678 74 11
> http://fisica.ciencias.uchile.cl/~emenendez/
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum

-- 
---------------------------------------------------------------------
Prof Nicola Marzari   Department of Materials Science and Engineering
13-5066   MIT   77 Massachusetts Avenue   Cambridge MA 02139-4307 USA
tel 617.4522758  fax 617.2586534  marzari at mit.edu  http://nnn.mit.edu



More information about the Pw_forum mailing list