<html><head></head><body><div class="ydp941b1e78yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:16px;"><div dir="ltr" data-setdir="false"><div><span style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;">Bonjour, </span><br style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;"><br style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;"><span style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;">Je suis étudiant en master 2 de géologie et dans le cadre de mon stage de recherche de fin d'étude j'étude les transitions de certaines phases via des méthodes de calculs ab-initio. Actuellement je m'intéresse à la Reidite (ZrSiO4-tétragonal). Plus particulièrement, je suis en train d'effectuer des calculs de fréquences phononiques en calculant des matrices dynamiques et de constantes de forces pour plusieurs pressions (0 à 300 kbar avec un step de 20kbar). Cependant, je rencontre le problème suivant : </span><br style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;"><span style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;"><br></span></div><div><span style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;">Pour les pressions ; 6</span><span style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 16px;"><span style="color: rgb(0, 0, 0);">0kbar / 80kbar / 160kbar / 180kbar / 220kbar / 240kbar / 280kbar / 300kbar, l'erreur suivante "q-mesh breaks symmetry" s'affiche dans mes fichiers *ph.out et le calcul s'arrête (voir la première image). <br><br><div class="ydpa0c5b736img-preview-wrapper">Pour les autres pressions tout fonctionne très bien ...<br></div><br>Nous avons exploré plusieurs pistes mais sans succès (toujours la même erreur) : <br><br>- Utilisation de différentes versions de q.e (6.2 / 6.6 / 5.1.1 / 6.4.1 ).<br>- Utilisation de plusieurs grilles de q-points (333 / 332 / 444 / 442 / 555 / 552).<br>- Utiliser ibrav = 7 en reconstituant celldm(1) et celldm(3) au lieu de ibrav =0. <br><br>L'incompréhension réside principalement dans le fait que les calculs fonctionnent très bien pour certaines pressions et pas pour certaines. <br>Je mets dans le mail mon script "job.reidite-60kbarT_ib7" (celui qui calcule à P=60kbar), il contient à la fois les informations pour ibrav = 0 et ibrav = 7. <br><br><div class="ydpa0c5b736img-preview-wrapper">Je m'excuse d'avance si certains points ne sont pas clairs (c'est la première fois que je publie et que je travaille dans ce domaine). Je reste entièrement à disposition si vous avez la moindre question. <br></div><br>Je vous remercie de l'attention que vous porterez à mon problème. <br><br>Cordialement, <br><br>Pierre.</span></span></div><br></div></div></body></html>