From owner-chemistry@ccl.net Wed Jul 6 01:38:01 2016 From: "Sharvan Kumar shravan.aug12#gmail.com" To: CCL Subject: CCL: Help Regarding NICS-XY scan plot Message-Id: <-52279-160706013539-26518-WKsu4QYBbsLppWiWiPNiFQ[]server.ccl.net> X-Original-From: "Sharvan Kumar" Date: Wed, 6 Jul 2016 01:35:38 -0400 Sent to CCL by: "Sharvan Kumar" [shravan.aug12 a gmail.com] Dear Sir, I am sharvan, a PhD student from JNU, New Delhi, India. I did the NICS-XY of calculation s-Indacene(DOI: 10.1039/c6sc00950f) and other reported molecules NICS in this paper by aroma software, when I tried to plot it, I have seen the plots are looking same in some what but not exactly matching. The difference is that they have shown + and - value and the highest peaks are at 0, -1 and 1 but I am not finding such. In my case it starts from 0. Please tell me what should I have to do to get the same plot as reported so I can repeat and do the same for my molecule. Thanks regards, Sharvan From owner-chemistry@ccl.net Wed Jul 6 12:02:00 2016 From: "Pierre Archirel pierre.archirel*_*u-psud.fr" To: CCL Subject: CCL:G: G gaussian calculation with freq=anharm Message-Id: <-52280-160706120036-32513-opXYHdAC9/cnXQr+42wzcA[-]server.ccl.net> X-Original-From: "Pierre Archirel" Date: Wed, 6 Jul 2016 12:00:34 -0400 Sent to CCL by: "Pierre Archirel" [pierre.archirel[]u-psud.fr] Dear colleagues, I am using gaussian 09 RevD01 with freq=anharm on a rather large molecule (36 atoms) and the calculation has been now running for 150 hours on an 8 core machine. I am wondering if the total cpu time can be evaluated. I have already got 100 DFT energies, is it possible to know the total number of DFT calculations, to be done? Is it related to the 3n-6 degrees of freedom of the molecule? Many thanks in advnce, Pierre Archirel LCP, Orsay, France pierre.archirel//u-psud.fr From owner-chemistry@ccl.net Wed Jul 6 13:50:01 2016 From: "Brian Skinn bskinn,+,alum.mit.edu" To: CCL Subject: CCL:G: G gaussian calculation with freq=anharm Message-Id: <-52281-160706134839-5145-gZsqE3CAJG5HlSDknGdafA-*-server.ccl.net> X-Original-From: Brian Skinn Content-Type: multipart/alternative; boundary=001a1141afe2a7aefd0536fb2d82 Date: Wed, 6 Jul 2016 13:48:14 -0400 MIME-Version: 1.0 Sent to CCL by: Brian Skinn [bskinn!A!alum.mit.edu] --001a1141afe2a7aefd0536fb2d82 Content-Type: text/plain; charset=UTF-8 Pierre, Assuming you're using the VPT2 method built into G09, a total of 6N-11 Hessian computations are required for a full anharmonic treatment of a nonlinear system. This would be 205 Hessians for your 36-atom molecule. The 6N-11 comes about due to the use of a two-sided finite difference approach to find the derivatives of the force constants with respect to displacements along each of the normal modes, plus one calculation for the unperturbed geometry: 2*(3N-6)+1 = 6N-12+1 = 6N-11. -Brian On Wed, Jul 6, 2016 at 12:00 PM, Pierre Archirel pierre.archirel*_*u-psud.fr wrote: > > Sent to CCL by: "Pierre Archirel" [pierre.archirel[]u-psud.fr] > Dear colleagues, > I am using gaussian 09 RevD01 with freq=anharm on a rather large molecule > (36 atoms) and the calculation has been now running for 150 hours on an 8 > core machine. I am wondering if the total cpu time can be evaluated. I have > already got 100 DFT energies, is it possible to know the total number of > DFT calculations, to be done? Is it related to the 3n-6 degrees of freedom > of the molecule? > Many thanks in advnce, > Pierre Archirel > LCP, Orsay, France > pierre.archirel//u-psud.fr> > > --001a1141afe2a7aefd0536fb2d82 Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable
Pierre,

Assuming you'= ;re using the VPT2 method built into G09, a total of 6N-11 Hessian computat= ions are required for a full anharmonic treatment of a nonlinear system. Th= is would be 205 Hessians for your 36-atom molecule.

The 6N-11 comes about due to the use of a two-sided finite difference app= roach to find the derivatives of the force constants with respect to displa= cements along each of the normal modes, plus one calculation for the unpert= urbed geometry: 2*(3N-6)+1 =3D 6N-12+1 =3D 6N-11.

=
-Brian


On Wed, Jul 6, 2016 at 12:00 PM, Pierre Archirel p= ierre.archirel*_*u-psud.fr <owner-c= hemistry : ccl.net> wrote:

Sent to CCL by: "Pierre=C2=A0 Archirel" [pierre.archirel[]u-psud.fr] Dear colleagues,
I am using gaussian 09 RevD01 with freq=3Danharm on a rather large molecule=
(36 atoms) and the calculation has been now running for 150 hours on an 8 c= ore machine. I am wondering if the total cpu time can be evaluated. I have = already got 100 DFT energies, is it possible to know the total number of DF= T calculations, to be done? Is it related to the 3n-6 degrees of freedom of= the molecule?
Many thanks in advnce,
Pierre Archirel
LCP, Orsay, France
pierre.archirel//u-psud.fr



-=3D This is automatically added to each message by the mailing script =3D-=
E-mail to subscribers: CHEMISTRY : ccl.n= et or use:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/ccl/s= end_ccl_message

E-mail to administrators: CHEM= ISTRY-REQUEST : ccl.net or use
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/ccl/s= end_ccl_message

Subscribe/Unsubscribe:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/chemistry/sub_un= sub.shtml

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

Job: http://www.ccl.net/jobs
Conferences: http://server.ccl.net/chemist= ry/announcements/conferences/

Search Messages: http://www.ccl.net/chemistry/sear= chccl/index.shtml
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/spammers.txt

RTFI: http://www.ccl.net/chemistry/aboutccl/ins= tructions/



--001a1141afe2a7aefd0536fb2d82-- From owner-chemistry@ccl.net Wed Jul 6 14:25:01 2016 From: "Guzman, Francisco fg58[#]njit.edu" To: CCL Subject: CCL:G: G gaussian calculation with freq=anharm Message-Id: <-52282-160706135151-5935-UTywMtgBiJPETAY3q384VA ~ server.ccl.net> X-Original-From: "Guzman, Francisco" Content-Type: multipart/alternative; boundary=001a1141d9b611c9660536fb39aa Date: Wed, 6 Jul 2016 13:51:45 -0400 MIME-Version: 1.0 Sent to CCL by: "Guzman, Francisco" [fg58!^!njit.edu] --001a1141d9b611c9660536fb39aa Content-Type: text/plain; charset=UTF-8 Pierre, What method/basis set? What other options did you run the job with? Are you optimizing as well or did you already do this? The number of cycles depends on some of your options used. I think default number of cycles is 128 or 512 if you used scf=qc. If you have the Firsch and Foresman book (Exploring Chemistry with Electronic structure) I recall seeing some info on this in the introduction. On Wed, Jul 6, 2016 at 12:00 PM, Pierre Archirel pierre.archirel*_*u-psud.fr wrote: > > Sent to CCL by: "Pierre Archirel" [pierre.archirel[]u-psud.fr] > Dear colleagues, > I am using gaussian 09 RevD01 with freq=anharm on a rather large molecule > (36 atoms) and the calculation has been now running for 150 hours on an 8 > core machine. I am wondering if the total cpu time can be evaluated. I have > already got 100 DFT energies, is it possible to know the total number of > DFT calculations, to be done? Is it related to the 3n-6 degrees of freedom > of the molecule? > Many thanks in advnce, > Pierre Archirel > LCP, Orsay, France > pierre.archirel//u-psud.fr> > > -- Francisco Guzman PhD Candidate Otto H. York Department of Chemical, Biological and Pharmaceutical Engineering New Jersey Institute of Technology Email: fg58%njit.edu, guzman.research%gmail.com Cell: 347-806-3977 --001a1141d9b611c9660536fb39aa Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable
Pierre,

What method/basis set? What= other options did you run the job with? Are you optimizing as well or did = you already do this?=C2=A0 The number of cycles depends on some of your opt= ions used. I think default number of cycles is 128 or 512 if you used scf= =3Dqc.

If you have the Firsch and Foresman book (Exploring Che= mistry with Electronic structure) I recall seeing some info on this in the = introduction.

On Wed, Jul 6, 2016 at 12:00 PM, Pierre Archirel pierre.archirel*_*u-psud.fr <owner-chemistry%ccl.net> wrote:

Sent to CCL by: "Pierre=C2=A0 Archirel" [pierre.archirel[]
u-psud.fr] Dear colleagues,
I am using gaussian 09 RevD01 with freq=3Danharm on a rather large molecule=
(36 atoms) and the calculation has been now running for 150 hours on an 8 c= ore machine. I am wondering if the total cpu time can be evaluated. I have = already got 100 DFT energies, is it possible to know the total number of DF= T calculations, to be done? Is it related to the 3n-6 degrees of freedom of= the molecule?
Many thanks in advnce,
Pierre Archirel
LCP, Orsay, France
pierre.archirel//u-psud.fr



-=3D This is automatically added to each message by the mailing script =3D-=
E-mail to subscribers: CHEMISTRY%ccl.n= et or use:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/ccl/s= end_ccl_message

E-mail to administrators: CHEM= ISTRY-REQUEST%ccl.net or use
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/ccl/s= end_ccl_message

Subscribe/Unsubscribe:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/chemistry/sub_un= sub.shtml

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

Job: http://www.ccl.net/jobs
Conferences: http://server.ccl.net/chemist= ry/announcements/conferences/

Search Messages: http://www.ccl.net/chemistry/sear= chccl/index.shtml
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/spammers.txt

RTFI: http://www.ccl.net/chemistry/aboutccl/ins= tructions/





--
Francisco Guzman
= PhD Candidate
Otto H. York Department of Chemical, Biological and = Pharmaceutical Engineering
New Jersey Institute of Technology
Email:=C2=A0fg58%njit.edu,=C2=A0guzman.research%gmail.com
Cell:=C2=A0347-806-3977
--001a1141d9b611c9660536fb39aa-- From owner-chemistry@ccl.net Wed Jul 6 14:59:01 2016 From: "Close, David M. CLOSED .. mail.etsu.edu" To: CCL Subject: CCL:G: G gaussian calculation with freq=anharm Message-Id: <-52283-160706145501-16483-om35ALVA4IZhmzWWEZWnJg%x%server.ccl.net> X-Original-From: "Close, David M." Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="us-ascii" Date: Wed, 6 Jul 2016 18:54:52 +0000 MIME-Version: 1.0 Sent to CCL by: "Close, David M." [CLOSED%a%mail.etsu.edu] Pierre: There is not information here to guess what's going on. It would be best to animate the steps to see what the problem is. Often after 15-20 iterations you are close to convergence. But then something happens. For example a small oscillation sets in involving two nearby energies that could result from a small rotation back and forth. Regards, Dave Close -----Original Message----- > From: owner-chemistry+closed==etsu.edu:ccl.net [mailto:owner-chemistry+closed==etsu.edu:ccl.net] On Behalf Of Pierre Archirel pierre.archirel*_*u-psud.fr Sent: Wednesday, July 06, 2016 12:01 PM To: Close, David M. Subject: CCL:G: G gaussian calculation with freq=anharm Sent to CCL by: "Pierre Archirel" [pierre.archirel[]u-psud.fr] Dear colleagues, I am using gaussian 09 RevD01 with freq=anharm on a rather large molecule (36 atoms) and the calculation has been now running for 150 hours on an 8 core machine. I am wondering if the total cpu time can be evaluated. I have already got 100 DFT energies, is it possible to know the total number of DFT calculations, to be done? Is it related to the 3n-6 degrees of freedom of the molecule? Many thanks in advnce, Pierre Archirel LCP, Orsay, France pierre.archirel//u-psud.frhttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt From owner-chemistry@ccl.net Wed Jul 6 18:31:01 2016 From: "Jens Spanget-Larsen spanget:+:ruc.dk" To: CCL Subject: CCL:G: SV: CCL:G: G gaussian calculation with freq=anharm Message-Id: <-52284-160706181100-7880-3hQ1sqbLGWEkS64JajNzDQ]|[server.ccl.net> X-Original-From: Jens Spanget-Larsen Content-Language: da-DK Content-Type: multipart/alternative; boundary="_000_A94E15A372E6194CA8719D62642F67448BA2FB5DMBX1adrucdk_" Date: Wed, 6 Jul 2016 22:10:54 +0000 MIME-Version: 1.0 Sent to CCL by: Jens Spanget-Larsen [spanget[a]ruc.dk] --_000_A94E15A372E6194CA8719D62642F67448BA2FB5DMBX1adrucdk_ Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hello Pierre! Gaussian freq=3Danharmo calculations are simply extremely time-consuming. A= few years ago we performed anharmonic analyses for a large number of syste= ms using B3LYP DFT and the modest basis set 6-31G* (see below). The largest= molecule we treated was usnic acid with 41 atoms; the calculation ran for = more than 54 days on a 4 core machine. It completed successfully. Yours, Jens >--< J. Spanget-Larsen, B. K. V. Hansen, P. E. Hansen: "OH stretching frequencie= s in systems with intramolecular hydrogen bonds. Harmonic and anharmonic an= alyses", Chem. Phys. 389, 107-115 (2011). DOI: 10.1016/j.chemphys.2011.09.0= 11. ------------------------------------------------------ JENS SPANGET-LARSEN Office: +45 4674 2710 Dept. of Science Fax: +45 4674 3011 Roskilde University Mobile: +45 2320 6246 P.O.Box 260 E-Mail: spanget|a|ruc.dk DK-4000 Roskilde, Denmark http://www.ruc.dk/~spanget ------------------------------------------------------ ________________________________________ Fra: owner-chemistry+spanget=3D=3Druc.dk|a|ccl.net [owner-chemistry+spanget= =3D=3Druc.dk|a|ccl.net] på vegne af Pierre Archirel pierre.archirel*_*u-= psud.fr [owner-chemistry|a|ccl.net] Sendt: 6. juli 2016 18:00 Til: Jens Spanget-Larsen Emne: CCL:G: G gaussian calculation with freq=3Danharm Sent to CCL by: "Pierre Archirel" [pierre.archirel[]u-psud.fr] Dear colleagues, I am using gaussian 09 RevD01 with freq=3Danharm on a rather large molecule (36 atoms) and the calculation has been now running for 150 hours on an 8 c= ore machine. I am wondering if the total cpu time can be evaluated. I have = already got 100 DFT energies, is it possible to know the total number of DF= T calculations, to be done? Is it related to the 3n-6 degrees of freedom of= the molecule? Many thanks in advnce, Pierre Archirel LCP, Orsay, France pierre.archirel//u-psud.fr -=3D This is automatically added to each message by the mailing script =3D-http://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt--_000_A94E15A372E6194CA8719D62642F67448BA2FB5DMBX1adrucdk_ Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable


Hello Pierre!

 

Gaussian freq=3Danharmo calculatio= ns are simply extremely time-consuming. A few years ago we performed anharm= onic analyses for a large number of systems using B3LYP DFT and the mo= dest basis set 6-31G* (see below). The largest molecule we treated was usnic acid with 41 atoms; the calculation ran for = more than 54 days on a 4 core machine. It completed successfully.

 

Yours, Jens >--<

 

J. Spanget-Larsen, B. K. V. Hansen= , P. E. Hansen: "OH stretching frequencies in systems with intramolecu= lar hydrogen bonds. Harmonic and anharmonic analyses", Chem. Phys. 389, 107-115 (2011). DOI: 10.1016/j.chemphys.2011.09.011.


  ------------------------------= ------------------------
  JENS SPANGET-LARSEN         = Office:      +45 4674 2710
  Dept. of Science        &nbs= p;   Fax:         +45= 4674 3011
  Roskilde University         = Mobile:      +45 2320 6246
  P.O.Box 260         &nb= sp;       E-Mail:     spanget|a|ruc.dk
  DK-4000 Roskilde, Denmark   http://www.ruc.dk/~spanget
  ------------------------------------------------------



________________________________________
Fra: owner-chemistry+spanget=3D=3Druc.dk|a|ccl.net [owner-chemistry+s= panget=3D=3Druc.dk|a|ccl.net] p&#229; vegne af Pierre Archirel pierre.arc= hirel*_*u-psud.fr [owner-chemistry|a|ccl.net]
Sendt: 6. juli 2016 18:00
Til: Jens Spanget-Larsen
Emne: CCL:G: G gaussian calculation with freq=3Danharm

Sent to CCL by: "Pierre  Archirel" [pierre.archirel[]u-psud.= fr]
Dear colleagues,
I am using gaussian 09 RevD01 with freq=3Danharm on a rather large molecule=
(36 atoms) and the calculation has been now running for 150 hours on an 8 c= ore machine. I am wondering if the total cpu time can be evaluated. I have = already got 100 DFT energies, is it possible to know the total number of DF= T calculations, to be done? Is it related to the 3n-6 degrees of freedom of the molecule?
Many thanks in advnce,
Pierre Archirel
LCP, Orsay, France
pierre.archirel//u-psud.fr



-=3D This is automatically added to each message by the mailing script =3D-=       http://www.ccl.net/cgi-bin/ccl/send_ccl_mess= age
      http://www.ccl.net/cgi-bin/ccl/send_ccl_mess= age

Subscribe/Unsubscribe:
     =

Job: http://www.ccl.net/jobs      

--_000_A94E15A372E6194CA8719D62642F67448BA2FB5DMBX1adrucdk_--