[QE-users] using ESM with boundary condition "bc1" for charged slab

[Ghosh, Prasenjit]

Thanks a lot Paolo for the clarification.

With regards,
Prasenjit

----- Original Message -----
From: "Paolo Giannozzi" <paolo.giannozzi at uniud.it>
To: "users" <users at lists.quantum-espresso.org>, "pghosh" <pghosh at iiserpune.ac.in>
Sent: Wednesday, February 14, 2024 9:46:36 PM
Subject: Re: [QE-users] using ESM with boundary condition "bc1" for charged slab

I inquired with the author of ESM and this is what he told me:
"Since BC1 has no screening medium in the system, we need to keep charge 
neutral in the DFT cell. Otherwise, the total energy is diverging. In 
the earlier versions, the manual included a note for such charge 
neutrality condition, but for added safety, a routine to check this has 
been incorporated into the program.

In an ESM-RISM calculation, even when using BC1, there is no need to 
impose the charge neutrality condition in the DFT cell, as the RISM part 
automatically generates a compensating charge to maintain the charge 
neutrality condition throughout the whole system."

Paolo

On 2/9/24 06:07, Ghosh, Prasenjit wrote:
> 	
> Non si ricevono spesso messaggi di posta elettronica da 
> pghosh at iiserpune.ac.in. Informazioni sul perché è importante 
> <https://secure-web.cisco.com/12AJuZ9ZsbH7OSSDaFdwYXUBzeEOCvjPihIc7H0WS6WY6Hx9lMBQkozsbXYvV0X-7Cjhn8W42VD0os00rHvSbmmrEi9OT0_ofj4puyNmirk20xRGg1Qm1khWm51fv0ofpaqoZddH-87afA2Pgm1MrQKy1gNaLOuKYcIYZntQvwvU0hMs5hfaS4PrRCLsZDsp9OCaIxO4vaVNWDkKq1Dp_SrxG2-wQtaH0xwsvLXOMFszwHv0gXQVQf0dCWZ2em23xQ-oZfJ4PKQBYOVwX1oRrNWqqqRlrpNBUua6pZDVF-W7DhIWS7oogHz5S4tLXZFMy/https%3A%2F%2Faka.ms%2FLearnAboutSenderIdentification>
> 	
> 
> Dear all,
> 
> I am performing a calculation for a negatively charged slab using
> assume_isolated = 'esm'
> esm_bc = 'bc1'
> 
> Till qe-7.0 the scf calculation proceeds. However, from version 7.1 the 
> calculation stops with the following error:
> "cannot have charge, when ESM-BC1".
> 
> This is the consequence of a check by the "esm_check" subroutine. What I 
> fail to understand is that this check is also present in qe-6.8, 
> however, there the calculation doesn't stop.
> Am I missing something?
> Also why is this stopped from version 7.1. Are there any issues in the 
> previous versions regarding the implementations of the same? I thought 
> that this is an effective way to perform charged slab calculations.
> 
> Any help in clarifying this will be highly appreciated.
> 
> With regards,
> Prasenjit
> 
> 
> _______________________________________________
> The Quantum ESPRESSO community stands by the Ukrainian
> people and expresses its concerns about the devastating
> effects that the Russian military offensive has on their
> country and on the free and peaceful scientific, cultural,
> and economic cooperation amongst peoples
> _______________________________________________
> Quantum ESPRESSO is supported by MaX (http://secure-web.cisco.com/1wg-Sn_ZUMtTyRMZsJKyKizaSqviZGjuik7-lGVJ130wHC5m66P1pvQ-ZUMUrYeXMA3qovdkUqGYUVjeTA3Mtg8obh8En3oN8CaZYpId-5ugTnWdq_I3lPum3Z5dVkPCsx--ir1MrwJldToxCEhBi8JOASttShakfuwRQsrHgG1vXd8luY3SSSwp3WjTDfNaWMAcQmKCxEWh5k67PcPvMMSqvvWl5Ky4RllnaHPcD1sveD9rSkht8dgXkN0Pths4Vvy0lzUAxnVmnPiJUJrWnnwkgrRp5-bL5nzsis57bBSuGRHvuKSB2q76txY6-ehxL/http%3A%2F%2Fwww.max-centre.eu)
> users mailing list users at lists.quantum-espresso.org
> https://lists.quantum-espresso.org/mailman/listinfo/users

-- 
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
Univ. Udine, via delle Scienze 208, 33100 Udine Italy, +39-0432-558216