<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="">
I am sorry, I may not be familiar with the terminology. You want to run a simulation with constant applied potential? Do you plan to do a molecular dynamics simulation or just total energy/geometry relaxations? With Environ you can perform simulations with
constant electrode charge, which can be converted into a constant potential simulations. These simulations could include a reasonable description of the solvent interface, as well as a reasonable description of the electrochemical diffuse layer. You could
vary the charge in the substrate (i.e. tot_charge in the pw input file) so as to obtain the potential drop that you want, or if you are interested in computing interfacial energies as a function of different potentials, you could use an approach based on Legendre
transforms such as the one described in N. Hörmann et al, J. Chem. Phys. 150, 041730 (2019).
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Oliviero<br class="">
<div class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Jul 22, 2020, at 9:53 AM, Antonio Crepaldi <<a href="mailto:antonio.crepaldi.lanza@gmail.com" class="">antonio.crepaldi.lanza@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">Dear Oliviero,
<div class=""><br class="">
</div>
<div class="">Thank you very much for your detailed explanation!</div>
<div class="">I would just like to ask you if it is possible to make a constant mu calculation with environ. I've read in the manual that this requires ESM, but I don't know if this is compatible with the environ setup.</div>
<div class="">Thank you again.</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Antonio</div>
</div>
<br class="">
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">Il giorno mar 21 lug 2020 alle ore 17:54 Andreussi, Oliviero <<a href="mailto:Oliviero.Andreussi@unt.edu" class="">Oliviero.Andreussi@unt.edu</a>> ha scritto:<br class="">
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div style="word-wrap:break-word;line-break:after-white-space" class="">I am sorry my answer can be misleading. The wave function and the density also depend on the environment, but the point is that PP does not recompute them, it just reads them from the saved
wfc files and plot them. PP instead needs to recompute the potentials to plot them and this is where the missing environ makes the results meaningless.
<div class=""><br class="">
</div>
<div class="">I hope this helps,</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Oliviero<br class="">
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Jul 21, 2020, at 10:47 AM, Andreussi, Oliviero <<a href="mailto:Oliviero.Andreussi@unt.edu" target="_blank" class="">Oliviero.Andreussi@unt.edu</a>> wrote:</div>
<br class="">
<div class="">
<div style="word-wrap:break-word;line-break:after-white-space" class="">Dear Antonio,
<div class=""><br class="">
</div>
<div class="">Thanks for reporting this issue and for trying Environ for your simulations. I am afraid the documentation of Environ needs to be improved on the point that you report, as this is not the first time I see this question. </div>
<div class=""><br class="">
</div>
<div class="">Environ is not coupled with the PP code. If PP is used to plot quantities that do not depend on the environment, such as the wavefunction, the density, etc., it should work fine. However, you cannot use PP to plot quantities that depend on the
environment, namely the electrostatic potential. For this purpose, you can turn on the verbosity of Environ using the verbose keyword. If you set verbose = 2 (or higher) you should get some .cube files that you can use to visualize properties of the simulation
cell, in particular the electrostatic potential. Cubefiles are read by many programs and can be processed relatively easily with scripts.</div>
<div class=""><br class="">
</div>
<div class="">Don’t hesitate to let me know if you have any problem with running Environ. You may also consider joining the google group created to support Environ users, there are not too many discussions right now, but it would be nice to build a collection
of known problems, solutions and suggestions. <a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fforum%2F%3Fhl%3Den%23!forum%2Fquantum-environ-users&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966371043&sdata=lPqnsBoTuLd%2FyGG3dFHx2gIRWosfJ9AsDh1VoCdGgxc%3D&reserved=0" originalsrc="https://groups.google.com/forum/?hl=en#!forum/quantum-environ-users" shash="N/VY5rJn0Q5NLpC21NhNwQkzLLqlY8X1ZM+4TDNzDrH6pwxlHWHKGF/CilZsH9A3LQ98kbGm98y4j4hXIDl7n4Hq5LC+RvmuhRVxu+TXEERGTDBTOlHUhC+ZzA+cmQonZ+l0QYlkxoQqJHQlWvL3U0UguoV/O9JZ9nxwq85nLiI=" target="_blank" class="">quantum-environ-users</a> (<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fforum%2F%3Fhl%3Den%23!forum%2Fquantum-environ-users&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966381038&sdata=2D%2Bjkelhj5dpTHSEtSr6PDkrKTPIInPQMvsJG2EznKw%3D&reserved=0" originalsrc="https://groups.google.com/forum/?hl=en#!forum/quantum-environ-users" shash="orQRXP7gqqF/3mqFcbIMx5EXLfocEeBqgqKngf6FwGcjqZLeTpsSwgMItobJvXT4OtBuBZWqs8EpG1NAhRuXe3gqbsrw5caD9DrbuWYkaHizZAp/9yXNZBpk9nmVaBgpHFWbsfjZWL+XOOa32jnU534p+Q9mzYntqBigv6qlqWY=" target="_blank" class="">https://groups.google.com/forum/?hl=en#!forum/quantum-environ-users</a>)</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Oliviero</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div dir="auto" style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space" class="">
<div dir="auto" style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space" class="">
<div dir="auto" style="word-wrap:break-word;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;text-decoration:none" class="">
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" target="_blank" class="">Email: oliviero.andreussi@unt.edu</a> </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;text-decoration:none" class="">
Phone: +1-(940)-369-5316<br class="">
Skype: olivieroandreussi<br class="">
Web: <a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.materialab.org%2F&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966381038&sdata=tK05iKxaG0dKOHXN4IUxF8twFc0f1Bjy4hriZTv9Bkk%3D&reserved=0" originalsrc="https://www.materialab.org/" shash="Wpvc9jtJnkWPqIhhsYQQ4Y051DnDx7gPLZftvuyDVxqfhsg/LeRqkrohZWoWAJfI0K1vhN56V83MENiqb7IhZgiS+rUhiQO+VlXnbPZa41DKeEQC8JT0k8q39y6jxVdv4aA1RVAL24T2kuUmk6A1l66Fqu1pEUBoBXfyZjpMUDU=" target="_blank" class="">https://www.materialab.org</a></div>
</div>
</div>
</div>
</div>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Jul 21, 2020, at 8:42 AM, Antonio Crepaldi <<a href="mailto:antonio.crepaldi.lanza@gmail.com" target="_blank" class="">antonio.crepaldi.lanza@gmail.com</a>> wrote:</div>
<br class="">
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">Dear QE users and developers,
<div class=""><br class="">
</div>
<div class="">I'm currently trying to simulate a charged 2D slab embedded in water. I'm using two planar countercharges as explained in Example 5 of Environ. The computation does not have any convergence problem, however, when I plot the total electrostatic
potential (plot_num=11 in pp.x) I do not get a flat potential in the vacuum region screened by the countercharges. I would expect to find a flat potential in that region and be able to use it as a reference for further calculations.</div>
<div class="">Here is my environ input</div>
<div class=""><br class="">
</div>
<div class="">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
&ENVIRON<br class="">
verbose = 0<br class="">
environ_thr = 1.d0<br class="">
environ_type = 'input'<br class="">
env_electrostatic = .true.<br class="">
env_static_permittivity = 80<br class="">
env_surface_tension = 0.D0<br class="">
env_pressure = 0.D0<br class="">
env_external_charges = 2<br class="">
/<br class="">
&BOUNDARY<br class="">
solvent_mode = 'full'<br class="">
/<br class="">
&ELECTROSTATIC<br class="">
pbc_correction = 'parabolic'<br class="">
pbc_dim = 2<br class="">
pbc_axis = 3<br class="">
tol = 5.D-13<br class="">
!inner_tol = 5.D-18<br class="">
/<br class="">
EXTERNAL_CHARGES (bohr)<br class="">
2. 0.0 0.0 -14.29 1.0 2 3<br class="">
2. 0.0 0.0 14.29 1.0 2 3</blockquote>
<div class=""><br class="">
</div>
<div class=""> Thank you very much for your help.</div>
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Antonio Crepaldi-Lanza</div>
<div class=""><br class="">
</div>
<div class="">--</div>
<div class="">Department of Physics "Aldo Moro"</div>
<div class="">University of Bari</div>
</div>
</div>
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966391035&sdata=38TOc%2F%2B4xmKTHySJFL6wZY0Wx8ABT7KimJzm2WrR0uk%3D&reserved=0" originalsrc="http://www.max-centre.eu/quantum-espresso" shash="Ct+A9sQ3BgA6FM/IGYjOINqfhfjWXwgr641tOE6Cqa+6UPuhz2Yod+3vMq+bYW6pXy2IHWGaRW7EqsGnY32XCxWdYpiBqNjdfIAxsbS8G0+dTJU20Sd7Hc2yRZUdojxyYFGjSK2lisIaQazFZs4XkKvtOyWApar6RSYrWbS5Pfo=" target="_blank" class="">https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7Cbf10c4b1071745838eb008d82d7c03d9%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637309357952927330&sdata=pvaXJ16pf8n%2Bmah6qe4wX7cgT9o7bvGnO4f6FrrYukM%3D&reserved=0</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">
users@lists.quantum-espresso.org</a><br class="">
<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966401027&sdata=%2BjuDjqfCUjrWeNawFWCl0wBj%2FfkH80kHuoiwVc%2BiDCI%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/mailman/listinfo/users" shash="GqzLWYAqwfDkZ9L8ibZ3WdMtC1Wf1QIu9AWvmjFgL5AlxLF9JxFpcsptFcRADus7mTyK1VF8Xlu5IOf7h1uxqAvgkPKfFsdUkxvA73UChdn2ca/SIDMEsl2HCIqdOiRx839svy8lTQ2joQ1uDb+I8jk0lZx7T2hj4m8eDbj5sbA=" target="_blank" class="">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7Cbf10c4b1071745838eb008d82d7c03d9%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637309357952937325&sdata=j2PM%2FJTHGhBVzsNqYw5ZmvWAj3XuuafBVypy2DG5am4%3D&reserved=0</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966401027&sdata=Rhy3j%2FoToWrjRIt7xt0Ja5f5%2FH9fUicdVqM404nd5kw%3D&reserved=0" originalsrc="http://www.max-centre.eu/quantum-espresso" shash="Rle8G1IqbGEJSEChrRK6tTWdEyHI+N+LyzQJk4AYzKRNYJx3YBRe0f2Lhvj0uMcSleNhx1cTUGlDTJKc2TshvfeNIQG03/MBW+pKlmgnMJi8FSwjQ9qxeYesAVt6F7wydWeGInXmCl8i9aI3G3hPOHGEPvmWKmQSmLSc9smiRjY=" target="_blank" class="">https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C3c06ec5dcb694d44a47808d82d8d85ee%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637309433158048698&sdata=oxwo%2B%2F%2BVi7N5kgQ7evuEJYeZMkZwz8s8wZIrm4Q%2Fcg0%3D&reserved=0</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">
users@lists.quantum-espresso.org</a><br class="">
<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966411028&sdata=BzQcluqBdgXXs2eUv0wYecEZYvalAcpZuehXN35Iff0%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/mailman/listinfo/users" shash="qYiDdxNZ0p+gXPT74wjMg56I/d8bUISzW6DQRtjbASEyA9tZIqJWFOjK7JNcQqMpu3ilUtK6OpkP7UPSMNTlmINbqaf/DxvhDtMioMNLPPjE2FzLtj8L4aboitSk7/FYPx+aYeBrwg2E4wyqSrLI97cv62OjNKlXBECx2hUXQiQ=" target="_blank" class="">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C3c06ec5dcb694d44a47808d82d8d85ee%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637309433158048698&sdata=X30Qq0KT%2BDtSZHCZrVo5%2BhCFPqGdBnGJ%2FFfE7%2FFQBDI%3D&reserved=0</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966421023&sdata=%2Fz6Ixmb7IoAl3UOx4LJUDCo62SJ6UgvOLMV2Y%2Bg15vs%3D&reserved=0" originalsrc="http://www.max-centre.eu/quantum-espresso" shash="w3Mfx8fK7P+SALlAjrm6sxBuxap27NC/3S6vcj0pjzaXCgAYRSunpUWoZiPk/nixs5cxdcUIGNdbSz7GFG5nL0TYLaklOydZrN5OJz3s7KFt2NH2j1Z2u+A1zcuTxBAiK9iwR3rQ1L+oOn9DPH1cY08e+djNX1CodZqu9aVXvJo=" rel="noreferrer" target="_blank" class="">www.max-centre.eu/quantum-espresso</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">
users@lists.quantum-espresso.org</a><br class="">
<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966421023&sdata=TtqbnKvUN6OnmJiTbZDscdstkhYEjHW%2BL%2B1iNEEUY5M%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/mailman/listinfo/users" shash="kNFep3iHpQMHTUAta+HjFx+9nBEZjJB6hkJzr/FQqCxPs4RmyUEYIQENf2UyrCVuhYrtue4bQK6TWplgkxyBjWGIKRxyj0ZiE8IEDUC519N5PeprXWQXkdbVUh1wssEo9p8BrVq+cTuOJ4IFDsI6O/r33wiZGLmjGNcXhhZ9Xxw=" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote>
</div>
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966461010&sdata=Cb6uC1pSzkL4ZO1bWOXWZd5A%2FjSjvD7mpoKCI53XecI%3D&reserved=0" class="">https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2Fquantum-espresso&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966461010&sdata=Cb6uC1pSzkL4ZO1bWOXWZd5A%2FjSjvD7mpoKCI53XecI%3D&reserved=0</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966461010&sdata=cilJg%2FMI8GIFB50fz%2BoUfKc0TIhYvtHrel%2BMYIO%2BlIs%3D&reserved=0" class="">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C883c9e7711f7412393b708d82e4f31a4%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637310264966461010&sdata=cilJg%2FMI8GIFB50fz%2BoUfKc0TIhYvtHrel%2BMYIO%2BlIs%3D&reserved=0</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</body>
</html>