<div>Yes, when anharmonicity get stronger, the perturbative approach of d3q eventually breaks down. It is always possible to use stochastic (sscha) or molecular dynamics methods. Each method has trade-offs, in terms of: computational cost, ability to cope with strong anharmonicity, treatment of quantum ions/very low temperatures, and last but not least in my opinion, giving an insight to the physical processes.</div><br><div>regards</div><br><br><div><signature id="initial">Dr. Lorenzo Paulatto<br>IdR @ IMPMC - CNRS UMR 7590 & Sorbonne Université<br>phone: +33 (0)1 442 79822 / skype: paulatz<br><a href="https://link.getmailspring.com/link/6A135B8E-AE49-4ACA-A7D4-E5CDAF809154@getmailspring.com/0?redirect=http%3A%2F%2Fwww.impmc.upmc.fr%2F~paulatto%2F&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">http://www.impmc.upmc.fr/~paulatto/</a> - <a href="https://link.getmailspring.com/link/6A135B8E-AE49-4ACA-A7D4-E5CDAF809154@getmailspring.com/1?redirect=https%3A%2F%2Fanharmonic.github.io%2F&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">https://anharmonic.github.io/</a><br>23-24/423 B115, 4 place Jussieu 75252 Paris CX 05</signature></div><div class="gmail_quote_attribution">On Mar 31 2021, at 4:35 am, Terry Frankcombe <t.frankcombe@adfa.edu.au> wrote:</div><blockquote><div><div>While we are talking about interesting ways to include anharmonic</div><div>effects, it's worth taking a look at SSCHA: http://sscha.eu/</div><br><div>Not necessarily computationally inexpensive.</div><br><br><div>On Wed, 2021-03-31 at 07:20 +0900, Kazume NISHIDATE wrote:</div><div>> > 2021/03/31 1:00、Lorenzo Paulatto <paulatz@gmail.com>のメール:</div><div>> ></div><div>> > Hello you can take into account anharmonicity with the d3q code which is included in QE,</div><div>></div><div>></div><div>> This may be the most sophisticated way to account for the</div><div>> anharmonicity from first principle. I'm impress by the ingeniousness.</div><div>> Thank you for your notice.</div><div>></div><div>> btw-</div><div>> The methods described in my previous message are heavily relying upon</div><div>> a brute force manner. The only advantage of the method should be the</div><div>> automatic inclusion of higher order terms beyond the third-order, but</div><div>> of course its resolution, especially when it was obtained by MD, would</div><div>> be very low.</div><div>></div><div>></div><div>></div><div>></div><div>> best regards</div><div>> kazume NISHIDATE</div><div>> 敬具　西館数芽</div><div>></div><div>> nisidate@iwate-u.ac.jp</div><div>> kazume.nishidate@gmail.com</div><div>></div><div>></div><div>></div><div>></div><div>></div><div>></div><div>></div><div>></div><div>></div><div>></div><div>> _______________________________________________</div><div>> Quantum ESPRESSO is supported by MaX (www.max-centre.eu)</div><div>> users mailing list users@lists.quantum-espresso.org</div><div>> https://lists.quantum-espresso.org/mailman/listinfo/users</div><br><br><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)</div><div>users mailing list users@lists.quantum-espresso.org</div><div>https://lists.quantum-espresso.org/mailman/listinfo/users</div></div></blockquote><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/6A135B8E-AE49-4ACA-A7D4-E5CDAF809154@getmailspring.com?me=ca1a67e4&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">