<br clear="all"><div> Dear Prof. Stefano Baroni and Vega Lew<br><br><br> Many many thanks for your detailed answers and helpful comments. This is scf energy of each relaxation steps :<br></div>
<div> <br> -411.0652023864<br>-411.1576447124<br>-411.1877929276<br>-411.1989785058<br>-411.1027361193**<br>-411.2006610355<br>-411.2078087316<br>-411.2087637085<br>-411.2094335080<br>-411.2112728535<br>-411.2137042660<br>
-411.2145098860<br>-411.2163180058<br>-411.2200044391<br>-411.2186337506**<br>-411.2238781353<br>-411.2272606876<br>-411.2335499997<br>-411.2410057850<br>-411.2486666714<br>-411.2432716037**<br>-411.2545344840<br>-411.2632404207<br>
-411.2686094326<br>-411.2697596541<br>-411.2688830229**<br>-411.2723156669<br>-411.2728025269<br>-411.2731592912<br>-411.2735337540<br>-411.2741229900<br>-411.2749931725<br>-411.2755105410<br>-411.2761617487<br>-411.2768678588<br>
-411.2772028256<br>-411.2775503379<br>-411.2785154793<br>-411.2823072885<br>-411.2298751795**<br>-411.2825004743<br>-411.2845206434<br>-411.2875009929<br>-411.2801470356**<br>-411.2954982819<br>-411.3096061836<br>-411.3466223605<br>
-411.367570164<br>-411.3724319076<br>-411.3795692495<br>-411.3812435563<br>-411.3823959522<br>-411.3829995766<br>-411.3833540501<br>-411.3835954004<br>-411.3839027669<br>-411.3841960285<br>-411.3843908495<br>-411.3844986818<br>
-411.3845647753<br>-411.3846244537<br>-411.3846644112<br>-411.3846945470<br>-411.3847346519 <br>-411.3848191008<br>-411.3848652595<br>-411.3849367730<br>-411.3849811533<br>-411.3851262954 <br>-411.3534987061**<br>-411.3850020755<br>
-411.3851753523<br>-411.3852665074<br>-411.3853358233<br>-411.3855056618<br>-411.3858383834<br>-411.3860833880<br>-411.3866016967<br>-411.3874393558<br>-411.3878819548<br>-411.3881741927<br>-411.3884732905<br>-411.3886017423<br>
-411.3886998952<br>-411.3887918295<br>-411.3888586159 <br>-411.3888992755<br>-411.3889304888<br>-411.3889622530<br>-411.3889906826 <br>-411.3890146960<br>-411.3890447538<br>-411.3890806506 <br>-411.3891138607<br>-411.3891306978<br>
-411.3891510321<br>-411.3891853920 <br>-411.3892253795<br>-411.3892680354<br>-411.3893152439<br>-411.3893766070<br>-411.3894890010<br>-411.3895516220<br>-411.3896539375<br>-411.3896995617<br>-411.3897936720<br>-411.3898406152<br>
-411.3898796419<br>-411.3899004418 <br></div>
<div> </div>
<div> ** shows the step in which the energy rises. I don't think that these points indicate a TS or something like that since one may see such points through any optimization procedure; am I correct ? <br><br> Also, here is the input : <br>
<br> &CONTROL<br> calculation = "relax",<br> pseudo_dir = "/home/soft/qe4.2/espresso-4.2/pseudo",<br> outdir = "/home/AB/tmp",<br> restart_mode = 'from_scratch'<br> etot_conv_thr= 1.0D-4,<br>
forc_conv_thr= 1.0D-3,<br> nstep=120<br> /<br>&SYSTEM<br> ibrav = 4,<br> a = 4.864990484,<br> b = 4.864990484,<br> c = 16.418933296,<br> cosab = -0.5,<br> cosac = 1.0,<br>
cosbc = 1.0,<br> nat = 26,<br> ntyp = 3,<br> ecutwfc = 40.D0,<br> ecutrho = 480.D0,<br> occupations = 'smearing'<br> smearing ='mp', <br> degauss = 0.03,<br> nspin = 2,<br> starting_magnetization(1)= 0.01,<br>
starting_magnetization(2)= 0.01,<br> starting_magnetization(3)= 0.01,<br> london=.true.,<br> nosym=.true.<br> /<br>&ELECTRONS<br> conv_thr = 1.D-6,<br> mixing_beta = 0.7D0,<br> diagonalization = "david",<br>
/<br>&IONS<br> ion_dynamics="bfgs"<br>/<br>ATOMIC_SPECIES<br>C 12.0107 C.pbe-rrkjus.UPF<br>Ni 58.69 Ni.pbe-nd-rrkjus.UPF<br>O 15.9994 O.pbe-rrkjus.UPF<br>ATOMIC_POSITIONS {angstrom}<br>C 0.000000000 0.000000000 0.000000000 0 0 0<br>
C 2.461165257 -0.001796837 0.002770193<br>C -1.232804188 2.130698110 0.002678921<br>C 1.230011567 2.131364695 0.003376901<br>C -0.003333507 1.422214727 0.001510708<br>C -1.234100323 3.553382228 0.001755289<br>
C 2.458567713 1.421054385 0.003703007<br>C 1.228614904 3.554143930 0.002513562<br>C -0.235163108 0.160826352 3.247131267<br>C 2.244119337 0.161060766 3.248352597<br>C -1.457317029 2.298122349 3.233987380<br>
C 1.004488573 2.308122685 3.246477338<br>C 1.003490453 0.877294779 3.260744638<br>C -0.228726740 3.007907781 3.234369538<br>C 3.466400555 0.879138405 3.235278922<br>C 2.237111888 3.008093688 3.233133671<br>
C 1.014869884 -0.620629524 6.400327421<br>C 3.476740887 -0.544001537 6.494151711<br>C -0.280354238 1.625208172 6.378065136<br>C 2.312052832 1.626645989 6.400900301<br>Ni 1.023326058 0.872911705 5.337347244<br>
O 1.023326000 0.872911000 8.700000000<br>O 2.012052000 1.326645000 8.700000000<br>C -0.206008752 3.013580161 6.506783380<br>C 3.478548719 0.866737695 6.505990211<br>C 2.235793764 3.015530163 6.513155404 <br>
K_POINTS {automatic}<br>4 4 1 1 1 1<br><br><div>
<div><i><font color="#000066" face="comic sans ms,sans-serif">-----------------</font></i></div>
<div><i><font color="#000066" face="comic sans ms,sans-serif">Masoud Nahali</font></i></div>
<div><font size="1"><font face="times new roman,serif"><font size="2"><i><font size="1"><font face="times new roman,serif"><font size="2"><i>PhD Student of Physical Chemistry<br>Sharif University of Technology</i> <br></font></font></font> </i></font></font></font></div>
</div> <br></div>
<div><br></div><br>
<br> Stefano Baroni wrote :<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
On Sep 29, 2010, at 5:32 PM, vega lew wrote:<br>
<br>
> Through the relaxation of a oxygen molecule on a surface by QE4.2, it dissociated and adsorbed on the surface. The scf energy of each step has decreased in compared with the previous step until the relaxation completed. Now, can one say that there is not any TS for O2 dissociation here ?<br>
><br>
>> Probably not. The potential energy surface might be not very corrugated and the energy barrier for the dissociation might be not very large, so that the dissociation of O2 could happen in a relax calculation. I think there might be a very minute TS in the dissociation process.<br>
<br>
here, I beg to differ. If a state, say B, can be reached from another state, say A, following a path along which the energy always decreases, this means *by definition* that there is *no* energy barrier between A and B. I admit that if the energy barrier is so "minute" as to be smaller that the numerical accuracy of the minimization algorithm, it can escape proper detection. (notice that this can be said of most, if not all, numerical computations of any quantity)<br>
<br>
><br>
> Or is the reaction spontaneous ?<br>
<br>
It is spontaneous<br>
<br>
>> If the total energy of system goes down continuously and the adsorption energy of dissociated state is lower than the one of molecular state, you could say the reaction is favorable. Whether it is spontaneous, it is depend on how large of the energy barrier.<br>
<br>
in the reported conditions, the energy barrier is zero<br>
<br>
>> If there was a TS in this procedure (dissociation and adsorption of O2) would the optimization algorithm of QE4.2 stop the work or continue ?<br>
<br>
In the presence of an energy barrier, the system would stay stuck in the starting configuration<br>
<br>
Stefano<br>
<br>
> If there is only one bond broken in the process you considered, you could try constrained optimization technique. For you information, please refer to B.Hammer (Hammer, B.; Jacobsen, K. W.; Norskov, J. K. Phys. Rev. Lett. 1992, 69, 1971) or PJ Hu's work. (Alavi, A.; Hu, P.; Deutsch, T.; Silvestrelli, P. L.; Hutter, J. Phys. Rev. Lett. 1998, 80, 3650). Or recent developed method by ZP Liu (J. Chem. Theory Comput. 2010, 6, 1136?1144 )<br>
<br>
If I have not misunderstood some important point, no constrained minimization or other advanced techniques are needed here.<br>
<br>
Stefano B<br>
<br></blockquote></div><br> <br><div><br></div>