[CP2K-user] [CP2K:19852] cp2k input for semiemperical

Saurabh Negi saurabhnegi95628 at gmail.com
Tue Jan 30 12:38:48 UTC 2024


My metadynamics is not running properly 
Please check my input file

&GLOBAL
  PROJECT 1a
  RUN_TYPE MD
  PRINT_LEVEL MEDIUM
&END GLOBAL

#&EXT_RESTART
#RESTART_FILE_NAME 1a-1.restart
#&END EXT_RESTART

&MOTION
  &MD
    ENSEMBLE NVT # ENSEMBLE NVE
    STEPS 200000
    TIMESTEP 1.0
    TEMPERATURE 295.0
    &THERMOSTAT
      &NOSE
        LENGTH 3 # 3 is the default
        MTS 2 # 2 is the default
        TIMECON 300.0
        YOSHIDA 3 # 3 is the default
      &END NOSE
    &END THERMOSTAT
  &END MD
 
  &FREE_ENERGY
    &METADYN
      DO_HILLS
      NT_HILLS 100
      WW 3.0e-3


      &METAVAR
        SCALE 0.3
        COLVAR 1

        &WALL
          POSITION [angstrom] 4.0
          TYPE QUADRATIC
          &QUADRATIC
            DIRECTION WALL_PLUS
            K [kcalmol] 150
          &END QUADRATIC
&END WALL
      &END METAVAR

      &PRINT
        &COLVAR
           COMMON_ITERATION_LEVELS 3
           &EACH
             MD 1
           &END
        &END
      &END
    &END METADYN
  &END


  &PRINT
    &RESTART_HISTORY
      &EACH
        MD 100
      &END
    &END RESTART_HISTORY
    &VELOCITIES
      &EACH
        MD 1
      &END
    &END
  # Normal restart file
  &RESTART
    &EACH
      MD 1
    &END
  &END RESTART
  &END PRINT
&END MOTION

&FORCE_EVAL
  METHOD Quickstep              ! Electronic structure method (DFT,...)
  &DFT
    CHARGE = 0
    MULTIPLICITY = 2
    LSD
    BASIS_SET_FILE_NAME  BASIS_SET.txt
    POTENTIAL_FILE_NAME  POTENTIAL.txt
    # if WFN restart is required
    WFN_RESTART_FILE_NAME 1a-RESTART.wfn
    &MGRID
     #CUTOFF 1000
      CUTOFF 300
      NGRIDS 4
      REL_CUTOFF 60
    &END MGRID
    &QS
      METHOD PM6
      EPS_DEFAULT 1.0E-10 # 1.0E-10 is the default
    &END QS
   &SCF
     SCF_GUESS RESTART
     EPS_SCF 1.0E-5 # 1.0E-5 is the default
     MAX_SCF 600 # 50 is the default
     &OT
     MINIMIZER CG
     PRECONDITIONER FULL_ALL
     LINESEARCH 3PNT
     &END
     &OUTER_SCF
       EPS_SCF 1.0E-5 # 1.0E-5 is the default
       MAX_SCF 20 # Keep > 100 when OT is on
     &END
   &END SCF  
   &XC
     &XC_FUNCTIONAL BLYP
     &END XC_FUNCTIONAL
     &vdW_POTENTIAL
       DISPERSION_FUNCTIONAL PAIR_POTENTIAL # POTENTIAL_TYPE is alias to 
DISPERSION_FUNCTIONAL
       &PAIR_POTENTIAL
         TYPE DFTD3
         PARAMETER_FILE_NAME dftd3.dat
         REFERENCE_FUNCTIONAL PBE
       &END PAIR_POTENTIAL
    &END vdW_POTENTIAL
   &END XC
   &POISSON                    ! Solver requested for non periodic 
calculations
      PERIODIC XYZ          ! Type of solver
   &END POISSON
  &END DFT  
  &SUBSYS
    &CELL
  A     20.0000000     0.00000000     0.00000000
  B     0.00000000     20.0000000     0.00000000
  C     0.00000000     0.00000000     20.00000000

      PERIODIC XYZ              ! Non periodic calculations. That's why the 
POISSON section is needed 
    &END CELL
    &TOPOLOGY                    ! Section used to center the atomic 
coordinates in the given box. Useful for big molecules
      COORD_FILE_FORMAT xyz
      COORD_FILE_NAME  1a.xyz
    &END
    &KIND H
      ELEMENT H
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q1
    &END KIND    
    &KIND O
      ELEMENT O
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q6
    &END KIND    
    &KIND C
      ELEMENT C
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q4
    &END KIND   
      &KIND Cl
      ELEMENT Cl
      BASIS_SET DZVP-GTH-BLYP
      POTENTIAL GTH-BLYP-q7
    &END KIND  

     &COLVAR
      &DISTANCE
        ATOMS 2 5
      &END DISTANCE
     &END COLVAR
  &END SUBSYS
&END FORCE_EVAL


-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/d6bda0e3-9068-4e04-9802-2bd1ce3e5aa5n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240130/b33de460/attachment.htm>


More information about the CP2K-user mailing list