<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2800.1170" name=GENERATOR></HEAD>
<BODY>
<DIV>Dear CPMD users,</DIV>
<DIV>I want to computer the binding energy of ammonium-benzene complex. It is well know </DIV>
<DIV>this complex contains three stable structures (we denote them as 1, 2 and 3, 1 means one </DIV>
<DIV>hydrogen atom of ammonium toward benzene, 2 means two hydrogen atoms of ammonium</DIV>
<DIV>toward benzene, 3 means three hydrogen atoms of ammonium toward benzene). Previous </DIV>
<DIV>B3LYP and MP2 calculation using large basis sets shown 2 is most stable and CPMD/70Ry</DIV>
<DIV>calculation also proven this. The calculated binding energy of 2 is also similar. But for 1 and 3,</DIV>
<DIV>CPMD and Gaussian program give different results. In B3LYP and MP2 calculation, 3 is more </DIV>
<DIV>stable than 1. While in CPMD calculation, 1 is more stable than 3. Also in CPMD calculation, </DIV>
<DIV>with the increasing of cell size, the deviation from experimental results increasing. I have tried </DIV>
<DIV>to increasing the convergence of optimization. But it seems little effect. I wonder if you can </DIV>
<DIV>give me some advice about my calculation. </DIV>
<DIV>Input file of 1:</DIV>
<DIV><FONT size=2></FONT> </DIV>
<DIV><FONT size=2>&CPMD</FONT></DIV>
<DIV><FONT size=2>RESTART WAVEFUNCTION LATEST</FONT></DIV>
<DIV><FONT size=2>OPTIMIZE GEOMETRY</FONT></DIV>
<DIV><FONT size=2>CONVERGENCE</FONT></DIV>
<DIV><FONT size=2>1.D-5 5.D-5</FONT></DIV>
<DIV><FONT size=2> EMASS</FONT></DIV>
<DIV><FONT size=2> 500</FONT></DIV>
<DIV><FONT size=2> MIRROR</FONT></DIV>
<DIV><FONT size=2>&END</FONT></DIV>
<DIV><FONT size=2> &SYSTEM</FONT></DIV>
<DIV><FONT size=2> SYMMETRY</FONT></DIV>
<DIV><FONT size=2> 1</FONT></DIV>
<DIV><FONT size=2> CELL</FONT></DIV>
<DIV><FONT
size=2> 15.1178 1.0 1.0 0 0 0</FONT></DIV>
<DIV><FONT size=2> CUTOFF</FONT></DIV>
<DIV><FONT size=2> 70.0</FONT></DIV>
<DIV><FONT size=2> CHARGE</FONT></DIV>
<DIV><FONT size=2> 1</FONT></DIV>
<DIV><FONT size=2>&END</FONT></DIV>
<DIV><FONT size=2>&ATOMS</FONT></DIV>
<DIV><FONT size=2>*C_MT_GIA_BLYP KLEINMAN-BYLANDER</FONT></DIV>
<DIV><FONT size=2>LMAX=P</FONT></DIV>
<DIV><FONT size=2>6</FONT></DIV>
<DIV><FONT size=2> 0.0023 2.6441 -1.2327</FONT></DIV>
<DIV><FONT size=2> 2.2960 1.3201 -1.2362</FONT></DIV>
<DIV><FONT size=2> 2.2955 -1.3326 -1.2336</FONT></DIV>
<DIV><FONT size=2> 0.0009 -2.6548 -1.2281</FONT></DIV>
<DIV><FONT size=2>-2.2930 -1.3312 -1.2271</FONT></DIV>
<DIV><FONT size=2>-2.2922 1.3216 -1.2297</FONT></DIV>
<DIV><FONT size=2>*N_MT_GIA_BLYP KLEINMAN-BYLANDER</FONT></DIV>
<DIV><FONT size=2> LMAX=P</FONT></DIV>
<DIV><FONT size=2> 1</FONT></DIV>
<DIV><FONT size=2>-0.0069 0.0238 4.7434</FONT></DIV>
<DIV><FONT size=2>*H_MT_BLYP KLEINMAN-BYLANDER</FONT></DIV>
<DIV><FONT size=2> LMAX=S</FONT></DIV>
<DIV><FONT size=2>10</FONT></DIV>
<DIV><FONT size=2> 0.0032 4.7040 -1.2833</FONT></DIV>
<DIV><FONT size=2> 4.0772 2.3470 -1.2920</FONT></DIV>
<DIV><FONT size=2> 4.0760 -2.3608 -1.2869</FONT></DIV>
<DIV><FONT size=2> 0.0004 -4.7150 -1.2769</FONT></DIV>
<DIV><FONT size=2>-4.0740 -2.3586 -1.2757</FONT></DIV>
<DIV><FONT size=2>-4.0728 2.3496 -1.2809</FONT></DIV>
<DIV><FONT size=2>-0.0029 -0.0039 2.7647</FONT></DIV>
<DIV><FONT size=2>-1.6017 0.9524 5.3801</FONT></DIV>
<DIV><FONT size=2> 1.5842 0.9530 5.3884</FONT></DIV>
<DIV><FONT size=2>-0.0085 -1.8058 5.4252</FONT></DIV>
<DIV><FONT size=2>&END</FONT></DIV>
<DIV><FONT size=2>&DFT</FONT></DIV>
<DIV><FONT size=2> NEWCODE</FONT></DIV>
<DIV><FONT size=2> FUNCTIONAL BLYP</FONT></DIV>
<DIV><FONT size=2>&END</FONT></DIV>
<DIV><FONT size=2></FONT> </DIV>
<DIV><FONT size=2>Input file for 2:</FONT></DIV>
<DIV><FONT size=2></FONT> </DIV>
<DIV><FONT
size=2> &ATOMS
<BR> *C_MT_GIA_BLYP KLEINMAN-BYLANDER
<BR> LMAX=P
<BR> 6
<BR> 0.0001 2.6507
-1.1977
<BR> -0.0005 -2.6471
-1.1761
<BR> -2.2965 1.3281
-1.1968
<BR> 2.2964 1.3277
-1.1951
<BR> 2.2963 -1.3242
-1.1845
<BR> -2.2970 -1.3237
-1.1861
<BR>
<BR> *N_MT_GIA_BLYP KLEINMAN-BYLANDER
<BR>
LMAX=P
<BR>
1
<BR> 0.0008 -0.0084
4.6011
<BR>
<BR> *H_MT_BLYP
KLEINMAN-BYLANDER
<BR>
LMAX=S
<BR> 10
<BR> 0.0003 4.7116
-1.2507
<BR> -0.0006 -4.7083
-1.2129
<BR> -4.0767 2.3580
-1.2364
<BR> 4.0768 2.3573
-1.2334
<BR> 4.0767 -2.3541
-1.2144
<BR> -4.0774 -2.3534
-1.2175
<BR> -0.0017 1.5778
3.4438
<BR> 0.0049 -1.5530
3.3854
<BR> -1.5944 -0.0313
5.7201
<BR> 1.5946 -0.0247
5.7222
<BR>
<BR> &END
</FONT></DIV>
<DIV><FONT size=2>Input file for 3:</FONT></DIV>
<DIV><FONT size=2>*C_MT_GIA_BLYP KLEINMAN-BYLANDER
<BR>LMAX=P
<BR>6
<BR> 0.0000 1.2865
2.6672
<BR> 2.2953 1.2470
1.3452
<BR> 2.2956 1.1696
-1.3075
<BR> 0.0000 1.1273
-2.6313
<BR>-2.2956 1.1696
-1.3075
<BR>-2.2953 1.2470
1.3452
<BR>
<BR>*N_MT_GIA_BLYP KLEINMAN-BYLANDER <BR>
LMAX=P
<BR> 1
<BR>0.0000 -4.6951
-0.0804
<BR>
<BR>*H_MT_BLYP KLEINMAN-BYLANDER
<BR>
LMAX=S
<BR>10
<BR> 0.0000 1.4037
4.7256
<BR> 4.0760 1.3412
2.3727
<BR> 4.0774 1.1969
-2.3365
<BR> 0.0000 1.1327
-4.6935
<BR>-4.0774 1.1969
-2.3365
<BR>-4.0760 1.3412
2.3727
<BR> 0.0000 -6.6392
-0.2002
<BR>-1.5821 -4.0530
0.8709
<BR> 1.5821 -4.0530
0.8709
<BR> 0.0000 -3.8717 -1.8560
<BR>&END
</FONT></DIV></BODY></HTML>