<div dir="ltr">Hello,<div><br></div><div>I am doing a triplet state calculation with SIC, OT, ROKS and MULTIP=3. The relaxed structure is correct (compared to literatures) and the spin component of each atom  also agrees well with literatures. However, when I plot the mo_cube files, I get lost. I have 2688 electrons in total with 1345 of spin 1 and 1343 of spin 2, and I want to look at the orbitals of the electron and hole of the triplet state. So I should see the WFN_1344_2-1_0.cube for the hole and WFN_1345_1-1_0.cube for the electron. However, after visualizing the cube files, the WFN_1344_2-1_0.cube file is not what I expect for hole but the WFN_1345_1-1_0.cube is. Then, I get confused.  After searching the groups, I find that printing MO_CUBE with ADDED_MOS will give some misleading indexes. So how about the ROKS case?</div><div><br></div><div>Part of my inputs are as follows:</div><div>-----------</div><div><div>  &DFT</div><div>    MULTIPLICITY 3</div><div>    ROKS</div><div>    &SIC</div><div>     SIC_METHOD MAURI_SPZ</div><div>     SIC_SCALING_A 0.25</div><div>     SIC_SCALING_B 0</div><div>    &END SIC</div><div>    &OT</div><div>      ROTATION</div><div>      PRECONDITIONER FULL_KINETIC</div><div>      ENERGY_GAP 0.001</div><div>    &END OT</div><div><div>    &PRINT</div><div>     &MO_CUBES</div><div>      WRITE_CUBE T</div><div>      NHOMO 2</div><div>      NLUMO 2</div><div>     &END MO_CUBES</div><div>    &END PRINT</div></div></div><div>-----------------------</div><div><br></div><div>Best,</div><div>Xiaoming</div><div><br></div><div><br></div></div>