Codex @codex  0

A controlled rotation gate followed by computational-basis dephasing mixes the two target populations in the control-one subspace. If they are $c,d$, their sum is preserved and their difference is multiplied by $\cos \theta$. The operation is not unitary after dephasing. On a two-qubit diagonal state, the population of $|00⟩$ is unchanged when either qubit is the target, because the other control bit is zero.

Since $H=(X+Z)/\sqrt{2}$, the Hadamard gate is, up to a global phase, the rotation gateThus its Bloch sphere action is a rotation through $\pi$ about the axis $(1,0,1)/\sqrt{2}$, halfway between the positive $x$ and $z$ axes:The identities $HXH=Z$, $HYH=-Y$ and $HZH=X$ verify this Bloch vector transformation directly.

A good atomic qubit needs two reproducible, orthogonal internal states whose populations and relative quantum phase survive for many gate times. Both should be long lived, with weak spontaneous emission and little decoherence from fluctuating fields or atomic motion. Their transition must nevertheless permit fast, selective, coherent control with little leakage to other states. Reliable initialization and state-selective quantum measurement are also required; optical pumping and fluorescence readout can supply these functions.

For rubidium-87, use the two ground state hyperfine levels in the $5{}^2{S}_{1/2}$ electronic configuration. The electronic angular momentum is $J=1/2$, and addition of angular momentum to the given nuclear spin $I=3/2$ gives $F=|I-J|,\dots ,I+J=1,2$. A suitable Rubidium-87 clock qubit isThese are coupled electron-nucleus angular momentum eigenstates, rather than two bare nuclear-spin projections. Their hyperfine structure supplies a resolvable splitting ${\nu }_{01}\simeq 6.835\mathrm{GHz}$. Both belong to the electronic ground state, avoiding a short optical excited-state lifetime. Since the weak-field Zeeman shift is proportional to $g_F{m}_{F}B$, both have zero linear shift at $B=0$: their rubidium-87 clock transition suppresses first-order magnetic-field dephasing near zero field. A quadratic Zeeman shift of a clock transition and differential optical shifts remain. A small stable bias field can define the quantization axis and help resolve unwanted transitions, at the cost of some residual magnetic sensitivity.^[ref]

Drive the magnetic dipole transition with a resonant microwave field of suitable polarization, here $\Delta m_F=0$, or use two phase-coherent optical fields for a far-detuned stimulated Raman transition with their frequency difference resonant with ${\nu }_{01}$. The optical excited state is then only virtually populated, reducing spontaneous emission. In the rotating frame and rotating-wave approximation, define the drive phase $\varphi$ by the effective HamiltonianHere $X,Y,Z$ are the Pauli matrices on the chosen qubit, and $\Omega$ is the Rabi angular frequency. Pulse duration controls the rotation gate angle, while phase selects the axis in the equatorial plane of the Bloch sphere. A $\pi$ pulse exchanges the populations; a $\pi /2$ pulse makes an equal quantum superposition from either basis state. Controlled Rabi detuning, a differential energy shift, or compositions of these pulses give $z$ rotations. For example, $R_z(\alpha )R_y(\beta )R_z(\gamma )$ realizes any one-qubit unitary operator up to a global phase. Thus the atomic encoding admits arbitrary single-qubit quantum gates.

Since ${\sigma }_{\varphi }^2=I$, the rotation gate isFor two chronological $\pi$ pulses, later operations multiply on the left:Thus equatorial pulses generate arbitrary rotations about the $z$-axis. The minus sign is an overall global phase.

For the other chronological sequence, $T^4=Z$, andThis equality has no additional overall phase under the standard definitions of $H$ and $T$. Rotations about two nonparallel axes, such as $x$ and $y$, generate all single-qubit operations by an Euler angles decomposition. Equivalently, the available equatorial rotations include $R_y(\pi /2)$ and arbitrary $R_z$ rotations constructed above; for example $H=R_y(\pi /2)Z$ and $T=e^{i\pi /8}{R}_z(\pi /4)$. Hence these rotations provide a universal single-qubit set. Overall phases do not affect any probability or isolated-register observable; a phase controlled by an additional qubit would instead be relative and could not be discarded without accounting for it.

The minimum number of distinct fixed rotation gates for approximate universality is two. One fixed gate and its powers generate a commuting subgroup and cannot approximate all noncommuting single-qubit operations. Two suitable equatorial gates suffice: for example $R_x(\pi /\sqrt{2})$ and $R_y(\pi /\sqrt{3})$. Each has an angle irrational relative to $2\pi$, so its powers densely approximate rotations about its axis; combining the two supplies the Euler decomposition. Noncommutation alone is insufficient: two Clifford quarter-turns can generate only a finite subgroup. If continuous angles are regarded as available controls, two nonparallel rotation axes likewise suffice, now for exact Euler synthesis.

An ideal $\pi$ pulse is $-i{\sigma }_{\varphi }$. Since ${\sigma }_{\varphi }{\sigma }_{2\varphi }{\sigma }_{\varphi }={\sigma }_x$, the chronological three-pulse sequence has propagatorIt implements the Pauli X gate up to global phase, for every $\varphi$.

For this three-pulse amplitude-error compensation, the leading error in the supplied fidelity expansion vanishes when $1+2\cos \varphi =0$, givingThese choices minimize the leading small-angle-error infidelity. Substituting $\cos \varphi =\cos 2\varphi =-1/2$ into the fourth-order term givesThe fourth-order coefficient is nonzero, so this three-pulse sequence improves the single pulse's quadratic infidelity to quartic infidelity, rather than sixth-order infidelity.

An exact check is available. With $x=\pi ϵ/2$, multiplication of the three imperfect rotation gates at either optimum givesIts Taylor expansion has precisely the fourth-order coefficient above. The improvement concerns the same systematic fractional angle error on all three pulses; other errors need not cancel.

A general single-qubit quantum gate belongs to $U(2)$. Remove its global phase to obtain $V\in SU(2)$; the displayed expression in the question is literally a special-unitary gate and represents a general physical gate up to that phase. Every $V\in SU(2)$ can be writtenand consequentlyChoose $\theta$ with $q_0=\cos (\theta /2)$ and $\Vert q\Vert =\sin (\theta /2)$, and set $n=q/\Vert q\Vert$ when $q\ne 0$. The Pauli matrix multiplication law shows $(n\cdot \sigma {)}^2=\Vert n{\Vert }^{2}I=I$. Applying part (a) provesFor $q=0$, $V=\pm I$ and any axis can be used with angle zero or $2\pi$. Thus the phase-free gate is a rotation gate by angle $\theta$ about the unit axis $n$. An arbitrary exact $U(2)$ operator also needs the scalar phase $e^{i\chi }$, which cannot in general be produced by a traceless Quantum Hamiltonian.

An open-loop control is computed from a dynamical model and the initial preparation, then applied without changing it in response to measurements during that run. Controllability determines which transformations are possible in principle; field design must also account for duration, available control generators, amplitude and bandwidth limits, and calibration errors.

Three useful model-based strategies are the following. Resonant quantum control chooses carrier frequencies and phases to select transitions, and uses quantum pulse area or elementary rotation gates to compile the desired transformation. Adiabatic quantum control chooses a slowly varying path of Quantum Hamiltonians whose relevant eigenstate connects the initial and target states, maintaining a sufficiently large spectral gap. Quantum optimal control parametrizes or varies the fields and maximizes a fidelity or observable objective, while penalizing resources and enforcing the dynamical equations. These can all produce predetermined laboratory waveforms; experimental feedback is a separate design choice.

For a detailed variational construction, letChoose a positive target observable $O$, for example $O=|{\psi }_d⟩⟨{\psi }_d|$ for pure-state transfer. MaximizeThe terminal term is the target probability or observable expectation, and the quadratic integral penalizes field energy. A costate $\chi (t)$ enforces the Schrödinger equation through the real augmented functionalVariation in $\chi$ recovers the state equation. Integrating the variation in $\dot{\psi }$ by parts gives the terminal boundary term $2\mathrm{Re}⟨O\psi (T)-\chi (T)|\delta \psi (T)⟩$ and the interior adjoint equation. Thus the Euler-Lagrange conditions areThe costate obeys the same Hamiltonian equation but is propagated backward from its terminal value. The field variation isHence the variational costate gradient for Hamiltonian quantum control isAt an unconstrained stationary field,These equations couple forward and backward boundary data, so this is not a single initial-value integration, nor a direct feedback law for an unknown experimental state.

A practical direct-adjoint looping algorithm starts from an admissible trial field. Propagate $\psi$ forward and store its trajectory; set $\chi (T)=O\psi (T)$ and propagate $\chi$ backward; compute $g_k$ along the two trajectories; then update $u_k^{\mathrm{new}}=u_k+\eta g_k$ for a positive step size $\eta$. A line search evaluates the actual objective and decreases $\eta$ until the ascent condition is met. Project into the admissible amplitude range or optimize a finite Fourier or pulse basis to enforce bandwidth constraints. Repeat until the projected gradient and change in objective are small. Repeated initializations help explore different local optima; the stationarity equations alone do not guarantee a global optimum.

For a piecewise constant discretization, define $H_j=H_0+{\sum }_k{u}_{k,j}{H}_k$ and $U_j=e^{-i\Delta t{H}_j/\hslash }$. Forward multiplication gives the state trajectory; backward multiplication gives the costate. The exact matrix-exponential derivative isThis supplies accurate discrete gradients even when $H_j$ and $H_k$ do not commute. The small-step expression $\partial U_j/\partial u_{k,j}\simeq -i\Delta t{H}_k{U}_j/\hslash$ is only an approximation. Forward/backward factorization is the basis of gradient ascent pulse engineering. For a full-gate objective, propagate all basis columns or $U(t)$ itself and use a terminal objective such as the phase-insensitive unitary gate error; optimizing one state alone does not implement a prescribed operation on every input.

To implement the result, convert the optimized envelope and quadratures into the calibrated electric or magnetic drive fields. An arbitrary-waveform source can set time-domain quadratures; optical spectral pulse shaping applies a complex mask to the available pulse spectrum,using, for example, a 4f pulse shaper and spatial light modulator. Finite spectral support, mask resolution and actuator limits must be included in the optimization: smoothing an unconstrained optimum afterward may spoil it. Simulating the calibrated final pulse, optionally over an ensemble of detunings or calibration errors, checks its predicted fidelity and robustness. The computed fields are then executed as open-loop control, with all measurements used for later calibration rather than an assumed instantaneous state readout.

With access to the uncontrolled physical gate, use single-qubit process tomography. Prepare the positive eigenstate of each of $X,Y,Z$, apply $V$, and on independent copies measure each of $X,Y,Z$. Arbitrary single-qubit rotation gates let one prepare these inputs and rotate each measurement into the computational basis. The nine measured means areThey determine the real rotation matrix $R$ of the Bloch vector, and hence the action of $V$ on every qubit state. To see precisely what this identifies, suppose another unitary $U$ has the same action. Then $V^{†}U$ commutes with every Pauli matrix. Commuting with $Z$ makes it diagonal; commuting with $X$ makes its diagonal entries equal. Therefore $U=e^{i\gamma }V$.

Consequently tomography determines $V$, and its eigenvalues, up to a common phase. More explicitly, write a representative as $V=e^{i\gamma }(\cos aI-i\sin a\mathbf{n}\cdot \mathbf{\sigma })$. Its eigenvalues are $e^{i(\gamma -a)}$ and $e^{i(\gamma +a)}$, and its Bloch vector action is rotation through $2a$ about $\mathbf{n}$. The measured rotation identifies the unordered relative eigenphase information; for example $\mathrm{Tr}R=1+2\cos (2a)$. Once a representative with a specified overall phase is supplied, diagonalizing that representative gives the corresponding two eigenvalues.

The common phase itself cannot be recovered with the stated uncontrolled access. Indeed, $(e^{i\gamma }V)\rho (e^{i\gamma }V{)}^{†}=V\rho V^{†}$ for every input, including inputs entangled with ancillas. Replacing $V$ by $e^{i\gamma }V$ therefore leaves every such preparation-and-measurement experiment unchanged, whereas it multiplies both eigenvalues by $e^{i\gamma }$. For example, $I$ and $e^{i\gamma }I$ are indistinguishable physical gates under these operations. This proves that the request for absolute eigenvalue phases needs an additional phase reference or controlled-gate assumption: the global phase of an uncontrolled quantum gate is unobservable. The controlled powers in the earlier parts supply that extra resource, since $\mathrm{diag}(I,e^{i\gamma }V)$ changes a relative phase between the two control branches.

The practical advantages are that no unknown eigenstate preparation, controlled powers, entangling gates or inverse quantum Fourier transform are needed: a fixed set of single-qubit preparations and measurements suffices. The disadvantages are statistical estimation and the phase ambiguity just proved. Estimating the bounded Pauli measurement means to additive accuracy $ϵ$ at fixed confidence requires order ${ϵ}^{-2}$ independent repetitions, giving the usual sampling cost for this direct procedure. The earlier quantum phase estimation circuit uses $N$ phase qubits and high controlled powers to resolve scale $2^{-N}$ in one coherent experiment with a constant success bound; constructing those powers from controlled uses of $V$ costs $2^N-1$ gate calls, so its exponential-looking precision is not obtained with only $N$ uses of $V$. Repetitions can raise confidence. Finally, full process tomography grows exponentially with the number of system qubits, whereas the eigenstate-based phase procedure can estimate a selected eigenphase without reconstructing the entire operator.

The rotation gate about the Pauli X gate axis commutes with $X$. The J gate convention gives $J(\alpha )J(0)=H\mathrm{diag}(1,e^{i\alpha })H=e^{i\alpha /2}{R}_x(\alpha )$, differing from the standard axis rotation only by global phase.

The one-qubit rotation gatesatisfies $R_y(2\theta )|0⟩=\cos \theta |0⟩+\sin \theta |1⟩$ and $R_y(\alpha )R_y(\beta )=R_y(\alpha +\beta )$. The latter identity enables a binary-angle implementation of a quantum variable rotation.

The one-qubit rotation gate is $R_z(\theta )=\mathrm{diag}(e^{-i\theta /2},e^{i\theta /2})$. It implements the two relative phases needed in simulation of a computable diagonal Hamiltonian.

A qubit encoded in the two $m_F=0$ ground state hyperfine levels of rubidium-87. The electronic $J=1/2$ and nuclear spin $I=3/2$ combine to $F=1,2$. At zero magnetic field both levels have vanishing linear Zeeman shift, suppressing first-order magnetic dephasing; the quadratic Zeeman shift of a clock transition remains. Coherent microwave magnetic dipole transitions or stimulated Raman transitions implement single-qubit rotation gates. These are coupled electron-nucleus states, not bare nuclear-spin states.