From owner-chemistry@ccl.net Mon Feb 23 01:33:01 2009 From: "Leyla Ramin l.ramin- -aeromech.usyd.edu.au" To: CCL Subject: CCL: GAMESS - "THE WORK ARRAY FOR CM2 IS NOT SUFFICIENT" Message-Id: <-38701-090223012942-10918-+0+zdxx0K8deVJ+wpDogQQ]|[server.ccl.net> X-Original-From: "Leyla Ramin" Content-class: urn:content-classes:message Content-Type: multipart/alternative; boundary="----_=_NextPart_001_01C9957F.F54E3F5B" Date: Mon, 23 Feb 2009 17:26:34 +1100 MIME-Version: 1.0 Sent to CCL by: "Leyla Ramin" [l.ramin^-^aeromech.usyd.edu.au] This is a multi-part message in MIME format. ------_=_NextPart_001_01C9957F.F54E3F5B Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hello, =20 I am writing an MD code for lennard-jones type of molecules using = neighbour list and link cell method, However the total running average = of energy in the system is not conserved and it goes down by time and I = can't find any reason for that. Do you have any suggestion? =20 Leyla Ramin =20 PhD Candidate School of Mechanical Engineering University of Sydney Sydney, NSW 2006, Australia E-mail: lram5228_+_usyd.edu.au =20 =20 ------_=_NextPart_001_01C9957F.F54E3F5B Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable =0A= =0A= =0A= =0A=
=0A=
Hello,
=0A=
 
=0A=
I am writing an MD code for = lennard-jones type of molecules using neighbour list and link cell = method, However the total running average of energy in the system is not = conserved and it goes down by time and I can't find any reason for that. = Do you have any suggestion?
=0A=
 
=0A=
=0A=
=0A=
Leyla = Ramin
=0A=
 
=0A=
 PhD = Candidate
=0A=
 School of = Mechanical Engineering
=0A=
 University of = Sydney
=0A=
 Sydney, NSW = 2006, Australia
=0A=
=0A=

 
------_=_NextPart_001_01C9957F.F54E3F5B-- From owner-chemistry@ccl.net Mon Feb 23 04:29:01 2009 From: "Leyla Ramin lram5228*_*usyd.edu.au" To: CCL Subject: CCL: Energy in MD Message-Id: <-38702-090223013228-11244-dR2HXvC7Tru17AJ86AI1rA-x-server.ccl.net> X-Original-From: "Leyla Ramin" Date: Mon, 23 Feb 2009 01:32:23 -0500 Sent to CCL by: "Leyla Ramin" [lram5228\a/usyd.edu.au] Hello, I am writing an MD code for lennard-jones type of molecules using neighbour list and link cell method, However the total running average of energy in the system is not conserved and it goes down by the time and I can't find any reason for that. Do you have any suggestion? Regards Leyla Ramin From owner-chemistry@ccl.net Mon Feb 23 05:43:01 2009 From: "neeraj misra neerajmisra(_)hotmail.com" To: CCL Subject: CCL: interaction study Message-Id: <-38703-090223054138-32521-5c75bLucz3YvDPJO1fT9qA^^^server.ccl.net> X-Original-From: "neeraj misra" Date: Mon, 23 Feb 2009 05:41:35 -0500 Sent to CCL by: "neeraj misra" [neerajmisra:+:hotmail.com] I WOULD BE THANKFUL IF ANY OF YOU COULD SUGGEST ME SOME SOFTWARE/METHOD WHICH CAN HELP IN STUDYING THE INTERACTION OF A NANOPARTICLE WITH BIOLOGICAL MATTER. From owner-chemistry@ccl.net Mon Feb 23 07:10:00 2009 From: "Victor Bernshtein chr21vb-x-technion.ac.il" To: CCL Subject: CCL: Energy in MD Message-Id: <-38704-090223061330-15667-soKY9vKel0mcYy0+6Wh7Hw~!~server.ccl.net> X-Original-From: Victor Bernshtein Content-Type: text/plain; charset="us-ascii"; format=flowed Date: Mon, 23 Feb 2009 12:37:09 +0200 Mime-Version: 1.0 Sent to CCL by: Victor Bernshtein [chr21vb-.-technion.ac.il] Hi, In my experience, energy conservation problem usually indicates on a large integration time-step or/and bugs in forces calculations. I would check derivatives. Remember that these are partial (over x, y, and z) derivatives. Good luck, Victor At 08:32 AM 2/23/2009, you wrote: >Sent to CCL by: "Leyla Ramin" [lram5228\a/usyd.edu.au] >Hello, > >I am writing an MD code for lennard-jones type of molecules using >neighbour list and link cell method, However the total running >average of energy in the system is not conserved and it goes down by >the time and I can't find any reason for that. Do you have any suggestion? > >Regards >Leyla Ramin From owner-chemistry@ccl.net Mon Feb 23 07:45:00 2009 From: "Lukasz Cwiklik cwiklik*gmail.com" To: CCL Subject: CCL: Energy in MD Message-Id: <-38705-090223073734-13407-FKO5Y/sSYtKfNpHpx/lWTw:server.ccl.net> X-Original-From: Lukasz Cwiklik Content-Transfer-Encoding: quoted-printable Content-Type: text/plain; charset=ISO-8859-1 Date: Mon, 23 Feb 2009 14:37:20 +0200 MIME-Version: 1.0 Sent to CCL by: Lukasz Cwiklik [cwiklik~!~gmail.com] On Mon, Feb 23, 2009 at 8:32 AM, Leyla Ramin lram5228*_*usyd.edu.au wrote: > Sent to CCL by: "Leyla Ramin" [lram5228\a/usyd.edu.au] > I am writing an MD code for lennard-jones type of molecules using neighbo= ur list and link cell method, However the total running average of energy i= n the system is not conserved and it goes down by the time and I can't find= any reason for that. Do you have any suggestion? Hello, An easy test, see if conservation of energy changes with shorter timestep. Another option, if cutoff is used, see if cutoff value changes the conserved quantity. Best, Lukasz --=20 Lukasz Cwiklik http://cwiklik.wordpress.com From owner-chemistry@ccl.net Mon Feb 23 08:22:00 2009 From: "Konrad Hinsen hinsen|-|cnrs-orleans.fr" To: CCL Subject: CCL: Energy in MD Message-Id: <-38706-090223081604-11794-Vn635zLIjODI6fSDeFSUIQ]~[server.ccl.net> X-Original-From: "Konrad Hinsen" Date: Mon, 23 Feb 2009 08:16:00 -0500 Sent to CCL by: "Konrad Hinsen" [hinsen]|[cnrs-orleans.fr] On Feb 23, 2009, at 7:32, Leyla Ramin lram5228*_*usyd.edu.au wrote: > I am writing an MD code for lennard-jones type of molecules using neighbour list and link cell > method, However the total running average of energy in the system is not conserved and it goes > down by the time and I can't find any reason for that. Do you have any suggestion? First, check that the forces are calculated correctly. A simple yet effective way is to compare the forces to numerical gradients of the potential energy function calculated with a small step size. By decreasing the step size, the numerical derivatives must converge to the separately calculated forces if those are correct. Once you have confidence in the forces, run your MD simulation with a very small step size, say 1/10 of what you would normally use. If the energy is still not constant, something is probably wrong with your integrator. If you get a constant energy for a very small step size but not for a normal one, check the Lennard-Jones parameters, in particular the cutoff and its implementation through the link-cell method. Another way to investigate cutoff problems is to compare, for a small systems, each energy and force term to the corresponding one calculated with infinite cutoff. -- --------------------------------------------------------------------- Konrad Hinsen Centre de Biophysique Molculaire, CNRS Orlans Synchrotron Soleil - Division Expriences Saint Aubin - BP 48 91192 Gif sur Yvette Cedex, France Tel. +33-1 69 35 97 15 E-Mail: hinsen * cnrs-orleans.fr Web: http://dirac.cnrs-orleans.fr/~hinsen/ --------------------------------------------------------------------- From owner-chemistry@ccl.net Mon Feb 23 08:56:01 2009 From: "Vincent Leroux vincent.leroux() loria.fr" To: CCL Subject: CCL: Energy in MD Message-Id: <-38707-090223082213-14497-4PzlZAuaHsQwQECSAKDaxA||server.ccl.net> X-Original-From: Vincent Leroux Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=ISO-8859-1; format=flowed Date: Mon, 23 Feb 2009 14:22:06 +0100 MIME-Version: 1.0 Sent to CCL by: Vincent Leroux [vincent.leroux . loria.fr] Hi Leyla, The reason might be that a mechanics system being allowed to lower its energy will do so until it reaches a minimum state... Did you implement the thermostat, barostat and piston in your MD code? VL Leyla Ramin lram5228*_*usyd.edu.au a écrit : > Sent to CCL by: "Leyla Ramin" [lram5228\a/usyd.edu.au] > Hello, > > I am writing an MD code for lennard-jones type of molecules using neighbour list and link cell method, However the total running average of energy in the system is not conserved and it goes down by the time and I can't find any reason for that. Do you have any suggestion? > > Regards > Leyla Ramin > > From owner-chemistry@ccl.net Mon Feb 23 09:32:00 2009 From: "Ting Li Ting.Li##chem.kuleuven.be" To: CCL Subject: CCL: Energy in MD Message-Id: <-38708-090223072543-8207-dX/Vu6ELg3ec1eddPKRB7w * server.ccl.net> X-Original-From: Ting Li Content-Transfer-Encoding: 7bit Content-Type: text/plain; charset=ISO-8859-1; format=flowed Date: Mon, 23 Feb 2009 12:25:51 +0100 MIME-Version: 1.0 Sent to CCL by: Ting Li [Ting.Li%x%chem.kuleuven.be] Hi, Leyla, Slow energy drifting is one of typical issue in the implementation of MD code. There are many reasons, including too large timestep, inappropriate intergrator, boundary problem, potential truncation and energy conservation characteristic of the implemented algorithm. To understand all of these, you have to consult the simulation books flying out there. And you didn't mention what ensemble you are simulating, so it's very difficult to judge the source of issue. My suggestion is instead of handling a complex system, simply use the LJ fluid model to verify your code. You may choose NVE ensemble to check the energy conservation, because this will isolate the thermostat/barostat problem out. I also suggest you follow the case studies in Smit and Frenkel's book. Their book is really good for simulation developers. Best regards, Ting Li --------------------------------------------------------- Molecular and Nanomaterials, Department of Chemistry Katholieke Universiteit Leuven (KULeuven) Celestijnenlaan 200F, B-3001 Heverlee, Belgium Email: Ting.Li**chem.kuleuven.be Tel: 32-16-327614 or 327481 http://simuworld.linkka.com --------------------------------------------------------- Leyla Ramin lram5228*_*usyd.edu.au wrote: > Sent to CCL by: "Leyla Ramin" [lram5228\a/usyd.edu.au] > Hello, > > I am writing an MD code for lennard-jones type of molecules using neighbour list and link cell method, However the total running average of energy in the system is not conserved and it goes down by the time and I can't find any reason for that. Do you have any suggestion? > > Regards > Leyla Ramin> > > > > Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm From owner-chemistry@ccl.net Mon Feb 23 10:07:01 2009 From: "Jamie Platts Platts-x-cardiff.ac.uk" To: CCL Subject: CCL: Any molecular mechanics model good at pi-pi stacking calculation? Message-Id: <-38709-090223081501-11291-FTXi51Bp2IOE0vnDhT/gXA() server.ccl.net> X-Original-From: Jamie Platts Content-Type: multipart/alternative; boundary="=_alternative 004127DE80257566_=" Date: Mon, 23 Feb 2009 11:51:38 +0000 MIME-Version: 1.0 Sent to CCL by: Jamie Platts [Platts_._cardiff.ac.uk] This is a multipart message in MIME format. --=_alternative 004127DE80257566_= Content-Type: text/plain; charset="US-ASCII" Dear Alex For semi-empirical, standard methods are unlikely to work well at all. However, there are a number of efforts to include the missing dispersion effects - see for example: SCC-DFTB-D: Elstner, M.; Hobza, P.; Frauenheim, T.; Suhai, S.; Kaxiras, E. J. Chem. Phys., 2001, 114, 5149. AM1-D and PM3-D: J.P. McNamara, I.H. Hillier, PCCP, 2007, 9, 2362. OMx-D: T. Tuttle, W. Thiel, PCCP, 2008, 10, 2159. 200 atoms is probably not beyond the scope of modern DFT methods, especially if you use the resolution of identity (RI) or density fitting (DF) approximations which can speed things up dramatically. There are many dispersion-corrected DFT methods now available in a number of packages. Some DFT methods can also give good descriptions of stacking, for instance the M05- and M06- familes from Truhlar's group, Becke's half-and-half (BHandH) - even LDA isn't bad! Most molecular mechanics forcefields should include non-bonded terms that represent pi-pi stacking, although details will depend on how close your system is to those used to parameterize the forcefield. Hope this helps, Jamie > From: "Alex Liu ylvn8_-_umkc.edu" To: "Platts, Jamie " Date: 22/02/2009 05:03 Subject: CCL: Any molecular mechanics model good at pi-pi stacking calculation? Sent to CCL by: "Alex Liu" [ylvn8,,umkc.edu] The molecules I am working on have over 200 atoms. I am wondering is there any MM or semiempirical package that is able to handle pi-pi stacking engeries relatively accurately? Alex ylvn8{=}umkc.eduhttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt--=_alternative 004127DE80257566_= Content-Type: text/html; charset="US-ASCII"
Dear Alex

For semi-empirical, standard methods are unlikely to work well at all. However, there are a number of efforts to include the missing dispersion effects - see for example:
SCC-DFTB-D: Elstner, M.; Hobza, P.; Frauenheim, T.; Suhai, S.; Kaxiras, E. J. Chem. Phys., 2001, 114, 5149.
AM1-D and PM3-D: J.P. McNamara, I.H. Hillier, PCCP, 2007, 9, 2362.
OMx-D: T. Tuttle, W. Thiel, PCCP, 2008, 10, 2159.

200 atoms is probably not beyond the scope of modern DFT methods, especially if you use the resolution of identity (RI) or density fitting (DF) approximations which can speed things up dramatically. There are many dispersion-corrected DFT methods now available in a number of packages. Some DFT methods can also give good descriptions of stacking, for instance the M05- and M06- familes from Truhlar's group, Becke's half-and-half (BHandH) - even LDA isn't bad!

Most molecular mechanics forcefields should include non-bonded terms that represent pi-pi stacking, although details will depend on how close your system is to those used to parameterize the forcefield.

Hope this helps,
Jamie



From: "Alex Liu ylvn8_-_umkc.edu" <owner-chemistry+*+ccl.net>
To: "Platts, Jamie " <platts+*+cardiff.ac.uk>
Date: 22/02/2009 05:03
Subject: CCL: Any molecular mechanics model good at pi-pi stacking calculation?






Sent to CCL by: "Alex  Liu" [ylvn8,,umkc.edu]
The molecules I am working on have over 200 atoms. I am wondering is there any MM or semiempirical package that is able to handle pi-pi stacking engeries relatively accurately?

Alex
ylvn8{=}umkc.edu


     
http://www.ccl.net/cgi-bin/ccl/send_ccl_message
     
http://www.ccl.net/cgi-bin/ccl/send_ccl_message
     
http://www.ccl.net/chemistry/sub_unsub.shtml

Before posting, check wait time at:
http://www.ccl.net

Job:
http://www.ccl.net/jobs
Conferences:
http://server.ccl.net/chemistry/announcements/conferences/

Search Messages:
http://www.ccl.net/htdig  (login: ccl, Password: search)
     
http://www.ccl.net/spammers.txt

RTFI:
http://www.ccl.net/chemistry/aboutccl/instructions/




--=_alternative 004127DE80257566_=-- From owner-chemistry@ccl.net Mon Feb 23 10:51:01 2009 From: "Hui Yang Hui.Yang::utoledo.edu" To: CCL Subject: CCL: GAMESS - "THE WORK ARRAY FOR CM2 IS NOT SUFFICIENT" Message-Id: <-38710-090223094406-16549-y+gCr+eSyyz9uA+o1zoVZA[a]server.ccl.net> X-Original-From: "Hui Yang" Content-Type: multipart/alternative; boundary="----=_NextPart_000_001F_01C9959B.3F02FE00" Date: Mon, 23 Feb 2009 09:43:56 -0500 MIME-Version: 1.0 Sent to CCL by: "Hui Yang" [Hui.Yang:_:utoledo.edu] This is a multi-part message in MIME format. ------=_NextPart_000_001F_01C9959B.3F02FE00 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hi, Many thanks to Jiabo and Chris. Now we've updated to the latest version = of GAMESSPLUS and the problem is solved. Thank you again for your time. Hui Yang Ph.D student University of Toledo ----- Original Message -----=20 From: Christopher Cramer cramer]^[umn.edu=20 To: Yang, Hui =20 Sent: Sunday, February 08, 2009 10:03 AM Subject: CCL: GAMESS - "THE WORK ARRAY FOR CM2 IS NOT SUFFICIENT" Jiabo coded CM2 in GAMESSPLUS, so his memory is very good. Relevant to the current version: In the latest version of GAMESSPLUS (v2008-2) the MAX size of the CM2 work array is defined by the formula: MXSATM*(MXSATM+1)+3*MXSBAS*(MXSBAS+1)/2+2*MXSBAS*MXSBAS+1 where MXSATM is the maximum number of atoms and MXSBAS is the maximum = number of contracted basis functions permissible for solvation calculations = using SMx models in GAMESSPLUS. C CURRENT GAMESSPLUS LIMITATIONS FOR SOLVATION MODULE: C ON TOTAL NUMBER OF ATOMS IS MXSATM=3D250 C ON TOTAL NUMBER OF CONTRACTED BASIS FUNCTIONS IS MXSBAS=3D2500. Thus MAX=3D250*(250+1)+3*2500*(2500+1)/2+2*2500*2500+1=3D =3D21,941,501>10,282,931 In the latest version of GAMESSPLUS (v2008-2) the limits should = suffice. http://comp.chem.umn.edu/gamessplus/ Perhaps Hui uses an older version of GAMESSPLUS with smaller limits on MXSATM and MXSBAS. Chris On Feb 7, 2009, at 3:09 PM, Jiabo Li jiaboli ~~ yahoo.com wrote: Hi Hui, I assume that you compiled GAMESS with GAMESSPLUS for SMx = solvation model. The message you got indicates that the scratch space = needed for your calculation is larger than the fixed scratch space for = SMx calculations. The only way to solve this problem is to modify the = scratch space parameter the GAMESSPLUS source code and recompile it.=20 Here is the steps that you can try. 1. Find the source file smx.src in the GAMESSPLUS package. 2. Do a global change of MAXWK=3D5500000 into = MAXWK=3D110000000 in smx.src 3. Recompile the smx and relink gamess to create a new gamess = executable. Please let me know if this solves your problem or not. Regards, Jiabo --- On Fri, 2/6/09, Yang, Hui hyang3{=3D}UTNet.UToledo.Edu = wrote: From: Yang, Hui hyang3{=3D}UTNet.UToledo.Edu = Subject: CCL: GAMESS - "THE WORK ARRAY FOR CM2 IS NOT = SUFFICIENT" To: "Li, Jiabo " Date: Friday, February 6, 2009, 12:54 PM Hi, I am trying to use GAMESS to calculate the solvation free = energy of a biological system which has 584 electrons. Here is the error. ... DIRECT SCF CALCULATION, SCHWRZ=3DT FDIFF=3DF THE WORK ARRAY FOR CM2 IS NOT SUFFICIENT MAX =3D 5600000 REQUIRMENT =3D 10282931 EXECUTION OF GAMESS TERMINATED ABNORMALLY AT Fri Feb 6 = 15:08:45 2009 STEP CPU TIME =3D .95 TOTAL CPU TIME =3D 674.0 ( = 11.2 MIN) TOTAL WALL CLOCK TIME=3D 675.2 SECONDS, CPU = UTILIZATION IS 99.83% 22548143 WORDS OF DYNAMIC MEMORY USED *** ERROR TERMINATION *** IN COMPUTE PROCESS 0 ... And here is the header to my input file: $CONTRL SCFTYP=3DRHF RUNTYP=3DENERGY ICHARG=3D0 MAXIT=3D500 NPRINT=3D-5 EXETYP=3DRUN $END $SYSTEM TIMLIM=3D6000 MWORDS=3D900 $END $SCF DIRSCF=3D.TRUE. FDIFF=3D.FALSE. $END $BASIS GBASIS=3DN31 NGAUSS=3D6 NDFUNC=3D1 DIFFSP=3D.TRUE. POLAR=3DPOPLE $END $GUESS GUESS=3DHUCKEL $END $CM2 ISCRF=3D2 ICMD=3D8C1 IAQU=3D1 Dielec=3D78.0 $END $DATA HF/6-31+G*C1 C 6.0 45.984800 111.010300 25.847100 ...... Does anyone know what does "THE WORK ARRAY FOR CM2 IS NOT = SUFFICIENT" mean? And what is "MAX =3D 5600000"? Thank you for your time. Hui =20 -- Christopher J. Cramer University of Minnesota Department of Chemistry 207 Pleasant St. SE Minneapolis, MN 55455-0431 -------------------------- Phone: (612) 624-0859 || FAX: (612) 626-2006 Mobile: (952) 297-2575 cramer a umn.edu http://pollux.chem.umn.edu/~cramer (website includes information about the textbook "Essentials of Computational Chemistry: Theories and Models, 2nd Edition") ------=_NextPart_000_001F_01C9959B.3F02FE00 Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable
Hi,
 
Many thanks to Jiabo and Chris. = Now we've=20 updated to the latest version of GAMESSPLUS and the problem is=20 solved.
 
Thank you again for your = time.
 
Hui Yang
Ph.D student
University of Toledo
----- Original Message -----
From:=20 Christopher Cramer = cramer]^[umn.edu=20
Sent: Sunday, February 08, 2009 = 10:03=20 AM
Subject: CCL: GAMESS - "THE = WORK ARRAY=20 FOR CM2 IS NOT SUFFICIENT"

Jiabo coded CM2 in GAMESSPLUS, so his memory is very = good.

Relevant to the current version:

In the latest version of
GAMESSPLUS (v2008-2) the MAX size of = the CM2=20 work array is defined by
the=20 = formula:
MXSATM*(MXSATM+1)+3*MXSBAS*(MXSBAS+1)/2+2*MXSBAS*MXSBAS+1
= where=20 MXSATM is the maximum number of atoms and MXSBAS is the maximum = number
of=20 contracted basis functions permissible for solvation calculations = using
SMx=20 models in GAMESSPLUS.
C CURRENT GAMESSPLUS LIMITATIONS FOR = SOLVATION=20 MODULE:
C  ON TOTAL NUMBER OF ATOMS IS MXSATM=3D250
C =  ON TOTAL=20 NUMBER OF CONTRACTED BASIS FUNCTIONS IS MXSBAS=3D2500.

Thus=20 = MAX=3D250*(250+1)+3*2500*(2500+1)/2+2*2500*2500+1=3D
=3D21,941,501>= 10,282,931

In=20 the latest version of GAMESSPLUS (v2008-2) the limits should = suffice.
http://comp.chem.umn.edu/ga= messplus/

Perhaps=20 Hui uses an older version of GAMESSPLUS with smaller
limits on = MXSATM=20 and
MXSBAS.

Chris

On Feb 7, 2009, at 3:09 PM, Jiabo Li jiaboli ~~ yahoo.com = wrote:

--


Christopher J. Cramer

University of Minnesota

Department of Chemistry

207 Pleasant St. SE

Minneapolis, MN 55455-0431

--------------------------

Phone:  (612) 624-0859 || = FAX:  (612) = 626-2006

Mobile: (952) 297-2575

cramer a = umn.edu

http://pollux.chem.umn.edu/~c= ramer

(website includes information about the textbook = "Essentials

    of Computational=20 Chemistry:  Theories and Models, = 2nd=20 Edition")



------=_NextPart_000_001F_01C9959B.3F02FE00-- From owner-chemistry@ccl.net Mon Feb 23 11:38:01 2009 From: "Richard Leo Wood woodx278]*[umn.edu" To: CCL Subject: CCL: Any molecular mechanics model good at pi-pi stacking calculation? Message-Id: <-38711-090223113629-881-bP80qyjeyaNQQXdqqsnzsg||server.ccl.net> X-Original-From: "Richard Leo Wood" Date: Mon, 23 Feb 2009 11:36:26 -0500 Sent to CCL by: "Richard Leo Wood" [woodx278=-=umn.edu] I think a lot of this depends upon the type of molecules that you have. One should do a validation study before choosing a particular functional to use for your set of molecules. Be aware that if you have sulfur containing molecules, straight DFTB is not a good method, but its "followup", DFTB+ works very well. Richard > From: owner-chemistry(!)ccl.net > To: rwoodphd(!)msn.com > Subject: CCL: Any molecular mechanics model good at pi-pi stacking calculation? > Date: Mon, 23 Feb 2009 11:51:38 +0000 > > Dear Alex > > For semi-empirical, standard methods are unlikely to work well at all. > However, there are a number of efforts to include the missing dispersion > effects - see for example: > SCC-DFTB-D: Elstner, M.; Hobza, P.; Frauenheim, T.; Suhai, S.; Kaxiras, E. > J. Chem. Phys., 2001, 114, 5149. > AM1-D and PM3-D: J.P. McNamara, I.H. Hillier, PCCP, 2007, 9, 2362. > OMx-D: T. Tuttle, W. Thiel, PCCP, 2008, 10, 2159. > > 200 atoms is probably not beyond the scope of modern DFT methods, > especially if you use the resolution of identity (RI) or density fitting > (DF) approximations which can speed things up dramatically. There are many > dispersion-corrected DFT methods now available in a number of packages. > Some DFT methods can also give good descriptions of stacking, for instance > the M05- and M06- familes from Truhlar's group, Becke's half-and-half > (BHandH) - even LDA isn't bad! > > Most molecular mechanics forcefields should include non-bonded terms that > represent pi-pi stacking, although details will depend on how close your > system is to those used to parameterize the forcefield. > > Hope this helps, > Jamie >
Hi Hui,
 
I assume that you compiled GAMESS with GAMESSPLUS for SMx = solvation model. The message you got indicates that the = scratch space=20 needed for your calculation is larger than the fixed scratch = space for=20 SMx calculations. The only way to solve this problem is to = modify the=20 scratch space parameter the GAMESSPLUS source code and = recompile it.=20
 
Here is the steps that you can try.
 
1. Find the source file smx.src in the GAMESSPLUS = package.
2. Do a global change of MAXWK=3D5500000 into = MAXWK=3D110000000 in=20 smx.src
3. Recompile the smx and relink gamess to create a new = gamess=20 executable.
 
Please let me know if this solves your problem or = not.
 
Regards,
 
Jiabo

--- On Fri, 2/6/09, Yang, Hui=20 hyang3{=3D}UTNet.UToledo.Edu=20 <owner-chemistry=3D-=3Dccl.net> = wrote:
From:=20 Yang, Hui hyang3{=3D}UTNet.UToledo.Edu=20 <owner-chemistry=3D-=3Dccl.net>
Subject: CCL: = GAMESS - "THE=20 WORK ARRAY FOR CM2 IS NOT SUFFICIENT"
To: "Li, Jiabo "=20 <jiaboli=3D-=3Dyahoo.com>
Date: Friday, February 6, = 2009, 12:54=20 PM

Hi,

I am trying to use GAMESS to = calculate=20 the solvation free energy of a biological system which has = 584=20 electrons.

Here is the = error.

...
 DIRECT SCF=20 CALCULATION, SCHWRZ=3DT   FDIFF=3DF
  THE = WORK ARRAY=20 FOR CM2 IS NOT SUFFICIENT
  MAX =3D  5600000 = REQUIRMENT =3D=20 10282931
 EXECUTION OF GAMESS TERMINATED ABNORMALLY = AT Fri=20 Feb  6 15:08:45 2009
 STEP CPU TIME=20 =3D      .95 TOTAL CPU TIME=20 =3D      674.0 (   11.2=20 MIN)
 TOTAL WALL CLOCK = TIME=3D     =20 675.2 SECONDS, CPU UTILIZATION IS  = 99.83%
 22548143=20 WORDS OF DYNAMIC MEMORY USED
 *** ERROR TERMINATION = *** IN=20 COMPUTE PROCESS     0
...

And = here is=20 the header to my input file:
 $CONTRL SCFTYP=3DRHF=20 RUNTYP=3DENERGY=20 = ICHARG=3D0
        =20 MAXIT=3D500 NPRINT=3D-5 EXETYP=3DRUN $END
 $SYSTEM = TIMLIM=3D6000=20 MWORDS=3D900 $END
 $SCF DIRSCF=3D.TRUE. = FDIFF=3D.FALSE.=20 $END
 $BASIS  GBASIS=3DN31 NGAUSS=3D6=20 = NDFUNC=3D1
        =20 DIFFSP=3D.TRUE. POLAR=3DPOPLE $END
 $GUESS  = GUESS=3DHUCKEL=20 $END
 $CM2 ISCRF=3D2=20 ICMD=3D8C1
      IAQU=3D1=20 = Dielec=3D78.0
 $END
 $DATA
 HF/6-31+G*C1
 = ;C=20 6.0 45.984800 111.010300 = 25.847100
 ......

Does=20 anyone know what does "THE WORK ARRAY FOR CM2 IS NOT = SUFFICIENT"=20 mean? And what is "MAX =3D  5600000"?

Thank you = for your=20 = time.
Hui