Dear all, How can I print out the reduced density matrix with spin Alpha and Beta seperately in Molpro 2002.6? I know using the keyword "Load" to get charge density matrix, but it don't divide two different spin. Thanks. BingwuWang Chemistry Department Peking University, Beijing, China email: bingwuwang at vip.sina.com