Dear all, I tried to generate PDOS by psodxml program. There are several points that I can not understand. First, How to find the orbital that I am interested in. In the m_orbital_chooser.f90 there are explanation, what does orbid stand for? I tried but always failed, because I donot know the orbid meaning. ! 1. Want only s-orbitals ! wantit = ( orbid%l == 0 ) ! 2. Want only n=3 orbitals ! wantit = ( orbid%n == 3 ) ! 2. Want 3p orbitals ! wantit = ( ( orbid%n == 3 ) .and. (orbid%l == 0 ) ) ! 3. Want Oxygen orbitals wantit = ( orbid%species == "O-wc" ) !wantit = .true. ! wantit = ( orbid%species == "O" ) Second, when I choose the last " wantit = ( orbid%species == "O" )" for the example "h2o_dos.PDOS" ,why I can get this 2d orbital. I get this "
Orbital:(n.l.m.z.index.atom_index.species): 2 0 0 1 1 1 O ------------------------> 2s Orbital:(n.l.m.z.index.atom_index.species): 2 0 0 2 2 1 O ------------------------> 2s* Orbital:(n.l.m.z.index.atom_index.species): 2 1 -1 1 3 1 O ------------------------> 2px Orbital:(n.l.m.z.index.atom_index.species): 2 1 0 1 4 1 O ------------------------> 2py Orbital:(n.l.m.z.index.atom_index.species): 2 1 1 1 5 1 O ------------------------> 2pz Orbital:(n.l.m.z.index.atom_index.species): 2 1 -1 2 6 1 O ------------------------> 2px* Orbital:(n.l.m.z.index.atom_index.species): 2 1 0 2 7 1 O ------------------------> 2py* Orbital:(n.l.m.z.index.atom_index.species): 2 1 1 2 8 1 O ------------------------> 2pz* Orbital:(n.l.m.z.index.atom_index.species): 2 2 -2 1 9 1 O ------------------------> 2dxy Orbital:(n.l.m.z.index.atom_index.species): 2 2 -1 1 10 1 O ------------------------> 2yz Orbital:(n.l.m.z.index.atom_index.species): 2 2 0 1 11 1 O ------------------------> 2zx Orbital:(n.l.m.z.index.atom_index.species): 2 2 1 1 12 1 O ------------------------> 2 z2 Orbital:(n.l.m.z.index.atom_index.species): 2 2 2 1 13 1 O ------------------------> 2 (x2-y2) " I must misunderstanding something, so could anyone tell me what is the problem,please? I will be very appreciated. Thank you very much. Regards, Gao Min