Hi Mauro, To add to what I replied before, I realized that if the implementation of ESM in QE does not include confinement, you could you the Environ one. Thinking at this conversation I realized that, while ESM and Environ are not compatible for what concerns electrostatics, it is probably possible to use the confining potential of Environ in a normal QE calculation with the ESM turned on. I don’t think anyone has tried it, but I don’t see any reason why it should not work. Just to give you some more options.
Best, Oliviero Oliviero Andreussi -- Assistant Professor Department of Physics University of North Texas Email: [email protected]<mailto:[email protected]> Phone: +1-(940)-369-5316 Skype: olivieroandreussi Web: https://www.materialab.org On Jun 28, 2020, at 2:49 PM, Mauro Sgroi <[email protected]<mailto:[email protected]>> wrote: Dear Minoru, Thanks a lo for the clarification. Kind regards, Mauro. Il sab 27 giu 2020, 08:04 大谷実 <[email protected]<mailto:[email protected]>> ha scritto: Dear Mauro, There are possibly two barriers to an MD calculation with ESM. One is the barrier for electrons, and the other is atoms. The former prevents electrons from leaking into the vacuum region when the surface is negatively charged. The latter prevents atoms from going (evaporating) into the vacuum region. We usually use a trigonometric function for the barrier of electrons and a parabolic function for atoms. You can find the original function forms in the following reference: Surf. Sci. 601 5237 (2007). I am not sure such barriers are introduced to QE. At least I haven't implemented it. Best regards, Minoru -------------------------------------------------------------------------------------------------------- National Institute of Advanced Industrial Science and Technology, Research Centre for Computational Design of Advanced Functional Materials E-mail: [email protected]<mailto:[email protected]> tel : +81-29-861-5202 -------------------------------------------------------------------------------------------------------- ________________________________________ From: users <[email protected]<mailto:[email protected]>> on behalf of Mauro Sgroi <[email protected]<mailto:[email protected]>> Sent: Saturday, June 27, 2020 7:57 AM To: Quantum ESPRESSO users Forum Subject: [QE-users] Surface-model slab: ESM and definition of a liquid region inside the vacuum Dear all, I'm trying to simulate the interface between a charged doped graphene layer and a liquid electrolyte. I would like to treat the solvent explicitly. So I should place some molecules of the solvent in the vacuum region of my slab model. I've seen that with the OpenMX code it is possible to use the ESM method and define a potential to confine the liquid in a certain portion of the space between periodic images of the substrate. This prevents the liquid filling all the vacuum space of the supercell during MD. Here is an explicative figure: http://www.openmx-square.org/openmx_man3.7/node118.html<https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.openmx-square.org%2Fopenmx_man3.7%2Fnode118.html&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003170363&sdata=ynvYs1MyEcIn8aXpMQCCQY0t2vVOo3WNosDaVnY1k48%3D&reserved=0><https://jpn01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.openmx-square.org%2Fopenmx_man3.7%2Fnode118.html&data=02%7C01%7Cminoru.otani%40aist.go.jp%7C9d3031d5a9ae40fcea8a08d81a2467b2%7C18a7fec8652f409b8369272d9ce80620%7C0%7C0%7C637288090956080025&sdata=ngwDC21Ut6Pdxjmy6lWI2DMkeXHeJRBHk1EtqTtLAjw%3D&reserved=0<https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.openmx-square.org%2Fopenmx_man3.7%2Fnode118.html&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003180354&sdata=4ryVUpUdmwE5R6M3zQaXIWCrcbdsR0J3xvpw1bf0GyY%3D&reserved=0>> Is something similar implemented in Quantum Espresso? Or is it possible to obtain it using constraints on the coordinates? Thanks a lot in advance and best regards, Mauro Sgroi. Centro Ricerche FIAT. _______________________________________________ Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso<https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003180354&sdata=4C0hanUtewLe60vWGbAuvEPlkLZkoWe3wjtsBuKQHsg%3D&reserved=0>) users mailing list [email protected]<mailto:[email protected]> https://lists.quantum-espresso.org/mailman/listinfo/users<https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003190349&sdata=O4IMkcCOwrxCxqtLe1GUNszyF1weqtIVIlIaU%2ByxUnQ%3D&reserved=0> _______________________________________________ Quantum ESPRESSO is supported by MaX (https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003220334&sdata=dbdCUDtHSq%2FbxqSzSCnaEKU4KQhQV9ldplsOTopJBJ0%3D&reserved=0) users mailing list [email protected]<mailto:[email protected]> https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C90aab5056260404119ea08d81b9c6f82%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637289706003220334&sdata=3eQ752zSjynnRO4zvKFc65ZtauirtxDoukKAmLSgBAw%3D&reserved=0
_______________________________________________ Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso) users mailing list [email protected] https://lists.quantum-espresso.org/mailman/listinfo/users
