<div dir="ltr">Dear Giuseppe,<div><br></div><div>thanks a lot!</div><div><br></div><div>Paolo</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno lun 28 gen 2019 alle ore 11:05 Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
Dear Paolo<br>
It's a miracle that your calculation went somewhere... ;-D<br>
You are using a plane-wave cutoff of 1.47 Ry and a density cutoff of <br>
18.37 Ry. I would expect higher values, and it is always a good idea <br>
to test the convergence of your pseudopotentials.<br>
As a suggestion, don't overuse scientific notation, unless you can't <br>
help. It generally leads to "distraction errors" such as these.<br>
HTH<br>
Giuseppe<br>
<br>
Paolo Costa <<a href="mailto:paolo.costa85@gmail.com" target="_blank">paolo.costa85@gmail.com</a>> ha scritto:<br>
<br>
> Dear QE users,<br>
><br>
> I am trying to optimize (relax) Pd dimer (Pd2) as cluster placed on cubic<br>
> lattice.<br>
> However I got a value of Pd-Pd bond length which is far too long (3.36 A)<br>
> compare to the experimental one (about 2.6 A).<br>
> Am I doing something wrong? Here below the input file:<br>
><br>
> &CONTROL<br>
> calculation = "relax"<br>
> prefix ='Pd2'<br>
> forc_conv_thr = 1.94469e-03<br>
> max_seconds = 1.34369e+14<br>
> nstep = 800<br>
> pseudo_dir = "/home/pcosta/pseudo"<br>
> verbosity ='high'<br>
> disk_io ='none'<br>
> tprnfor = .TRUE.<br>
> tstress = .TRUE.<br>
> /<br>
><br>
> &SYSTEM<br>
> a = 4.50000e+01<br>
> degauss = 1.00000e-02<br>
> ecutrho = 1.83747e+01<br>
> ecutwfc = 1.46997e+00<br>
> lda_plus_u = .FALSE.<br>
> ibrav = 1<br>
> nat = 2<br>
> nspin = 2<br>
> ntyp = 1<br>
> occupations = "smearing"<br>
> smearing = "gaussian"<br>
> starting_magnetization(1) = 2.00000e-01<br>
> assume_isolated = 'makov-payne'<br>
> vdw_corr = 'grimme-d3'<br>
> /<br>
><br>
> &ELECTRONS<br>
> conv_thr = 1.00000e-08<br>
> diagonalization = "david"<br>
> electron_maxstep = 800<br>
> mixing_beta = 4.00000e-01<br>
> startingpot = "atomic"<br>
> startingwfc = "atomic+random"<br>
> /<br>
><br>
> &IONS<br>
> ion_dynamics = "bfgs"<br>
> trust_radius_min=1e-5<br>
> /<br>
><br>
> &CELL<br>
> cell_dofree = "all"<br>
> cell_dynamics = "bfgs"<br>
> press_conv_thr = 5.00000e-01<br>
> /<br>
><br>
> K_POINTS {automatic}<br>
> 3 3 1 0 0 0<br>
><br>
> ATOMIC_SPECIES<br>
> Pd 106.42000 Pd.pbe-n-kjpaw_psl.1.0.0.UPF<br>
><br>
> ATOMIC_POSITIONS {angstrom}<br>
> Pd 16.194860 19.259990 18.659720<br>
> Pd 18.251640 18.665930 18.594650<br>
><br>
> Thanks.<br>
> --<br>
> Paolo Costa, Ph.D.<br>
> Postdoctoral Researcher<br>
> Department of Chemistry and Biomolecular Sciences<br>
> University of Ottawa<br>
> 10 Marie Curie, Ottawa, ON K1N 6N5, Canada<br>
> Room number: DRO 326 (D'Iorio Hall)<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_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>