<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Dear Alex,
<div class=""><br class="">
</div>
<div class="">As Paolo replied, jellium is unavoidable in periodic simulations of charged systems. If you are using periodic conditions to study an isolated system or other partially periodic systems, then there are ways to avoid the jellium and correct for
the pbc artifacts. In your molecular cation optimization, you can have results equivalent to free boundary conditions if you specify assume_isolated = ‘martyna-tuckermann’ and choose a cell size twice as large as your cation. If you plan to study charged defects
in periodic structures (crystals, liquids, amorphous), then there are very few things you can do apart from increasing the cell size. </div>
<div class=""><br class="">
</div>
<div class="">Calculations on cations in QE do work fine, forces are usually ok and geometry optimizations converge, unless there are other physical/numerical problems in your input. Mind that anions, instead, may pose convergence problems in QE due to the
well-known problems of DFT to localize the extra electron. This issue is not seen in localized basis set simulations, as the electron has no way to delocalize, but in plane-wave setups calculations on anions may not converge, unless one uses some tricks. </div>
<div class=""><br class="">
</div>
<div class="">What do you mean with the sentence “without jellium the optimized structure is very close to the real one”? Without jellium means for a charge neutral molecule? The real one means the one of Gaussian?</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; text-decoration: none;">
Oliviero Andreussi<br class="">
--<br class="">
Assistant Professor<br class="">
Department of Physics<br class="">
University of North Texas<br class="">
<a href="mailto:oliviero.andreussi@unt.edu" class="">Email: oliviero.andreussi@unt.edu</a> <br class="">
Skype: olivieroandreussi<br class="">
Web: <a href="https://sites.google.com/site/olivieroandreussi" class="">https://sites.google.com/site/olivieroandreussi</a></div>
</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Apr 25, 2018, at 6:03 AM, Aleksandra Oranskaia <<a href="mailto:aleksandra.oranskaia@kaust.edu.sa" class="">aleksandra.oranskaia@kaust.edu.sa</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Hi dear users and developers of QE,</div>
<div class=""><br class="">
</div>
<div class="">If one is interested in studying charged state defects -- is it correct to run supercell optimizations with a compensating jellium background ? </div>
<div class="">Let’s say, if one is interested in +1 point defect state -- is it correct to add tot_charge = +1 ?</div>
<div class=""><br class="">
</div>
<div class="">I quickly checked if the gradients for molecular cation optimization are reasonable with such an approach. For this I took small cation with known structure and:</div>
<div class="">- gaussian09-d.01 optimization with 6-11G**/B3LYP, +1 charge, 1 spin multiplicity gave perfectly matching structure with the real one </div>
<div class="">- QE-6.2.1 optimization with vdw-df2/USPP and tot_charge = +1 destroyed the cation, it dissociated; without “jellium” the optimized structure is very close to the real one.</div>
<div class=""><br class="">
</div>
<div class="">I would greatly appreciate your advice on how to obtain realistic structures of the supercells with charged defects,</div>
<div class=""><br class="">
</div>
<br class="">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">
Thanks in advance,<br class="">
Alex.<br class="">
___<br class="">
Aleksandra Oranskaia (M.Sc.)<br class="">
ChemS PhD student, KAUST<br class="">
Phone: +966 50 1335254<br class="">
<br class="">
<br class="">
</div>
<div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">
<br class="">
</div>
<br class="Apple-interchange-newline">
</div>
<br class="Apple-interchange-newline">
</div>
<br class="">
</div>
<br class="">
<div class="">
<hr class="">
</div>
<font face="Arial" size="1" class="">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this
message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font>_______________________________________________<br class="">
users mailing list<br class="">
<a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class="">
https://lists.quantum-espresso.org/mailman/listinfo/users</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>