Quantum computers have an enormous impact on quantum chemical calculations. Approaches to calculate the energies of atoms and molecules on quantum computers by utilizing quantum phase estimation (QPE) and the variational quantum eigensolver (VQE) have been well documented, and dozens of methodological improvements to decrease computational costs and to mitigate errors have been reported until recently. However, the possible methodological implementation of observables on quantum computers such as calculating the spin quantum numbers of arbitrary wave functions, which is a crucial issue in quantum chemistry, has been discussed less. Here, we propose a quantum circuit to simulate the time evolution of wave functions under an S2 operator, exp(−iS2t)|Ψ〉, and integrate it into the QPE circuit enabling us to determine the spin quantum number of the arbitrary wave functions. We demonstrate that the spin quantum numbers of up to three spins can be determined by only one qubit measurement in QPE.