From owner-chemistry@ccl.net Thu Oct 28 00:36:01 2010 From: "xunlei ding dingxunlei-,-gmail.com" To: CCL Subject: CCL:G: freezing angle Error in internal coordinate system Message-Id: <-43024-101028003445-17635-FJWZI1aTdbEjVGPzaO0pJQ:server.ccl.net> X-Original-From: xunlei ding Content-Type: multipart/alternative; boundary=0016e6570eae1b853c0493a5dc8e Date: Thu, 28 Oct 2010 12:34:38 +0800 MIME-Version: 1.0 Sent to CCL by: xunlei ding [dingxunlei!=!gmail.com] --0016e6570eae1b853c0493a5dc8e Content-Type: text/plain; charset=ISO-8859-1 Hi All, If doing a full optimization, I think the best way is to use opt(cartesian). But for partical optimization, this does not work. Sometimes this error may be solved by using G03 instead of G09, but I don't think this is a good way. Ding 2010/9/27 Decai Yu decaiyu23]![yahoo.com > Dear All, > > I am trying to do a partial geometry optimization. > In my structures, I have a number of atoms fixed and I also have one angle > fixed. > When I changed the value of the fixed angle from 135 to 180, the smaller > angle runs are finished without error. > But for angle value larger than 175, it can only run one step (a full SCF > relaxation). > > The program ends with the following error: > > GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad > Berny optimization. > NTrRot= 421 NTRed= 574 NAtoms= 51 NSkip= 421 IsLin=F > Error in internal coordinate system. > > I would appreciate any help or advice. > > Best Regards, > > Decai > > > > > --0016e6570eae1b853c0493a5dc8e Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
Hi All,
=A0
If doing a full optimization, I think the best way=A0 is to use opt(ca= rtesian).
But for partical optimization, this does not work.
Sometimes this error may be solved by using G03 instead of G09, but I = don't think this is a good way.
=A0
Ding


=A0
2010/9/27 Decai Yu decaiyu23]![yahoo.com <owner-chemistry{=}ccl.net>
Dear All,
=A0
I am trying to do a partial=A0 geometry optimization.
In my structures, I have a number of atoms fixed and I also have one a= ngle fixed.
When I changed the value of the fixed angle from 135 to 180, the small= er angle runs are finished without error.
But for angle value larger than 175, it can only run one step (a full = SCF relaxation).
=A0
The program=A0ends with the following error:
=A0
=A0GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGra= dGrad
=A0Berny optimization.
=A0NTrRot=3D=A0=A0 421 NTRed=3D=A0=A0 57= 4 NAtoms=3D=A0=A0=A0 51 NSkip=3D=A0=A0 421 IsLin=3DF
=A0Error in interna= l coordinate system.

I would appreciate any help or advice.
=A0
Best Regards,
=A0
Decai
=A0


--0016e6570eae1b853c0493a5dc8e-- From owner-chemistry@ccl.net Thu Oct 28 04:05:00 2010 From: "Lars Goerigk lars.goerigk#,#uni-muenster.de" To: CCL Subject: CCL:G: Double hybrid functionals for excited states with Gaussian09 Message-Id: <-43025-101028040035-16019-/vR4robuXdGjEygBcI0W0Q^^server.ccl.net> X-Original-From: Lars Goerigk Content-Transfer-Encoding: 7bit Content-Type: Text/Plain; charset="iso-8859-15" Date: Thu, 28 Oct 2010 10:00:27 +0200 MIME-Version: 1.0 Sent to CCL by: Lars Goerigk [lars.goerigk _ uni-muenster.de] TD-B2PLYP works as follows: First, you perform a ground state B2LYP calculation. The perturbative correction is not needed if you are not interested in any ground state properties, because you are just interested in the resulting KS-orbitals. Then a standard TD-DFT treatment (i.e. TD-B2LYP) is carried out, which gives you excitation energies and transition moments of TD-BHLYP quality. The B2LYP orbitals and TD-B2LYP excitation energies and amplitudes are then used for a CIS(D) type treatment to correct the TD-B2LYP energies. I do not know if such a treatment is possible in Gaussian by adjusting some IOPs. I assume that the CIS(D) routine will crash if TD-DFT results are used as an input. However, the routine may not fail with an CIS like input. The Tamm-Dancoff- Approximation is basically CIS. I do not know if Gaussian does TDA-DFT. Basically you can then carry out the CIS(D) correction based on your TDA-B2LYP result. This gives you TDA-B2PLYP energies. This was also discussed in the 2007 Grimme/Neese paper. TDA-DFT usually gives you very similar energies to TD-DFT. But again, I do not know whether Gaussian can handle this with IOPs, but it's worth testing it. Maybe the Gaussian people can help you with the IOPs if you contact them; maybe it's even possible to do full TD-B2PLYP. Another option would be ORCA, which is freely available. It does TDA-B2PLYP, but not full TD-B2PLYP. I hope these information are somehow useful for you. Cheers, Lars On Wednesday 27 October 2010 05:08:27 pm you wrote: > Sent to CCL by: "Lara Martinez Fernandez" [lara.martinez*o*uam.es] > I am trying to run some TD-DFT calculation (for excited states) using > double hybrid functionals with Gaussian09. I have some doubts about how > TD-B2PLYP,TD-B2GP-PLYP work. If I dont misunderstood, first I have to run > a B2LYP calculation for the ground state and then use this orbitals and > amplitudes to run the CIS(D) correction for the rest of excited states. I > dont know if this procedure is possible with Gaussian09. Do anyone know > how to do that? I guess there must be any keyword to read the B2LYP (of > the ground state) orbitals and not optimise them in the CIS(D) > calculation. > > > Thanks in advance, Lara> From owner-chemistry@ccl.net Thu Oct 28 06:59:01 2010 From: "makowskm!^!chemia.uj.edu.pl" To: CCL Subject: CCL:G: MP2 Calculation Message-Id: <-43026-101028063634-27985-nzeXCRrDdzyEbVeo7mHnaw]|[server.ccl.net> X-Original-From: makowskm#chemia.uj.edu.pl Content-Transfer-Encoding: 8bit Content-Type: text/plain;charset=iso-8859-2 Date: Thu, 28 Oct 2010 12:36:26 +0200 (CEST) MIME-Version: 1.0 Sent to CCL by: makowskm^^chemia.uj.edu.pl > > Sent to CCL by: "Leonardo Moreira da Costa" [leomdcosta a yahoo.com.br] > I have a very simple doubt. > > It is said in the Gaussian help (C:\G03W\help\k_maxdisk.htm), about > maxDisk, that MP2 energies and gradients obey MaxDisk, which must be at > least 2ON2. > > It is said that, but they do not specify what are the ''O'' and ''N'' > parameters of the formula. > > Anyone can, please, explain me what are these two parameters? I have > looked for it in many resources, but I did not obatin a certain meaning > (number of the electrons, orbitals, atoms, ...). O denotes the number of occupied orbitals, N refers to the number of basis functions. Regards, Marcin From owner-chemistry@ccl.net Thu Oct 28 10:14:01 2010 From: "Emilio Xavier Esposito emilio.esposito=-=gmail.com" To: CCL Subject: CCL: The Anaheim ACS meeting abstract deadline: November 2nd, 2010 Message-Id: <-43027-101028101232-22052-ND5QMa9A0JucFBORaK/eFA .. server.ccl.net> X-Original-From: Emilio Xavier Esposito Content-Type: text/plain; charset=ISO-8859-1 Date: Thu, 28 Oct 2010 10:12:03 -0400 MIME-Version: 1.0 Sent to CCL by: Emilio Xavier Esposito [emilio.esposito!^!gmail.com] Hi The abstract submission deadline for the Anaheim ACS meeting is Tuesday, November 2nd, 2010. Have you submitted your abstract? Symposia include: - Computers in Chemical & Pharmaceutical Research Award Symposium Honoring Thom Dunning - Chemical Computing Group Excellence Award (Graduate students; open submission) - Hewlett-Packard Outstanding Junior Faculty Award (Open submission) - Thomas Kuhn Paradigm Shift Award (Open submission) - Challenges in Industrial Computational Methods (Invitation only) - Docking & Scoring (Invitation only) - Drug Discovery - Molecular Mechanics - Quantum Chemistry - The COMP Poster Session Abstracts can be submitted using the Program and Abstract Creation System (PACS) at http://abstracts.acs.org . Please feel free to contact us if you have questions about the COMP technical programs for the upcoming meetings or if you are interested in organizing a symposium within the COMP Division for the Fall 2011 or Spring 2012 National Meeting. We have a great technical program planned for the Anaheim meeting and we hope that you are able to attend. Sincerely Emilio Xavier Esposito, Jeffry D. Madura, and Scott Wildman The COMP Programming Board COMP Division of the ACS Follow COMP on twitter! http://twitter.com/acscomp From owner-chemistry@ccl.net Thu Oct 28 14:19:00 2010 From: "lamees hegazy lameesshams(0)yahoo.com" To: CCL Subject: CCL: capping residues for charmm force field Message-Id: <-43028-101028123028-15148-ibKYqJuzNb2h97BwLoQa4g,+,server.ccl.net> X-Original-From: lamees hegazy Content-Type: multipart/alternative; boundary="0-240090501-1288283423=:44810" Date: Thu, 28 Oct 2010 09:30:23 -0700 (PDT) MIME-Version: 1.0 Sent to CCL by: lamees hegazy [lameesshams%%yahoo.com] --0-240090501-1288283423=:44810 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable I am working on developing parameters for a ligand connected to my protein.= =A0The ligand is=A0oxalyl-aspartyl anhydride=A0=A0 O-CO-CO-O-CO-ASP. As the= first step for optimization, I wil cut between O and CO-=A0ASP and add cap= ping residue NH-CH3 so=A0 I will parameterize this ligand O-CO-CO-O-NH-CH3.= I woud like to do the parameterization accurately and I appreciate your su= ggessions.=0A=0A=0A --0-240090501-1288283423=:44810 Content-Type: text/html; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable <= /tr>
I am working on developing parameters fo= r a ligand connected to my protein. The ligand is oxalyl-aspartyl= anhydride   O-CO-CO-O-CO-ASP. As the first step for optimization= , I wil cut between O and CO- ASP and add capping residue NH-CH3 so&nb= sp; I will parameterize this ligand O-CO-CO-O-NH-CH3. I woud like to do the= parameterization accurately and I appreciate your suggessions.

=0A=0A --0-240090501-1288283423=:44810-- From owner-chemistry@ccl.net Thu Oct 28 14:54:00 2010 From: "Hans De Winter hans.dewinter/a\silicos.com" To: CCL Subject: CCL:G: PIRAMID 1.0.0: Shape-based alignment of molecules Message-Id: <-43029-101028090330-31168-HeHgdGz0wQpta+b9wjWuIQ^server.ccl.net> X-Original-From: Hans De Winter Content-Type: multipart/alternative; boundary=Apple-Mail-11-726154964 Date: Thu, 28 Oct 2010 15:03:20 +0200 Mime-Version: 1.0 (Apple Message framework v936) Sent to CCL by: Hans De Winter [hans.dewinter_._silicos.com] --Apple-Mail-11-726154964 Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes Content-Transfer-Encoding: 7bit Dear, Silicos NV is pleased to announced today the release of its new open source software tool called PIRAMID under the General Public License (GPL) v2. PIRAMID has been written on top of Open Babel's C++ API. PIRAMID is a shape-based alignment tool for the alignment of pairs of molecules. It is based on Gaussian volumes for the description molecular shape as was originally introduced by Grant and Pickup ( J. Phys. Chem. 1995, 99, 3503-3510). The software is downloadable from Silicos website ( http://www.silicos.be ) or from the Open Babel related projects webpage ( http://openbabel.org/wiki/Related_Projects ). Please contact Silicos at info[at]silicos.com for more information. Kind regards, Hans --Apple-Mail-11-726154964 Content-Type: text/html; charset=US-ASCII Content-Transfer-Encoding: 7bit Dear,

Silicos NV is pleased to announced today the release of its new open source software tool called PIRAMID under the General Public License (GPL) v2. PIRAMID has been written on top of Open Babel's C++ API.

PIRAMID is a shape-based alignment tool for the alignment of pairs of molecules. It is based on Gaussian volumes for the description molecular shape as was originally introduced by Grant and Pickup ( J. Phys. Chem. 1995, 99, 3503-3510).

The software is downloadable from Silicos website ( http://www.silicos.be ) or from the Open Babel related projects webpage ( http://openbabel.org/wiki/Related_Projects ).

Please contact Silicos at info[at]silicos.com for more information.

Kind regards,
Hans

--Apple-Mail-11-726154964-- From owner-chemistry@ccl.net Thu Oct 28 19:28:00 2010 From: "Adam Sykes admin(0)quanano.co.uk" To: CCL Subject: CCL: capping residues for charmm force field Message-Id: <-43030-101028180347-22636-w9w9wAyG2ckNsbf69kEQBQ-$-server.ccl.net> X-Original-From: Adam Sykes Content-Type: multipart/alternative; boundary=00221534d727d27f3f0493b4835b Date: Thu, 28 Oct 2010 23:03:59 +0100 MIME-Version: 1.0 Sent to CCL by: Adam Sykes [admin]~[quanano.co.uk] --00221534d727d27f3f0493b4835b Content-Type: text/plain; charset=ISO-8859-1 I would suggest taking a look at http://mackerell.umaryland.edu/Empirical_FF_Dev.html as a starting point. Best wishes, Adam Dr. Adam Sykes MChem., PhD. *Managing Director* *QuaNano Ltd* e: adam.sykes!=!quanano.co.uk m: +44 (0)7855 379 397 t: +44 (0)151 324 2866 w: http://www.quanano.co.uk *From:* owner-chemistry+adam.sykes==quanano.co.uk!=!ccl.net [mailto: owner-chemistry+adam.sykes ==quanano.co.uk!=! ccl.net] *On Behalf Of *lamees hegazy lameesshams(0)yahoo.com *Sent:* 28 October 2010 17:30 *To:* Sykes, Adam C *Subject:* CCL: capping residues for charmm force field I am working on developing parameters for a ligand connected to my protein. The ligand is oxalyl-aspartyl anhydride O-CO-CO-O-CO-ASP. As the first step for optimization, I wil cut between O and CO- ASP and add capping residue NH-CH3 so I will parameterize this ligand O-CO-CO-O-NH-CH3. I woud like to do the parameterization accurately and I appreciate your suggessions. --00221534d727d27f3f0493b4835b Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable

I would suggest taking a look at http://mackerell.umaryland.edu/Empiri= cal_FF_Dev.html as a starting point.

=A0

Best wishes,


Adam

=A0

Dr. Adam Sykes MChem., PhD.

Managing Director

QuaNano Ltd

e:=A0=A0 adam.sy= kes!=!quanano.co.uk

m: +44 (0)7855 379 397

t:=A0=A0 +44 (0)151 324 2866

w:=A0 http://www.quana= no.co.uk

=A0

From: owner-chemistry+adam.sykes=3D=3Dquanano.co= .uk!=!ccl.net [mailto:owner-chemistry+ada= m.sykes=3D=3Dquanano.co.uk!=!ccl.net] On Behalf Of lamees hegazy lameesshams(0)yahoo.com
Sent: 28 October 2010 17:30
To: Sykes, Adam C
Subject: CCL: capping residues for charmm force field

=A0

I am working on developing parameters for a ligand connected to my protein.=A0The ligand is=A0oxalyl-aspartyl anhydride=A0= =A0 O-CO-CO-O-CO-ASP. As the first step for optimization, I wil cut between O= and CO-=A0ASP and add capping residue NH-CH3 so=A0 I will parameterize this ligand O-CO-CO-O-NH-CH3. I woud like to do the parameterization accuratel= y and I appreciate your suggessions.

=A0

--00221534d727d27f3f0493b4835b--