<div dir="ltr">Hi, Dear CP2K users,<br><div>    </div><div>      When I simulate  proton
permeation of graphene,  I can't keep H-positive 1 charge in the process.  The only keep the system 1charge,but the charge of H  keep changing.</div><div>  </div><div><br></div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT gre-O</div><div>  RUN_TYPE BAND </div><div>  PRINT_LEVEL low</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD QS              ! Electronic structure method (DFT,...)</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME  BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME  GTH_POTENTIALS</div><div><br></div><div>    CHARGE 1</div><div>    MULTIPLICITY  1</div><div><br></div><div>    &POISSON                    ! Solver requested for non periodic calculations</div><div>      PERIODIC XYZ</div><div>    &END POISSON</div><div><br></div><div>   &PRINT </div><div>      &HIRSHFELD</div><div>        SELF_CONSISTENT T</div><div>        SHAPE_FUNCTION DENSITY</div><div>      &END HIRSHFELD</div><div>    &END PRINT</div><div><br></div><div>    &SCF                        ! Parameters controlling the convergence of the scf. This section should not be changed. </div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 9.9999999999999998E-07</div><div>      MAX_SCF 30</div><div>      &OT ON                  </div><div>          MINIMIZER DIIS            </div><div>      &END OT</div><div>      &OUTER_SCF</div><div><span style="white-space:pre">     </span>EPS_SCF 9.9999999999999998E-07          </div><div><span style="white-space:pre">        </span>MAX_SCF 50         </div><div>      &END OUTER_SCF</div><div>    &END SCF</div><div><br></div><div>    &QS</div><div>      EPS_DEFAULT 9.9999999999999998E-013</div><div>      WF_INTERPOLATION PS</div><div>      EXTRAPOLATION_ORDER 3</div><div>      METHOD GPW</div><div>    &END QS</div><div><br></div><div>    &MGRID</div><div>      NGRIDS 5</div><div>      CUTOFF 400</div><div>      REL_CUTOFF 60</div><div>    &END MGRID   </div><div>  </div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE </div><div>      &END XC_FUNCTIONAL</div><div>      &vdW_POTENTIAL</div><div>         DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>            TYPE DFTD3</div><div>            REFERENCE_FUNCTIONAL PBE</div><div>            PARAMETER_FILE_NAME dftd3.dat         </div><div>       &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div><br></div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 15 15 15</div><div>      PERIODIC XYZ</div><div>    &END CELL</div><div><br></div><div>    &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules</div><div>      &CENTER_COORDINATES</div><div>      &END</div><div>      COORD_FILE_FORMAT xyz</div><div>      COORD_FILE_NAME  is.xyz</div><div>    &END</div><div><br></div><div>    &KIND H<br>      ELEMENT H</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END KIND</div><div>    &KIND C</div><div>      ELEMENT C</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>  &BAND</div><div>    NPROC_REP  12</div><div>    BAND_TYPE CI-NEB</div><div>    NUMBER_OF_REPLICA 18 </div><div>    K_SPRING 0.02</div><div>    &CONVERGENCE_CONTROL</div><div>      MAX_FORCE 0.0010</div><div>      RMS_FORCE 0.0050</div><div>      MAX_DR 0.002</div><div>      RMS_DR 0.005</div><div>    &END</div><div>    ROTATE_FRAMES TRUE</div><div>    ALIGN_FRAMES TRUE</div><div>    &CI_NEB</div><div>      NSTEPS_IT 2</div><div>    &END</div><div>    &OPTIMIZE_BAND</div><div>      OPT_TYPE DIIS</div><div>      OPTIMIZE_END_POINTS F</div><div>      &DIIS</div><div>        MAX_STEPS 1000</div><div>        N_DIIS 3</div><div>      &END</div><div>    &END</div><div>    &PROGRAM_RUN_INFO</div><div>    &END</div><div>    &CONVERGENCE_INFO</div><div>    &END</div><div><br></div><div>    &REPLICA</div><div>      COORD_FILE_NAME is.xyz</div><div>    &END</div><div>    &REPLICA</div><div>      COORD_FILE_NAME fs.xyz</div><div>    &END</div><div>  &END BAND</div><div>&END MOTION</div></div><div><br></div><div>  </div><div>Best,</div><div> Yan Huang</div></div>