<div dir="ltr"><div><div>Hi, All:<br><br></div><div>I am very grateful somebody can help me the following problem.<br></div><div><br></div>I am trying to use QE to do Si  crystal lattice optimization. please see the input of my file.  However, after calculations, I sue grep -n "volume" Si.rx.out , I got results that always show the volume of the unit cell doesn't change. <br>[hongtang@login3 results]$ grep -n "volume" Si.rx.out<br>39:     unit-cell volume          =     265.3020 (a.u.)^3<br>1812:     new unit-cell volume =    265.30200 a.u.^3 (    39.31370 Ang^3 )<br>3146:     new unit-cell volume =    265.30200 a.u.^3 (    39.31370 Ang^3 )<br>4462:     new unit-cell volume =    265.30200 a.u.^3 (    39.31370 Ang^3 )<br>5776:     new unit-cell volume =    265.30200 a.u.^3 (    39.31370 Ang^3 )<br>5807:     unit-cell volume          =     265.3020 (a.u.)^3<br>________________________________________________<br><br></div>input  <br>&CONTROL<br>  calculation = "vc-relax",<br>  pseudo_dir  = "/home/hongtang/tmqe6.0/pseudo",<br>  prefix      = "Si",<br>  outdir      = "/home/hongtang/tmqe6.0/tempdir",<br>  dt          = 100,<br>  verbosity   = 'high',<br>  tstress     = .true.<br>  tprnfor     = .true. <br>  nstep       = 900,<br>  etot_conv_thr = 1.D-5,<br>  forc_conv_thr  = 1.D-4,<br>/<br>&SYSTEM<br>  ibrav       = 2,<br>  celldm(1)   = 10.2D0,<br>  nat         = 2,<br>  ntyp        = 1,<br>  ecutwfc     = 60.D0,<br>  nosym = .true.<br>/<br>&ELECTRONS<br>  conv_thr   = 1.D-7,<br>  mixing_beta = 0.7D0,<br>  electron_maxstep  = 300<br>/<br>&IONS<br>   bfgs_ndim      = 1,<br>   trust_radius_max = 1.0D0,<br>   trust_radius_ini  = 0.6D0,<br>   trust_radius_min = 0.1D0,<br>/<br>&CELL<br>  cell_dynamics = 'bfgs',<br>/<br>ATOMIC_SPECIES<br>Si 28.09 Si.tpss-mt.UPF   <br>ATOMIC_POSITIONS<br>Si 0.00 0.00 0.00<br>Si 0.25 0.26 0.24<br>K_POINTS {automatic}<br>7 7 7 0 0 0 <br><br><br><br><br><br></div>