<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi all,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I would like to calculate the derivative of the KS orbitals with respect to k for the conduction states. The derivative of the occupied orbitals can be calculated using the linear response routine cgsolve_all.f90 in LR_Modules, which solves the Equation 30
of the paper <i>Rev. Mod. Phys., 73, 515-562</i>. For the conduction states, I refer to Eq. 25 of the same paper and adapt to the k derivative (with v the velocity operator):</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
( H - e_n ) dpsi_n/dk = -( 1 - |psi_n><psi_n| ) v |psi_n></div>
<div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I replaced d0psi of cgsolve_all with the right-hand side vector of the above equation and remove the alpha*Pv term. I managed to get the CG iterations converged but with much more steps (~100) than the default implementation for the occupied states. To check
whether the obtained first order wave functions are correct or not, I calculated some properties (orbital magnetization in my case) which only needs occupied orbitals. However, the method using the above equation produced different results compared with the
QE default implementation which agrees with literature. So, here is my question: is there anything wrong to use the above equation to calculate the wave function derivative? </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,</div>
<div id="Signature">
<div>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">--------------------------------------------</div>
<div style="font-family:Tahoma; font-size:13px">Xiaoming Wang
<div>Postdoctoral Research Associate</div>
<div>Department of Physics and Astronomy</div>
<div>University of Toledo</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>