[CP2K-user] How to limit the size of QM region in adbf-QMMM simulation?
Kalven Davis
mxh... at 163.com
Mon Jul 19 15:12:56 UTC 2021
I am using ADBF-qmmm to simulate the interfacial reaction of carbonyl
compounds. In the core region, a carbonyl compound and a hydrated ion are
selected. Because the core region involves protons, with the continuous
transfer of protons, the QM region becomes larger , and the amount of
calculation is increasing. I set the maximum number of atoms in QM area to
200, but the calculation stopped midway with an error: add_new_label tried
to add more atoms than allowed by &FORCE_MIXING&MAX_N_QM!How can I limit
the size of QM area? Thank you for any possible help.
FORCE_MIXING parameter:
&FORCE_MIXING
R_CORE 0.0 0.0
R_QM 3.0 3.5
R_BUF 3.0 3.5
QM_KIND_ELEMENT_MAPPING O O !Mapping from elements to QM_KINDs for
adaptively included atoms.
QM_KIND_ELEMENT_MAPPING H H
QM_KIND_ELEMENT_MAPPING C C
QM_KIND_ELEMENT_MAPPING N N
MAX_N_QM 200
MOMENTUM_CONSERVATION_TYPE EQUAL_A
MOMENTUM_CONSERVATION_REGION QM
QM_EXTENDED_SEED_IS_ONLY_CORE_LIST .TRUE.
&END FORCE_MIXING
&WALLS
TYPE QUADRATIC
WALL_SKIN 0.5 0.5 0.5
K [kcalmol] 100.0
&END WALLS
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20210719/9371ae23/attachment.htm>
More information about the CP2K-user
mailing list