aiida-quantumespresso icon indicating copy to clipboard operation
aiida-quantumespresso copied to clipboard

Ph Symmetry labels not correct when Raman active

Open mkotiuga opened this issue 1 year ago • 1 comments

The symmetry labels are parsed as the last field on lines with '-->' in the ph output:

https://github.com/aiidateam/aiida-quantumespresso/blob/5d4a7d9405b757e2ecf65a72c1ee92aa2fb36a39/src/aiida_quantumespresso/parsers/parse_raw/ph.py#L165

In cases where the mode is Raman active (thus the last field is "R" or "I+R" - if both Raman and IR active), the R/I+R label is kept and not the point symmetry label - Only the IR only case is treated. Note the "R" and "I" labels are only relevant for gamma point calcuations. https://github.com/aiidateam/aiida-quantumespresso/blob/5d4a7d9405b757e2ecf65a72c1ee92aa2fb36a39/src/aiida_quantumespresso/parsers/parse_raw/ph.py#L169

I suggest is just better to keep everything after the '-->' as it is all information related to the symmetry (labels in different conventions) so different users may have different preferences. i.e. something like:

symm_label = line.split('-->')[-1] replacing L169-172

running the following pw.in and ph.in should produce a system with R and I+R labels at gamma (the SSSP - with the updated Barium should be OK) ph.in.txt pw.in.txt

mkotiuga avatar Jul 21 '23 16:07 mkotiuga

Thanks @mkotiuga! I'll try to fix this asap.

mbercx avatar Jul 24 '23 14:07 mbercx