Bilaterally Symmetrical Transformation between Independent Operators and Rotations

This report describes an approach for representation of quantum operators through rotations and rotation through quantum operators. The approach of the proposed method transforms rotation in a kind of a unitary matrix that corresponds to the rotation. Operations with qubits are very similar to the rotation, but with an added phase coefficient. This fact is used to create a process for rotation between unitary matrices. This approach could be used to modifying the controls to apply in a different basis.


Introduction
The fundamental concept for realization of the transitions is based on the fact that the quantum operation is always just a unitary matrix, which may be a linear interpolation between the matrices: ( ) . The operations with single qubits are very similar to the rotations, but with an added coefficient of the phase.This fact will be used to create a method for transformation of rotational into qubit operations.
Rotations of the eigenvector of qubit could be viewed as operations in a single Bloch sphere [1]- [3].A more formal approach is to provide two or three fixed orthogonal axes of rotation.This allows arbitrary rotations in three or less steps, for example, by using the angle of the construction of Euler.Logical qubits, composed of two or more physicalturns, are possible physically different binding mechanisms to control different axes of rotation of the field unit.For example, logical qubits formed in double quantum dots use local gradients for generating rotations around an axis, and switched interactions generate rotations around an orthogonal axis [4] [5].In this paper, cases of intermediate situations between these two extremes are investigated, where rotations can be performed around an arbitrary axis constrained to lie in one plane.The main motivation for our work is the example of single exchange of logical Qubit in triple quantum dot [6] [7].We will build on the known results [8] that the internal connections between the physical spins in one exchange-only qubit can be used to generate a consistent set of rotations in the XZ basis of qubit [9].
In the development process of a circuit, quantum simulator [10]- [15] was required to seek solutions to approximate the results of the quantum operations, in order to animate what occurs without any sharp jumps.More generally, a transition was necessary between the two operations without any jumps.This paper will show that such rotations in a single plane reduce the total number of steps required for Qubit operations.For example, a transformation of a state that converts specific baseline to a certain final state can be carried out in one stage, random rotation single-Qubit can be done in two stages.We offer constructive solutions to these problems.Given that the space of the unit matrix is very similar to the space of rotations, it is useful to turn into a kind of rotation of the unit matrix which corresponds to the rotation in some useful way.The simplest method of achieving this is to provide such a rotation as a quaternion, and then replace the quaternion components with the relevant times of Pauli matrices: If the input rotation is a vector v, where the direction of v is the axis of rotation and the length v is the number of rotations in radians, then the unitary matrix will turn this rotation in: where ˆxyz vσ gives the vector of Pauli in v.The problem with the linear transformation is that the intermediate matrices may not be valid operations.The linear transformation tends to create a matrix entries, which are too close to zero, i.e. the resulting matrices will shrink the values instead of retaining their length (which is not a desired effect, as the whole point of using unitary matrices is to preserve the length).Actually, the goal is to be made a transformation without leaving the space of the unitary matrices.A compact way for the parameterization of the space of the unitary matrices is: The above equation includes four constants and three variables.The constants are the identity matrix (I), the square root of −1 (i), the constant of Euler (e) and the vector of Pauli matrices xyz σ .The three variables are the angle φ , the angle θ and the single vector v .Each of the variables plays a different role.φ is a global phase coefficient.It's what distinguishes the group of unitary matrices U(2) from "the special unitary group SU(2)".v and θ correspond to a rotation.v is like an axis to rotate around, and θ is how much to rotate around the said axis.What does it mean that v and θ are like a rotation?It becomes a bit clearer when the compact parameterization from above is expanded.Through incorporation of the Pauli matrices and splitting of v in x y z , the following is obtained: This is not yet clear enough.The following part is omitted: the equation to convert from an axis-angle style rotation to a single quaternion style rotation: The resemblance is visible.Without paying attention to the mysterious division by 2 of the angles, the Pauli matrices actually play the role of quaternion constants: i, j and k.If each of the Pauli matrices is multiplied by i, the following is obtained: This, in turn, looks very similar to the way in which the quaternions are defined: . Why is this similarity with the rotations important?Because it will be used for linear transformation.There are already existing methods for smooth linear transformation between quaternions and they will be applied in order to be handle the rotation part of the unitary operation.Then the remaining phase part simply must be interpolated between two angles.
The main contribution of this paper is the proposed method of mapping rotations in unitary matrices.

Converting Rotation into Qubit Operations
The general comparison has got many good features, but also two flaws that, if possible, should be avoided.The first flaw is that the angles are divided by 2. Therefore, a rotation of 360˚ does not return to the starting position, but to I − .If we use this mapping, the rotation should be of 720˚, which is confusing.The second flaw is that there is no rotation, which corresponds to Pauli matrices.Rotation of 180˚ around the X axis gives x iσ instead of x σ .This means that when using these rotations to define operations for quantum computer, a NOT exit can- not be even made!It is possible these flaws to be avoided, but this will be at the expense of something else.

Phase Correction
Both of the aforementioned flaws are caused by one problem: the global phase coefficient is wrong.In order to fix it, we have to add a specially prepared coefficient of the counter-phase.Currently a half-turn around the X axis gives x iσ .To turn it into x σ , a phase correction factor of -i will be needed.Similarly, a coefficient with phase correction of −1 for one full turn is also needed.After three half rotations around the X axis (or one half turn around the negative X axis) a coefficient with phase correction of i is needed.The coefficients with phase correction are the same for the Y and Z axes with respect to the amount of turning.Given the above information, it's clear that is needed a correction coefficient such as , where s is either +1, or −1, depending on the direction of rotation.Choosing s is a delicate moment, which will be discussed later.For now, it will be used in the mapping as taken for granted: Because the mapping formula already has two coefficients involving half-angles being multiplied together, there are possibilities for their simplification.Let us first reveal e ix in the trigonometric functions: ( ) Then it will be divided, while the trigonometric multiplications are simple, but the elements are not yet grouped based on the matrix coefficients: Let's make the trigonometric coefficients more similar by subtracting the sign coefficients in relation to the parity (i.e.

( ) ( )
( ) ( ) ( ) Double-Angle formulas expressing trigonometric functions of an angle 2x in terms of functions of an angle x, ( ) Simplification of the coefficients: Subtraction of is from the component on the right side: The coefficient of s is moved into the trigonometric functions, so they can be merged in the exponential functions: This is already simple enough.The fact that the angles are no longer being divided by two 2, shows, that the problem with 720˚ is solved.Evaluating the mapping at ± should be canceled in the results by selecting s in an appropriate manner.How should be selected whether s must be +1 or −1?A naive solution would be to always use s = 1 or −1, but that would cause (for example) ,0,0 despite starting from equivalent rotations.Negating a vector has to negate s, because, otherwise, repeatedly doing and undoing a rotation would give matrices that were not inverses of each other.This in turn will cause an increase in the phase coefficient instead of its change together with the rotations.The reversal of s without introducing a discontinuity in the phase correction requires v to be an integer from half of the turns, as this is the only time when π π e e i i θ θ − = .However, s must be reversed for all pairs v vs v − − , including those, which are not half rotations, and the axis, can be rotated so as to correspond to reversals.Thus, regardless of the actions, there always will be a discontinuity in the phase correction angle (Figure 1).
The plane, in which the discontinuity occurs, can be rotated so that no rotational axes are used on it.Note that, even if you avoid rotation axes along the bad plane, you will still run into problems when combining multiple rotations.This is a consequence of wanting half turns to correspond to the Pauli matrices.There's no way to avoid accumulating a global phase factor via multiple rotations because, although rotating a half turn around the X then Y then Z axis does not give a final rotation, and x y z iI I . In the end, despite the elimination of the two flaws from the original mapping (the rotation at 720˚ and the inability to obtain the Pauli matrices) now there is a discontinuity, where the sign of the result is switched when the rotational axis is swiveled, and there is a problem with the phase accumulation at multiple rotations.Whether or not these disadvantages are preferable over the initial ones depends on the particular application.

Implementation
Here is python code implementing the final mapping from above: A notable problem in the above code is the division of the rotation angle in order to compute the unit vector along the axis of rotation.This will cause problems at small rotations.Thus, the part e isθ θ − must be rewritten so it does not involve a division.To avoid the division can be used the sine function Let's first reveal the exponential functions into trigonometric, then they should be simplified and the half-angle identities should be used and they should be simplified again: Although the sine also has a division by zero, it can be used the fact that sin x acts as x near zero in order to eliminate the problem.Upon approaching to zero, the series x + starts being rounded to 1.Given the stable sine function, and the arbitrarily chosen plane-of-terribleness we can write the function rotation_to_matrix: defquantum_ rota-tion_to_matrix (x=0, y=0, z=0): It returns a unitary matrix that corresponds, in a useful, but not unique way, to a rotation around the axis <x, y, z> by sqrt(x^2 + y^2 + z^2) radians.
# Discontinuance of the phase correction in an inconvenient plane s = math.copysign(  The above method will be applied at the method for full transformation, when a spherical transformation is being carried out.In order to make an angular interpolation the obvious shall be carried out: the difference between the two angles is learned, care should be taken to recourse to a roundabout way and then a linear transformation should be made.To take correctly the sign of the difference is a difficult task, but it is already explained.When everything is put together, we obtain:

Demonstration
Below is given a demonstration of the described above method with the developed by the author of the article quantum simulator (Figure 2).Matrices for start and end can be entered in the text boxes and (after the entered matrices are adjusted, they should be unitary) a continuous transition between the two matrices is shown.It is difficult to be checked ostensibly, whether the intermediate matrices are unitary, but it can be seen that the movement is smooth and the colored area remains relatively constant.

Conclusion
Single qubit operations are a lot like rotations, but with an added phase coefficient.This fact can be used to create a method for linear transformation between 2 unitary matrices.The method described above, works, but is not optimal.For example, it does not ensure a constant angular speed.Also, in some cases it doesn't take the shortest possible path.The common method for mapping rotations into unitary matrices is smooth, but cannot generate the Pauli matrices and requires a 720˚ turn to get back to the starting point.By applying a phase correction, we can fix those issues, but we are forced to introduce a phase discontinuity with the respect to the axis of rotation.
are the correct results, with the exception of the detail that s

Figure 1 .
Figure 1.Discontinuity in phase correction angle.Note: It not is necessary to have only one discontinuity.There may be several.

Figure 2 .
Figure 2. Linear transformation.(Note: The input correction is done by doing a singular value decomposition and omitting the non-unitary factor.This turns out to be really effective.)