[CP2K-user] [CP2K:19164] regarding cp2k error CP2K requires molecules to be contiguous and we have detected a non contiguous one
Anmol
anmoljindal0 at gmail.com
Tue Aug 1 08:04:17 UTC 2023
Dear Cp2k community
i was running an NVT simulation for water after 3000 steps it is showing me
the error
CP2K requires molecules to be contiguous and we have detected a non *
* [ABORT] contiguous one!! In particular a molecule defined from index
( *
* \___/ 4 ) to ( 41 ) contains
other *
* | molecules, not connected! Too late at this stage everything
should *
* O/| be already ordered! If you have not yet employed the
REORDERING *
* /| | keyword, please do so.It may help to fix this issue.
*
* / \
topology_generate_util.F:326 *
how to resolve this
below i have attached input file
##############
!Specify Run Type and Print level
&GLOBAL
PROJECT water
RUN_TYPE MD
PRINT_LEVEL MEDIUM
&END GLOBAL
!Force avail section
&FORCE_EVAL
METHOD Quickstep
!STRESS_TENSOR ANALYTICAL
&PRINT
&FORCES ON
&END FORCES
&END PRINT
&DFT
BASIS_SET_FILE_NAME GTH_BASIS_SET
POTENTIAL_FILE_NAME POTENTIAL
CHARGE 0
MULTIPLICITY 1
&XC
&XC_FUNCTIONAL
&MGGA_XC_B97M_V
&END MGGA_XC_B97M_V
&END XC_FUNCTIONAL
&VDW_POTENTIAL
&NON_LOCAL
TYPE RVV10
&END NON_LOCAL
&END VDW_POTENTIAL
&END XC
&QS
EPS_DEFAULT 1.0E-07
&END QS
&MGRID
CUTOFF 400
NGRIDS 4
REL_CUTOFF 70
&END MGRID
&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-07
MAX_SCF 20
&OT ON
MINIMIZER DIIS
PRECONDITIONER FULL_ALL
&END OT
&OUTER_SCF
MAX_SCF 100
&END OUTER_SCF
&END SCF
&END
&SUBSYS
&CELL
ABC [angstrom] 10.00 10.00 10.00
&END CELL
&TOPOLOGY
COORD_FILE_NAME new_water.pdb
COORD_FILE_FORMAT PDB
&END
&KIND H
BASIS_SET TZV2P-GTH-q1
POTENTIAL GTH-BLYP-q1
&END KIND
&KIND O
BASIS_SET TZV2P-GTH-q6
POTENTIAL GTH-PBE-q6
&END KIND
&END SUBSYS
&END FORCE_EVAL
&MOTION
&GEO_OPT
OPTIMIZER LBFGS
MAX_ITER 100
MAX_FORCE 1.0E-03
RMS_DR 1.0E-03
RMS_FORCE 1.0E-03
MAX_DR 1.0E-03
&LBFGS
&END
&END
&MD
ENSEMBLE NVT
TEMPERATURE [K] 300
TIMESTEP [fs] 0.5
STEPS 30000
&THERMOSTAT
REGION MASSIVE
TYPE GLE
&GLE
NDIM 5
A_SCALE [ps^-1] 1.00
A_LIST 1.859575861256e+2 2.726385349840e-1 1.152610045461e+1
-3.641457826260e+1 2.317337581602e+2
A_LIST -2.780952471206e-1 8.595159180871e-5 7.218904801765e-1
-1.984453934386e-1 4.240925758342e-1
A_LIST -1.482580813121e+1 -7.218904801765e-1 1.359090212128e+0
5.149889628035e+0 -9.994926845099e+0
A_LIST -1.037218912688e+1 1.984453934386e-1 -5.149889628035e+0
2.666191089117e+1 1.150771549531e+1
A_LIST 2.180134636042e+2 -4.240925758342e-1 9.994926845099e+0
-1.150771549531e+1 3.095839456559e+2
&END GLE
&END THERMOSTAT
&END
&PRINT
&TRAJECTORY
ADD_LAST NUMERIC
FORMAT XYZ
&EACH
MD 1
&END EACH
&END TRAJECTORY
&RESTART
ADD_LAST NUMERIC
&EACH
MD 1
&END EACH
&END RESTART
&RESTART_HISTORY ON
&END
&FORCES
FILENAME forces_nvt
&EACH
MD 1
&END EACH
&END FORCES
&STRESS
FILENAME stress_nvt
&EACH
MD 1
&END EACH
&END STRESS
&END PRINT
&END
#####################
--
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/46c77221-04e6-4cd9-ad0b-d46d3a275d5cn%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230801/2aa178c7/attachment.htm>
More information about the CP2K-user
mailing list