Dear Tereza,<br>May be you have a symmetry-breaking error. OPTG move
the starting geommetry along the energy gradient. If you specify
symmetry to start and the atoms begin to move, the new set of integrals
have to be calculated at a lower symmetry and the program can't figure
out what to do. use "nosym" in geometry card, if the optimization
procedure work well you will get a nearly symmetric structure...<br><br><br><div class="gmail_quote">On Mon, Nov 2, 2009 at 2:39 PM, tereza sedivcova <span dir="ltr"><<a href="mailto:tereza.sedivcova@gmail.com">tereza.sedivcova@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Tim Lillestolen and others MOLPRO users,<br>
did you solve your problem?<br>
Because I have the same now, and I can not see a bug (bugs) in my Z-matrix.<br>
<br>
I am interested in FSO3 radical in its C3v symmery, however even if I tryed it<br>
for Cs or C1 symmetry I obtain the same error.<br>
The inputs for all three geometries are below.<br>
<br>
Thank you very very much for your help,<br>
Tereza<br>
**************************************************************<br>
Laboratory of High Resolution Molecular Spectroscopy<br>
Department of Analytical Chemistry<br>
Institute of Chemical Technology, Prague<br>
**************************************************************<br>
<br>
input:<br>
***C3v symmetry point group<br>
<br>
!<br>
memory,500,m<br>
gprint,orbital,civector<br>
geometry={<br>
Angstrom;<br>
f;<br>
s,f,fs;<br>
o1,s,so,f,oso;<br>
o2,s,so,f,oso,o1,D;<br>
o3,s,so,f,oso,o2,D;<br>
}<br>
<br>
fs = 1.542;<br>
so = 1.443;<br>
oso = 105.404;<br>
D = 120.000;<br>
<br>
basis=vdz;<br>
int;<br>
{uhf; occ, 17,8; wf,49,2,1;}<br>
<br>
{multi; occ, 19,9; closed,15,6; wf,49,2,1;state,1;<br>
maxiter,100;}<br>
<br>
{optg;}<br>
forces;<br>
<br>
put, molden, 1.molden;<br>
<br>
&&&&&&&&&&&&&&&&&&&&&&&&&&&&<br>
***Cs<br>
memory,500,m<br>
gprint,orbital,civector<br>
geometry={<br>
Angstrom;<br>
s;<br>
f,s,fs;<br>
o1,s,so1,f,oso1;<br>
o2,s,so1,f,oso1,o1,D1;<br>
o3,s,so2,f,oso2,o2,D2;<br>
}<br>
<br>
fs = 1.542;<br>
so1 = 1.443;<br>
so2 = 1.443;<br>
oso1 = 105.404;<br>
oso2 = 105.404;<br>
D1 = 120.000;<br>
D2 = 120.000;<br>
<br>
basis=vdz;<br>
int;<br>
{uhf; occ, 17,8; wf,49,2,1;}<br>
<br>
{multi; occ, 19,9; closed,15,6; wf,49,2,1;state,1;<br>
maxiter,100;}<br>
<br>
{optg;}<br>
forces;<br>
<br>
put, molden, 1.molden;<br>
<br>
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&<br>
***,C1<br>
<br>
memory,500,m<br>
gprint,orbital,civector<br>
geometry={<br>
Angstrom;<br>
s;<br>
f,s,fs;<br>
o1,s,so1,f,oso1;<br>
o2,s,so2,f,oso2,o1,D1;<br>
o3,s,so3,f,oso3,o2,D2;<br>
}<br>
<br>
fs = 1.542;<br>
so1 = 1.443;<br>
so2 = 1.443;<br>
so3 = 1.443;<br>
oso1 = 105.404;<br>
oso2 = 105.404;<br>
oso3 = 105.404;<br>
D1 = 120.000;<br>
D2 = 120.000;<br>
<br>
basis=vdz;<br>
int;<br>
{uhf; occ, 17,8; wf,49,2,1;}<br>
<br>
{multi; occ, 19,9; closed,15,6; wf,49,2,1;state,1;<br>
maxiter,100;}<br>
<br>
{optg;}<br>
forces;<br>
<br>
put, molden, 1.molden;<br>
<br>
<br>
<br>
<br>
<br>
On Wed, Nov 5, 2008 at 12:24 PM, lillestolen<br>
<<a href="mailto:timothy.lillestolen@nottingham.ac.uk">timothy.lillestolen@nottingham.ac.uk</a>> wrote:<br>
> Dear Molpro users,<br>
><br>
> I feel like I'm probably missing something simple here, but I still<br>
> can't figure out what I am doing wrong. I am trying to optimise the<br>
> geometry of carbonic acid (H2CO3) with C2V symmetry at the<br>
> B3LYP/aug-cc-pVDZ level of theory. I am including the input file at<br>
> the end of the message. The first energy evaluation works just fine,<br>
> and the program detects the correct C2V symmetry. However, after the<br>
> first displacements are calculated I get the following error:<br>
><br>
><br>
> Optimization point 1<br>
><br>
> Variable Last Current<br>
> Next Gradient Hessian<br>
><br>
> E(KS-SCF000) / Hartree 0.00000000 -264.87588720 0.00000000<br>
> OC1 / ANGSTROM 0.00000000 1.22200000 1.28411786<br>
> -0.25660683 4.00134221<br>
> OC2 / ANGSTROM 0.00000000 1.22200000 1.25951398<br>
> -0.32956964 8.18713291<br>
> OCO / DEGREE 0.00000000 120.00000000<br>
> 119.07590136 0.00134251 0.00164275<br>
> DIH4 / DEGREE 0.00000000 180.00000000<br>
> 180.00000000 0.00000000 0.00000000<br>
> HO / ANGSTROM 0.00000000 0.88000000 0.94600782<br>
> -0.41770827 6.28691438<br>
> HOC / DEGREE 0.00000000 120.00000000<br>
> 117.35334871 0.00102677 0.00046113<br>
> DIH5 / DEGREE 0.00000000 0.00000000<br>
> 0.00000000 0.00000000 0.00000000<br>
> DIH6 / DEGREE 0.00000000 180.00000000<br>
> 180.00000000 0.00000000 0.00000000<br>
> Convergence: 0.00000000 (line search)<br>
> 0.19172405 0.32724871 (total)<br>
> The request for symmetry elements X,Y, cannot be honoured<br>
> ? Error<br>
> ? Symmetry problem<br>
> ? The problem occurs in zmatrix.f:zmat_orient<br>
><br>
> ERROR EXIT<br>
> CURRENT STACK: MAIN<br>
><br>
><br>
> **********************************************************************************************************************************<br>
> The request for symmetry elements X,Y, cannot be honoured<br>
><br>
> Prototype coordinates<br>
> 1 2 3<br>
> 4 5 6<br>
> 1 0.0000000 0.0000000 0.0000000 0.0000000<br>
> 0.0000000 0.0000000<br>
> 2 -0.0240505 1.2534949 1.1355597 -2.4015615<br>
> 0.1491318 0.3359890<br>
> 3 -0.0252998 2.0378084 -2.1038457 0.0864578<br>
> -3.5947566 3.5720959<br>
> ? Error<br>
> ? Symmetry problem<br>
> ? The problem occurs in zmatrix.f:zmat_orient<br>
><br>
> RECURSIVE CALL OF FEHLER. EMERGENCY STOP.<br>
> FILE INFORMATION MAY BE LOST<br>
><br>
> RECURSIVE CALL OF ERROR HANDLER, FILE INFORMATION MAY BE LOST<br>
><br>
> EMERGENCY STOP<br>
><br>
><br>
> I get the same message if I try and do the optimisation with hartree<br>
> fock (replacing ks,b3lyp with hf). Am I doing something incorrectly<br>
> with my z-matrix? Or is there another option that I am missing out on?<br>
> I really need to get the final geometry with the correct symmetry as I<br>
> am calculating some atomic properties using these geometries and I need<br>
> symmetric atoms. Much thanks in advance for any advice/help,<br>
><br>
> best regards,<br>
><br>
> Tim Lillestolen<br>
> Wheatley Group<br>
> University of Nottingham<br>
><br>
><br>
> input file:<br>
><br>
> ***,h2co3<br>
><br>
> basis=avdz<br>
><br>
> geometry={<br>
> angstrom;<br>
> c;<br>
> o1, c,oc1;<br>
> o2, c,oc2,o1,oco;<br>
> o3, c,oc2,o1,oco,o2,dih4;<br>
> h1,o2,ho,c ,hoc,o3,dih5;<br>
> h2,o1,ho,c ,hoc,o2,dih6;<br>
> }<br>
><br>
> oc1 = 1.222;<br>
> oc2 = 1.222;<br>
> ho = 0.88;<br>
><br>
><br>
> oco= 120.000;<br>
> hoc= 120.000;<br>
><br>
> dih4= 180.000;<br>
> dih5= 0.000;<br>
> dih6= 180.000;<br>
><br>
> ks,b3lyp;<br>
> optg,symmetry=auto,coord=zmat;<br>
><br>
><br>
> This message has been checked for viruses but the contents of an attachment<br>
> may still contain software viruses, which could damage your computer system:<br>
> you are advised to perform your own checks. Email communications with the<br>
> University of Nottingham may be monitored as permitted by UK legislation.<br>
><br>
> _______________________________________________<br>
> Molpro-user mailing list<br>
> <a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
> <a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
><br>
_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
</blockquote></div><br>