<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">from: <a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idm6452736" class="">http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idm6452736</a><div class=""><br class=""></div><div class="">nstep -> number of molecular-dynamics or structural optimization steps</div><div class="">performed in this run</div><div class=""><br class=""></div><div class="">So if you set nstep=1 just one optimisation step is carried out.</div><div class=""><br class=""></div><div class="">Giovanni</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 7 Oct 2016, at 05:25, quantum <<a href="mailto:quantum@mail.ncyu.edu.tw" class="">quantum@mail.ncyu.edu.tw</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class=""><br class="">Dear QE Friends:<br class=""><br class="">I am trying to optimize a crystal by allowing the atomic positions of hydrogens to relax. All other atoms of the crystal are held fixed. The QE input file is listed by the end of this message.<br class=""><br class="">QE starts with the self-consistent calculation using the initial guess of atomic positions. After the first scf cycle, QE prints out a new set of atomic positions, and ends without a new scf cycle with the new set of atomic positions. (Please see the QE output attached below.)<br class=""><br class="">The EQ output suggests that the forces acting on all hydrogen atoms are greater than "forc_conv_thr = 1.0D-4." <br class=""><br class="">Why does the simulation end with only one scf cycle?<br class=""><br class="">Thank you in advance.<br class=""><br class=""><br class="">**** QE output file (partial) starts here ******************<br class=""><br class="">! total energy = -1426.31830754 Ry<br class=""> Harris-Foulkes estimate = -1426.31830738 Ry<br class=""> estimated scf accuracy < 0.00000009 Ry<br class=""><br class=""> The total energy is the sum of the following terms:<br class=""><br class=""> one-electron contribution = -615.01611222 Ry<br class=""> hartree contribution = 407.92628821 Ry<br class=""> xc contribution = -287.90966771 Ry<br class=""> ewald contribution = -702.94572199 Ry<br class=""> Dispersion Correction = -0.28187072 Ry<br class=""> one-center paw contrib. = -228.09122311 Ry<br class=""><br class=""> convergence has been achieved in 9 iterations<br class=""><br class=""> Forces acting on atoms (Ry/au):<br class=""><br class=""> atom 1 type 1 force = 0.00000104 -0.00894702 -0.00000051<br class=""> atom 2 type 1 force = -0.00000045 -0.00894443 -0.00000014<br class=""> ........<br class=""> atom 33 type 4 force = -0.01684128 -0.06271162 0.16881369<br class=""> atom 34 type 4 force = -0.01684137 -0.06271191 0.16881443<br class=""> atom 35 type 4 force = -0.01684128 0.06271162 0.16881368<br class=""> atom 36 type 4 force = -0.01684137 0.06271191 0.16881444<br class=""> atom 37 type 4 force = 0.01684128 -0.06271162 -0.16881368<br class=""> atom 38 type 4 force = 0.01684137 -0.06271191 -0.16881444<br class=""> atom 39 type 4 force = 0.01684128 0.06271162 -0.16881369<br class=""> atom 40 type 4 force = 0.01684137 0.06271191 -0.16881443<br class=""> atom 41 type 4 force = 0.05926419 0.21504962 -0.01148228<br class=""> atom 42 type 4 force = 0.05926450 0.21505075 -0.01148234<br class=""> atom 43 type 4 force = 0.05926419 -0.21504962 -0.01148228<br class=""> atom 44 type 4 force = 0.05926450 -0.21505072 -0.01148234<br class=""> atom 45 type 4 force = -0.05926419 0.21504962 0.01148228<br class=""> atom 46 type 4 force = -0.05926450 0.21505072 0.01148234<br class=""> atom 47 type 4 force = -0.05926419 -0.21504962 0.01148228<br class=""> atom 48 type 4 force = -0.05926450 -0.21505075 0.01148234<br class=""><br class=""> ...........<br class=""><br class=""> Total force = 2.613470 Total SCF correction = 0.000031<br class=""><br class=""> Total Dispersion Force = 0.007543<br class=""><br class=""> ...........<br class=""><br class=""> BFGS Geometry Optimization<br class=""><br class=""> number of scf cycles = 1<br class=""> number of bfgs steps = 0<br class=""><br class=""> energy new = -1426.3183075352 Ry<br class=""><br class=""> new trust radius = 0.5000000000 bohr<br class=""> new conv_thr = 0.0000001000 Ry<br class=""><br class=""><br class="">**** QE input file starts here ***********************************<br class="">&control<br class=""> title = 'Gypsum (CaSO4.2H20 Z=4)'<br class=""> prefix = 'gypsum-SSSP_acc_PBE-exp0.56G'<br class=""> outdir = '/work2/u2equa02/tempdir/gypsum-SSSP_acc_PBE-exp0.56G'<br class=""> pseudo_dir = '/home/u2equa02/quantumhome/pseudo'<br class=""> calculation = 'relax'<br class=""> restart_mode = 'from_scratch'<br class=""> disk_io = 'low'<br class=""> verbosity = 'high'<br class=""> tstress = .true.<br class=""> tprnfor = .true.<br class=""> nstep = 1<br class=""> etot_conv_thr = 1.0D-5<br class=""> forc_conv_thr = 1.0D-4<br class=""> iprint = 1<br class="">/<br class="">&system<br class=""> ibrav = 0<br class=""> ntyp = 4<br class=""> nat = 48<br class=""> ecutwfc = 100<br class=""> ecutrho = 800<br class=""> vdw_corr = 'grimme-d2'<br class="">/<br class="">&electrons<br class=""> electron_maxstep = 200<br class=""> conv_thr = 1.0D-7<br class=""> diago_thr_init = 1.0D-4<br class=""> startingpot = 'atomic'<br class=""> startingwfc = 'atomic+random'<br class=""> mixing_mode = 'plain'<br class=""> mixing_beta = 0.5<br class=""> mixing_ndim= 8<br class=""> diagonalization = 'david'<br class="">/<br class="">&ions<br class=""> ion_dynamics = 'bfgs'<br class="">/<br class=""><br class="">ATOMIC_SPECIES<br class=""> Ca 40.0780000000 Ca_pbe_v1.uspp.F.UPF<br class=""> S 32.0650000000 S_pbe_v1.2.uspp.F.UPF<br class=""> O 15.9994000000 O.pbe-n-kjpaw_psl.0.1.UPF<br class=""> H 1.0079000000 H.pbe-rrkjus_psl.0.1.UPF<br class=""><br class="">ATOMIC_POSITIONS crystal<br class="">Ca 0.00000000000000 0.17110000000000 0.25000000000000 0 0 0<br class="">Ca 0.50000000000000 0.67110000000000 0.25000000000000 0 0 0<br class="">Ca 0.00000000000000 0.82890000000000 0.75000000000000 0 0 0<br class="">Ca 0.50000000000000 0.32890000000000 0.75000000000000 0 0 0<br class=""> S 0.00000000000000 0.32880000000000 0.75000000000000 0 0 0<br class=""> S 0.50000000000000 0.82880000000000 0.75000000000000 0 0 0<br class=""> S 0.00000000000000 0.67120000000000 0.25000000000000 0 0 0<br class=""> S 0.50000000000000 0.17120000000000 0.25000000000000 0 0 0<br class=""> O 0.08200000000000 0.27600000000000 0.59100000000000 0 0 0<br class=""> O 0.58200000000000 0.77600000000000 0.59100000000000 0 0 0<br class=""> O 0.08200000000000 0.72400000000000 0.09100000000000 0 0 0<br class=""> O 0.58200000000000 0.22400000000000 0.09100000000000 0 0 0<br class=""> O 0.91800000000000 0.27600000000000 0.90900000000000 0 0 0<br class=""> O 0.41800000000000 0.77600000000000 0.90900000000000 0 0 0<br class=""> O 0.91800000000000 0.72400000000000 0.40900000000000 0 0 0<br class=""> O 0.41800000000000 0.22400000000000 0.40900000000000 0 0 0<br class=""> O 0.20300000000000 0.38700000000000 0.91700000000000 0 0 0<br class=""> O 0.70300000000000 0.88700000000000 0.91700000000000 0 0 0<br class=""> O 0.20300000000000 0.61300000000000 0.41700000000000 0 0 0<br class=""> O 0.70300000000000 0.11300000000000 0.41700000000000 0 0 0<br class=""> O 0.79700000000000 0.38700000000000 0.58300000000000 0 0 0<br class=""> O 0.29700000000000 0.88700000000000 0.58300000000000 0 0 0<br class=""> O 0.79700000000000 0.61300000000000 0.08300000000000 0 0 0<br class=""> O 0.29700000000000 0.11300000000000 0.08300000000000 0 0 0<br class=""> O 0.78700000000000 0.06900000000000 0.91800000000000 0 0 0<br class=""> O 0.28700000000000 0.56900000000000 0.91800000000000 0 0 0<br class=""> O 0.78700000000000 0.93100000000000 0.41800000000000 0 0 0<br class=""> O 0.28700000000000 0.43100000000000 0.41800000000000 0 0 0<br class=""> O 0.21300000000000 0.06900000000000 0.58200000000000 0 0 0<br class=""> O 0.71300000000000 0.56900000000000 0.58200000000000 0 0 0<br class=""> O 0.21300000000000 0.93100000000000 0.08200000000000 0 0 0<br class=""> O 0.71300000000000 0.43100000000000 0.08200000000000 0 0 0<br class=""> H 0.74200000000000 0.08700000000000 0.76600000000000 1 1 1<br class=""> H 0.24200000000000 0.58700000000000 0.76600000000000 1 1 1<br class=""> H 0.74200000000000 0.91300000000000 0.26600000000000 1 1 1<br class=""> H 0.24200000000000 0.41300000000000 0.26600000000000 1 1 1<br class=""> H 0.25800000000000 0.08700000000000 0.73400000000000 1 1 1<br class=""> H 0.75800000000000 0.58700000000000 0.73400000000000 1 1 1<br class=""> H 0.25800000000000 0.91300000000000 0.23400000000000 1 1 1<br class=""> H 0.75800000000000 0.41300000000000 0.23400000000000 1 1 1<br class=""> H 0.75600000000000 0.02000000000000 0.92300000000000 1 1 1<br class=""> H 0.25600000000000 0.52000000000000 0.92300000000000 1 1 1<br class=""> H 0.75600000000000 0.98000000000000 0.42300000000000 1 1 1<br class=""> H 0.25600000000000 0.48000000000000 0.42300000000000 1 1 1<br class=""> H 0.24400000000000 0.02000000000000 0.57700000000000 1 1 1<br class=""> H 0.74400000000000 0.52000000000000 0.57700000000000 1 1 1<br class=""> H 0.24400000000000 0.98000000000000 0.07700000000000 1 1 1<br class=""> H 0.74400000000000 0.48000000000000 0.07700000000000 1 1 1<br class=""><br class="">K_POINTS automatic<br class="">12 6 18 0 0 0<br class=""><br class="">CELL_PARAMETERS bohr<br class=""> 11.83535471815098 0.00000000000000 0.00000000000000<br class=""> 0.00000000000000 28.65391722678003 0.00000000000000<br class=""> -4.39283675093634 0.00000000000000 9.75628273508813<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">Tsung-Lung Li, Ph. D.<br class="">Professor<br class="">Department of Applied Physics<br class="">National Chia-Yi University<br class="">300 Hsueh-Fu Road, Chiayi 60004, Taiwan<br class="">Phone: 886-5-2717904. FAX: 886-5-2717909.<br class="">E-mail:<a href="mailto:quantum@mail.ncyu.edu.tw" class="">quantum@mail.ncyu.edu.tw</a> <br class="">URL:<a href="http://web.ncyu.edu.tw/~quantum" class="">http://web.ncyu.edu.tw/~quantum</a><br class=""><br class="">_______________________________________________<br class="">Pw_forum mailing list<br class=""><a href="mailto:Pw_forum@pwscf.org" class="">Pw_forum@pwscf.org</a><br class="">http://pwscf.org/mailman/listinfo/pw_forum<br class=""><br class=""></div></div></blockquote></div><br class=""><div class="">
-- <br class=""><br class="">Giovanni Cantele, PhD<br class="">CNR-SPIN<br class="">c/o Dipartimento di Fisica<br class="">Universita' di Napoli "Federico II"<br class="">Complesso Universitario M. S. Angelo - Ed. 6<br class="">Via Cintia, I-80126, Napoli, Italy<br class="">e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it" class="">giovanni.cantele@spin.cnr.it</a><br class="">Phone: +39 081 676910<br class="">Skype contact: giocan74<br class=""><br class="">ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" class="">http://www.researcherid.com/rid/A-1951-2009</a><br class="">Web page: <a href="http://people.na.infn.it/~cantele" class="">http://people.na.infn.it/~cantele</a><br class="">
</div>
<br class=""></div></body></html>