[CP2K-user] [CP2K:21098] Job failed for large structure but ran perfectly for smaller structure

Abrar Fahim Navid fahimnavid1999 at gmail.com
Thu Jan 30 18:07:54 UTC 2025


Dear CP2K community,

I am currently running an md calculation using nose hoover thermostat. The 
calculation runs perfectly for a supercell (96 atoms). But when I use a 
slab with 348 atoms using the same setup, it fails. I used 3072 cpu's here. 
Could you please enlighten me if I am doing anything wrong or I should 
increase the node number more? 

Best,
Abrar Navid

*My input file:*

&FORCE_EVAL

  METHOD SIRIUS             ! Using SIRIUS method

  STRESS_TENSOR ANALYTICAL

  &PW_DFT

    &CONTROL

        PROCESSING_UNIT cpu

         STD_EVP_SOLVER_NAME LAPACK

         GEN_EVP_SOLVER_NAME LAPACK

    &END CONTROL

 


   &PARAMETERS

      ELECTRONIC_STRUCTURE_METHOD pseudopotential

      SMEARING_WIDTH 0.01

      USE_SYMMETRY false

      NUM_MAG_DIMS 0

!      GK_CUTOFF 5                      !ecutwfc

 !     PW_CUTOFF 20

      GK_CUTOFF 7 !ecutwfc

      PW_CUTOFF 15 !ecutrho

      ENERGY_TOL 1e-10

!      POTENTIAL_TOL 1e-8

!      NUM_DFT_ITER 100 !??

      NGRIDK 3 3 3 !kpoint

      SHIFTK 0 0 0   

    &END PARAMETERS

    &ITERATIVE_SOLVER

       ENERGY_TOLERANCE 1e-5

!       NUM_STEPS 20 !??

       SUBSPACE_SIZE 4

       TYPE DAVIDSON

       CONVERGE_BY_ENERGY 1

    &END ITERATIVE_SOLVER

       &MIXER

       BETA 0.3

       TYPE BROYDEN2

       MAX_HISTORY 8

    &END MIXER


  &END PW_DFT


  &DFT

    BASIS_SET_FILE_NAME BASIS_MOLOPT  ! DZVP-MOLOPT basis sets

    POTENTIAL_FILE_NAME GTH_POTENTIALS

    CHARGE 0

    &SCF

EPS_SCF 1.0E-4 ! etot_conv_thr = 1.0e-4

MAX_SCF 100   !n=400

    &END SCF

    &MGRID

      CUTOFF 300 !Energy cutoff 300

    &END MGRID

    &QS

     METHOD GPW           ! GPW method combining Gaussian basis sets and 
plane waves

    &END QS

    &XC

&XC_FUNCTIONAL

&GGA_X_PBE

&END GGA_X_PBE

&GGA_C_PBE

                &END GGA_C_PBE

! &LIBXC

!               FUNCTIONAL XC_LDA_X

!        &END LIBXC

!        &LIBXC

!               FUNCTIONAL XC_LDA_C_PZ

!        &END LIBXC

! 

        &END XC_FUNCTIONAL

    &END XC

  &END DFT


  &SUBSYS

    &CELL

    A 12.278427332000  0.000000000000  0.000000000000

    B 0.000000000000  12.278427332000   0.000000000000

    C 0.000000000000  0.000000000000  63.345622170000

    &END CELL


    &COORD

      SCALED

[coordinate of atoms;]


     &END COORD


    &KIND C

      POTENTIAL upf C.json

    &END KIND


    &KIND N

      POTENTIAL upf N.json

    &END KIND

    &KIND H

      POTENTIAL upf H.json

    &END KIND


    &KIND I

      POTENTIAL upf I.json 

    &END KIND

    &KIND Pb

      POTENTIAL upf Pb.json

    &END KIND

    &KIND Br

      POTENTIAL upf Br.json

    &END KIND    


  &END SUBSYS

&END FORCE_EVAL


&MOTION

  &GEO_OPT

      OPTIMIZER BFGS

      MAX_FORCE 1.0E-3 !forc_conv_thr = 1.0e-3

 !     MAX_ITER 400

  &END GEO_OPT

  &MD

    ENSEMBLE NPT_F           ! Molecular dynamics with NPT_F ensemble

    TIMESTEP 0.1             ! Timestep in fs should be 0.5

    STEPS 100

    TEMPERATURE 300          ! Initial temperature in Kelvin

    &THERMOSTAT

      TYPE NOSE

      &NOSE

        TIMECON 50           ! Thermostat time constant in fs

        LENGTH 3             ! Three chains for Nose-Hoover thermostat

      &END NOSE

    &END THERMOSTAT

    &BAROSTAT

!      TYPE MTK               ! Martyna-Tuckerman-Klein barostat

      PRESSURE 1.0           ! Target pressure in atm

      TIMECON 50             ! Barostat time constant in fs

    &END BAROSTAT

  &END MD

&END MOTION

&GLOBAL

  PROJECT MAPI

  RUN_TYPE MD               ! Molecular dynamics run

  PRINT_LEVEL MEDIUM

&END GLOBAL

-- 
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 visit https://groups.google.com/d/msgid/cp2k/5a25b394-a726-47e0-a91d-1fc51ca946f5n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250130/aa33f0bc/attachment-0001.htm>


More information about the CP2K-user mailing list