<html><head></head><body><div class="ydpf2ed9143yahoo-style-wrap" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 13px;"><div></div>
<div><br></div><div dir="ltr" data-setdir="false"><div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">Dear Dr. Tamas </div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">your ideas are very helpful. Your are right.</div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">i have just noticed that it is a surface reaction + diffusion of dissociated products (SO & O) to most stable sites ( one reaction + 2 diffusion step). Dissociation occurs in image 2 immediately , image 3,4,5 shows diffusion of SO to most stable site (a bit far site , 3 Angstrom), (image 5 is rotation of SO), image 5 shows diffusion of "O" and again rotation of "SO" . </div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">So how to break down this steps into simpler steps than can be easily handle by NEB ? </div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">how to deal many rotations of <span style="color: rgb(0, 0, 0);">"SO" molecule</span> as it takes most of the images ? </div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;"><br clear="none"></div><div dir="ltr" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">attached is axsf output file for the neb, please view it with xcrysden </div></div><br></div>
</div><div id="yahoo_quoted_4350263946" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Monday, October 26, 2020, 07:20:48 PM GMT+4, Omer Mutasim <omermutasim@ymail.com> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv3904053764"><div><div class="yiv3904053764ydp40198310yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div></div>
<div><br clear="none"></div><div dir="ltr">Dear Dr. Tamas </div><div dir="ltr">your ideas are very helpful. Your are right.</div><div dir="ltr">i have just noticed that it is a surface reaction + diffusion of dissociated products (SO & O) to most stable sites ( one reaction + 2 diffusion step). Dissociation occurs in image 2 immediately , image 3,4,5 shows diffusion of SO to most stable site (a bit far site , 3 Angstrom), (image 5 is rotation of SO), image 5 shows diffusion of "O" and again rotation of "SO" . </div><div dir="ltr">So how to break down this steps into simpler steps than can be easily handle by NEB ? </div><div dir="ltr">how to deal many rotations of <span><span style="color:rgb(0, 0, 0);font-family:Helvetica Neue, Helvetica, Arial, sans-serif;">"SO" molecule</span></span> as it takes most of the images ? </div><div dir="ltr"><br clear="none"></div><div dir="ltr">attached is axsf output file for the neb, please view it with xcrysden </div>
</div><div class="yiv3904053764yqt9764352963" id="yiv3904053764yqt11135"><div class="yiv3904053764yahoo_quoted" id="yiv3904053764yahoo_quoted_4285867970">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Monday, October 26, 2020, 06:31:22 PM GMT+4, Tamas Karpati <tkarpati@gmail.com> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div dir="ltr">please note that in case your preoptimized first and/or last structures<br clear="none"></div><div dir="ltr">are not the direct reactant and product structures but a R+diffusion<br clear="none"></div><div dir="ltr">or P+diffusion step results, your MEP would (and it does) look<br clear="none"></div><div dir="ltr">as if you have modelled a two (or even more) steps "reaction" even if<br clear="none"></div><div dir="ltr">just one of them is actually "chemistry".<br clear="none"></div><div dir="ltr"><br clear="none"></div><div dir="ltr">On Mon, Oct 26, 2020 at 3:11 PM Omer Mutasim <<a rel="nofollow" shape="rect" ymailto="mailto:omermutasim@ymail.com" target="_blank" href="mailto:omermutasim@ymail.com">omermutasim@ymail.com</a>> wrote:<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> yes, there is a barrier for the reverse reaction.<br clear="none"></div><div dir="ltr">> i have check the initial & final structure again , it was relaxed until force is less than 0.003.<br clear="none"></div><div dir="ltr">> I do also agree with you that i should use 1 neb with barrier, and this is what i'm doing exactly.<br clear="none"></div><div dir="ltr">> but for this dissociation reaction step (SO2 = SO + O), i think this is the simplest elementary step i can get form SO2 molecule , it can't be broken down into a simpler elementary reaction steps.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> One this i should mention is that : for the initial structure (SO2*) , I didn't consider the most stable adsorption site for SO2* (E_ads=-0.3 eV), there is neighboring site that is a bit less stable (E_ads = -0.2 eV) which i've used for NEB. I didn't consider the former (most stable site) because the distance between molecule and surface is 3.5 A , , however for the less stable site , the distance is 1.5 A, so i thought it is not proper to consider this physisorbed state in NEB, Please correct me if i'm wrong.<br clear="none"></div><div dir="ltr">> thanks a lot for your help.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> Regards<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> On Monday, October 26, 2020, 05:00:26 PM GMT+4, Tamas Karpati <<a rel="nofollow" shape="rect" ymailto="mailto:tkarpati@gmail.com" target="_blank" href="mailto:tkarpati@gmail.com">tkarpati@gmail.com</a>> wrote:<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> Dear Omer,<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> I'd like to underline what Antoine has said and suggest that<br clear="none"></div><div dir="ltr">> both your first and last structures are saddle points (of order K and L,<br clear="none"></div><div dir="ltr">> respectively), rather than minima. Also I agree that your reaction<br clear="none"></div><div dir="ltr">> is not barrierless -in accordance with chemical intuition.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> To make it simple, I recommend<br clear="none"></div><div dir="ltr">> - make your 2nd image the first,<br clear="none"></div><div dir="ltr">> - make your 6th image the last,<br clear="none"></div><div dir="ltr">> - do use more images (even for a single step reaction,<br clear="none"></div><div dir="ltr">> but your MEP indicates two TS-es meaning a 2 step reaction<br clear="none"></div><div dir="ltr">> which -on more elaboration- may turn out to be an N step one;<br clear="none"></div><div dir="ltr">> only you need more points to see its real E-profile).<br clear="none"></div><div dir="ltr">> - rerun your NEB job.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> In fact you should check for imaginary second derivatives by the Phonon code<br clear="none"></div><div dir="ltr">> to ensure about each minima/maxima/TS being what they look in such a MEP.<br clear="none"></div><div dir="ltr">> Of course, it is painfully slow and many just skip this step.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> Other ideas:<br clear="none"></div><div dir="ltr">> - check the geometry of your 3 local minima and 2 TS-es<br clear="none"></div><div dir="ltr">> to see if they correspond to chemically rational structures.<br clear="none"></div><div dir="ltr">> - if they look so, decide whether the second minimum is really<br clear="none"></div><div dir="ltr">> your product, and if it is so then this should be your last image<br clear="none"></div><div dir="ltr">> (and not the 6th as I said above) for rerunning your NEB.<br clear="none"></div><div dir="ltr">> - of course, you can model a multistep reaction by a single NEB job,<br clear="none"></div><div dir="ltr">> choose the approach that best serves your postprocessing task.<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> Hope this helps,<br clear="none"></div><div dir="ltr">> t<br clear="none"></div><div dir="ltr">><br clear="none"></div><div dir="ltr">> On Mon, Oct 26, 2020 at 1:40 PM Antoine Jay <<a rel="nofollow" shape="rect" ymailto="mailto:ajay@laas.fr" target="_blank" href="mailto:ajay@laas.fr">ajay@laas.fr</a>> wrote:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > There is an energy barrier:<br clear="none"></div><div dir="ltr">> > the one between your intermediate minima and your final state.<br clear="none"></div><div dir="ltr">> > There is no barrier between initial and intermediate minima.<br clear="none"></div><div dir="ltr">> > You should wonder why you have an intermediate minima that is lower in energy (<0.4eV) than the final inserted molecule, this is why I was asking if it was enough relaxed.<br clear="none"></div><div dir="ltr">> > Maybe the first exothermic reaction gives enough energy for the second...<br clear="none"></div><div dir="ltr">> > But for sure, when you have such a multi-barriers reaction, a 7 images neb is not enough.<br clear="none"></div><div dir="ltr">> > If you need accurate results, it is better to have 1 neb per barrier, as you have 1 CI per path.<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Regards,<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Antoine Jay<br clear="none"></div><div dir="ltr">> > LAAS-CNRS<br clear="none"></div><div dir="ltr">> > Toulouse, France<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Le Lundi, Octobre 26, 2020 09:10 CET, Omer Mutasim <<a rel="nofollow" shape="rect" ymailto="mailto:omermutasim@ymail.com" target="_blank" href="mailto:omermutasim@ymail.com">omermutasim@ymail.com</a>> a écrit:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Dear Dr. Jay<br clear="none"></div><div dir="ltr">> > I have relaxed the initial and final structures before neb.<br clear="none"></div><div dir="ltr">> > Regarding simulation box, i'm using sqrt(3)*sqrt(3) supercell, the other five reaction steps converged well.<br clear="none"></div><div dir="ltr">> > However, i have seen in the literature that similar catalyst resulted in such barrier-less dissociation.<br clear="none"></div><div dir="ltr">> > So my question goes like : with this oscillated MEP , can i conclude it is barrier-less reaction ? or it is even necessary for barrier-less step to have no oscillation ?<br clear="none"></div><div dir="ltr">> > does changing the adsorption site of the reactant (SO2) to less stable site might solve the issue ?<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Regards<br clear="none"></div><div dir="ltr">> > On Monday, October 26, 2020, 11:25:36 AM GMT+4, Antoine Jay <<a rel="nofollow" shape="rect" ymailto="mailto:ajay@laas.fr" target="_blank" href="mailto:ajay@laas.fr">ajay@laas.fr</a>> wrote:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Dear Omer,<br clear="none"></div><div dir="ltr">> > I think your initial and final minima have not been well relaxed.<br clear="none"></div><div dir="ltr">> > When you fix the initial and final structures in a neb you must have relaxed them before, otherwise, you will have negative energy barriers.<br clear="none"></div><div dir="ltr">> > Moreover, you may have rotation of molecules that return local minima if your simulation box is too small.<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Regards,<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Antoine Jay<br clear="none"></div><div dir="ltr">> > LAAS-CNRS<br clear="none"></div><div dir="ltr">> > Toulouse, France<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Le Lundi, Octobre 26, 2020 06:39 CET, Omer Mutasim <<a rel="nofollow" shape="rect" ymailto="mailto:omermutasim@ymail.com" target="_blank" href="mailto:omermutasim@ymail.com">omermutasim@ymail.com</a>> a écrit:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Thanks a lot Dr. Tamas & Dr. Jay. , it is very efficient procedure, it worked for me now for all reaction steps. cheers<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > However for one elementary step , particularly SO2 dissociation ( SO2 = SO+O ) i got the following activation barrier, (it hasn't finished yet, but expected to remain around these values since it doesn't change much):<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > ------------------------------ iteration 297 ------------------------------<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > activation energy (->) = 0.000000 eV<br clear="none"></div><div dir="ltr">> > activation energy (<-) = 0.308512 eV<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > 1 -92402.4972907 0.036606 T<br clear="none"></div><div dir="ltr">> > 2 -92402.8008646 0.020347 F<br clear="none"></div><div dir="ltr">> > 3 -92402.6789202 0.048720 F<br clear="none"></div><div dir="ltr">> > 4 -92403.2726990 0.102631 F<br clear="none"></div><div dir="ltr">> > 5 -92403.0642888 0.050277 F<br clear="none"></div><div dir="ltr">> > 6 -92403.2377599 0.067121 F<br clear="none"></div><div dir="ltr">> > 7 -92402.8058032 0.029355 T<br clear="none"></div><div dir="ltr">> > activation energy (->) = 0.000000 eV<br clear="none"></div><div dir="ltr">> > activation energy (<-) = 0.308512 eV<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Attached is the MEP curve. As you see in MEP graph , there is oscillation in energies.<br clear="none"></div><div dir="ltr">> > is it normal to get this oscillated MEP curve for such barrier-less reaction step ? if not, how to get rid of this oscillations ?<br clear="none"></div><div dir="ltr">> > does using "CI" can increase this barrier a bit ?<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Thanks in advance<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Regards<br clear="none"></div><div dir="ltr">> > On Wednesday, October 21, 2020, 11:04:14 PM GMT+4, Omer Mutasim <<a rel="nofollow" shape="rect" ymailto="mailto:omermutasim@ymail.com" target="_blank" href="mailto:omermutasim@ymail.com">omermutasim@ymail.com</a>> wrote:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Very helpful ideas.<br clear="none"></div><div dir="ltr">> > But after pre-converging with inexpensive parameters, i will get first & last image that are different than my actual images with higher parameters ( k-pointss, cutoff,..)<br clear="none"></div><div dir="ltr">> > So then how i can use this pre-converged path for my actual settings?<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Sent from Yahoo Mail for iPhone<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > On Wednesday, October 21, 2020, 3:06 PM, Tamas Karpati <<a rel="nofollow" shape="rect" ymailto="mailto:tkarpati@gmail.com" target="_blank" href="mailto:tkarpati@gmail.com">tkarpati@gmail.com</a>> wrote:<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Dear Omar,<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > Hope it helps, just some ideas:<br clear="none"></div><div dir="ltr">> > - I could tell more if you would attach the whole input file (ie. the<br clear="none"></div><div dir="ltr">> > structures).<br clear="none"></div><div dir="ltr">> > - Without knowing the structures only I can give some hints:<br clear="none"></div><div dir="ltr">> > -- Try using smaller PW basis and lower ecutwfc, ecutrho to speed up<br clear="none"></div><div dir="ltr">> > your simulation.<br clear="none"></div><div dir="ltr">> > -- When you obtain something more reliable result, you can change<br clear="none"></div><div dir="ltr">> > back to the higher basis.<br clear="none"></div><div dir="ltr">> > -- Try leaving opt_scheme at its default value.<br clear="none"></div><div dir="ltr">> > -- For such a reaction (dissociation of such a polarized molecule) you should<br clear="none"></div><div dir="ltr">> > expect a barrier, therefore CI_scheme should be anything except for no-CI.<br clear="none"></div><div dir="ltr">> > -- The best is if you can specify the CI manually in the<br clear="none"></div><div dir="ltr">> > CLIMBING_IMAGES section<br clear="none"></div><div dir="ltr">> > (choose the CI_scheme accordingly).<br clear="none"></div><div dir="ltr">> > Bests,<br clear="none"></div><div dir="ltr">> > t<br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > On Tue, Oct 20, 2020 at 6:53 PM Omer Mutasim <<a rel="nofollow" shape="rect" ymailto="mailto:omermutasim@ymail.com" target="_blank" href="mailto:omermutasim@ymail.com">omermutasim@ymail.com</a>> wrote:<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > Dear All<br clear="none"></div><div dir="ltr">> > > I'm doning NEB for dissociation reaction of SO2 to SO +O. But it is not converging for more than a week, and the path length is increasing.<br clear="none"></div><div dir="ltr">> > > Please tell me what is wrong in my input file:<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > below is the input & output files:<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > Input file:<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > BEGIN<br clear="none"></div><div dir="ltr">> > > BEGIN_PATH_INPUT<br clear="none"></div><div dir="ltr">> > > &PATH<br clear="none"></div><div dir="ltr">> > > restart_mode = 'restart'<br clear="none"></div><div dir="ltr">> > > string_method = 'neb',<br clear="none"></div><div dir="ltr">> > > nstep_path = 800,<br clear="none"></div><div dir="ltr">> > > ds = 1.D0,<br clear="none"></div><div dir="ltr">> > > opt_scheme = "broyden",<br clear="none"></div><div dir="ltr">> > > num_of_images = 7,<br clear="none"></div><div dir="ltr">> > > CI_scheme = 'no-CI',<br clear="none"></div><div dir="ltr">> > > path_thr = 0.05D0,<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > /<br clear="none"></div><div dir="ltr">> > > END_PATH_INPUT<br clear="none"></div><div dir="ltr">> > > BEGIN_ENGINE_INPUT<br clear="none"></div><div dir="ltr">> > > &CONTROL<br clear="none"></div><div dir="ltr">> > > calculation = "relax"<br clear="none"></div><div dir="ltr">> > > prefix = 'SO2_neb'<br clear="none"></div><div dir="ltr">> > > outdir = './outdir'<br clear="none"></div><div dir="ltr">> > > pseudo_dir = '/home/yQE-test/pseudo/'<br clear="none"></div><div dir="ltr">> > > restart_mode = 'from_scratch'<br clear="none"></div><div dir="ltr">> > > forc_conv_thr = 1.0e-03<br clear="none"></div><div dir="ltr">> > > etot_conv_thr = 1e-04<br clear="none"></div><div dir="ltr">> > > nstep = 200<br clear="none"></div><div dir="ltr">> > > !tefield = .TRUE<br clear="none"></div><div dir="ltr">> > > !dipfield = .TRUE<br clear="none"></div><div dir="ltr">> > > /<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > &SYSTEM<br clear="none"></div><div dir="ltr">> > > ibrav = 0<br clear="none"></div><div dir="ltr">> > > ecutrho = 270<br clear="none"></div><div dir="ltr">> > > ecutwfc = 45<br clear="none"></div><div dir="ltr">> > > nat = 111<br clear="none"></div><div dir="ltr">> > > ntyp = 4<br clear="none"></div><div dir="ltr">> > > occupations='smearing',smearing='gaussian',degauss=0.005<br clear="none"></div><div dir="ltr">> > > vdw_corr = 'DFT-D2'<br clear="none"></div><div dir="ltr">> > > !edir = 3 , emaxpos = 0.6808, eopreg = 0.08 , eamp = 0.001,<br clear="none"></div><div dir="ltr">> > > nspin = 2<br clear="none"></div><div dir="ltr">> > > starting_magnetization(1)= 0.01<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > /<br clear="none"></div><div dir="ltr">> > > &ELECTRONS<br clear="none"></div><div dir="ltr">> > > conv_thr = 1e-06<br clear="none"></div><div dir="ltr">> > > electron_maxstep = 200<br clear="none"></div><div dir="ltr">> > > mixing_mode ='local-TF'<br clear="none"></div><div dir="ltr">> > > mixing_beta = 0.3<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > /<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > &IONS<br clear="none"></div><div dir="ltr">> > > /<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > K_POINTS {automatic}<br clear="none"></div><div dir="ltr">> > > 3 3 1 0 0 1<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ATOMIC_SPECIES<br clear="none"></div><div dir="ltr">> > > Ni 58.69340 Ni.pbe-n-rrkjus_psl.0.1.UPF<br clear="none"></div><div dir="ltr">> > > P 30.97376 P.pbe-n-rrkjus_psl.1.0.0.UPF<br clear="none"></div><div dir="ltr">> > > S 32.065 S.pbe-n-rrkjus_psl.1.0.0.UPF<br clear="none"></div><div dir="ltr">> > > O 15.9999 O.pbe-n-rrkjus_psl.1.0.0.UPF<br clear="none"></div><div dir="ltr">> > > CELL_PARAMETERS {angstrom}<br clear="none"></div><div dir="ltr">> > > 11.765383541833 0.0000000000 0.0000000000<br clear="none"></div><div dir="ltr">> > > -5.88269177091652 10.1891210324947 0.0000000000<br clear="none"></div><div dir="ltr">> > > 0.0000000000 0.0000000000 30.9938690567585<br clear="none"></div><div dir="ltr">> > > BEGIN_POSITIONS<br clear="none"></div><div dir="ltr">> > > FIRST_IMAGE<br clear="none"></div><div dir="ltr">> > > ATOMIC_POSITIONS (angstrom)<br clear="none"></div><div dir="ltr">> > > S -1.181561037 6.155418563 12.124345096<br clear="none"></div><div dir="ltr">> > > O -1.100425541 4.672437254 11.356300976<br clear="none"></div><div dir="ltr">> > > O 0.190308001 6.839217965 11.448732238<br clear="none"></div><div dir="ltr">> > > Ni -2.738525121 4.763450297 0.239145520<br clear="none"></div><div dir="ltr">> > > Ni 3.139579474 1.358483744 0.232252034<br clear="none"></div><div dir="ltr">> > > Ni 3.135766403 8.150575392 0.235327906<br clear="none"></div><div dir="ltr">> > > Ni -4.673593720 8.104467836 1.780118367<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > output file:<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > Program NEB v.6.4.1 starts on 16Oct2020 at 11:35:32<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > This program is part of the open-source Quantum ESPRESSO suite<br clear="none"></div><div dir="ltr">> > > for quantum simulation of materials; please cite<br clear="none"></div><div dir="ltr">> > > "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br clear="none"></div><div dir="ltr">> > > "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);<br clear="none"></div><div dir="ltr">> > > URL <a rel="nofollow" shape="rect" target="_blank" href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>",<br clear="none"></div><div dir="ltr">> > > in publications or presentations arising from this work. More details at<br clear="none"></div><div dir="ltr">> > > <a rel="nofollow" shape="rect" target="_blank" href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</a><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > Parallel version (MPI), running on 80 processors<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > MPI processes distributed on 5 nodes<br clear="none"></div><div dir="ltr">> > > R & G space division: proc/nbgrp/npool/nimage = 80<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > parsing_file_name: input.in<br clear="none"></div><div dir="ltr">> > > Reading input from pw_1.in<br clear="none"></div><div dir="ltr">> > > Message from routine read_upf::<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > initial path length = 11.3145 bohr<br clear="none"></div><div dir="ltr">> > > initial inter-image distance = 1.8857 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > string_method = neb<br clear="none"></div><div dir="ltr">> > > restart_mode = from_scratch<br clear="none"></div><div dir="ltr">> > > opt_scheme = broyden<br clear="none"></div><div dir="ltr">> > > num_of_images = 7<br clear="none"></div><div dir="ltr">> > > nstep_path = 800<br clear="none"></div><div dir="ltr">> > > CI_scheme = no-CI<br clear="none"></div><div dir="ltr">> > > first_last_opt = F<br clear="none"></div><div dir="ltr">> > > use_freezing = F<br clear="none"></div><div dir="ltr">> > > ds = 1.0000 a.u.<br clear="none"></div><div dir="ltr">> > > k_max = 0.1000 a.u.<br clear="none"></div><div dir="ltr">> > > k_min = 0.1000 a.u.<br clear="none"></div><div dir="ltr">> > > suggested k_max = 0.6169 a.u.<br clear="none"></div><div dir="ltr">> > > suggested k_min = 0.6169 a.u.<br clear="none"></div><div dir="ltr">> > > path_thr = 0.0500 eV / A<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 1 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 6.2 self-consistency for image 1<br clear="none"></div><div dir="ltr">> > > tcpu = 3675.5 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 7662.5 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 11422.7 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 15346.3 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 19108.7 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > > tcpu = 22571.1 self-consistency for image 7<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 70.216194 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 71.022062 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92468.8536637 23.505267 F<br clear="none"></div><div dir="ltr">> > > 3 -92442.9691259 150.213122 F<br clear="none"></div><div dir="ltr">> > > 4 -92406.7311409 330.353055 F<br clear="none"></div><div dir="ltr">> > > 5 -92431.0052901 207.333777 F<br clear="none"></div><div dir="ltr">> > > 6 -92469.0661237 51.663167 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 11.314 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 1.886 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 2 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 26119.7 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 28731.5 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 31027.4 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 34094.2 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 36988.0 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 22.531451 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 23.337319 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92469.5101428 22.300995 F<br clear="none"></div><div dir="ltr">> > > 3 -92454.4158842 70.627594 F<br clear="none"></div><div dir="ltr">> > > 4 -92461.1206812 34.307062 F<br clear="none"></div><div dir="ltr">> > > 5 -92464.4669859 46.783708 F<br clear="none"></div><div dir="ltr">> > > 6 -92471.4896125 37.765708 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 11.384 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 1.897 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 3 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 39172.0 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 41888.0 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 44777.8 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 47642.0 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 50615.2 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 13.435341 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 14.241209 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92471.6434742 16.119604 F<br clear="none"></div><div dir="ltr">> > > 3 -92463.5119937 28.367753 F<br clear="none"></div><div dir="ltr">> > > 4 -92468.1466546 16.740841 F<br clear="none"></div><div dir="ltr">> > > 5 -92472.7705146 11.019872 F<br clear="none"></div><div dir="ltr">> > > 6 -92475.3040517 10.662908 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 11.502 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 1.917 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 4 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 53323.8 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 56077.9 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 59014.9 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 61990.6 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 64608.8 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 6.530687 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 7.336554 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92474.0378392 7.910468 F<br clear="none"></div><div dir="ltr">> > > 3 -92470.4166483 13.061889 F<br clear="none"></div><div dir="ltr">> > > 4 -92471.2528453 9.923078 F<br clear="none"></div><div dir="ltr">> > > 5 -92474.2165523 4.209611 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.2787664 3.450159 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 11.724 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 1.954 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 5 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 67273.9 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 70152.2 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 73153.1 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 76203.5 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 78824.8 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 3.710859 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 4.516727 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92475.5387437 3.615458 F<br clear="none"></div><div dir="ltr">> > > 3 -92473.5622985 4.695416 F<br clear="none"></div><div dir="ltr">> > > 4 -92473.2364760 5.734324 F<br clear="none"></div><div dir="ltr">> > > 5 -92474.4833513 2.877608 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.4425910 2.876700 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 12.038 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 2.006 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 6 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 81355.0 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 84402.5 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 87564.5 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 90568.2 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 93110.5 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 2.560838 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 3.366706 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.0900010 1.040106 F<br clear="none"></div><div dir="ltr">> > > 3 -92474.7832671 2.988289 F<br clear="none"></div><div dir="ltr">> > > 4 -92474.3864972 2.085630 F<br clear="none"></div><div dir="ltr">> > > 5 -92474.8266397 1.998808 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.6377600 0.667994 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 12.364 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 2.061 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 7 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 95248.7 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 98189.3 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 101337.9 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 104423.7 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 107076.7 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 2.125802 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 2.931670 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.0736630 1.319140 F<br clear="none"></div><div dir="ltr">> > > 3 -92475.4151167 1.955048 F<br clear="none"></div><div dir="ltr">> > > 4 -92474.8215329 1.921925 F<br clear="none"></div><div dir="ltr">> > > 5 -92475.0627346 2.135695 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.7117640 0.696381 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 12.868 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 2.145 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 8 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 108885.3 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 111194.4 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 113961.2 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 116506.3 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 118361.2 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 2.073805 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 2.879673 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.1622863 0.861666 F<br clear="none"></div><div dir="ltr">> > > 3 -92475.4162307 2.192183 F<br clear="none"></div><div dir="ltr">> > > 4 -92474.8735300 1.631538 F<br clear="none"></div><div dir="ltr">> > > 5 -92475.0684015 1.824977 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.7113576 0.606060 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 12.633 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 2.105 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > .<br clear="none"></div><div dir="ltr">> > > reading file 'SO2_neb.path'<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > string_method = neb<br clear="none"></div><div dir="ltr">> > > restart_mode = restart<br clear="none"></div><div dir="ltr">> > > opt_scheme = broyden<br clear="none"></div><div dir="ltr">> > > num_of_images = 7<br clear="none"></div><div dir="ltr">> > > nstep_path = 800<br clear="none"></div><div dir="ltr">> > > CI_scheme = no-CI<br clear="none"></div><div dir="ltr">> > > first_last_opt = F<br clear="none"></div><div dir="ltr">> > > use_freezing = F<br clear="none"></div><div dir="ltr">> > > ds = 1.0000 a.u.<br clear="none"></div><div dir="ltr">> > > k_max = 0.1000 a.u.<br clear="none"></div><div dir="ltr">> > > k_min = 0.1000 a.u.<br clear="none"></div><div dir="ltr">> > > suggested k_max = 0.6169 a.u.<br clear="none"></div><div dir="ltr">> > > suggested k_min = 0.6169 a.u.<br clear="none"></div><div dir="ltr">> > > path_thr = 0.0500 eV / A<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 26 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 6.2 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 3713.3 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 7137.5 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 10796.2 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 14447.3 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 0.640765 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 1.446632 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.3065704 0.250516 F<br clear="none"></div><div dir="ltr">> > > 3 -92477.0673165 0.278078 F<br clear="none"></div><div dir="ltr">> > > 4 -92476.7787332 0.431860 F<br clear="none"></div><div dir="ltr">> > > 5 -92476.4800774 0.453182 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.5576488 0.296200 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 17.545 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 2.924 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 27 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 18237.5 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 20736.0 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 23008.1 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 25721.8 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 28310.3 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 0.643277 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 1.449145 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.3040583 0.235899 F<br clear="none"></div><div dir="ltr">> > > 3 -92477.0080434 0.474599 F<br clear="none"></div><div dir="ltr">> > > 4 -92476.8143156 0.678632 F<br clear="none"></div><div dir="ltr">> > > 5 -92476.4592005 0.597470 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.4827638 0.224064 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 18.367 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 3.061 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 28 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 30382.0 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 32498.8 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 34597.5 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 37250.7 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 39649.7 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 0.651733 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 1.457601 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.2956021 0.236096 F<br clear="none"></div><div dir="ltr">> > > 3 -92476.9388664 0.680956 F<br clear="none"></div><div dir="ltr">> > > 4 -92476.8025379 0.874373 F<br clear="none"></div><div dir="ltr">> > > 5 -92476.3933083 0.734403 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.4272820 0.239132 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 19.115 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 3.186 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 29 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 41622.7 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 43787.2 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 45892.1 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 48482.6 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 50617.1 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 0.661553 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 1.467420 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.2857825 0.249692 F<br clear="none"></div><div dir="ltr">> > > 3 -92476.8823826 0.778237 F<br clear="none"></div><div dir="ltr">> > > 4 -92476.7843580 1.002202 F<br clear="none"></div><div dir="ltr">> > > 5 -92476.3323697 0.748960 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.3885082 0.238984 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 19.742 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 3.290 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 30 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 52474.3 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 54679.9 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 57012.5 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 59877.3 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > > tcpu = 62490.2 self-consistency for image 6<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > activation energy (->) = 0.704760 eV<br clear="none"></div><div dir="ltr">> > > activation energy (<-) = 1.510628 eV<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > image energy (eV) error (eV/A) frozen<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > 1 -92476.9473351 0.023792 T<br clear="none"></div><div dir="ltr">> > > 2 -92476.2668397 0.342411 F<br clear="none"></div><div dir="ltr">> > > 3 -92476.7810889 0.907920 F<br clear="none"></div><div dir="ltr">> > > 4 -92476.7414553 1.153276 F<br clear="none"></div><div dir="ltr">> > > 5 -92476.2425749 0.898295 F<br clear="none"></div><div dir="ltr">> > > 6 -92476.3370447 0.341313 F<br clear="none"></div><div dir="ltr">> > > 7 -92477.7532028 0.024858 T<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > path length = 20.745 bohr<br clear="none"></div><div dir="ltr">> > > inter-image distance = 3.457 bohr<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > ------------------------------ iteration 31 ------------------------------<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > tcpu = 64431.3 self-consistency for image 2<br clear="none"></div><div dir="ltr">> > > tcpu = 66430.0 self-consistency for image 3<br clear="none"></div><div dir="ltr">> > > tcpu = 68491.4 self-consistency for image 4<br clear="none"></div><div dir="ltr">> > > tcpu = 70987.8 self-consistency for image 5<br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > ><br clear="none"></div><div dir="ltr">> > > _______________________________________________<br clear="none"></div><div dir="ltr">> > > Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> > > users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> > > <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">> > _______________________________________________<br clear="none"></div><div dir="ltr">> > Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> > users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> > <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > _______________________________________________<br clear="none"></div><div dir="ltr">> > Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> > users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> > <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> ><br clear="none"></div><div dir="ltr">> > _______________________________________________<br clear="none"></div><div dir="ltr">> > Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> > users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> > <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">> _______________________________________________<br clear="none"></div><div dir="ltr">> Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">> _______________________________________________<br clear="none"></div><div dir="ltr">> Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">> users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr">> <a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br clear="none"></div><div dir="ltr">_______________________________________________<br clear="none"></div><div dir="ltr">Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none"></div><div dir="ltr">users mailing list <a rel="nofollow" shape="rect" ymailto="mailto:users@lists.quantum-espresso.org" target="_blank" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br clear="none"></div><div dir="ltr"><a rel="nofollow" shape="rect" target="_blank" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></div>
</div>
</div></div></div></div></div>
</div>
</div></body></html>