Phonons, DFPT, electron-phonon, electric-field response, mechanical response…
Moderators: mverstra, joaocarloscabreu
-
disae
- Posts: 2
- Joined: Tue Oct 18, 2016 4:37 pm
Post
by disae » Fri Nov 18, 2016 5:00 pm
Dear everyone,
I'm a new abinit user and I'm trying to calculate the phonon DOS for a NiTi compound, but I'm not getting the same result as reported in the literature.
I just modified the input files in the tutorials to match my system. The general shape of the PDOS looks ok, but the frequencies seem to be shifted up by ~1 THz, and that branch near gamma looks bad.
Can someone help me to understand what's wrong with my model?
I'm using the norm-conserving pseudopotentials available in the abinit webpage.
Abinit output:
- My calculation
- pdos.png (41.03 KiB) Viewed 5195 times
Reported in literature:
- Reported in literature
- niti_aus.jpg (32.6 KiB) Viewed 5195 times
anaddb.in- anaddb input
- (4.3 KiB) Downloaded 413 times
Thanks
-
ilukacevic
- Posts: 271
- Joined: Sat Jan 16, 2010 12:05 pm
- Location: Dept. of Physics, University J. J. Strossmayer, Osijek, Croatia
-
Contact:
Post
by ilukacevic » Tue Nov 22, 2016 8:13 am
Dear disae,
1. First, this is not phonon DOS, but phonon dispersion curves.
2. The differences you see can be the result of many factors: approximation used in the calculation, unconverged calculation parameters, unrelaxed structure, and similar. There are many things you should check first, before trying to compare your results with some other. And you have to be specially careful if you are comparing to the experimental results. Frequency shifts of 1 THz are not unseen before. All in all, I would say that your results are good.
Best regards,
Igor