Hessian Matrix calculation

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Locked
Message
Author
ali
Newbie
Newbie
Posts: 2
Joined: Mon Nov 23, 2009 7:27 pm

Hessian Matrix calculation

#1 Post by ali » Tue Nov 24, 2009 3:03 am

Hi, everyone.
I will appreciate it, if anyone answer my question.

I want to calculate the hessian matrix for a material with cubic symmetry like sc, bcc, fcc materials by vasp.4.6:

one of row of the matrix is as below:

1X 1Y 1Z 2X
1X 0.238518 0.000000 0.000000 2.425284

Some of the parameter that i used in my INCAR file:

IBRION = 5
NFREE = 2
POTIM = 0.015


I used the KPOINTS:

Line-mode
cart
0 0 0 ! gamma
0 0 1 ! X

I have some problems with the result:

1 - The sign of the value of some elements like main diagonal (1X-1X) are sensitive to choose the KPOINTS,how should i select my KPOINTS if i want to obtain all diagonal elemnts positive?

2 - I don't know why elements such as (1X-1Y or 1X-2Z, etc)
are zero?

Thanks in advance.
Last edited by ali on Tue Nov 24, 2009 3:03 am, edited 1 time in total.

panda

Hessian Matrix calculation

#2 Post by panda » Tue Nov 24, 2009 11:50 pm

Can't you just calculate it numerically? I am confused as to what you are asking, but here is some information on how to calculate the Hessian

http://www.dreamincode.net/forums/showtopic53310.htm
Last edited by panda on Tue Nov 24, 2009 11:50 pm, edited 1 time in total.

physicalattraction
Newbie
Newbie
Posts: 44
Joined: Wed May 21, 2008 1:47 pm
License Nr.: 916
Location: Eindhoven, The Netherlands

Hessian Matrix calculation

#3 Post by physicalattraction » Wed Nov 25, 2009 7:44 am

1) For the calculation of the Hessian matrix, VASP moves the atoms a certain distance (POTIM) out of their equilibrium position, relaxes the electronic system, and calculates the electronic forces on the atoms. During this electronic relaxation, integration over the full k-space is required. That is why I recommend you not to use the line mode in making a k-grid, which is more suited for band structure calculations, but using a Monkhorst-Pack grid. I don't know of which size, but you should be sure that the equilibrium is converged with respect to the number of k-points before calculating the Hessian matrix.
2) The 1X-2Z element is the force on atom 2 in the z-direction if atom 1 is moved in the x-direction. You could imagine that these forces are nearly zero, depending on your system, but the fact that is is exactly zero is worrisome. Maybe it's due to the line mode KPOINTS file.
Last edited by physicalattraction on Wed Nov 25, 2009 7:44 am, edited 1 time in total.

panda

Hessian Matrix calculation

#4 Post by panda » Wed Nov 25, 2009 7:08 pm

Oh, sorry, for some reason I was thinking they wanted to do it by hand to compare with VASP, I understand the question now, thanks for clarifying.
Last edited by panda on Wed Nov 25, 2009 7:08 pm, edited 1 time in total.

Locked