<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hello,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I have a molecular crystal (1 unit cell = 56 atoms). I use SG15 ONCV 1.2 pseudopotential and Grimme D2 correction.<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Using vc-relax with PBE, the structure converges, no problem.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
But with PBE0, density and pressure oscillate without dampening (see figure below - Pressure remains negative, and not close to zero). Is there any particular strategy with PBE0?</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<ul>
<li><span>vc-relax with PBE and then starting from the result and use PBE0?</span></li><li><span>relax with PBE0 and then vc-relax?</span></li><li><span>any simulation parameter to change (ion dynamics and/or cell dynamics - currently BFGS)?</span></li><li><span>any strategy with PBE0?</span></li></ul>
<div></div>
<div><img size="28643" contenttype="image/png" style="max-width: 100%; user-select: none;" unselectable="on" tabindex="-1" data-outlook-trace="F:1|T:1" src="cid:ff10511e-c1ff-4370-a85d-f491329e2243"><br>
</div>
<div>Thank you,</div>
<div>Xavier<br>
</div>
<div><br>
</div>
</div>
</body>
</html>