(fwd) CCSD(T) geometry optimization with xyz coordinates

Ludger Palm Ludger.Palm at lrz-muenchen.de
Thu Jul 18 13:32:09 BST 2002


Thanks a lot for the quick reply.

Does this mean that geometry optimization for e.g. CCSD is only possible 
using z-matrices?

I tried a nonlinear molecule, methane:

***,methan
memory,20,m
geomtyp=xyz
geometry={angstrom
5
methan
  H    0.000000    0.000000   -1.083658
  C    0.000000    0.000000    0.000000
  H    0.000000    1.021683    0.361219
  H    0.884803   -0.510841    0.361219
  H   -0.884803   -0.510841    0.361219
}
basis=cc-pVTZ
hf
ccsd
optg

and Molpro gave me the message

 Permutation not found: symmetry=           1   i=          28  
label=2px     

bflab:                                                                        
                                                                                
                                    32py    32pz            32py   
32pz        
     33d0    33d2+   33d1-                   51s     51s     51s    
51s     52p
 y    52pz    52py    52pz    52py    52pz    53d0    53d2+   53d1-  
53d0    53
 d2+   53d1-   54f1-   54f0    54f3-   54f2+ 

 ERROR EXIT

"H. -J. Werner" wrote:
> 
> This problem occurs since the program tries to do a numerical optimization
> in 3N coordinates (wich is of course nonsense in this case). It increments
> an x coordinate by some small value, and then the symmetry is reduced, and
> the error occurs. It would certainly work without symmetry, but this would
> be a big wast of time. You should use a z-matrix input which contains in this
> case just one parameter (the bond distance) to be optimized.
> H.-J. Werner
> ----- Forwarded message from Ludger Palm <Ludger.Palm at lrz-muenchen.de> -----
> 
> From: Ludger Palm <Ludger.Palm at lrz-muenchen.de>
> To: molpro-user at stchem.bham.ac.uk
> Subject: CCSD(T) geometry optimization with xyz coordinates
> Date: Thu, 18 Jul 2002 14:02:03 +0200
> Message-Id: <3D36AE3B.20D74571 at lrz-muenchen.de>
> X-Mailer: Mozilla 4.77 [de] (X11; U; Linux 2.4.6 i686)
> Organization: Leibniz-Rechenzentrum
> 
> Dear Molpro users,
> 
> I try to do a geometry optimization of a molecule where the
> geometry is given in xyz coordinates. Here is a little example:
> 
> ***,co
> memory,20,m
> geomtyp=xyz
> geometry={angstrom
> 2
> carbonmonoxide
>  O   0.000000000   0.000000000  -0.569524962
>  C   0.000000000   0.000000000   0.569524962
> }
> basis=cc-pVTZ
> hf
> ccsd
> optg
> 
> Molpro stops the execution during the ccsd geometry optimization
> with the message
> 
>  Permutation not found: symmetry=           1   i=           5
> label=2px
>  bflab:                                  12pz    12pz    12pz    13d0
> 13d2+
>    13d0    13d2+   14f0    14f2+   21s     21s     21s     21s
> 22pz    22pz
>     22pz    23d0    23d2+   23d0    23d2+   24f0    24f2+
> 
>  ERROR EXIT
> 
> What's going on?
> Should Molpro be able to do such a geometry optimization?
> 
> --
> With best regards
> 
> Ludger Palm
> 
>  Dr. Ludger Palm                           | Phone: +(0)89-289-28792
>  Leibniz-Rechenzentrum der                 | Fax:   +(0)89-2809460
>  Bayerischen Akademie der Wissenschaften   | Email: palm at lrz.de
>  - High Performance Computing Group -      |
>  Barer Strasse 21                          |
>  80333 Muenchen                            | WWW: http://www.lrz.de
>  Germany                                   | Hotline: +(0)89-289-28800
> 
> ----- End forwarded message -----
> 
> --
> Prof. Hans-Joachim Werner
> Institute for Theoretical Chemistry
> University of Stuttgart
> Pfaffenwaldring 55
> D-70569 Stuttgart, Germany
> Tel.: (0049) 711 / 685 4400
> Fax.: (0049) 711 / 685 4442
> e-mail: werner at theochem.uni-stuttgart.de

-- 
With best regards
Mit freundlichen Gruessen

Ludger Palm


 Dr. Ludger Palm                           | Phone: +(0)89-289-28792
 Leibniz-Rechenzentrum der                 | Fax:   +(0)89-2809460     
 Bayerischen Akademie der Wissenschaften   | Email: palm at lrz.de       
 - High Performance Computing Group -      |                           
 Barer Strasse 21                          |                           
 80333 Muenchen                            | WWW: http://www.lrz.de    
 Germany                                   | Hotline: +(0)89-289-28800



More information about the Molpro-user mailing list