<div dir="ltr">Hi all,<div><br></div><div>We're trying to simulate an n-octane crystal (triclinic, coordinates and cell vectors from literature) -- just a static cell + geometry optimization. Here's the issue: when simulating a single cell (1 1 1), the crystal cohesive energy is far above the experiment -- an order of magnitude larger in absolute value. As a vacuum ref we're simulating an isolated octane molecule. In all cases, a NON_LOCAL vdw energy correction is used.</div><div><br></div><div>When simulating a total of eight cells (2 2 2), and the simulation is still running, so the results are preliminary, things appear to be much more reasonable. The input file is below (I can provide the coordinates as well, if needed). Can you just by simple inspection tell me if anything is obviously wrong?</div><div><br></div><div>Thanks a lot,</div><div><br></div><div>Alex</div><div><br></div><div>**********************************</div><div><div>&GLOBAL</div><div>  PROJECT opt</div><div>  PRINT_LEVEL MEDIUM </div><div>  RUN_TYPE CELL_OPT</div><div>&END GLOBAL</div><div> &MOTION</div><div>  &GEO_OPT</div><div>    OPTIMIZER BFGS </div><div>  &END</div><div>  &CELL_OPT</div><div>    EXTERNAL_PRESSURE [bar] 1.0</div><div>    OPTIMIZER BFGS</div><div>  &END</div><div>&END MOTION</div><div> </div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME BASIS_SET2</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS2</div><div>    &MGRID</div><div>      NGRIDS 5</div><div>      CUTOFF 500</div><div>    &END MGRID</div><div>    &QS</div><div>      METHOD GPW</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS RESTART </div><div>      EPS_SCF 5.0E-6</div><div>      MAX_SCF 300</div><div>      ADDED_MOS  100 </div><div>      &MIXING</div><div>          METHOD BROYDEN_MIXING</div><div>          ALPHA    0.1</div><div>          NBROYDEN   8</div><div>      &END MIXING</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>      &VDW_POTENTIAL NON_LOCAL</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div> <span class="Apple-tab-span" style="white-space:pre">  </span>ABC<span class="Apple-tab-span" style="white-space:pre"> </span>4.160 4.750 11.000</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>ANGLES<span class="Apple-tab-span" style="white-space:pre">      </span>94.80  84.50 105.10</div><div>        PERIODIC XYZ</div><div><span class="Apple-tab-span" style="white-space:pre">  </span>SYMMETRY TRICLINIC</div><div>        MULTIPLE_UNIT_CELL 2 2 2 </div><div>    &END</div><div>    &TOPOLOGY</div><div>       COORD_FILE_NAME noct_2mol.pdb</div><div>       COORD_FILE_FORMAT PDB</div><div>       MULTIPLE_UNIT_CELL  2 2 2 </div><div>    &END</div><div>    &KIND C</div><div>      BASIS_SET  DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL  GTH-PBE-q4</div><div>    &END KIND</div><div>    &KIND H</div><div>      BASIS_SET  DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL  GTH-PBE-q1</div><div>    &END KIND</div><div>&END SUBSYS</div><div>&END</div><div><br></div></div></div>