<p>Hello everyone,</p><p>I am studying the reduction mechanism of my anion on a Na surface using CP2K. However, I am encountering an issue where no reduction occurs on the Na surface, while the same setup shows anion reduction on a Li surface. The CP2K simulation does not produce any errors, making it challenging to diagnose the problem.</p><p>My workflow involves:</p><ol><li>Minimization using LAMMPS</li><li>Geometry optimization</li><li>AIMD (Ab Initio Molecular Dynamics)</li></ol><p>I have attached my AIMD script for reference. Could someone please help me identify what might be causing the lack of anion reduction on the Na surface? Any insights or suggestions would be greatly appreciated.</p><p>Thank you in advance for your assistance!<br />AIMD input file:<br />&GLOBAL<br />  PROJECT 100_IL<br />  RUN_TYPE MD<br />  IOLEVEL LOW<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  METHOD QS<br />  &DFT<br />    BASIS_SET_FILE_NAME ./BASIS_MOLOPT<br />    POTENTIAL_FILE_NAME ./GTH_POTENTIALS<br />    CHARGE 0<br />    LSD<br />    &QS<br />      EPS_DEFAULT 1.0E-10<br />    &END QS<br />    &MGRID<br />      CUTOFF 600<br />      NGRIDS 4<br />      REL_CUTOFF 60<br />    &END MGRID<br />    &SCF                              <br />      SCF_GUESS ATOMIC<br />      MAX_SCF 200<br />      EPS_SCF 1.0E-4<br />      &OT<br />        PRECONDITIONER FULL_SINGLE_INVERSE<br />        MINIMIZER DIIS<br />      &END OT<br />      &OUTER_SCF<br />        MAX_SCF 200<br />        EPS_SCF 1.0E-4<br />      &END OUTER_SCF<br />      &PRINT<br />        &RESTART OFF<br />        &END<br />      &END PRINT<br />    &END SCF<br />    &XC<br />      &XC_FUNCTIONAL PBE<br />      &END XC_FUNCTIONAL<br />      &VDW_POTENTIAL<br />        POTENTIAL_TYPE PAIR_POTENTIAL <br />        &PAIR_POTENTIAL<br />          PARAMETER_FILE_NAME dftd3.dat<br />          TYPE DFTD3<br />          REFERENCE_FUNCTIONAL PBE<br />          R_CUTOFF [angstrom] 16<br />        &END PAIR_POTENTIAL<br />      &END VDW_POTENTIAL<br />    &END XC<br />  &END DFT<br /><br />  &SUBSYS<br />    &CELL<br />      ABC [angstrom] 12.5000 12.5000 30.000<br />    &END CELL<br />    &TOPOLOGY<br />      COORD_FILE_NAME EM_TFSI-pos-1.xyz<br />      COORD_FILE_FORMAT XYZ<br />    &END TOPOLOGY<br /><br />    &KIND H<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q1<br />    &END KIND<br />    &KIND F<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q7<br />    &END KIND<br />    &KIND O<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q6<br />    &END KIND<br />    &KIND C<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q4<br />    &END KIND<br />    &KIND S<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q6<br />    &END KIND<br />    &KIND N<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q5<br />    &END KIND<br />    &KIND Si<br />      BASIS_SET TZVP-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q4<br />    &END KIND<br />    &KIND Na<br />      BASIS_SET DZVP-MOLOPT-SR-GTH<br />      POTENTIAL GTH-PBE-q9<br />    &END KIND<br />  &END SUBSYS<br />&END FORCE_EVAL<br /><br />&MOTION<br />  &GEO_OPT<br />    OPTIMIZER LBFGS<br />    MAX_ITER 100<br />    MAX_DR [bohr] 0.003<br />    &BFGS<br />    &END BFGS<br />  &END GEO_OPT<br /><br />  &MD<br />    ENSEMBLE NVT<br />    TEMPERATURE [K] 350<br />    TIMESTEP [fs] 1<br />    STEPS 30000<br />    &THERMOSTAT<br />      TYPE NOSE<br />      &NOSE<br />        TIMECON 50<br />      &END NOSE<br />    &END THERMOSTAT<br />  &END MD<br /><br />  &PRINT<br />    &TRAJECTORY<br />      &EACH<br />        MD 1<br />      &END EACH<br />    &END TRAJECTORY<br />    &VELOCITIES OFF<br />    &END VELOCITIES<br />    &FORCES OFF<br />    &END FORCES<br />    &RESTART_HISTORY<br />      &EACH<br />        MD 5000<br />      &END EACH<br />    &END RESTART_HISTORY<br />    &RESTART<br />      BACKUP_COPIES 2<br />      &EACH<br />        MD 1<br />      &END EACH<br />    &END RESTART<br />  &END PRINT<br /><br />  &CONSTRAINT<br />    &FIXED_ATOMS<br />      COMPONENTS_TO_FIX XYZ<br />      LIST 1 19 37 7 25 43 13 31 49 2 20 38 8 26 44 14 32 50 3 21 39 9 27 45 15 33 51<br />    &END FIXED_ATOMS<br />  &END CONSTRAINT<br />&END MOTION<br /></p><p><br /></p><p><br /></p><p><br /></p><p>Regards,</p><p>Sreehari</p>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e8e6b622-d45d-44be-b241-e6525bfb86dbn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/e8e6b622-d45d-44be-b241-e6525bfb86dbn%40googlegroups.com</a>.<br />