Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

Questions will be moved to this forum when we consider them out of scope for support from our side: for instance when we do not have the necessary expertise to come up with an answer.
Another user still might, though!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
manoj_wijesingha
Newbie
Newbie
Posts: 19
Joined: Wed Sep 21, 2022 6:52 pm

Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#1 Post by manoj_wijesingha » Mon Jan 30, 2023 7:06 pm

Hi,
I need to get the vibrational frequency calculation for the Metal doped Graphene. After the calculations, there are four imaginary frequencies, and some imaginary frequencies are higher (1.028868 cm-1, 1.334421 cm-1, 3.595799 cm-1, 125.990998 cm-1).

Then I wanted to run for H2O molecule as a test frequency calculation. The results are the same as the Metal doped Graphene. I could run the vibrational frequency calculation for H2O molecule after the relax calculation. But there are three imaginary frequencies for this H2O molecule. I don’t understand why is showing the imaginary frequencies for the optimized minimum structure.
Could you please explain me and figure out and what should I do for the frequency calculation.

I have attached my files of H2O.
Thank you.

INCAR:
System = FREQ
ISTART = 0
ICHARG = 2
PREC=Accurate
EDIFF=1E-6

ENMAX = 450
ISMEAR = 0
SIGMA = 0.05
LREAL=Auto

IBRION = 5
ISIF = 2
NSW = 1
EDIFFG = -0.01
ALGO=Fast
NFREE = 2
POTIM=0.005

KPOINTS:
KPoints-point
0
Monkhorst Pack
3 3 1
0 0 0

POSCAR
H2O _2
1.0000000000000000
8.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 8.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 8.0000000000000000
O H
1 2
Direct
-0.0000340711360937 0.0000000000000000 0.0000000000000000
0.0745271855680487 0.9040070284914102 0.0000000000000000
0.0745271855680487 0.0959929715085898 0.0000000000000000


Imaginary frequencies of first frequency calculation (H2O):
5.875007 cm-1, 148.058529 cm-1, 202.491147 cm-1


Then I increased EDIFF to 1E-8. After that the imaginary frequency (H2O) are;
3.087036 cm-1, 53.908671 cm-1, 77.342731 cm-1

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#2 Post by fabien_tran1 » Mon Jan 30, 2023 8:41 pm

Hi,

Could you also provide the OUTCAR of your calculation and mention which PAW dataset (POTCAR) you used? Have you used Phonopy for calculating the vibrational frequencies?

manoj_wijesingha
Newbie
Newbie
Posts: 19
Joined: Wed Sep 21, 2022 6:52 pm

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#3 Post by manoj_wijesingha » Mon Jan 30, 2023 9:16 pm

Hi,
Thank you so much and really appreciate your response. I have attached the OUTCAR files and POTCAR files. I used potpaw_PBE.54 dataset for the POTCAR. I did not use Phonopy for the frequency calculation.
Could you please assist me for the calculation of vibrational frequency.

Thank you.
Files.zip
You do not have the required permissions to view the files attached to this post.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#4 Post by fabien_tran1 » Mon Jan 30, 2023 10:38 pm

For the moment, I can mention two things:

1) In principle, the calculation of phonon frequencies with finite differences should be done with a supercell that is large enough to describe properly long-wavelength phonons. With a single cell, only the frequencies at the GAMMA point may be correct, while at other points the values may be (wrongly) imaginary. See, e.g., https://iopscience.iop.org/article/10.1 ... 075/ac78b3

2) The calculation of phonon frequencies requires rather accurate settings. As you observed, there is a rather large difference between EDIFF=1E-6 and 1E-8.

manoj_wijesingha
Newbie
Newbie
Posts: 19
Joined: Wed Sep 21, 2022 6:52 pm

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#5 Post by manoj_wijesingha » Tue Jan 31, 2023 12:31 am

Hi,
Thank you for your informative response. So, do you suggest Phonophy for the vibrational frequency calculations and is this better for the large super cell (5x5x1) calculation?
Thank you.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#6 Post by fabien_tran1 » Tue Jan 31, 2023 7:49 am

Not necessarily. For the moment, the important point is to consider a supercell. See also Phonons_from_finite_differences.

manoj_wijesingha
Newbie
Newbie
Posts: 19
Joined: Wed Sep 21, 2022 6:52 pm

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#7 Post by manoj_wijesingha » Tue Jan 31, 2023 2:42 pm

Hi,
Thank you for your guidance. If the supercell is large, how we can calculate the correct frequencies and how we can identify the structures are minima or transition structures? Could you please explain what should I do for the frequency calculation for the super cell 5x5x1.
This test H2O molecule is very simple calculation but it gives imaginary frequencies. Could you please suggest for me the process to compute frequency calculations to identify the minima and the transition states.

Thank you.

fabien_tran1
Global Moderator
Global Moderator
Posts: 316
Joined: Mon Sep 13, 2021 11:02 am

Re: Vibrational Frequency Calculation: Imaginary Frequencies from Optimized Minimum structures (H2O)

#8 Post by fabien_tran1 » Tue Jan 31, 2023 3:03 pm

The topic has been moved to the forum "From Users for Users". Actually, now this is more about learning (from literature and documentation) than solving technical problems (.e.g., bugs).

Post Reply