I am trying to do calculations on adsorptions of singlet and triplet oxygen molecule on a Cu(100) surface? Apart from the ISPIN tag in INCAR, does VASP allow you to specify either triplet or singlet states? So far, my results do not change when I change the values of ISPIN. If no, any other software can do that?