<div dir="ltr"><div><div><div><div>Dear Reinaldo, are you sure that these values of ecutwfc = 40 and ecutrho = 200 describe your system well? Probabily you should increase these values, in order to describe the transition metal atom (Cu) and the total energy (forces and pressure depends of the total energy). Try to converge the total energy of the system respect to ecutwfc value performing self-consistent calculation for different vaues of ecutwfc. Furthermore, for US pseudopotential the ecutrho value is tipically 8 to 12 times ecutwfc (see QE manual).<br><br></div>Best<br><br></div>Dr. Arles V. Gil Rebaza<br></div>Instituto de Física La Plata<br></div>La Plata - Argentina. <br></div><div class="gmail_extra"><br><div class="gmail_quote">2018-06-01 11:20 GMT-03:00 Reinaldo Pis Diez <span dir="ltr"><<a href="mailto:reinaldo.pisdiez@gmail.com" target="_blank">reinaldo.pisdiez@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear folks,<br>
<br>
I usually help some experimentalist fellows from the Center I belong by calculating geometries, cell parameters, vibrational frequencies, electronic transitions, etc, of transition metal complexes. As the cif file is often available, I use it to create the input file for PWScf using Virtual Nano Lab from Quantum Wise.<br>
<br>
Although convergence of the vc-relax job is usually achieved, some cases are more problematic. Job ends with a high pressure and with atomic forces far from thresholds when 'bfgs' is used both for &ions and &cell namelists. Using 'damp' does not help to get geometry and cell parameters optimized. Thus my (very general) question is: does anybody know an strategy to systematically get convergence when doing vc-relax? A link to a tutorial or a document is welcome.<br>
<br>
For completeness, the input of one of the problematic cases is appended below. By the way, I'm using "Program PWSCF v.6.1 (gpu tag v1.0)" compiled with PGI compilers and MKL libraries on a CentOS box with a Titan X (Pascal) gpu.<br>
<br>
Thanks in advance<br>
<br>
Reinaldo<br>
Center of Inorganic Chemistry<br>
Natl Univ of La Plata, Argentina<br>
<br>
Input follows....<br>
<br>
&CONTROL<br>
title = 'Cu-ovatpnh2'<br>
calculation = 'vc-relax'<br>
restart_mode = 'from_scratch'<br>
outdir = './'<br>
pseudo_dir = './'<br>
prefix = 'Cu-ovatpnh2'<br>
etot_conv_thr = 0.0001<br>
forc_conv_thr = 0.0004<br>
tstress = .true.<br>
tprnfor = .true.<br>
nstep = 200<br>
/<br>
<br>
&SYSTEM<br>
ibrav = 14<br>
celldm(1) = 9.9320792573d0<br>
celldm(2) = 2.0161515879d0<br>
celldm(3) = 2.1156030541d0<br>
celldm(4) = -0.1705259449d0<br>
celldm(5) = -0.0219091103d0<br>
celldm(6) = -0.2238523262d0<br>
nat = 59<br>
ntyp = 6<br>
ecutwfc = 40<br>
ecutrho = 200<br>
input_dft = 'pbe'<br>
occupations = 'smearing'<br>
smearing = 'mv'<br>
degauss = 0.005<br>
nspin = 2<br>
starting_magnetization(1) = 0.7<br>
/<br>
<br>
&ELECTRONS<br>
conv_thr = 1d-06<br>
mixing_beta = 0.6<br>
mixing_mode = 'plain'<br>
mixing_ndim = 8<br>
diagonalization = 'david'<br>
/<br>
<br>
&IONS<br>
ion_dynamics = 'bfgs'<br>
/<br>
<br>
&CELL<br>
cell_dynamics = 'bfgs'<br>
press = 0.5<br>
/<br>
<br>
ATOMIC_SPECIES<br>
Cu 63.546000d0 cu_pbe_v1.2.uspp.F.UPF<br>
C 12.010700d0 c_pbe_v1.2.uspp.F.UPF<br>
H 1.007940d0 h_pbe_v1.4.uspp.F.UPF<br>
N 14.006700d0 n_pbe_v1.2.uspp.F.UPF<br>
O 15.999400d0 o_pbe_v1.2.uspp.F.UPF<br>
S 32.065000d0 s_pbe_v1.4.uspp.F.UPF<br>
<br>
ATOMIC_POSITIONS (crystal)<br>
C 0.082231085 0.204570350 0.616545553<br>
C 0.917768915 0.795429650 0.383454447<br>
C 0.172973917 0.203448044 0.733466608<br>
C 0.827026083 0.796551956 0.266533392<br>
H 0.367673209 0.258122966 0.770942028<br>
H 0.632326791 0.741877034 0.229057972<br>
C 0.987334077 0.125087827 0.797766960<br>
C 0.012665923 0.874912173 0.202233040<br>
H 1.019591574 0.112838755 0.892010655<br>
H -0.019591574 0.887161245 0.107989345<br>
C 0.756454381 0.067986708 0.729401520<br>
C 0.243545619 0.932013292 0.270598480<br>
H 0.580696185 0.004985342 0.755622580<br>
H 0.419303815 0.995014658 0.244377420<br>
C 0.224564542 0.266351629 0.518206870<br>
C 0.775435458 0.733648371 0.481793130<br>
H 0.388829444 0.346391746 0.560989945<br>
H 0.611170556 0.653608254 0.439010055<br>
H 0.304022379 0.192372172 0.457292153<br>
H 0.695977621 0.807627828 0.542707847<br>
C 0.956901991 0.245893003 0.340167231<br>
C 0.043098009 0.754106997 0.659832769<br>
H 0.022416710 0.153595990 0.313033710<br>
H 0.977583290 0.846404010 0.686966290<br>
C 0.771131988 0.269306261 0.255617666<br>
C 0.228868012 0.730693739 0.744382334<br>
C 0.652850645 0.378765293 0.282708985<br>
C 0.347149355 0.621234707 0.717291015<br>
C 0.462994923 0.392175655 0.192940332<br>
C 0.537005077 0.607824345 0.807059668<br>
C 0.403148927 0.302768506 0.084001267<br>
C 0.596851073 0.697231494 0.915998733<br>
H 0.259814627 0.314740185 0.016757517<br>
H 0.740185373 0.685259815 0.983242483<br>
C 0.525750373 0.196269824 0.058490144<br>
C 0.474249627 0.803730176 0.941509856<br>
H 0.478830729 0.128376775 0.971151210<br>
H 0.521169271 0.871623225 0.028848790<br>
C 0.705969389 0.179966480 0.142934709<br>
C 0.294030611 0.820033520 0.857065291<br>
H 0.804698818 0.098652977 0.124221693<br>
H 0.195301182 0.901347023 0.875778307<br>
C 0.184454672 0.523754793 0.135830920<br>
C 0.815545328 0.476245207 0.864169080<br>
H 0.287576875 0.544363666 0.053314825<br>
H 0.712423125 0.455636334 0.946685175<br>
H 0.118308380 0.611071476 0.178999457<br>
H 0.881691620 0.388928524 0.821000543<br>
H 0.013394359 0.440437861 0.110093147<br>
H 0.986605641 0.559562139 0.889906853<br>
N 0.053284394 0.319764373 0.444382321<br>
N 0.946715606 0.680235627 0.555617679<br>
O 0.705759527 0.466685324 0.381694456<br>
O 0.294240473 0.533314676 0.618305544<br>
O 0.356827705 0.498376849 0.225196110<br>
O 0.643172295 0.501623151 0.774803890<br>
S 0.766003561 0.110202377 0.587178900<br>
S 0.233996439 0.889797623 0.412821100<br>
Cu -0.000000000 0.500000000 0.500000000<br>
<br>
K_POINTS {automatic}<br>
3 3 3 0 0 0<br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso<wbr>.org/mailman/listinfo/users</a></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Arles V. <br></div></div></div>
</div>