phonons: getting a low enough force

Phonons, DFPT, electron-phonon, electric-field response, mechanical response…

Moderators: mverstra, joaocarloscabreu

Locked
brehmj
Posts: 35
Joined: Thu Jan 20, 2011 3:18 pm

phonons: getting a low enough force

Post by brehmj » Tue Jul 30, 2013 10:57 pm

I am stumped on what to do about the following:
1. I relax a structure of say 20 atoms to < 1x10-6 Ha/bohr (tolmxf 1.0d-6)

2. I use this relaxed structures acell and xred as inputs to the beginning of a phonon calculation that
has tolvrs 10d-18 .

3. when I get the forces back after the completion of the tolvrs,
they are bigger than 1x10-6 ha/bohr -- sometimes by 1-2 orders of magnitude.

My questions:
1. why is this?
2. what can I do to improve the forces for the tolvrs part?

tnx,
jb

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: phonons: getting a low enough force

Post by jzwanzig » Sun Aug 18, 2013 1:24 pm

Are you using the same ecutsm in your GS run for the phonons, that you used in the structural relaxation step? You should, it makes a big difference.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

Locked