<html><head>
<meta name="Generator" content="OpenText GroupWise Web">
</head><body><!DOCTYPE html><html><head><title></title><style>body{font-size: 14px;font-family: Helvetica,sans-serif;}</style></head><body style="overflow: auto;"><div class="GroupWiseMessageBody" id="GroupWiseSection_1774174461514_922B92C00DB9000097756E6937383562_"><div>Hi Lorenzo,</div><p>I believe the only way to do this is via wannierization. But let me know if you find other ways.</p><p>Cheers,<br>Lukas</p></div><div class="GroupWiseMessageBody" id="GroupWiseSection_1774174461498_922B92C00DB9000097756E6937383562_"><span class="GroupwiseReplyHeader"><br>>>> On 03/17/2026 at 12:43 AM, Lorenzo Sponza <lorenzo.sponza@onera.fr> <lorenzo.sponza@onera.fr> wrote:<br></span><p>Dear QE users,</p><p>Is it possible to get the matrix element < PSI | H_{SO} | PSI > or its modulus squared out of a Quantum ESPRESSO calculation or using one of its native post-processing tools? Here PSI is a well defined Bloch state of a given band index and k point and H_{SO} is the spin-orbit component of the DFT Hamiltonian. </p><p>Thanks in advance to whomever will answer. </p><p>Cheers</p><div id="signature">--<br><div class="pre" style="margin: 0; padding: 0; font-family: monospace;">Dr. Lorenzo Sponza<br><br>Chargé de Recherche au CNRS<br>Université Paris-Saclay, ONERA, CNRS Laboratoire d'étude des microstructures (LEM)<br>29 Avenue de la division Leclerc, 92322 Châtillon<br>Tel: +33146734464<br><br>ETSF Research Team Leader<br>European Theoretical Spectroscopy Facility<br><a href="https://sites.unimi.it/etsf/" target="_blank" rel="noopener noreferrer">https://sites.unimi.it/etsf/</a></div></div></div></body></html></body></html>