Geometric Algebra Jordan–Wigner Transformation for Quantum Simulation

Quantum simulation qubit models of electronic Hamiltonians rely on specific transformations in order to take into account the fermionic permutation properties of electrons. These transformations (principally the Jordan–Wigner transformation (JWT) and the Bravyi–Kitaev transformation) correspond in a quantum circuit to the introduction of a supplementary circuit level. In order to include the fermionic properties in a more straightforward way in quantum computations, we propose to use methods issued from Geometric Algebra (GA), which, due to its commutation properties, are well adapted for fermionic systems. First, we apply the Witt basis method in GA to reformulate the JWT in this framework and use this formulation to express various quantum gates. We then rewrite the general one and two-electron Hamiltonian and use it for building a quantum simulation circuit for the Hydrogen molecule. Finally, the quantum Ising Hamiltonian, widely used in quantum simulation, is reformulated in this framework.


Introduction
Analysis of large molecules is a great challenge in computational chemistry.Studying the Shrödinger equation in most cases does not lead to analytic solutions, and therefore requires one to rely on simulations.The major goal in computational chemistry is to find the state energy spectrum of many body interacting fermionic systems and to investigate its wavefunctions and structural properties.Classical chemical simulations permit one to tackle efficiently physical, chemical, and structural properties of molecules.Due to the increase in computational power over the past decades, direct simulations have been made possible; for example, Density Functional Theory and Molecular Dynamics methods have enabled calculations for systems of more than a thousand atoms [1,2], but these methods suffer from an exponentially growing computational time.Some alternative methods have been developed in order to ease the computations for complex molecules, such as quantum Monte Carlo [3], tensor networks [4], post Hartree Fock methods [5], Green's function methods [6], many-body perturbation theory [7], etc.However, each method suffers from its own limitations; for example, up to now, tensor networks have been demonstrated to be efficient only for low entangled systems [8].
The general difficulty in simulating quantum mechanics with classical computers is linked to the exponential scaling of the Hilbert space dimension, with respect to the number of atoms in the molecule.In the 1980s, Richard Feynman [9] made the proposition to simulate quantum mechanics with a device capable of imitating the quantum rules rather than simulate it on a classical computer.This idea was taken further by David Deutsch and led to the concept of quantum computers [10].Quantum computers and quantum simulation are today a major research subject, involving private companies such as IBM [11], Google [12], Quantinuum [13], and Pasqal [14], along with universities all around the world.The seeked goal is to provide an advantage over classical computers [15] and to create an efficient tool for simulating quantum physical systems [16].
In quantum simulation, the molecular Hamiltonian under the Born-Oppenheimer approximation (BOA), where the atomic nuclei are considered fixed, is represented in the second-quantized form [17], using fermionic annihilation and creation operators.These operators anticommute for distinct electrons in a molecule; this means that they do not behave as qubits in a quantum circuit .Quantum circuits use the operations of matrix product to apply consecutive gates and Kronecker products to compose different qubits.This leads to the usual quantum circuit diagram where the consecutive boxes (from left to right) represent the applied quantum gates and the horizontal parallel wires represent the composed qubits behaving as commuting bosons.This fact makes the translation from the Hamiltonian form to quantum gates not straightforward.To practically obtain a quantum circuit from a Hamiltonian, one uses specific transformations such as the Jordan-Wigner transformation (JWT) [17,18] or the Bravyi-Kitaev transformation [19,20], which reestablishes the correct commutation properties.Then the Suzuki-Trotter approximation [17] is used to build the entire equivalent circuit.It is then possible to extract the useful information about the Hamiltonian, ground energies, for example, using different quantum algorithms.
Although JWT is effective in going from one formalism to another, it is commonly applied in an heuristic ad-hoc way.One could wonder whether alternative mathematical formalisms could make the application of this transformation more natural.We will show here that Geometric Algebra (GA) is an interesting candidate.
GA combines the algebra of quaternions, discovered by W.H. Hamilton, permitting one to describe geometric operations such as rotations and Grassman algebra based on the notion of exterior product, used to create objects of higher dimensions such as areas and volumes.The synthesis and generalization was done by Clifford [21] by combining quaternions and Grassmann's algebra into a single mathematical structure.Clifford himself named this structure "geometric algebra", although the term Clifford Algebra (CA) is more often used in the literature.
It must be emphasized that for historical reasons most physicists did not adopt GA for calculations on vector spaces, preferring Gibbs vector analysis, which almost completely replaced quaternion algebra by the end of the 19th century.However, in the second half of the 20th century, some physicists, lead by David Hestenes [22], used GA methods in classical mechanics, electromagnetism, quantum physics, and special and general relativity [23].
GA now also has many other applications in applied sciences outside physics, such as robotics [24] and computer graphics [21].More generally, GA has shown its potential for optimizing geometry-related problems [24].Current research has also used GA in quantum computing [25][26][27][28].Recently, a formalism for quantum computing was proposed built upon complex Clifford algebras using the Witt Basis (WB) [29].
In this paper, we provide a direct link between geometric algebra and quantum simulation, using both the WB and JWT.The goal here is to establish the equivalence between the algebra generated by the WB and the algebra generated by the creation and annihilation operators in the electronic Hamiltonians used in quantum simulation.
Section 2 recalls the basics of GA and then defines the concept of WB and of complex CA; it concludes with the spinor representation and the qubit analogy.The JWT is discussed in Section 3, and then the formulation of JWT in the framework of GA and WB is developed; the method is then applied to give the expressions of the usual 1-qubit and 2-qubit quantum gates and discusses the application on qubits in quantum circuits.Quantum simulation Hamiltonians in the BOA form are then expressed in the GA framework and a concrete example, the hydrogen molecule, is investigated in Section 4. In Section 5, an example of a general anisotropic XY Hamiltonian, which can be reduced to a quantum Ising Hamiltonian, commonly used as an ersatz in quantum simulation, is discussed in our framework.

Clifford Algebra (CA)
A CA is defined using a real vector space, V, with a quadratic form, ϕ, over a field [30].The algebra is produced by the quotient space T(V)/I ϕ (V), where T is the tensor product of vector spaces and For an orthonormal basis, {e 1 , . . ., e n }, spanning the vector space, V, the algebra is generated by the basis verifying e i e j + e j e i = 2ϕ(e i , e j ). ( An intuitive definition of a CA can be given using a finite dimensional vector space, V, equipped with the "geometric product" defined as the sum of an inner product and an outer product, which both have a geometrical meaning [21,22], as shown below.

Definition 1. Outer product
The outer product, noted by ∧, is an antisymmetric bilinear form giving the anticommutativity for the basis vectors: e i ∧ e j = −e j ∧ e i . ( The outer product of two vectors forms a new object, called a bivector.

Definition 2. Blade
The blades are generated by the outer product of vectors of V.The grade of the blade is the number of vectors composing the blade: is called an r-blade.

Definition 3. Inner product
The inner product is a bilinear symmetric non-degenerate form defined on V.
For an orthonormal basis, {e i }, one has: where p ≤ n.This characterizes the CA, noted Cℓ p,q , with p + q = n.For a general signature, p denotes the numbers of positive norm basis vectors and q the number of negative norm basis vectors.
The definition of the inner product is extended to all blades by the formulae where a 1 ∧ • • • ∧ a r and b 1 ∧ • • • ∧ b p are r-and p-blades, and u is a vector.Then, the inner product lowers the grade of a blade.
The subspace Span{e i |1 ≤ i ≤ p} forms an euclidean space, and Span{e i |p + 1 ≤ i ≤ n} forms an anti-euclidean space with signature −1.

Definition 4. Geometric product
The geometric product is the sum of the outer and inner products: In the case of a vector, the outer product is the antisymmetric part and the inner product is the symmetric part.
More generally, for an r-blade, A r , and an p-blade, B p , we have : A r .B p = (−1) r(p+1) B p .A r .
Any combination made of non-linearly independent vectors vanishes because of the anticommutativity of the outer product.In this way, Cℓ p,q can be represented as a direct sum of i-blade independent sub-spaces: The highest blade element has rank p + q = n and is unique and squares to −1; it is named the pseudoscalar.The rank n of the pseudo-scalar is interesting as it specifies the dimension of the vector space over which the algebra is defined.
Moreover, the number of independent blades of rank k ∈ {0, . . ., n} is given by the binomial coefficient Meaning that the generic CA Cℓ p,q , with p + q = n, globally contains 2 n independent blades.This last fact shows that one can look for a correspondence between the 2 n blades with a system of n Fermions (each Fermion possessing two degrees of freedom) or n qubits; this will be the object of the following paragraphs.

Mother Algebra and Witt Basis (WB)
The mother algebra is defined from the Euclidean vector space, V, and its dual V * [31]: The corresponding CA is noted as R n,n , called the mother algebra.If w i is a basis of V and w * i is a basis of V * , then R n,n is the algebra generated by the blades: The inner product on R n,n is defined as follows: Therefore, this set of vectors, {w i , w * i |1 ≤ i ≤ n}, forms a non-orthonormal basis of the mother algebra.This basis is called the Witt basis (WB).
We can use the following reversion/conjugation properties to transform a blade into its dual when using the WB:

Definition 5. Reversion/Hermitian conjugation
On a real Clifford space, the Hermitian conjugation is the anti-automorphism of C p,q , which maps vectors to their dual: The WB verifies the relations: w i w j + w j w i = 0 (15) An orthonormal basis can be defined from the WB as: giving e 2 i = 1 and φ2 i = −1.Therefore, the mother algebra can be decomposed as a tensor product of two CA, one Euclidean, and the other anti-Euclidean: where R n = Cℓ(e 1 , . . ., e n ) and R n = Cℓ( φ1 , . . ., φn ).Here, the e i are Euclidean and the φi are anti-Euclidean vectors.

Witt Basis as a Complex CA
The previous properties can be interpreted as a complex structure of the mother algebra [28,30].Indeed, by defining the vectors ϕ j = i φj , where i is the imaginary unit, we obtain the complex CA: with R ′ n = Cℓ(ϕ 1 , . . ., ϕ n ), where now the vectors ϕ i are Euclidean.In this case, the WB can be written as: As noted in [31], the anticommutation properties ( 14) to (19) of the WB are analogous to those of the creation operator, c † , and annihilation operator, c, used in quantum physics for fermions, used to write Hamiltonians in the second quantization form.
For this reason, it will be more convenient in the following to use the fermion notation:

Spinors for Quantum Computing
The WB can be used to build spinor spaces and also to represent qubit states.

Definition 6. Spinor space
A spinor space is a minimal left ideal of the CA generated by a hermitian idempotent element.It can be written using the mother algebra R n,n of the WB [31]: where I ∈ R n,n , such that I † = I and is an idempotent I = I 2 .
The resolution of the identity is obtained for the entire vector space by using the complete set of idempotent operators: This means that the entire algebra can be decomposed into a direct sum of 2 n spinor spaces.
Quantum computations with n qubits can then be described using these complex CA R n,n by means of the WB, as outlined in [29]; an n-qubit state vector becomes an element of a spinor.
As shown above in (34), the spinor space can be represented by one of the 2 n projectors I i 1 . . .I i t K k 1 . . .K k s .In quantum physics, the natural choice is the empty spinor space, where the projector I = I 1 . . .I n represents the vacuum state |0 . . .0⟩.
Using the creation operators, c † i , and the annihilation operators, c i , we have: In this way, the n-qubit state can be represented in this spinor formulation : In a more general formulation [32], one could equivalently consider the projector I i 1 . . .I i t K k 1 . . .K k s as the original spinor, corresponding now to a qubit state where i ℓ qubits are in state |0⟩ and k ℓ qubits are in state |1⟩.In this paper, we will adopt the representation of (36).

Second Quantization
Let us consider a system made of N fermions.The state of the system is described by the state vector where x j counts as the number of particles in the state j.Each x j is either 0 or 1 due to the Pauli exclusion principle.The annihilation and creation operators, c j and c † j , give the following relations when applied to |ψ⟩: and they correspond to removing or adding a particle in the state j.These operators are not hermitian, so they do not correspond to observable quantities.The number operator is hermitian and returns the numbers 0 for no-particle and 1 if the particle is present in the state j.
The creation and annihilation operators verify the following properties: The two first equations give zero because of the Pauli exclusion principle (i.e., x j ∈ {0, 1}).The last equation also signifies the closure that Moreover, the annihilation/creation operators acting on different states always anticommute ∀i ̸ = j, It is clear that one again finds the same commutation properties as in the case of the WB discussed above in the framework of GA, but here we started exclusively from the quantum theory; this analogy between the two methods justifies our approach.
This representation is called second quantization; it is useful for modeling the Hamiltonians of molecules, where the states represent the possible orbital states for electrons and their spins.This representation can then be used to simulate the molecule by translating this Hamiltonian into quantum gates and applying a quantum algorithm to extract the orbital energies.
Creation and annihilation operators act on fermion states, which is why they anticommute (see (45)), while quantum gates act on qubits, which are represented by wires in a quantum circuit and behave as bosons whose creation and annihilation operators, a † i and a j , commute.Therefore, fermion creation and annihilation operators and qubit gates do not correspond to the same mathematical structure; this is the reason why a transformation is required to pass from one framework to the other [17].

Jordan-Wigner Transformation (JWT)
The Jordan-Wigner transformation (JWT) uses the following operators [17,33]: where the σ + j and σ − j are spin operators related to the 1-qubit Pauli gates X j and Y j by The inverse transformation is given by: One verifies that for even products of these operators one finds again the usual 1-qubit projection operators: The inverse Jordan-Wigner transformation therefore enables one to translate the Hamiltonian of the second quantized form, written with the c † i and c i operators into a quantum circuit.

Permutation Properties of the JWT Operators
When using the WB formalism for the elementary 1-qubit operators c j , c † j , c j c † j , c † j c j , the permutation of the operators and the qubits states is ruled by the following property: where |λ j | is the number of Witt basis vectors contained in λ j , (i.e., Proof.The aim here is to exchange the lambdas and the kets thanks to the anticommutation properties of the WB (44).Then, we only have to prove that Depending of the values of λ j−1 and x j−1 , the product λ j−1 x j−1 takes the values presented in Table 1.
Table 1.Expression of the two-term product λ j−1 x j−1 for all possible factors.In the left-most column are the four possible left factors, λ j−1 , in the product, and in the upper-most line the two possible right factors, x j−1 , in the product.right factor left factor The WB vectors with different indices always anticommute; this means that, in the case of λ j ∈ {c j , c † j }, λ j and x j−1 anticommute if x j−1 = 1 and commute if x j−1 = 0.In the case of λ j ∈ {c † j c j , c j c † j }, λ j always commutes with x j−1 .This can be written as (55).The formula (54) can then be found by exchanging the λ j successively with all the previous |x k ⟩.
Let us recall that the operator n j given in (40) acts on the qubit state as The above operator can be written using the property (−1) ∑ j−1 k=1 n k = ∏ n j=1 (1 − 2c † j c j ).Then, using (54), we finally get

Reformulation of Quantum Gates
This result enables one to directly express quantum gates in the GA formalism using the JWT.The problem of sign change due to anticommutativity, which was thoroughly discussed in [29] for the translation of the tensor product of qubits in a GA framework, is here directly encoded inside the JWT.
The JWT applied to the basis matrices leads to: which can be easily inverted: Theorem 1. JWT for 1-qubit gates The Jordan-Wigner transformation using the Witt basis provides the following expressions for the Pauli 1-qubit gates: This means that, if ∀j, Λ j ∈ {X j , Y j , Z j }, then we have (66)

Two-Qubit Systems
Let us consider the case of a 2-qubit system; this means that we work in the CA R 2,2 .The WB elements reduce to (w * 1 , w 1 , w * 2 , w 2 ), here noted as (c 1 , c † 1 , c 2 , c † 2 ).According to the preceding spinor formulation, the qubit states are represented by the following multivectors: This shows that for 2-qubit states, the Kronecker product is equivalent to the geometric product in R 2,2 .In order to introduce the 1-qubit gates in a 2-qubit circuit, we use the JWT in R 2,2 .For the gates acting on the first qubit, the JWT simply returns: For the gates acting on the second qubit, however, the JWT adds a correction factor, to the X and Y gates: General 2-qubit gates can be written using the sums and products of the above 1qubit gates.For instance, for the controlled-NOT, controlled-Z, and SWAP gates, we have, respectively: giving in our formulation in R 2,2 :

Controlled Gates
In particular, one can derive the expression of controlled gates in R n,n .Let A k be an arbitrary 1-qubit quantum gate acting on the qubit k.The controlled gate C(A k ) j (where j is the control qubit and k the target) is given by:

Electronic Hamiltonians for Quantum Simulation
The aim of quantum simulation is to analyze the electronic structure and the associated energies of a molecule.The Born-Oppenheimer approximation (BOA) is commonly used where all nuclei of the molecule are assumed to stay at fixed positions since nuclei have a mass larger than electrons by a factor of 10 3 .Therefore, the dynamics of the molecule is described by the electronic Hamiltonian, H, which is the sum of a single electron Hamil-tonian, H (1) , and a two-electron Hamiltonian, H (2) .The wave function corresponding to the different energy levels, E (eigenvalues), of the molecule is an eigenfunction |ψ⟩ of the operator H given by the Shrödinger equation: The ground state, ψ g , corresponds to the lowest energy, E g .To describe electrons in the second quantization, we use a basis of N orbital states, noted as χ j .Each wave function, χ j , is a single electron wave function and is the tensor product of a spatial component, |ϕ⟩, and a spin component, |↑⟩ or |↓⟩.
In the second quantization formalism, the single-electron and two-electron parts of the BOA Hamiltonians are based only on the relative electron dynamics and are expressed as: The first coefficient, h ij , is the single-electron overlap integral; it consists of the kinetic energy and the potential energy for the one electron interaction from state i to state j.The coefficients h ii correspond to the unperturbed self-energy of the state, |χ i ⟩.
The second coefficient h ijkl is the two-electron overlap integral which corresponds to the different interaction energies of a system of 2 electrons initially in states i and j ending up in states k and ℓ.The two-electron interaction comprises : the Coulomb (i = ℓ, j = k) and exchange (i = k, j = ℓ) interactions which represent respectively the interactions between electrons and between electrons and nuclei, the double excitation operator (i ̸ = j ̸ = k ̸ = ℓ) corresponds to the excitation of both electrons in all four orbital exchanges.
The coefficients h ij and h ijkℓ for a given set of orbitals, { χ j }, can be computed classically by solving the Hartree-Fock equations [17].

Hamiltonian Operators in the WB
In order to simulate quantum systems using the WB, we express the different contributions to the molecular Hamiltonians; these are shown in Table 2 [20].

Operator
Second Quantized Form Energy number operator In quantum simulations, we are also interested in the exponential of these operators.We will compute them in the following sections.The decomposition of an exponential form is based on the following quantum gates: and The index corresponds to the position of the qubit affected by the operator.

Energy Number Operator
The energy number operator, H n i , gives the electron energy of the i-th state.It is the product of the number operator c † i c j by h ii , which is the mean value of H with respect to |χ i ⟩, as defined in Section 4.1.
The number operator is idempotent: Therefore, we have the exponential: The corresponding circuit is given in Figure 2.

Coulomb Operator
The Coulomb operators describe Coulomb interactions between two electrons in states i and j.They are written as: The product of the two number operators c † i c i and c † j c j is still idempotent: It is then possible to decompose the operator into the sum of involutive operators to obtain where Z j = c j c † j − c † j c j .Then, we have the exponential: besides, we have cos( and cos( ))(CNOT ij ) Thus, the exponential of the Coulomb operator is The corresponding circuit is given in Figure 3 i j R( Figure 3.Quantum circuit for exponential of the Coulomb operator.

Excitation Operator
The excitation operator has the form: Then, knowing that its exponential becomes In order to express the Hamiltonians using known gates, we will use the JWT.For the sin part, assuming i < j, we obtain For the cos part, we have which leads to the final form If j = i + 1 and h ij t = π 2 , this operator is the swap gate applied to the qubits i and j, and a phase factor on two states.More generally, it is related to the exponential of the swap gate, as shown on Figure 4.If j > i + 1, the phase inside the exponential can be transformed into its opposite depending on the state of the intermediate qubits due to the Z gates in the sin term of (99).
where the controlled T gates come from the phase factor on the states |10⟩⟨01| ij and |01⟩⟨10| ij .

Number-Excitation Operators
The number-excitation operators are: They obey the same kind of relations as the excitation operator: Similarly, the exponential gives: In this way, we obtain the gate of the excitation operator controlled by j: This last result is quite natural, since the number-excitation operator results from the combination of the number operator on qubit j and the excitation operator between qubits i and k.The number operator returns the number of qubits in the state j, which is either 0 or 1, so it is the same as controlling the excitation operator by qubit j. Figure 5 represents the corresponding circuit in the case where j = i + 1.

Double Excitation Operators
The double excitation operators are: They obey the same kind of relations as the excitation operator: Similarly, assuming i < j < k < ℓ, the exponential gives: If the condition i < j < k < ℓ is not verified, the Z gates do not simplify in the same way, but the global form of the operator remains the same; if h ijkℓ t = π 2 , similarly to the case of the simple excitation operator, we obtain a kind of 16 × 16 swap gate.Here again, it is possible to write the general form of (108) as a power of that matrix, with phase corrections applied to some states (either the 2 exchanged states or the 14 others).This matrix exchanges two composite sates (not two qubits).necessitating the intermediate stage employing the JWT (or the BKT).In our method, the fermionization obtained by the JWT is built in in the GA formalism.

GA Representation of the Quantum Ising and XY Hamiltonians
The quantum Ising Hamiltonian serves as a versatile and widely studied model in quantum simulation, providing valuable insights into the behavior of complex quantum systems and serving as a test bed for developing new simulation techniques and exploring emergent phenomena.While it may not capture all the complexities of real-world quantum systems, its simplicity, universality, and experimental feasibility make it an interesting tool for quantum simulation research [13,34,35].
Here we adopt a general approach, inspired from [34,35], considering a general anisotropic transverse Z field XY Hamiltonian, where the quantum Ising Hamiltonian with transverse Z field is a particular case (see hereafter).Both Hamiltonians include a term known as the transverse field, which introduces quantum mechanical effects.The general anisotropic transverse field XY Hamiltonian is given by: The first sum accounts for the degree of anisotropy interaction energy, and the second one for the transverse Z field.γ is the anisotropic coefficient, and taking γ = ±1 leads to the transverse Z field Ising Hamiltonian, while taking γ = 0 leads to the isotropic XX model.The λ factor is the strength of the transverse Z field.Using the reverse JWT, we can write this Hamiltonian in the GA formalism as: Then, by taking e −iHt , the first sum gives the excitation operators.The nλ term reduces to a global phase shift, which is meaningless.The last sum gives the number operators.
In the second sum (including the γ factor), we recognize an excitation operator that undergoes the transformation U → X i UX i .Indeed, we have: Then, the quantum circuit can be written by using: where the e −iH ex i,i+1 ∆t has been calculated in (99).
The corresponding circuit diagram is given in Figure 6.

Conclusions
In this article, we have shown that Geometric Algebra (GA) using the Witt basis (WB) approach is equivalent to the Jordan-Wigner Transformation (JWT) commonly used in quantum simulation.This shows that the formalism of GA can be used to perform operations on Hamiltonians expressed in the second quantized form and to translate them into quantum circuits.
The important issue of our method is to propose a direct computational method in quantum simulation using GA compared to other methods where the direct representations of the excitation operators using quantum circuits cannot be systematically derived, requiring work on each particular case.
GA provides a clear and intuitive framework for representing quantum states and operations and allows a geometric interpretation of quantum operations.This geometric perspective can provide deeper insights into the behavior of quantum systems and facilitate the design and analysis of quantum algorithms and circuits.Dirac notation is translated to GA in a natural way: all expressions are represented in a particular algebra and thus may be manipulated and implemented as algebra elements directly, without any need for matrix representation as is sometimes necessary within the Dirac notation formalism.
Further directions could be explored to improve the method, for example, by further investigating the relation between the geometric product of GA and the tensor product between quantum states in the context of quantum computing.Concerning the Witt Basis (WB) approach, the question rises: is there a more general physical meaning and why does the WB structure describe the fermion second quantization representation so well?
The advantages of GA in scaling compared to qubit methods may vary depending on the specific application and context.GA can lead to more efficient representations of quantum states and operations compared to traditional approaches, such as matrix representations.This reduced complexity can help mitigate the challenges associated with scaling quantum algorithms to larger systems, including computational overhead and resource constraints.Research in this area is ongoing, and further developments may uncover additional advantages of GA in quantum computing scalability.
While JWT has been widely used in quantum computing qubit gate models for simulating fermionic systems, it can introduce stability challenges due to its complexity and error-prone nature.Our approach naturally accommodates the representation of fermionic states and operators and thus offers a more direct, natural, and potentially more stable alternative for representing fermionic systems in quantum computing applications.Also, our method could be applied for a GA representation of the Bravyi-Kitaev transformation, which is also widely used in quantum simulation because it seems less consuming in quantum gate resources.
We plan in the future to investigate some algorithmic implementation of the method presented here, for example, for the determination of molecule ground states.

Figure 2 .
Figure 2. Quantum circuit for exponential of the number operator.

Figure 4 .
Figure 4.Quantum circuit for exponential of the Excitation operator assuming j = i + 1.

Figure 5 .
Figure 5.Quantum circuit for exponential of the excitation-number operator.

Figure 6 .
Figure 6.Quantum circuit for the quantum anisotropic Ising model simulation given in (117).