Hello developers and Daltonusers,
I want to repeat a simple scfcalculation on the Catom, but it is not clear to me, how one adjust the different Lvalues in Dalton.
The CAtom, with the following configuration has three different states:
(1s^2 2s^2) 2p^2:
^3P, ^1D, ^1S.
So in D2hsym, there are for the px,y,z orbs, b3u,b2u and b1u.
The groundstate is ^3P and, also according to chap.24 in the manual, I did the input as shown in the appendix. But how one has to treat the other possibillities is not clear to me ?
Knowing how to adjust different Lvalues, is also important if one wants to calculate different states to a certain configuration of a molecule.
Please help me with this topic.
Thank you!
Best greetings
Alfred
CAtom

 Posts: 34
 Joined: 11 Jan 2014, 13:59
 First name(s): Alfred
 Last name(s): Güthler
 Affiliation: privat
 Country: Germany
CAtom
 Attachments

 C_engy_3B1g_C_631gdp.out
 (33.56 KiB) Downloaded 17 times

 Posts: 600
 Joined: 15 Oct 2013, 05:37
 First name(s): Peter
 Middle name(s): Robert
 Last name(s): Taylor
 Affiliation: Tianjin University
 Country: China
Re: CAtom
For the 3P state your input is correct but you should recognize that this is only one component of a triply (spatially) degenerate state. The pz orbital that would result will not have the same spatial form as the singly occupied px and py. This may not matter to you, of course. To obtain three equivalent p orbitals a fractional occupation of 2/3 of an electron in each p component is required. This can be enforced in some openshell SCF codes (although a dwindling number as the years go by...) but not in Dalton. An alternative possible with some MCSCF codes (but not Dalton) is to calculate the average of three states, with occupations, pxpy, pxpz, and pypz. This is exactly equivalent to the openshell SCF calculation with fractional occupation numbers.
For the singlet states the situation is more complex. The 1S state has the configuration (for full spatial symmetry) px^2 + py^2 + pz^2 which obviously has spin multiplicity 1 and symmetry 1 in D2h. But the five components of 1D appear in symmetries 1 (twice), 4, 6, and 7. You can thus run a CASSCF and take out a root in symmetry 4 or 6 or 7 and get a component of a 1D state, but again the three p orbitals will not be equivalent. The components in symmetry 1 will also not have equivalent p orbitals, and both will lie below 1S. So you could try taking the third root in symmetry 1 out of a CASSCF calculation to get your 1S state.
None of this is terribly satisfactory from a symmetry perspective. In general, if one wants orbitals in calculations on atoms or diatomics or other highsymmetry systems the easiest tactic is stateaveraged MCSCF as I mentioned above, but this is not available in Dalton.
Best regards
Pete
For the singlet states the situation is more complex. The 1S state has the configuration (for full spatial symmetry) px^2 + py^2 + pz^2 which obviously has spin multiplicity 1 and symmetry 1 in D2h. But the five components of 1D appear in symmetries 1 (twice), 4, 6, and 7. You can thus run a CASSCF and take out a root in symmetry 4 or 6 or 7 and get a component of a 1D state, but again the three p orbitals will not be equivalent. The components in symmetry 1 will also not have equivalent p orbitals, and both will lie below 1S. So you could try taking the third root in symmetry 1 out of a CASSCF calculation to get your 1S state.
None of this is terribly satisfactory from a symmetry perspective. In general, if one wants orbitals in calculations on atoms or diatomics or other highsymmetry systems the easiest tactic is stateaveraged MCSCF as I mentioned above, but this is not available in Dalton.
Best regards
Pete

 Posts: 34
 Joined: 11 Jan 2014, 13:59
 First name(s): Alfred
 Last name(s): Güthler
 Affiliation: privat
 Country: Germany
Re: CAtom
Dear Mr. Pete Taylor,
thank you for your quick and interesting answer.
You say that I could use SAstate averaging, but one has to use different codes, to perform these kind calculations. But I have to understand throughly what e.g. SA3 means. Surely more than the averaging of three states and it is interesting to note, that you speak of three components. It's a bit confusing for me, but do you mean the averaging of the three states, namly the ^3P, ^1D, ^1S, or that of the 3 pstates ?
In any way, this topic does not concern that of this forum and therefore I will stop here.
I want to mention that according to Herzberg, there is a connection of atomic and molecular states, which is marvelous for me. In fact, in the case of the Catom, which has 6 electrons, there should be also a molecule  Li2  which has similar states.
Finally one last question regarding, which I beg to answer me is concerning : STATES. What do you mean by the term " states", electronic, spectroscopic or molecular states ?
Of course these terms are in conection, but it is not always clear to me, what is meant by "states".
Thank you !
Best regards
Alfred
P.S.:
There are (2S+1)(2L+1) functions which desribe Terms to a certain Electronenconfiguraion. Therefore in the case of ^1D (L=2,S=0) , there are 5 Eigenfunctions.
thank you for your quick and interesting answer.
You say that I could use SAstate averaging, but one has to use different codes, to perform these kind calculations. But I have to understand throughly what e.g. SA3 means. Surely more than the averaging of three states and it is interesting to note, that you speak of three components. It's a bit confusing for me, but do you mean the averaging of the three states, namly the ^3P, ^1D, ^1S, or that of the 3 pstates ?
In any way, this topic does not concern that of this forum and therefore I will stop here.
I want to mention that according to Herzberg, there is a connection of atomic and molecular states, which is marvelous for me. In fact, in the case of the Catom, which has 6 electrons, there should be also a molecule  Li2  which has similar states.
Finally one last question regarding, which I beg to answer me is concerning : STATES. What do you mean by the term " states", electronic, spectroscopic or molecular states ?
Of course these terms are in conection, but it is not always clear to me, what is meant by "states".
Thank you !
Best regards
Alfred
P.S.:
There are (2S+1)(2L+1) functions which desribe Terms to a certain Electronenconfiguraion. Therefore in the case of ^1D (L=2,S=0) , there are 5 Eigenfunctions.
Who is online
Users browsing this forum: No registered users and 3 guests