[molpro-user] Problem with partial optimization

rsb20 at uakron.edu rsb20 at uakron.edu
Fri May 22 13:17:38 BST 2009


Hi molpro users;
I would like to thank you, especially Berger and Peter, for the 
suggestion to fix problem with molpro. I am still facing problem with 
the partial optimization of methylamine. For this I need to constrain 
the dihedral angle. I have been able to do that in g03 but when I tried 
on molpro it is showing following errors:

Input
***,ch3nh2
r1=1.54, r2=1.069, r3=1.071, r4=1.07, r5=1.07, r6=1.07, a1=109.708, 
a2=108.138,
a3=109.87, a4=108.33, a5=108.95, d1=180, d2=60.0, d3=60.0, d4=180.0
geometry={angstrom,
           C;
           N, C, r1;
           H1, N, r2, C, a1;
          H2, N, r3, C, a2, H1, -d1;
          H3, C, r4, N, a3, H2, d2;
          H4, C, r5, N, a4, H2, -d3;
          H5, C, r6, N, a5, H2, -d4}
cartesian
basis=6-31g**
hf
mp2
optg
constraint,-170,deg,dihedral,atoms=[H2,N,C,H1];

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

The error in output :
?error
?internal coordinator/constraint is allowed only for slapof optimizer
?the problem occurs in opt_inp1
----------------------------------------------------------------------------

Then I tried to do alternatively as
input file
***,ch3nh2
geometry={
C     -0.000583670    0.001974822   -1.479670491
N      0.001542237   -0.000651370    1.430505570
H1    -0.005609368    1.900503337    2.113464323
H2     0.009738396   -1.924531645    2.058816552
H3     1.649583947   -0.941586470   -2.168982072
H4    -1.659433862   -0.963836160   -2.115230653
H5     -0.008755367    1.914969870   -2.134571020
}
basis=6-31g**
hf
mp2
optg
  optg,savexyz=ch3nh2.xyz
  coord,bmat                  !Optimization in natural internal coordinates
  method,diis                 !Optimization method: Geometry DIIS

  optg,coord=bmat,method=diis
constraint,-170,deg,dihedral,atoms=[H2,N,C,H1];

--=============================
The error in output is
  ? Error
  ? Too few geometry parameters given on zmatrix line
  ? The problem occurs in zmatrix

I would appreciate your help for the suggestion to fix this error.
Thank you
Ram




More information about the Molpro-user mailing list