From chemistry-request@ccl.net Fri Jul  8 07:34:45 2005
Received: from server.ccl.net (ccl [127.0.0.1])
	by server.ccl.net (8.13.1/8.13.1) with ESMTP id j68BYi8O026034
	for <chemistry _()_ ccl.net>; Fri, 8 Jul 2005 07:34:44 -0400
Received: (from apache@localhost)
	by server.ccl.net (8.13.1/8.13.1/Submit) id j68BYinl026033
	for chemistry _()_ ccl.net; Fri, 8 Jul 2005 07:34:44 -0400
Date: Fri, 8 Jul 2005 07:34:44 -0400
Message-Id: <200507081134.j68BYinl026033 _()_ server.ccl.net>
X-Authentication-Warning: server.ccl.net: apache set sender to chemistry-request _()_ ccl.net using -f
From: "Gustavo, L.C., Moura" <gustavo _()_ mercury.chem.pitt.edu>
To: chemistry _()_ ccl.net
X-Web-Message-Number: 050708073319-25962
Subject: W:Symmetry breaking during CASSCF 
X-Spam-Status: No, score=-2.8 required=5.0 tests=ALL_TRUSTED autolearn=failed 
	version=3.0.4
X-Spam-Checker-Version: SpamAssassin 3.0.4 (2005-06-05) on server.ccl.net

Dear CCL Readers,
     I have a planar conjugated organic molecule that, on paper, has D2h
symmetry. I am working with the hypothesis that the molecule is a
biradical singlet. I want to perform a CASSCF(2,2) on this molecule
to estimate the biradical character of this molecule. When I
calculate the molecule with the constraint that the molecule has D2h
symmetry, the calculation has severe convergence problems. If I
remove this constraint, the molecule distorts to a C2v symmetry and I
obtain a biradical character ~20%. I checked the orbitals being used
for the CASSCF calculation and they look right (one nonbonding and
another slightly antibonding).
     Is the situation described above correct? Can an organic molecule
have such a symmetry breaking as a result of a CASSCF calculation? Is
the molecule really a C2v system? Any suggestions are welcome.
     Thank you very much in advance for your help.
     Sincerely yours,
                                Gustavo L.C. Moura
                           gustavo _()_ mercury.chem.pitt.edu