<br><div>Deal all,</div><div><br></div><div>I'm doing a cell optimisation run of a YAG (Y3Al5O12) crystal. I wonder about the meaning of the tags EXTERNAL_PRESSURE and PRESSURE_TOLERANCE in the CELL_OPT section of the input. The default values are both 100 bar. Here are the few points I don't understand,</div><div><br></div><div>1. If the cell optimisation run is targeted at the input EXTERNAL_PRESSURE, shouldn't the default PRESSURE_TOLERANCE be given a much smaller value, like 5.0 or 1.0 bar?</div><div><br></div><div>2. Why 100 bar? 1 atm = 1.013 bar sounds more reasonable.</div><div><br></div><div>3. My simple understanding of cell optimisation is to find the minimum point of the energy-volume (E-V) curve. Then what's the role of external pressure?</div><div><br></div><div>Here is my input file.</div><div><br></div><div>&GLOBAL<br>  project  pristine-opt<br>  run_type CELL_OPT<br>  print_level MEDIUM<br>&END GLOBAL<br><br>&FORCE_EVAL<br>  METHOD QS<br><br>  &SUBSYS<br>    &CELL <br>      ABC  12.0003 12.0003 12.0003  <br>      PERIODIC XYZ <br>      MULTIPLE_UNIT_CELL  1   1   1<br>      SYMMETRY CUBIC <br>    &END CELL<br><br>    &TOPOLOGY<br>      MULTIPLE_UNIT_CELL  1   1   1<br>      COORD_FILE_FORMAT XYZ<br>      COORD_FILE_NAME yag_from_cif_conven_unit.xyz <br>    &END TOPOLOGY<br><br>    &KIND Y<br>      ELEMENT Y<br>      BASIS_SET DZV-GTH-PADE<br>      POTENTIAL GTH-PBE-q11<br>    &END KIND <br><br><br>    &KIND Al<br>      ELEMENT Al<br>      BASIS_SET DZVP-MOLOPT-SR-GTH <br>      POTENTIAL GTH-PBE-q3<br>    &END KIND <br><br>    &KIND O<br>      ELEMENT O<br>      BASIS_SET DZVP-GTH-PBE<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br><br>    &KIND Ce<br>      ELEMENT Ce<br>      BASIS_SET DZVP-MOLOPT-SR-GTH-q12<br>      POTENTIAL GTH-PBE-q12<br>    &END KIND<br><br>    &KIND Yb<br>      ELEMENT Yb<br>      BASIS_SET DZVP-MOLOPT-SR-GTH-q24<br>      POTENTIAL GTH-PBE-q24<br>    &END KIND<br><br>  &END SUBSYS<br><br>  &DFT<br>    BASIS_SET_FILE_NAME BASIS_SET_YAG<br>    POTENTIAL_FILE_NAME POTENTIAL_YAG<br><br>    &MGRID<br>      NGRIDS 4<br>      CUTOFF 600<br>      REL_CUTOFF 60<br>    &END MGRID<br><br>    &QS<br>      EPS_DEFAULT 1.0E-12<br>    &END QS<br><br>    &POISSON<br>      PERIODIC XYZ<br>    &END POISSON<br><br>    &SCF<br>      SCF_GUESS ATOMIC<br>      EPS_SCF 5.0E-7<br>      MAX_SCF 200 <br><br>      CHOLESKY INVERSE <br><br>      &OT ON<br>        MINIMIZER DIIS<br>      &END OT<br>      #&DIAGONALIZATION <br>      #  ALGORITHM STANDARD<br>      #&END DIAGONALIZATION <br><br>      # ADDED_MOS 30<br>      #<br>      # &SMEAR ON<br>      # METHOD FERMI_DIRAC<br>      # ELECTRONIC_TEMPERATURE [K] 300<br>      #&END SMEAR <br><br>      # &MIXING<br>      #  METHOD BROYDEN_MIXING<br>      #  ALPHA 0.4<br>      #  BETA 0.5<br>      #  NBROYDEN 8<br>      # &END MIXING <br>    &END SCF<br><br>    &XC<br>      &XC_FUNCTIONAL PBE<br>      &END XC_FUNCTIONAL <br>    &END XC<br><br>    &KPOINTS<br>      EPS_GEO 1.0E-5<br>      SCHEME MONKHORST 3 3 3<br>    &END KPOINTS<br>  &END DFT<br><br>  &PRINT<br>    &TOTAL_NUMBERS ON<br>    &END TOTAL_NUMBERS<br>  &END PRINT<br>  <br>  STRESS_TENSOR ANALYTICAL <br>&END FORCE_EVAL<br><br>&MOTION<br>  &GEO_OPT<br>    OPTIMIZER LBFGS<br>    TYPE MINIMIZATION<br>    MAX_DR   5.0E-4<br>    MAX_FORCE 5.0E-4<br>    RMS_DR  5.0E-4<br>    RMS_FORCE 5.0E-4<br>    MAX_ITER 500<br>  &END GEO_OPT<br><br>  &CELL_OPT<br>    TYPE DIRECT_CELL_OPT<br>    OPTIMIZER LBFGS<br>    KEEP_SYMMETRY<br>    MAX_DR 5.0E-4<br>    MAX_FORCE 1.0E-4<br>    RMS_DR 5.0E-4<br>    RMS_FORCE 5.0E-4<br>    # defaults for p and delta p [bar]<br>    EXTERNAL_PRESSURE   100.0 0.0 0.0    0.0 100.0 0.0   0.0 0.0 100.0<br>    PRESSURE_TOLERANCE  100.0<br>  &END CELL_OPT<br>&END MOTION <br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Thank you,</div><div>Huaiyang Sun<br></div><div><br></div>