<div dir="ltr"><div><div><div><div>To whom it may concern,<br></div><div><br></div>I am attempting to use the stress tensor to compute the gradient of my unit cell. <br><br></div><div>Looking in your source code I noticed in move_ions.f90 you have the following code around line 153:<br><br>IF ( lmovecell ) THEN<br>  CALL cell_force( fcell, - transpose(bg)/alat, sigma, omega, press )<br>END IF<br><br>At this point in your code you are in reciprocal space correct (since you are using bg which is your reciprocal lattice vectors)? If so I'm operating in direct space so I was wondering if you had any insight into how I would compute this same quantity in cartesian? Also when I turn on tstress = true is the stress tensor that is printed out in Cartesian space or in Crystal?<br></div><div><br></div>Best,<br></div>Jessica McKinley<br></div>UC Riverside<br></div>