in addition, how to write ORCAINFO?  Is my file right?

! RKS B3LYP/G SV(P) TightSCF Opt  or  RKS B3LYP/G SV(P) TightSCF Opt  


--- 11年11月7日,周一, Gerrit Groenhof <ggro...@gwdg.de> 写道:


发件人: Gerrit Groenhof <ggro...@gwdg.de>
主题: [gmx-users] Re 1. orca and qm/mm (xi zhao)
收件人: gmx-users@gromacs.org
日期: 2011年11月7日,周一,下午9:23


Try to remove these lines, or put something there. The input is ignored, but 
since strings are used as input (for use in multui-layer oniom), leaving blank 
causes an error.

Gerrit
On 7 Nov 2011, at 14:21, gmx-users-requ...@gromacs.org wrote:

> Send gmx-users mailing list submissions to
>     gmx-users@gromacs.org
> 
> To subscribe or unsubscribe via the World Wide Web, visit
>     http://lists.gromacs.org/mailman/listinfo/gmx-users
> or, via email, send a message with subject or body 'help' to
>     gmx-users-requ...@gromacs.org
> 
> You can reach the person managing the list at
>     gmx-users-ow...@gromacs.org
> 
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gmx-users digest..."
> 
> 
> Today's Topics:
> 
>   1. orca and qm/mm (xi zhao)
> 
> 
> ----------------------------------------------------------------------
> 
> Message: 1
> Date: Mon, 7 Nov 2011 20:54:04 +0800 (CST)
> From: xi zhao <zhaoxiitc2...@yahoo.com.cn>
> Subject: [gmx-users] orca and qm/mm
> To: gmx-users@gromacs.org
> Message-ID:
>     <1320670444.82465.yahoomailclas...@web15103.mail.cnb.yahoo.com>
> Content-Type: text/plain; charset="utf-8"
> 
> All users:
> 
>  
> According to http://wwwuser.gwdg.de/~ggroenh/qmmm.html#code, I want to build 
> qm/mm calculation by using gromacs 4.5.3 + orca 2.8.0. I set up the BASENAME, 
> ORCA_PATH and  BASENAME.ORCAINFO as told in the instruction. 
> BASENAME=pyp_qm
> here is the BASENAME.ORCAINFO file:
> ! RKS B3LYP/G SV(P) TightSCF Opt
>  
> here is the md file:
> integrator               = md
> tinit                    = 0
> dt                       = 0.001
> nsteps                   = 500
> nstcomm                  = 1
> comm_grps                = system
>  
> emtol                    = 100.0
> emstep                   = 0.001
> nstcgsteep               = 50
>  
> nstxout                  = 1
> nstvout                  = 1
> nstfout                  = 1
> nstlog                   = 1
> nstenergy                = 1
> nstxtcout                = 1
> xtc_grps                 = system
> energygrps               = QMatoms rest_Protein SOL
>  
> nstlist                  = 10
> ns_type                  = grid
> pbc                      = xyz
> rlist                    = 1.0
>  
> coulombtype              = cut-off
> rcoulomb                 = 1.4
> epsilon_r                = 1
> vdwtype                  = Cut-off
> rvdw                     = 1.4
>  
> tcoupl                   = berendsen
> tc-grps                  = rest_Protein SOL QMatoms
> tau_t                    = 0.1 0.1 0 ; QM atoms are uncoupled
> ref_t                    = 300 300 300
> Pcoupl                   = Berendsen
> pcoupltype               = isotropic
> tau_p                    = 1.0
> compressibility          = 4.5e-5
> ref_p                    = 1.0
>  
> free_energy              = no
> init_lambda              = 0
> delta_lambda             = 0
> QMMM                     = yes
> QMMM-grps                = QMatoms
> QMmethod                 =
> QMbasis                  =
> QMMMscheme               = normal
> QMcharge                 = -1
> CASelectrons             =
> CASorbitals              =
> SH                       =
>  
> gen_vel                  = no
> gen_temp                 = 300
> gen_seed                 = 173529
>  
> constraints              = all-bonds
> constraint_algorithm     = LINCS
> unconstrained_start      = yes
> shake_tol                = 0.0001
> lincs_order              = 4
> lincs_warnangle          = 30
> morse                    = no
>   According to the instruction “In the ORCAINFO-file the method, basis set 
>and all other ORCA-specific keywords must be given. (This also means that 
>QMmethod and QMbasis from the mdp-file are ignored).”, the QMmethod and 
>QMbasis are blank, 
>    But When grompp
> grompp_c -f pyp.mdp -c pyp.gro -p pyp.top -n pyp.ndx -o pyp_qm.tpr
> ……….
> Fatal error:
> Invalid QMMM input: 1 groups 0 basissets and 0 methods.
>  
> How to deal with it? Please help me!
> 
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: 
> http://lists.gromacs.org/pipermail/gmx-users/attachments/20111107/8c5c63c5/attachment-0004.html
> 
> ------------------------------
> 
> -- 
> gmx-users mailing list
> gmx-users@gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at 
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> 
> End of gmx-users Digest, Vol 91, Issue 36
> *****************************************

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to