No. of occ. orb. > total no. of orb.

Find answers or ask questions regarding Dalton calculations.
Please upload an output file showing the problem, if applicable.
(It is not necessary to upload input files, they can be found in the output file.)

Post Reply
pandey28
Posts: 5
Joined: 02 Feb 2020, 17:18
First name(s): Sidd
Last name(s): Pandey
Affiliation: Purdue Physics
Country: United States

No. of occ. orb. > total no. of orb.

Post by pandey28 » 23 Nov 2020, 05:12

I am trying to calculate the third-order response for a singly-ionized Xenon ion with a half-filled 4dz^2 orbital. I am using the "aug_cc_pvdz_pp" basis set since the manual says it can be used for Xenon. However, I keep getting the following error:

Node 0: --- SEVERE ERROR, PROGRAM WILL BE ABORTED ---
No. of occ. orb. > total no. of orb.


Any insights will be appreciated.
Last edited by pandey28 on 23 Nov 2020, 16:33, edited 1 time in total.

hjaaj
Posts: 395
Joined: 27 Jun 2013, 18:44
First name(s): Hans Jørgen
Middle name(s): Aagaard
Last name(s): Jensen
Affiliation: Universith of Southern Denmark
Country: Denmark

Re: No. of occ. orb. > total no. of orb.

Post by hjaaj » 23 Nov 2020, 08:30

Please upload your output file (input files are not needed, as they are listed in the output).
(Use "Attachment" tab in "full editor" mode)

pandey28
Posts: 5
Joined: 02 Feb 2020, 17:18
First name(s): Sidd
Last name(s): Pandey
Affiliation: Purdue Physics
Country: United States

Re: No. of occ. orb. > total no. of orb.

Post by pandey28 » 23 Nov 2020, 15:55

I have already found a resolution to this problem: I believe the number of orbitals in my CAS space simply exceeded the number of orbitals in my basis set. So I tried with a smaller CAS space and it worked.
Attachments
Xe+[4d].txt
(24.36 KiB) Downloaded 32 times

hjaaj
Posts: 395
Joined: 27 Jun 2013, 18:44
First name(s): Hans Jørgen
Middle name(s): Aagaard
Last name(s): Jensen
Affiliation: Universith of Southern Denmark
Country: Denmark

Re: No. of occ. orb. > total no. of orb.

Post by hjaaj » 23 Nov 2020, 19:52

Yes, you can see the negative numbers for "NSSH", the number of secondary orbitals in each symmetry.

The first basis functions in symmetry 8 are f orbitals, thus zero with our basis set.

The symmetries 1-8 corresponds to

Code: Select all

sym 1: s,  dx2-y2, dz2 
sym 2: px, fxxx-xyy, fxzz
sym 3: py, fxxy-fyyy, fyzz
sym 4: dxy
sym 5: pz, fxxz-fyyz, fzzz
sym 6: dxz
sym 7: dyz
sym 8: fxyz
(You can in fact see this from the listing of symmetry orbitals in the output.)

When I have done calculations on atoms I identified which orbitals have a degeneracy of 1, 3, 5, 7, ... ( s, p, d, f, ... type, respectively).

Two more comments: I would use MP2 natural orbital occupation numbers to select good CAS spaces. And you need bigger basis sets to get good numbers for cubic response, in particular you need f-functions or your results will be garbage, and probably also g-functions.
This is because you need the second order wave function for cubic response, and
two dipole operators working 1) on a p-orbital gives an f-orbital and 2) on a d-orbital gives a g-orbital !! Even if you do not need them for the excitation energies.

Post Reply

Who is online

Users browsing this forum: No registered users and 1 guest