[CP2K-user] MD with 1706 atoms - input suggestions?
MD Simulation
mdsimula... at gmail.com
Tue Feb 4 20:38:51 UTC 2020
Hello,
I have a hydrated organic compound, a total of 1706 atoms of C, H, O, N, F,
but mainly waters.
I would like to run molecular dynamics but SCF convergence is taking a very
long time.
I have tried LC_SCF but still a bit slow. I've attache my input and xyz
file. My big question right now....am I setting up the run correctly?
I have tried running on 240 CPUs and still slow. I am using CP2K 7.1
Any ideas would be appreciated. My input file is also copy/pasted below.
Thanks for the help!
Stacey
************************************************************************************
@SET COORD_FILENAME Hydrated_6FDA.xyz
&GLOBAL
PROJECT Hydrated_6FDA
RUN_TYPE MD
PRINT_LEVEL LOW
&END GLOBAL
&FORCE_EVAL
METHOD Quickstep
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME GTH_POTENTIALS
CHARGE 0
MULTIPLICITY 1
&QS
LS_SCF
EPS_DEFAULT 1.0E-10
&END QS
&MGRID
NGRIDS 5
CUTOFF 300.0
REL_CUTOFF 60.0
&END MGRID
&LS_SCF
MAX_SCF 100
EPS_FILTER 1.0E-6
EPS_SCF 1.0E-6
MU -0.05
MIXING_FRACTION 0.45
S_PRECONDITIONER NONE
PERFORM_MU_SCAN OFF
PURIFICATION_METHOD TRS4
&CURVY_STEPS
&END CURVY_STEPS
&END LS_SCF
&SCF
EPS_SCF 1.0E-6
MAX_SCF 200
&END SCF
&XC
&XC_FUNCTIONAL BLYP
&END XC_FUNCTIONAL
&VDW_POTENTIAL
POTENTIAL_TYPE PAIR_POTENTIAL
&PAIR_POTENTIAL
PARAMETER_FILE_NAME dftd3.dat
TYPE DFTD3
REFERENCE_FUNCTIONAL BLYP
&END
&END VDW_POTENTIAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 25.000 25.000 25.000
&END CELL
&TOPOLOGY
COORD_FILE_FORMAT XYZ
COORD_FILE_NAME ${COORD_FILENAME}
&END TOPOLOGY
&KIND H
ELEMENT H
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q1
&END KIND
&KIND C
ELEMENT C
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q4
&END KIND
&KIND N
ELEMENT N
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q5
&END KIND
&KIND O
ELEMENT O
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q6
&END KIND
&KIND F
ELEMENT F
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q7
&END KIND
&END SUBSYS
&END FORCE_EVAL
&MOTION
&MD
ENSEMBLE NVT
TIMESTEP [fs] 0.5
STEPS 10000
TEMPERATURE [K] 298.15
&THERMOSTAT
REGION MASSIVE
TYPE NOSE
&NOSE
TIMECON 100.0
&END NOSE
&END THERMOSTAT
&END MD
&PRINT
&TRAJECTORY
FORMAT PDB
&EACH
MD 20
&END EACH
&END TRAJECTORY
&END PRINT
&END MOTION
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200204/f91cf255/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Hydrated_6FDA.xyz
Type: chemical/x-xyz
Size: 97323 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200204/f91cf255/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LS-AIMD_input.inp
Type: chemical/x-gamess-input
Size: 1868 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200204/f91cf255/attachment.inp>
More information about the CP2K-user
mailing list