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: Moderator, Global 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: 370
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: 370
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: 370
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: 370
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).

nidhi_tiwari
Newbie
Newbie
Posts: 19
Joined: Sat Nov 25, 2023 12:29 pm

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

#9 Post by nidhi_tiwari » Mon Jun 03, 2024 4:45 pm

Dear manoj_wijesingha,

I have also obtained similar result (three imaginary values out of six) in the frequency calculation of a single water molecule in a 7.34X7.34X20 Angstrom supercell, 4X4X1 Gamma KPOINTS, and EDIFF=1E-6.

I am really curious to know why an optimized molecule could show imaginary frequencies? I have attached the files of my calculation. Is there any mistake? If not, how can we rely of imaginary frequency values to distinguish transition states?
You do not have the required permissions to view the files attached to this post.

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

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

#10 Post by fabien_tran1 » Mon Jun 10, 2024 9:18 am

Among the 9 printed frequency modes of an isolated water molecule, only the first three are physically meaningful, see the following:
tutorials/latest/molecules/part3/#molecules-e10
https://www.chem.purdue.edu/jmol/vibs/h2o.html

Besides, the proper way to do a calculation on an isolated molecule is to choose only one k-point and to ensure that the cell is large enough to avoid interaction between the periodic images of the molecule.

nidhi_tiwari
Newbie
Newbie
Posts: 19
Joined: Sat Nov 25, 2023 12:29 pm

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

#11 Post by nidhi_tiwari » Wed Jun 12, 2024 1:23 pm

Thank you for the answer.

alex
Hero Member
Hero Member
Posts: 571
Joined: Tue Nov 16, 2004 2:21 pm
License Nr.: 5-67
Location: Germany

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

#12 Post by alex » Thu Jun 13, 2024 6:01 am

Hello, Fabian is right.

I might add some remarks if consider doing thermodynamics with partition sums caculated from these frequency analysis:
you need to do projection of the rotational and translational modes of your system (in case its gas phase, non-linear: 3N-6 degrees of freedom should ring a bell).
Some reading with the proper math: Wilson, JR., J. C. Decius, Paul C. Cross: Molecular Vibrations.

Or just use phonopy, it should also do the trick.

Best regards,

alex

Post Reply