<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Tyler,<div><br></div><div>I took a look at your input and the error below is because you have an extra curly brace after your wf command. Worse, however, is that this command is floating around all by itself - it needs to be associated with a program directive, e.g., {hf;wf,25,2,1}. As I believe Fernando just posted, the ground states of the halogen oxides are doublet-Pi states, which are 2-fold degenerate in symmetry 2 and 3 in C2v (B1 and B2). So your wf directive for ClO should be wf,25,2,1 or wf,25,3,1. In the case of BrO since you're using a 10-electron pseudopotential: wf,33,2,1 or wf,33,3,1. </div><div><br></div><div>Also the basis sets you are using are already in the Molpro basis set library:</div><div><br></div><div>In the ClO case:</div><div><br></div><div>basis=avtz</div><div><br></div><div>In the BrO case:</div><div><br></div><div>basis={default,avtz,Br=avtz-pp}</div><div><br></div><div>Last, since these are open-shell species, you'll need to use UCCSD(T) or RCCSD(T) (I personally favor uccsd(t) ). Note that both of these still use ROHF orbitals.</div><div><br></div><div>hope that helps,</div><div><br></div><div>-Kirk</div><div><br><div><div>On Apr 2, 2009, at 2:16 PM, Tyler Thompson wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Kirk,<div><br></div><div>Thanks for your help and for getting back to me so quickly. I am attaching the inputs for ClO and BrO and the output file from ClO. </div><div><br></div><div>I have looked through the manual but the examples that were provided were still unclear and I still seemed to get errors.</div> <div><br></div><div>Once again, thank you for your help. No one else in my research group is working with molpro so I am learning it on my own. Also, besides the manual do you know of any other resources that might be helpful for writing the input files for molpro?</div> <div><br></div><div>Thanks,</div><div>Tyler<br><br><div class="gmail_quote">On Thu, Apr 2, 2009 at 1:33 PM, Kirk Peterson <span dir="ltr"><<a href="mailto:kipeters@wsu.edu">kipeters@wsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"> <div style="word-wrap:break-word">Tyler,<div><br></div><div>from this error it looks like you've got a different syntax problem. Could you send along your input? In general the wf command requires all the electrons (except for any that have been replaced by an ECP). You only need a closed or open directive if the occupation is ambiguous.</div> <div><br></div><div>-Kirk</div><div><br></div><div>PS - the manual has lots of examples :^)</div><div><br><div><div><div></div><div class="h5"><div>On Apr 2, 2009, at 11:56 AM, Tyler Thompson wrote:</div><br></div></div><blockquote type="cite"> <div><div></div><div class="h5">I recently switched from Gaussian to Molpro for a number of calculations with ClO and BrO but have questions regarding the definition of the wave function. I am trying to calculate the optimized geometry and thermo data for the molecules. <div> <br></div><div>For ClO, I defined the wave function as wf, 25,1,1 and 25,1,3 but received an error stating </div><div><br></div><div> ILLEGAL COMMAND: BEGINWF</div><div> ? Error</div><div> ? Input error</div><div> ? The problem occurs in wfudef</div> <div><br></div><div>I think the problem is with the number of electrons. Should this number refer to the total number of electrons in the molecule or only to the valance electrons? Also, do I need to define open or closed shells for this species? </div> <div><br></div><div>Any help would be appreciated along with any examples that you can provide that way help me better learn this program. </div><div><br></div><div>Thanks,</div><div>Tyler Thompson</div><div><br></div><div> M.S. Energy Resources Engineering</div><div>Stanford University</div><div><a href="mailto:tylert@stanford.edu" target="_blank">tylert@stanford.edu</a></div></div></div> _______________________________________________<br> Molpro-user mailing list<br><a href="mailto:Molpro-user@molpro.net" target="_blank">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></div></div></blockquote></div><br></div> <span><clo1.inp></span><span><clo1.out></span><span><bro1.inp></span></blockquote></div><br></div></body></html>