<div dir="ltr">Dear Fedor,<br><br>thank you for your message. I have found a (fairly OK) solution to applying higher fields that seems to work OK with graphene and is actually described in example 10 of PW but I think the sample does not apply the technique. I quote:<br><br>"To perform a calculations with an electric field, an estimate of<br>the optimized wavefunctions is needed to build the electric field<br>operator (See: I. Souza, J.Iniguez and D. Vanderbilt, PRB 69, 085106,<br>2004). Therefore when lelfield ==.true. a copy of the wavefunctions<br>is read from disk (i.e. startingwfc should be 'file')."<br><br><br>therefore one first runs the scf input file with everything present but the efield turned off (efield_cart(1,2,3)=0.0). After converging this one copies edits the efield (efield_cart(3)=0.02) and adds "startingwfc='file'); this has allowed me to use much higher fields than a "cold start" (startingwfc='random') and the behavior seems consistent with literature! It seems that in the zero-field case all efield related inputs have to be present, otherwise I ran into some I/O error ( Error in routine gk_sort (1): array gk out-of-bounds)<br><br>Hope this helps people encountering convergence problems in the future!<br><br>With best regards,<br>Chris <br><br>On Fri, Jan 5, 2018 at 12:40 AM, Fedor Goumans <<a href="mailto:goumans@scm.com">goumans@scm.com</a>> wrote:<br>><br>> Dear Christoph,<br>><br>> While I don’t have an answer to your specific QE question, the application of an electric field to a 2D material in a plane wave code may be problematic to begin with.<br>> We had someone try out band gap closure in the ‘MoWSeS’ materials under electric fields (which would amount to sawtooth potentials in PW codes if I understood correctly).<br>> In that case PW codes would just give wrong behavior, while localized AO codes with prober 2D pbc (and a homogenous electric field) would allow for polarization across the surface.<br>> See e.g.: <a href="https://www.scm.com/highlights/closing-band-gap-2d-semiconductors/">https://www.scm.com/highlights/closing-band-gap-2d-semiconductors/</a><br>><br>> Hope this helps,<br>> Best wishes,<br>> Fedor<br>><br>> On Jan 1, 2018, at 13:13, Christoph Wolf <wolf.christoph@qns.science> wrote:<br>><br>> Dear all,<br>><br>> a happy new year!<br>><br>> I have recently played with 2-D materials (yes, graphene) and electric fields. I encountered several situations where convergence is really hard to achieve and I was wondering which of the following parameters is usually considered to be helping with convergence<br>><br>> &control<br>> ....<br>>     lelfield=.true.<br>>     gdir=3<br>>     nberrycyc=6<br>>    nppstr=1<br>> /<br>><br>> &ELECTRONS<br>>  ....<br>>  mixing_beta=0.4<br>>  mixing_mode='TF'<br>>  efield_cart(3)=0.00275<br>> /<br>><br>><br>> the system is separated in z-direction by ~20A of vacuum. the cutoffs for (NC) PPs were converged in the field-free case.<br>><br>> What is a suggested "systematic" approach to get convergence in fields of this magnitude?<br>><br>> Your help is greatly appreciated!<br>><br>> With best wishes for the new year,<br>><br>> Chris<br>> --<br>> Postdoctoral Researcher<br>> Center for Quantum Nanoscience, Institute for Basic Science<br>> Ewha Womans University, Seoul, South Korea<br>><br>> _______________________________________________<br>> Pw_forum mailing list<br>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>> <a href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a><br>><br>><br>> Dr. T. P. M. (Fedor) Goumans<br>> Business Developer<br>> Software for Chemistry & Materials<br>> Vrije Universiteit, FEW, Theoretical Chemistry<br>> De Boelelaan 1083<br>> 1081 HV Amsterdam, The Netherlands<br>> T +31 20 598 7625<br>> <a href="https://www.scm.com">https://www.scm.com</a><br>> <a href="https://twitter.com/SCM_Amsterdam">https://twitter.com/SCM_Amsterdam</a><br>> <a href="https://www.linkedin.com/company/software-for-chemistry-&-materials">https://www.linkedin.com/company/software-for-chemistry-&-materials</a><br>><br>><br>><br><br><br><br>--<br>Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea</div>