[CP2K:7317] How to set the positions of QMMM cell explicitly in QMMM calculation

[ganta.pra... at gmail.com]

Dear Teo,

I have a system as shown below and the box is QMMM box with electron 
density plotted around QM atoms.  I have QMMM periodicity in XZ direction.

I dont want CP2K to center QM system, based on min-max of QM atoms (I fear 
if the QMMM grid is not aligned properly to accommodate periodicity - any 
suggestions?) and therefore I have used CENTER NEVER. 

But I would like to use INITIAL_TRANSLATION_VECTOR 
<https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/QMMM.html#INITIAL_TRANSLATION_VECTOR> to 
translate my system in Y direction. It does not seem to work. Please let me 
know if you have any suggestions. 

<https://lh3.googleusercontent.com/-1RjbJuUC3Yw/Wz8_stXf0VI/AAAAAAAAD3U/OYcbEcc3XH85Aen5os89VtdZH8euiafbgCLcBGAs/s1600/QMMM.png>




Have a nice weekend.

Thanks and Regards,
Prasanth. 





On Wednesday, January 13, 2016 at 10:59:40 AM UTC+1, Teo wrote:
>
> &CELL in QMMM specify the cell of the QM region
> &CELL in &SUBSYS specify the cell of the entire system
>
> The QM cell is always translated in order to have the QM system centred in 
> the QM cell: you should not care about where it is located.
>
> Teo
>
> On 13 Jan 2016, at 10:15, Geng Sun <sung... at gmail.com <javascript:>> 
> wrote:
>
> Hello  everyone,
>
>        I have a question about how to set the QMMM calculation,
>
>        Under the example folder from the cp2k distributions, one of the 
> input file looks like below : (only the QMMM and SUBSYS are shown)
>
>        In the input file, we have two CELL sections in QMMM and SUBSYS. 
>  It seems that the cell in QMMM is encompassed by cell in SUBSYS
>
>        Here is my question,  what's the relative positions of  cell/QMMM 
> in cell/SUBSYS ?
>
>        is  cell/QMMM movable in MD simulations ?  what's the meaning of 
> the keywors   INITIAL_TRANSLATION_VECTOR 
> <https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/QMMM.html#INITIAL_TRANSLATION_VECTOR>
>  ?
>
> Thanks 
>
> Geng Sun
>
>   &QMMM
>     MM_POTENTIAL_FILE_NAME ../MM_POTENTIAL
>     &CELL
>       ABC 6.0 6.0 6.0
>     &END CELL
>     ECOUPL GAUSS
>     NOCOMPATIBILITY
>     &INTERPOLATOR
>       EPS_R 1.0e-14
>       EPS_X 1.0e-14
>       MAXITER 100
>     &END INTERPOLATOR
>     &MM_KIND H
>       RADIUS 0.44
>     &END MM_KIND
>     &MM_KIND O
>       RADIUS 0.78
>     &END MM_KIND
>     &QM_KIND H
>       MM_INDEX 2
>       MM_INDEX 3
>     &END QM_KIND
>     &QM_KIND O
>       MM_INDEX 1
>     &END QM_KIND
>   &END QMMM
>   &SUBSYS
>     &CELL
>       ABC 24.955 24.955 24.955
>     &END CELL
>     &COORD
>    O     0.000000     0.000000     0.000000    H2O1
>    H     0.000000     0.000000     1.000000    H2O1
>    H     0.942809     0.000000    -0.333333    H2O1
>    O    -1.617979    -0.948062    -2.341650    H2O2
>    H    -2.529195    -1.296822    -2.122437    H2O2
>    H    -1.534288    -0.833088    -3.331486    H2O2
>    O    -1.447990     2.117783     1.555094    H2O3
>    H    -1.501128     2.645178     2.403050    H2O3
>    H    -2.090603     1.352766     1.597519    H2O3
>     &END COORD
>     &KIND H
>       BASIS_SET SZV-GTH
>       POTENTIAL GTH-PADE-q1
>     &END KIND
>     &KIND O
>       BASIS_SET SZV-GTH
>       POTENTIAL GTH-PADE-q6
>     &END KIND
>     &TOPOLOGY
>     &END TOPOLOGY
>   &END SUBSYS
>
> -- 
> 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+... at googlegroups.com <javascript:>.
> To post to this group, send email to cp... at googlegroups.com <javascript:>.
> Visit this group at https://groups.google.com/group/cp2k.
> For more options, visit https://groups.google.com/d/optout.
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180706/bdfb59c0/attachment.htm>