<div dir="ltr">Dear QE users and developers,<div><br></div><div>I am trying to use the real space force constants file obtained from q2r.x to do some other calculations outside of QE. For this purpose, I am trying to understand the format of the force constants file. I am attaching a sample force constants file obtained by performing ph.x and q2r.x on bulk silicon with a 4 x 4 x 4 q-mesh. </div>

<div><br></div><div>In the 22nd line in the file, it says 1 1 1 1. I believe the first two 1's denote the x-direction and the last two 1's denote atom numbers within the two-atom unit cell of silicon. After this it loops over all 64 different combinations of indices in a 4 x 4 x 4 mesh and lists the force constants in the 4th column. What do these 64 different combinations mean? Do they denote unit cells? </div>

<div><br></div><div>For instance, let's take the first line: 1 1 1 2.6997. 2.6997 is the force constant between which atoms? </div><div><br></div><div>Thanks for your help!</div><div>Sridhar</div><div>Purdue University</div>

</div>