A novel procedure for evaluating Wigner coupling coefficients and Racah recoupling coefficients for U(4) in two group–subgroup chains is presented. The canonical $$\rm U(4) \supset U(3)\supset U(2)\supset U(1)$$
U
(
4
)
⊃
U
(
3
)
⊃
U
(
2
)
⊃
U
(
1
)
coupling and recoupling coefficients are applicable to any system that possesses U(4) symmetry, while the physical $$\mathrm{U(4)} \supset \textrm{SU}_{\mathrm{S}}(2) \otimes \textrm{SU}_{\mathrm{T}}(2)$$
U
(
4
)
⊃
SU
S
(
2
)
⊗
SU
T
(
2
)
coupling coefficients are more specific to nuclear structure studies that utilize Wigner’s supermultiplet symmetry concept. The procedure that is proposed sidesteps the use of binomial coefficients and alternating sum series and consequently enables fast and accurate computation of any and all U(4)-underpinned features. The inner multiplicity of a (S, T) pair within a single $$\mathrm U(4)$$
U
(
4
)
irreducible representation is obtained from the dimension of the null space of the $$\mathrm{SU(2)}$$
SU
(
2
)
raising generators, while the resolution for the outer multiplicity follows from the work of Alex et al. on $$\mathrm U(N)$$
U
(
N
)
. It is anticipated that a C++ library will ultimately be available for determining generic coupling and recoupling coefficients associated with both the canonical and the physical group–subgroup chains of U(4).