We present an implementation of the spin-dependent unitary group approach to calculate spin densities for configuration interaction calculations in a basis of spin symmetry-adapted functions.Using S 2 eigenfunctions helps to reduce the size of configuration space and is beneficial in studies of the systems where selection of states of specific spin symmetry is crucial. To achieve this, we combine the method to calculate U(n) generator matrix elements developed by Downward and Robb (Theor. Chim. Acta 1977, 46, 129) with the approach of Battle and Gould to calculate U(2n) generator matrix elements (Chem. Phys. Lett. 1993, 201, 284). We also compare and contrast the spin density formulated in terms of the spin-independent unitary generators arising from the group theory formalism and equivalent formulation of the spin density representation in terms of the one-and two-electron charge densities. K E Y W O R D S configuration interaction, configuration state function, quantum chemistry, spin density, unitary group approach 1 | I N TR ODU C TI ONThe unitary group approach (UGA) in quantum chemistry was pioneered by Paldus in 1974, [1] with practical algorithms including the most famous graphical approach by Shavitt, [2,3] having been developed in the following years. It offered a simple and straightforward way of evaluating the spinless generator matrix elements (ME) using implicit Gelfand-Tsetlin basis [4,5] and was suitable for calculating many-electron wavefunction expectation values for spin-independent Hamiltonians. [6] The growing demand for relativisitc corrections and other spin dependent properties in wavefunction methods inspired the subsequent development of spin-dependent generator ME in 80-90s by the several groups of authors [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23] and the latest work in this field is actually quite recent. [24] In particular, a theory for spin dependent operators formulated entirely within the unitary group formalism has been developed by Gould, Chandler, Paldus, and Battle. [7][8][9][10][11][12][13][14][15] As we shall discuss, the final formulae are easy to implement and are based on spin-independent U(n) and well-known U(2) generator ME in terms of the Gelfand-Tsetlin implicit basis. As a special, zero-order case, this theory provides a rather simple way for calculating spin density. [13,15,25] An obvious approach to spin density calculation with configuration interaction (CI) wavefunction would be to use the precomputed spindependent generator ME. The general single-particle reduced density operator (in a 2 3 2 matrix form) is given by [13] : qðr; r 0 Þ lm 5 X n i;j51