<div dir="ltr"><div>The ESM-RISM extension of QE, written by Satomichi Nisihara, is available here:<br>   <a href="https://github.com/nisihara1/q-e">https://github.com/nisihara1/q-e</a><br></div>It is currently aligned to QE v.6.1<br><br>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 15, 2018 at 12:05 PM, Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Minoru<br>
I've read the PRB paper you suggested and I would be very interested in testing such a new 3D-RISM extension of the ESM method of yours. The paper mention a QE implementation, but I've not seen any documentation in QE 6.2.1 yet. Shall I find it either in the unstable version or in a next stable release?<br>
Thank you in advance<br>
Giuseppe<div class="HOEnZb"><div class="h5"><br>
<br>
Quoting 大谷?g <<a href="mailto:minoru.otani@aist.go.jp" target="_blank">minoru.otani@aist.go.jp</a>>:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Hong,<br>
<br>
In the boundary condition bc1, i.e. the open boundary condition, we need to prepare a neutral system. Because there are divergent terms in the total energy and the electrostatic potential. If your system is neutral, i.e. tot_charge=0, you can get an appropriate result with ESM, but your system is charged, i.e. tot_charge!=0, you can not get a physically meaningful result.<br>
<br>
In other boundary conditions such as bc2 and bc3, the ESM automatically introduce compensation charge. Thus the whole system is neutral. In other words, there are no divergent terms in the total energy and the electrostatic potential for bc2 and bc3.<br>
<br>
As you can understand from above description, the problem is the divergent term. If you can ensure that the additional charge induced in the system with bc1 will completely be screened by another medium, you can use the bc1 boundary condition. In this case, the total energy and the electrostatic potential are physically meaningful values. In other words, the divergent term from another medium eliminates that from the system.<br>
<br>
If you want to dig into the details of this, you can read the following paper: Phys Rev. B 96, 115429 (2017).<br>
<br>
In short, you can not use bc1 with a charged system.<br>
<br>
Best regard,<br>
Minoru<br>
<br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<br>
National Institute of Advanced Industrial Science and Technology,<br>
Research Centre for Computational Design of Advanced Functional Materials<br>
E-mail: <a href="mailto:minoru.otani@aist.go.jp" target="_blank">minoru.otani@aist.go.jp</a><br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<br>
<br>
<br>
______________________________<wbr>__________<br>
From: users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-e<wbr>spresso.org</a>> on behalf of Honghui Wu <<a href="mailto:hhwuaa@connect.ust.hk" target="_blank">hhwuaa@connect.ust.hk</a>><br>
Sent: Saturday, March 10, 2018 2:38:38 AM<br>
To: <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
Subject: [QE-users] ESM calculation with tot_charge in QE<br>
<br>
Dear all,<br>
<br>
Can we use the charged system (tot_charge= a finite value) with the<br>
vacuum-slab-vacuum boundary condition (esm_bc = 'bc1') in Quantum ESPRESSO?<br>
<br>
In the OpenMX website, I find the ESM developer writes 'ESM.switch =<br>
on1: Both ESM (I) and (II) are semi-infinite vacuum media. In this case,<br>
note that the total charge of a calculation system should be neutral.<br>
The keyword 'scf.system.charge' should be set to be zero.'.<br>
However, I did not find the similar description of the esm_bc= 'bc1' :<br>
Vacuum-slab-vacuum (open boundary conditions).  Does it mean we can use<br>
the charged system (tot_charge= a finite value) with the<br>
vacuum-slab-vacuum boundary condition in Quantum ESPRESSO?<br>
<br>
Thank you!<br>
<br>
Best regards,<br>
Hong<br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso<wbr>.org/mailman/listinfo/users</a><br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso<wbr>.org/mailman/listinfo/users</a><br>
</blockquote>
<br>
<br>
<br></div></div><span class="HOEnZb"><font color="#888888">
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) <a href="tel:%2B39%20373%207305625" value="+393737305625" target="_blank">+39 373 7305625</a><br>
Tel <a href="tel:%2B%2039%2006%2090672342" value="+390690672342" target="_blank">+ 39 06 90672342</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316" target="_blank">+39 06 90672316</a><br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></font></span><div class="HOEnZb"><div class="h5"><br>
<br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso<wbr>.org/mailman/listinfo/users</a></div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>