March 28, 2024, 05:12:51 PM
Forum Rules: Read This Before Posting


Topic: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS  (Read 3637 times)

0 Members and 1 Guest are viewing this topic.

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Hey! So, I'm not sure if this is supposed to happen, as my teacher hasn't replied to my queries, but my program has been running for well over 2 hours now, and yesterday I left it running for 11 hours and got nothing.

I'm supposed to do a vibrational analysis of the molecule, but for that I need the correct geometry. The teacher gave us a few specific functionals we should be using, and for the current optimization I'm running I'm using PBE exchange + PBE correlation.


Please help, I have to run 5 more files and at +12 hours per file I won't be done in time...

Thanks


P.D.
Here's the input file I used.

 $BASIS GBASIS=N21 NGAUSS=3 $END
 $CONTRL SCFTYP=UHF RUNTYP=OPTIMIZE DFTTYP=PBE MAXIT=99 $END
 $SYSTEM MWORDS=5 $END
 $STATPT OPTTOL=0.0001 NSTEP=20 HESS=CALC $END
 $FORCE METHOD=SEMINUM VIBANL=.FALSE. $END

 $DATA
Title
C1
C     6.0     3.13984     0.61327     0.30139
C     6.0     3.45596    -0.69886    -0.03575
C     6.0     2.43944    -1.60595    -0.32763
C     6.0     1.08918    -1.22237    -0.28903
C     6.0     0.72629     0.10535     0.05056
C     6.0     1.78990     0.99942     0.34164
C     6.0    -0.66926     0.46599     0.08019
C     6.0    -1.64470    -0.51708    -0.23204
C     6.0    -1.24719    -1.81496    -0.56145
C     6.0     0.10012    -2.16249    -0.58936
C     6.0    -1.15322     1.76017     0.40510
C     6.0    -2.51803     2.06763     0.42023
C     6.0    -3.43710     1.08145     0.10889
C     6.0    -3.01262    -0.19915    -0.21420
C     6.0     4.21383     1.61069     0.62039
H     1.0     4.49121    -1.02857    -0.07523
H     1.0     2.71398    -2.62615    -0.58873
H     1.0     1.57721     2.03109     0.60893
H     1.0    -1.98725    -2.57557    -0.80203
H     1.0     0.36862    -3.18388    -0.85096
H     1.0    -0.46536     2.56266     0.65544
H     1.0    -2.83244     3.07429     0.67576
H     1.0    -3.76043    -0.95148    -0.45378
O     8.0    -4.77724     1.32932     0.10809
H     1.0    -4.92104     2.25938     0.34780
H     1.0     5.21164     1.16702     0.53966
H     1.0     4.09254     1.98085     1.64346
H     1.0     4.16441     2.45366    -0.07636
 $END

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #1 on: May 16, 2021, 03:47:20 PM »
Have you see your compound ?

All hydrogens lies on a different plane than carbon atoms !!!

Besides - you do not need UHF if molecule is not radical (and I don't see any special multiplicity).

Also - molecule should be planar so application of Cs point group will reduce calculation time.
« Last Edit: May 16, 2021, 04:00:15 PM by Orcio_Dojek »

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #2 on: May 16, 2021, 04:03:33 PM »
Have you see your compound ?

All hydrogens lies on a different plane than carbon atoms !!!

Besides - you do not need UHF if molecule is not radical (and even for radicals sometimes ROHF is good enough).

Also - molecule should be planar so application of Cs point group will reduce calculation time.


I have seen the molecule, the only hydrogen atoms that seem to be off the plane are the ones in the methyl group, how would I align those?

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #3 on: May 16, 2021, 04:13:52 PM »
Quote from: roo_sado
$DATA
Title
C1
C     6.0     3.13984     0.61327     0.30139
C     6.0     3.45596    -0.69886    -0.03575
C     6.0     2.43944    -1.60595    -0.32763
C     6.0     1.08918    -1.22237    -0.28903
C     6.0     0.72629     0.10535     0.05056
C     6.0     1.78990     0.99942     0.34164
C     6.0    -0.66926     0.46599     0.08019
C     6.0    -1.64470    -0.51708    -0.23204
C     6.0    -1.24719    -1.81496    -0.56145
C     6.0     0.10012    -2.16249    -0.58936
C     6.0    -1.15322     1.76017     0.40510
C     6.0    -2.51803     2.06763     0.42023
C     6.0    -3.43710     1.08145     0.10889
C     6.0    -3.01262    -0.19915    -0.21420
C     6.0     4.21383     1.61069     0.62039
H     1.0     4.49121    -1.02857    -0.07523
H     1.0     2.71398    -2.62615    -0.58873
H     1.0     1.57721     2.03109     0.60893
H     1.0    -1.98725    -2.57557    -0.80203
H     1.0     0.36862    -3.18388    -0.85096
H     1.0    -0.46536     2.56266     0.65544
H     1.0    -2.83244     3.07429     0.67576
H     1.0    -3.76043    -0.95148    -0.45378
O     8.0    -4.77724     1.32932     0.10809
H     1.0    -4.92104     2.25938     0.34780
H     1.0     5.21164     1.16702     0.53966
H     1.0     4.09254     1.98085     1.64346
H     1.0     4.16441     2.45366    -0.07636
 $END
You need to correct every one hydrogen + oxygen atom, not just methyl group.

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #4 on: May 16, 2021, 04:24:22 PM »

You need to correct every one hydrogen + oxygen atom, not just methyl group.


Right, I'm sorry, this is how the molecule looks right now. Are you saying I should be building the molecule aligned relative to, say XY?

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #5 on: May 16, 2021, 04:31:31 PM »
Application of Cs point group will reduce time of calculation by twice.

Application of RHF instead of UHF will also reduce time of calculation.

Still it will take many hours, as computing IR spectrum is the most time consuming.

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #6 on: May 16, 2021, 04:40:20 PM »
Application of Cs point group will reduce time of calculation by twice.

Application of RHF instead of UHF will also reduce time of calculation.

Still it will take many hours, as computing IR spectrum is the most time consuming.

Okay, I'll try to change those. What do I need to run it with the CS point group? I never learned anything about this, we always used C1...

Also I really appreciate you taking the time to reply.

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #7 on: May 16, 2021, 04:44:04 PM »
Quote
$DATA
Title
C1
Instead "C1" type "Cs".

But you will get error if molecule do not fits Cs symmetry.
« Last Edit: May 16, 2021, 05:04:48 PM by Orcio_Dojek »

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #8 on: May 16, 2021, 05:12:43 PM »
Increasing Mwords will also help (1 mword = only 8 mb of allocated memory).

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #9 on: May 16, 2021, 05:18:50 PM »
Quote
$DATA
Title
C1
Instead "C1" type "Cs".

But you will get error if molecule do not fits Cs symmetry.
I did get an error :/

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #10 on: May 16, 2021, 05:23:55 PM »
Molecule is planar so one of the coordinate (x, y, z) should be equal 0 for every atom, except of two hydrogen atoms in methyl group.

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #11 on: May 16, 2021, 05:48:31 PM »
Molecule is planar so one of the coordinate (x, y, z) should be equal 0 for every atom, except of two hydrogen atoms in methyl group.

So now it looks like this, but it still gives me the same error

 $BASIS GBASIS=N21 NGAUSS=3 $END
 $CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE DFTTYP=PBEX MAXIT=99 $END
 $SYSTEM MWORDS=1000 $END
 $STATPT OPTTOL=0.0001 NSTEP=20 HESS=CALC $END
 $FORCE METHOD=SEMINUM VIBANL=.TRUE. $END

 $DATA
Title
CS
C     6.0   -12.73450    -2.10578     0.00000
C     6.0   -12.08405    -0.88474     0.00000
C     6.0   -10.68539    -0.85494     0.00000
C     6.0    -9.89978    -2.03709     0.00000
C     6.0   -10.60561    -3.26886     0.00000
C     6.0   -12.01003    -3.28890     0.00000
C     6.0    -8.45821    -2.05304     0.00000
C     6.0    -7.78171    -3.29905     0.00000
C     6.0    -8.51118    -4.49083     0.00000
C     6.0    -9.90283    -4.47600     0.00000
C     6.0    -7.64600    -0.88824     0.00000
C     6.0    -6.24323    -0.94145     0.00000
C     6.0    -5.61576    -2.18393     0.00000
C     6.0    -6.37855    -3.34957     0.00000
C     6.0    -5.42158     0.31362     0.00000
H     1.0    -6.05237     1.20859     0.00000
H     1.0    -4.78872     0.34983    -0.89255
H     1.0    -4.78872     0.34983     0.89255
H     1.0    -8.10848     0.09483     0.00000
H     1.0    -4.53066    -2.25440     0.00000
H     1.0    -5.86237    -4.30757     0.00000
H     1.0    -8.00113    -5.45189     0.00000
H     1.0   -10.43469    -5.42536     0.00000
H     1.0   -12.55161    -4.23193     0.00000
H     1.0   -12.63468     0.05029     0.00000
H     1.0   -10.21463     0.12391     0.00000
O     8.0   -14.09478    -2.18948     0.00000
H     1.0   -14.46100    -1.29001     0.00000               
 $END

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #12 on: May 17, 2021, 09:34:43 AM »
When using point group you need to leave an empty line between the symbol of point group (here "Cs") and coordinates of first atom (C).

Plus - you need to specify type of coordinates (coord=prinaxis) in $contr group.

Offline Orcio_87

  • Full Member
  • ****
  • Posts: 440
  • Mole Snacks: +39/-3
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #13 on: May 17, 2021, 10:17:23 AM »
As of using Cs point group....

I was wondering is the molecule in fact planar (H9 hydrogen lies on the same plane or above), but somehow - planar structure is slightly less stable than structure with H9 out of plane, so yes - Cs point group is even better for purpose of calculations.
« Last Edit: May 17, 2021, 12:58:02 PM by Orcio_Dojek »

Offline roo_sado

  • New Member
  • **
  • Posts: 7
  • Mole Snacks: +0/-0
Re: Geometry Optimization of oxi-methylphenantrene using DFT in GAMESS
« Reply #14 on: May 17, 2021, 02:04:13 PM »
Thanks a lot! I was able to finish by changing the HESS to GUESS. It still took about 8-10 hours of around constant 100% CPU and 100% memory allocated usage.

Sponsored Links