<div dir="ltr"><div dir="ltr">Dear QE users,<div><br></div><div>I am trying to optimize (relax) Pd dimer (Pd2) as cluster placed on cubic lattice.</div><div>However I got a value of Pd-Pd bond length which is far too long (3.36 A) compare to the experimental one (about 2.6 A).</div><div>Am I doing something wrong? Here below the input file:</div><div><div><br></div><div>&CONTROL</div><div>     calculation   = "relax"</div><div>prefix ='Pd2'</div><div>    forc_conv_thr =  1.94469e-03</div><div>    max_seconds   =  1.34369e+14</div><div>    nstep         = 800</div><div>    pseudo_dir    = "/home/pcosta/pseudo"</div><div><span style="white-space:pre">       </span>verbosity     ='high'</div><div>disk_io ='none'</div><div>    tprnfor       = .TRUE.</div><div>    tstress       = .TRUE.</div><div>/</div><div><br></div><div>&SYSTEM</div><div>    a     =  4.50000e+01</div><div>    degauss                   =  1.00000e-02</div><div>    ecutrho                   =  1.83747e+01</div><div>    ecutwfc                   =  1.46997e+00</div><div><span style="white-space:pre">  </span>lda_plus_u                = .FALSE.</div><div>    ibrav                     = 1</div><div>    nat                       = 2</div><div>    nspin                     = 2</div><div>    ntyp                      = 1</div><div>    occupations               = "smearing"</div><div>    smearing                  = "gaussian"</div><div>    starting_magnetization(1) =  2.00000e-01</div><div>    assume_isolated = 'makov-payne'</div><div><span style="white-space:pre"> </span>vdw_corr                  = 'grimme-d3'</div><div>/</div><div><br></div><div>&ELECTRONS</div><div>    conv_thr         =  1.00000e-08</div><div>    diagonalization  = "david"</div><div>    electron_maxstep = 800</div><div>    mixing_beta      =  4.00000e-01</div><div>    startingpot      = "atomic"</div><div>    startingwfc      = "atomic+random"</div><div>/</div><div><br></div><div>&IONS</div><div>    ion_dynamics = "bfgs"</div><div><span style="white-space:pre">   </span>trust_radius_min=1e-5</div><div>/</div><div><br></div><div>&CELL</div><div>    cell_dofree    = "all"</div><div>    cell_dynamics  = "bfgs"</div><div>    press_conv_thr =  5.00000e-01</div><div>/</div><div><br></div><div>K_POINTS {automatic}</div><div> 3  3  1  0 0 0</div><div>  </div><div>ATOMIC_SPECIES</div><div>Pd    106.42000  Pd.pbe-n-kjpaw_psl.1.0.0.UPF</div><div><br></div><div>ATOMIC_POSITIONS {angstrom}</div><div>Pd     16.194860  19.259990  18.659720</div><div>Pd     18.251640  18.665930  18.594650</div><div><br></div><div>Thanks.</div>-- <br><div dir="ltr" class="gmail-m_6880443592494725121gmail_signature"><div dir="ltr">Paolo Costa, Ph.D.<div>Postdoctoral Researcher</div><div>Department of Chemistry and Biomolecular Sciences</div><div>University of Ottawa</div><div>10 Marie Curie, Ottawa, ON K1N 6N5, Canada</div><div>Room number: DRO 326 (D'Iorio Hall)</div></div></div></div></div></div>