<div dir="ltr"><div>Thanks for the answer and the links. But if I understand correctly, these codes only give us the irreps, and I need the specific representation matrices calculated from the QE wave-functions. The best way to get these seem to be editing this routine. </div><div><br></div><div>I've found my mistake, and a bug (I'll report for the developers as well). My mistake was that I was not using the "which_irr_so(iclass)" to identify the class. I was assuming that "iclass" would already list the classes in the correct order.</div><div><br></div><div>Now, the bug is that "which_ir_so" is returning the class 9 two times here, but one of these should be 6 instead. This can be checked with the following code</div><div><br></div><div><div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:"Droid Sans Mono","monospace",monospace;font-size:14px;line-height:19px;white-space:pre-wrap"><div>  <span style="color:rgb(86,156,214)">PRINT</span> *, <span style="color:rgb(206,145,120)">"======================================================"</span></div><div>  <span style="color:rgb(86,156,214)">PRINT</span> *, <span style="color:rgb(206,145,120)">"which_ir_so:"</span>, (which_irr_so(iclass), iclass=<span style="color:rgb(181,206,168)">1</span>, nclass)</div><div>  <span style="color:rgb(86,156,214)">PRINT</span> *, <span style="color:rgb(206,145,120)">"classes:"</span>, (name_class_so(iclass), iclass=<span style="color:rgb(181,206,168)">1</span>, nclass)</div><div>  <span style="color:rgb(86,156,214)">DO</span> irap=<span style="color:rgb(181,206,168)">1</span>,nrap</div><div>   <span style="color:rgb(86,156,214)">PRINT</span> *, <span style="color:rgb(206,145,120)">">>"</span>, (char_mat_so(irap, which_irr_so(iclass)), iclass=<span style="color:rgb(181,206,168)">1</span>, nclass)</div><div>  <span style="color:rgb(86,156,214)">ENDDO</span></div><div>  <span style="color:rgb(86,156,214)">PRINT</span> *, <span style="color:rgb(206,145,120)">"======================================================"</span></div></div></div><div><br></div><div>Which prints</div><div><br></div><div>which_ir_so(iclass):           1           5           3           9           9           7           2           4           8<br>classes: E    -E   2C3  -2C3  3C2' s_h 2S3  -2S3  3s_v<br></div><div><br></div><div>So, the classes are ordered as {E, 3C2', 2C3, 3s_v, 3s_v, 2S3, -E, -2C3, -2S3}</div><div><br></div><div>Notice that 3s_v appears twice, and s_h does not appear in the list.</div><div><br></div>Best,<br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="font-family:arial"><div>--<br>Gerson J. Ferreira</div><div>Prof. Dr. @ InFis - UFU </div><div>----------------------------------------------</div><div style="font-family:Arial,Helvetica,sans-serif"><font face="arial"><a href="http://gjferreira.wordpress.com/" target="_blank">gjferreira.wordpress.com</a></font><br></div><div><font size="2" face="arial">Institute of Physics</font></div><div>Federal University of Uberlândia, Brazil</div><div style="font-family:Arial,Helvetica,sans-serif"><span style="font-family:arial">----------------------------------------------</span></div></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Mar 6, 2022 at 10:35 PM Hongyi Zhao <<a href="mailto:hongyi.zhao@gmail.com">hongyi.zhao@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Mon, Mar 7, 2022 at 2:39 AM Gerson J. Ferreira<br>
<<a href="mailto:gersonjferreira@ufu.br" target="_blank">gersonjferreira@ufu.br</a>> wrote:<br>
><br>
> Dear QE users,<br>
><br>
> I need the matrix representation for the symmetry operators in the basis of the QE bands, so I'm checking how to edit the file sym_band.f90 for this purpose.<br>
><br>
> The SUBROUTINE find_band_sym_so already calculates the trace as (loops are implied)<br>
>><br>
>> trace(iclass,igroup)=trace(iclass,igroup) + DOT_PRODUCT (evc(:,ibnd),evcr(:,ibnd))<br>
><br>
><br>
> So, at first I imagined that a simple change would allow me to get the full matrices as (loops are implied)<br>
>><br>
>> matrep(iclass,igroup,i,j) = DOT_PRODUCT (evc(:,ibnd),evcr(:,jbnd))<br>
><br>
><br>
> But I've noticed that the traces are "wrong". If I print as<br>
>><br>
>> PRINT *, 'Class:', name_class_so(iclass)<br>
>> PRINT *, 'Trace:', trace(iclass,igroup)<br>
><br>
> Both the trace and the matrices (matrep) above don't match the expected results. For instance, all double group bar-irreps are showing trace = 0.<br>
><br>
> In the second part of this subroutine, where the code identifies the symmetry representations, I don't understand some of the IFs there, and the meaning of the variable "shift". So I guess I'm misreading something.<br>
><br>
> Could someone help me understand what I am doing wrong?<br>
<br>
First, the method implemented in "./PP/src/sym_band.f90" is just a<br>
rough partial implementation, as indicated by the comment in the above<br>
source code file:<br>
<br>
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!<br>
SUBROUTINE find_band_sym_so (ik,evc,et,nsym,s,ft,d_spin,gk, &<br>
     invs,rap_et,times,ngroup,istart,accuracy)<br>
<br>
  !<br>
  !   This subroutine finds the irreducible representations of the<br>
  !   double group which give the transformation properties of the<br>
  !   spinor wavefunctions evc.<br>
  !   Presently it does NOT work at zone border if the space group of<br>
  !   the crystal has fractionary translations (non-symmorphic space groups).<br>
  !<br>
  !<br>
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!<br>
<br>
2. As an alternative, if I understand correctly, the following<br>
packages are currently available for this purpose, as discussed here<br>
[1]:<br>
<br>
A.  <a href="https://github.com/qeirreps/qeirreps" rel="noreferrer" target="_blank">https://github.com/qeirreps/qeirreps</a><br>
B. <a href="https://github.com/stepan-tsirkin/irrep" rel="noreferrer" target="_blank">https://github.com/stepan-tsirkin/irrep</a><br>
<br>
[1] <a href="https://github.com/goodluck1982/SpaceGroupIrep/issues/6" rel="noreferrer" target="_blank">https://github.com/goodluck1982/SpaceGroupIrep/issues/6</a><br>
<br>
Regards<br>
-- <br>
Assoc. Prof. Hongsheng Zhao <<a href="mailto:hongyi.zhao@gmail.com" target="_blank">hongyi.zhao@gmail.com</a>><br>
Theory and Simulation of Materials<br>
Hebei Vocational University of Technology and Engineering<br>
No. 473, Quannan West Street, Xindu District, Xingtai, Hebei province<br>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
</blockquote></div>