Implement alternative Qiskit gates for (3,1,p) QRAC basis rotation

[garrison]

What is the expected enhancement?

Currently, we use the the Qiskit r gate to return the basis rotation corresponding to the (3,1,p) QRAC:

prototype-qrao/qrao/encoding.py

Lines 83 to 90 in 38527fd

	if base == 0:
	circ.r(-BETA, -np.pi / 4, i)
	elif base == 1:
	circ.r(np.pi - BETA, np.pi / 4, i)
	elif base == 2:
	circ.r(np.pi + BETA, np.pi / 4, i)
	elif base == 3:
	circ.r(BETA, -np.pi / 4, i)

However, this obscures the connection to the math in the paper. It would be more clear and helpful to the user if we instead implemented these basis transformations explicitly as written in Appendix IV of the paper (that is, with the relevant combination of rx, rz, and pauli x, y, z gates).

I have previously considered that we ought to keep it as a single gate, as that is the most precise way to specify the operation. However, I've since learned abou the Optimize1qGates[Decomposition] passes, which should be able to transpile to the optimal operations on hardware even if given as multiple single-qubit gates.

However, this might be at odds with parameterizing the gates for use with the primitives (#21), so we might want to see how that goes first.