[CP2K-user] [CP2K:21159] Re: Job failed for large structure but ran perfectly for smaller structure
Abrar Fahim Navid
fahimnavid1999 at gmail.com
Sun Feb 16 02:25:03 UTC 2025
Hi,
Is there anyone who could help me in this issue? I am really in trouble and
can't get rid of the problem yet.
On Thursday, 30 January 2025 at 13:21:14 UTC-6 Abrar Fahim Navid wrote:
> 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/40862d7b-a39b-4d19-95de-7d2d35b699dan%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250215/4db41a59/attachment-0001.htm>
More information about the CP2K-user
mailing list