<div dir="ltr">Hi Conrad,<div><br></div><div>Thanks so much for the reply. By retaining the configuration, I do mean keeping the symmetry. I currently have a submitted job with the outer atoms of the configuration fixed, and I will see the result.</div><div><br></div><div>This issue with the field is not only found in the charge -1 case, but also in the charge 0 case. In that case, there is not an unpaired charge, but the field is once again concentrated on one of the molecules, but this time on the upper side molecule. </div><div><br></div><div>And unfortunately, the system I am using only has CP2K 2.2 working. I will try and find a way to use more recent versions.</div><div><br></div><div>I also tried submitting another job, but with a differet basis set. I have seen in the forum that the basis used can have huge effects on the final converged result. I am going to try with DVZP-MOLOPT-SR-GTH. Would this make sense for a system like this? </div><div><br></div><div>Thanks,</div><div><br></div><div>Gary</div><div><br><br>On Sunday, August 14, 2016 at 4:36:27 AM UTC-7, Conrad wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Gary,<div><br><div><font color="#0000ff">On Thursday, August 11, 2016 at 7:58:08 PM UTC+1, Gary wrote:</font><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><font color="#0000ff">Hi,</font><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">I am trying to calculate the electric field above a system of molecules. The molecules are arranged on the xy plane. I have done the same for single molecule case, and that is working fine.</font></div><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">My question is, how should I optimise the geometry? If it's a system of molecules, I want to retain the configuration, but relax the system? There are 4 molecules, and each molecule has 50 atoms, so there are a total of 200 atoms. Is this too many for geometry optimisation? </font></div></div></blockquote><div><font color="#0000ff"><br></font></div><div><font color="#000000">By retaining the configuration, do you mean the symmetry? </font></div><div><font color="#000000">Is it too many? No. CP2K scales to thousands of atoms. As a first approximation of computer time, budget for a few hundred self-consistent energy calculations' worth. </font></div><div><font color="#0000ff"> </font></div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font color="#0000ff"> </font></div></div></blockquote><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font color="#0000ff">I did a calculation of the EFIELD_CUBE without optimisation, and the results are a little strange. The field is asymmetrical and seems to completely concentrate around one of the 4 molecules. The magnitudei of the field is also very high. I have attached the input and output files. There is a total charge of -1. </font></div></div></blockquote><div><br></div><div>Have a look at the SPIN DENSITY cube. What you're doing in this input is effectively introducing an excess electron by having a non-zero overall charge. The SPIN DENSITY cube will show you where the unpaired charge (de)localises and may give you some clues. I've never spent any time looking at EFIELD data so I can't comment on this directly. </div><div> <br></div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font color="#0000ff">I also want to point out that the SCF run does not converge in the inner loop (MAX_SCF 800); I am not sure if this is causing non-sensical data.</font></div></div></blockquote><div>This is okay. You might even want to reduce MAX_SCF for the inner loop. When using OT, preconditioning is repeated before running over the inner loop. In principle, shorter inner loops can improve convergence. </div><div><br></div><div>Anyway, this calculation is converged, so no issue there.</div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">Thanks,</font></div><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">Gary</font></div></div></blockquote><div><br></div><div>Finally, I notice that you're using version 2.2 of the code. This is really old, and may have some legacy bugs that haven't been patched.<br><br>Best wishes,<br>Conrad </div></div></div></div></blockquote></div></div>