<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear PWscf users and developers,<br>
<br>
I have 2 questions and a comment on finite electric field calculations in PWscf and your help and thoughts are very appreciated.<br>
<br>
<br>
Questions:<br>
------------<br>
(1) PWscf refuses smearing in Berry phase calculations (lelfield=.true.). However, smearing could help in accelerating the convergence of some in-gap states for semiconductors that arise because of a defect or impurity. Is there a fundamental reason for why
smearing should not be used with Berry phase calculations, or is it possible to modify PWscf to accept smearing?<br>
<br>
<br>
(2) I recall Prof. Giannozzi in this mailing list mentioned that DFT+U is not tested with Berry phase calculations. The code does not complain about DFT+U with Berry phase and I have performed preliminary tests on metal oxides and found out that dielectric
constant predicted using DFT+U and (lelfield=.true.) are reasonable compared to results obtained with other codes using DFTP(+U). I wonder if somebody did more tests on DFT+U with Berry phase and would like to share their conclusion.<br>
<br>
<br>
Comment:<br>
------------<br>
(3) The total energies that PWscf outputs in finite electric field calculations are in fact electric enthalpies. One could check that the ( Harris-Foulkes estimate + (total dipole x filed)) is the what PWscf prints as the total energy. I think it would
be nice to designate this output as the electric enthalpy espeically that this will be consistent with what PWscf prints out for calculations under pressure. (Of course with the implicit warning that the enthalpy inherits the uncertainty in the polarization
quantum).<br>
<br>
<br>
Thank you,<br>
Mostafa Youssef <br>
MIT<br>
<br>
<br>
</div>
</body>
</html>