<div dir="ltr"><div dir="ltr">Dear QE users and developers,<div><br></div><div>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>Here is my environ input</div><div><br></div><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">&ENVIRON<br>   verbose = 0<br>   environ_thr = 1.d0<br>   environ_type = 'input'<br>   env_electrostatic = .true.<br>   env_static_permittivity = 80<br>   env_surface_tension = 0.D0<br>   env_pressure = 0.D0<br>   env_external_charges = 2<br>/<br>&BOUNDARY<br>   solvent_mode = 'full'<br>/<br>&ELECTROSTATIC<br>   pbc_correction = 'parabolic'<br>   pbc_dim = 2<br>   pbc_axis = 3<br>   tol = 5.D-13<br>   !inner_tol = 5.D-18<br>/<br>EXTERNAL_CHARGES (bohr)<br>2. 0.0 0.0 -14.29 1.0 2 3<br>2. 0.0 0.0  14.29 1.0 2 3</blockquote><div><br></div><div> Thank you very much for your help.</div></div><div>Best,</div><div><br></div><div>Antonio Crepaldi-Lanza</div><div><br></div><div>--</div><div>Department of Physics "Aldo Moro"</div><div>University of Bari</div></div></div>