<div dir="ltr">Dear Thermo_pw users,<div><br></div><div>I am new to therm_pw.I am trying to do lattice parameter calculation (and later elastic constants) for sI clathrate hydrate (46 water molecule+8 guest gases).The structure is a cubic.</div><div><br></div><div>The input file:</div><div><a href="http://scf.in">scf.in</a></div><div>&CONTROL<br>  title = 'SI_8CO2_Relax'<br>  pseudo_dir = '.'<br>  !outdir = './outdir'<br>  calculation = 'relax'<br>  nstep=50<br>  etot_conv_thr = 1e-4<br>  forc_conv_thr = 1e-4<br> max_seconds=82500<br>/<br><br>&SYSTEM<br>  A = 12.029999733000<br>  nat= 162<br>  ntyp= 3<br>  ibrav= 1<br>  ecutwfc= 55.0<br>  ecutrho = 440.0<br>!ecuthro=8-10*ecutwfc<br>!occupations='smearing'<br>!gaussian smearing<br>  smearing='gaussian'<br>!ordinary Gaussian spreading (Default)<br>  degauss= 0.022<br>  input_dft= 'vdw-df2'<br>/<br><br>&ELECTRONS<br>  mixing_beta = 0.7<br>  conv_thr =  1.0d-08<br>  diagonalization='david'<br>  electron_maxstep = 150<br>/<br><br>&IONS<br>/<br><br>&CELL<br>/<br><br>ATOMIC_SPECIES<br>C   12.011  C.pbe-n-kjpaw_psl.1.0.0.UPF<br>H    1.008  H.pbe-kjpaw_psl.1.0.0.UPF<br>O   15.999  O.pbe-n-kjpaw_psl.0.1.UPF<br><br>ATOMIC_POSITIONS (angstrom)</div><div>C             5.6914623648        6.3477925243        9.0385162417<br>C             6.0203188824        0.0140445035        6.0076929324<br>C             0.0103067680        6.0124985177        0.0271263591<br>C             0.5255489420        9.0628026901        5.6877293485<br>C            -0.5089082791        2.9502728029        6.3664117395<br>C             2.9707134143        0.1359390055        0.0057036813<br>C             9.0362791318        0.2927456487       -0.2246132818<br>C             6.3078182653        5.6844139347        2.9923046100<br>H             7.6681600210        8.8103809640        5.1812970528<br>H             7.9474183614        7.5109661777        6.0550379444<br>H             8.0588580146        9.3162094407        7.4168881942<br>H             8.0234619914        2.8265464617        7.5198466803<br>H            -0.9928948776        2.8873069490        1.7034973611<br>H             8.4213026994        5.2152049586        6.1694929950<br>H             8.5729430479        1.3775298181        4.1473499668<br></div><div>...</div><div>O            -0.0990445262        0.0367157095        2.9958286391</div>O             5.2088485235       -0.1956393873        5.1874258136<br>K_POINTS automatic<br>2 2 2  0 0 0<div><br><div>thermo_control<br></div><div>&INPUT_THERMO<br>  what='mur_lc',<br>  lmurn=.FALSE.,<br>  frozen_ions=.FALSE.<br> /<br>thermo_control (END)<br></div><div><br></div><div><br></div><div>I encountered the following problems:</div><div><br></div><div>1-The code did not recognize that the structure is cubic and assigin it to triclinic Bravais lattice; ibrav=   14</div><div><br></div><div>-----------</div><div>  Program THERMO_PW v.6.5 starts on 18Jul2021 at  9:31:41 <br><br>     This program is part of the open-source Quantum ESPRESSO suite<br>     for quantum simulation of materials; please cite<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);<br>          URL <a href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>", <br>     in publications or presentations arising from this work. More details at<br>     <a href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</a><br><br>     Parallel version (MPI), running on   168 processors<br><br>     MPI processes distributed on     6 nodes<br>     K-points division:     npool     =       2<br>     R & G space division:  proc/nbgrp/npool/nimage =      84<br>     Reading input from _temporary_1<br>Warning: card &CELL ignored<br>     IMPORTANT: XC functional enforced from input :<br>     Exchange-correlation= VDW-DF2<br>                           (   1   4  13   0   2   0   0)<br>     Any further DFT definition will be discarded<br>     Please, verify this is what you really want<br><br>               file C.pbe-n-kjpaw_psl.1.0.0.UPF: wavefunction(s)  2S 2P renormalized<br>               file H.pbe-kjpaw_psl.1.0.0.UPF: wavefunction(s)  1S renormalized<br>               file O.pbe-n-kjpaw_psl.0.1.UPF: wavefunction(s)  2P renormalized<br><br>   Info: using nr1, nr2, nr3 values from input<br><br>     C_1 (1) is incompatible with the simple cubic Bravais lattice<br>     It is compatible with the <br>     triclinic Bravais lattice; ibrav=   14<br><br>     You might want to change the Bravais lattice or to<br><br>     understand why the symmetries are wrong before continuing<br>     The point group or the Laue class are not used to reduce the number of <br>     computed tensor components<br><br><br>     Calculating the volume that minimizes the energy<br><br>     FFT mesh: (  160,  160,  160 )<br><br>     Bravais lattice:<br><br>     ibrav=  1: simple cubic<br>     Cell parameters:<br></div><div> alat=  22.733405 a.u.<br><br><br>     Starting primitive lattice vectors:<br>     crystal axes: (cart. coord. in units of alat)<br><br>               a(1) = (   1.000000   0.000000   0.000000 )  <br>               a(2) = (   0.000000   1.000000   0.000000 )  <br>               a(3) = (   0.000000   0.000000   1.000000 )  <br><br>     Starting reciprocal lattice vectors:<br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br><br>               b(1) = (  1.000000  0.000000  0.000000 )  <br>               b(2) = (  0.000000  1.000000  0.000000 )  <br>               b(3) = (  0.000000  0.000000  1.000000 )  <br></div><div><br></div><div>2-The maximum ionic relaxation steps did not exceed the default 20 steps although I have specified 50 in the input file, so some geometries were not relaxed.</div><div><br></div><div><br></div><div>Your help is appreciated</div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:10pt;background-color:rgb(253,253,253)">Regards<span style="font-size:10pt">,</span></div><div style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:10pt;background-color:rgb(253,253,253)"><div><span style="color:rgb(102,102,153);font-family:calibri;font-size:14pt;font-weight:900">Ahmed OMRAN</span><br><span style="color:rgb(124,148,157);font-family:calibri;font-size:16pt;font-weight:900"></span></div><div><strong style="white-space:pre-wrap;font-family:calibri;font-size:11pt">Researcher</strong></div><div><p class="MsoNormal" style="line-height:16.8667px;font-size:11pt;font-family:calibri,sans-serif"><span lang="FR" style="font-size:9.5pt;line-height:14.5667px;font-family:verdana,sans-serif"></span></p><span style="font-family:calibri;font-size:11pt"><strong>Laboratoire Catalyse et Spectrochemie (LCS)</strong></span></div><div><span style="font-family:calibri;font-size:11pt">France </span></div></div><div style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:10pt;background-color:rgb(253,253,253)"><table width="783" cellspacing="5" border="0" style="height:149px;font-size:13.3333px;width:783px"><tbody><tr style="height:149px"><td width="594" valign="middle" style="font-family:verdana,arial,helvetica,sans-serif;padding-left:5px;height:149px;width:645px"><br></td></tr></tbody></table></div></div></div></div></div></div></div></div>