Quantum Computingvia Entanglement in Geometric Algebra Approach ()
1. Introduction. What Is Entanglement in Conventional Approach and in the Geometric Algebra Approach
In conventional quantum mechanics entanglement means that multiple objects share a single quantum state [1] [2] . Entanglement is the amount by which multiple objects share a quantum state. They remain indeterminate until they are disentangled by a measurement.
The simplest quantum mechanical state, qubit, reads:
,
,
It has just two observable “things” after measurement, say “up” for
and “down” for
, with probabilities
and
. The qubit dimension two is the number of different observable things available after making a measurement on the particle.
In the case of two particles the vector space
is generalized to density matrix defined on tensor product
. The appropriateness of tensor products is that the tensor product itself captures all the ways that basic things can “interact” with each other. Quantum entanglement means that aspects of one particle of an entangled pair depend on aspects of the other particle.
To understand the quantum entanglement, it is important to, first, consider quantum superposition. Quantum superposition is the idea that particles exist in multiple states at once. When a measurement is performed, it is as if the particle selects one of the states in the superposition. The quantum entanglement emerges from the reality of quantum superposition.
Consider the decay of a particle called the pi meson. When this particle decays, it produces an electron and a positron that have opposite spin and are moving away from each other. If the electron spin is measured to be up, then the measured spin of the positron could only be down, and vice versa. Because of quantum mechanics, the spin of each particle is both part up and part down until it is measured. Only when the measurement occurs does the quantum state of the spin “collapse” into either up or down—instantaneously collapsing the other particle into the opposite spin.
That’s crucial words: “instantaneously collapsing”. Common wisdom considers that as the physically real transformation of the other particle’s unknown state into the known one.
The scheme suggested in the geometric algebra approach is based on the manipulation and transfer of quantum states as operators acting on observables. Wave functions act in that context on static
elements through measurements, creating “particles.” [3]
Normalized wave functions as elements of
are naturally mapped onto unit sphere
(see Figure 1), [4] [5] . Two-state system is then just a couple of points on
, say
where
in some bivector basis
, with multiplication rules
,
,
.
Figure 1. Geometric representation of a
element.
Then it follows that two wave functions of an arbitrary two-function system are, in any case, connected by the Clifford translation1:
The product of exponents
is trivial in the case
(the case of an unspecified imaginary unit in conventional quantum mechanics,) namely:
. In the general case we have more complicated result, see (1.2.2) in [6] :
(1.1)
where
and
are vectors dual to planes
and
matching orientation of
.
Take a general definition:
The bivector part of the product
belongs to the plane orthogonal to the associated vector as seen in Figure 2.
To get the module
and angle of rotation in the plane of the bivector part of
we need the scalar part:
and representation of
as linear combination of the basis bivectors
:
Figure 2. Constructing the bivector plane of the product.
We get:
where
Thus, the module
is:
If
we also have usual formulas valid for complex numbers when formal imaginary unit i is replaced by
(or
.) Particularly, angles are added in multiplication, subtracted in division, and de Moivre’s formula takes place:
If
then
for positive integer n.
In general case, when
, the angle of rotation
is:
and the plane of the rotation normal has directional cosines:
The result of Clifford translation (1.1) is an
element. From knowing Clifford translation connecting any two wave functions as points
it follows that the result of measurement of any observable C by wave function
, for example
, immediately gives the result of (not made) measurement by
:
When assuming that observables are also identified by points on
and thus are connected by the formulas similar to the above one we get that measurements of any amount of observables by arbitrary set of wave functions are simultaneously available.
This is a geometrically clear and unambiguous explanation of strict connectivity of the results of measurements instead of quite absurd “entanglement” in conventional quantum mechanics.
2. Maxwell Equation in Geometric Algebra
Without charges and currents, the Maxwell equation is:
(2.1)
The circular polarized electromagnetic waves are the only type of waves following from the solution of Maxwell equations in free space done in geometric algebra terms. Indeed, let’s take the electromagnetic field in the form:
(2.2)
requiring that it satisfies (2.1)
Element
in (2.2) is a constant element of geometric algebra
and
is unit value bivector of a plane S in three dimensions, generalization of the imaginary unit [7] . The exponent in (2.2) is the unit value element of
:
,
Solution of (2.1) should be the sum of a vector (electric field
) and bivector (magnetic field
):
with some initial conditions:
Substitution of (2.2) into Maxwell’s (2.1) will show us what the solution looks like.
The derivative by time gives
The geometric algebra product
is:
or
,
depending on do we write
or
. The result should be the same because
is a scalar.
Commutativity
is valid only if
. The following agreement takes place between orientation of
, orientation of
and direction of vector k. The vector
is orthogonal to the plane of
and its direction is defined by orientations of
and
. Rotation of right/left hand screw defined by orientation of
gives the movement of right/left hand screw. This is the direction of the vector
. That means that the matching between
and
should be
(see Figure 3).
Assume first that orientation is
. Then Maxwell equation becomes:
or
Figure 3. Two possible orientations of the wave vector.
Left hand side of equation is sum of vector and bivector, while right hand side is scalar
plus bivector
, plus pseudoscalar
, plus vector
. It follows that both and h lie on the plane of
and then:
,
Thus,
and we get equation
from which particularly follows
and
.
The result for the case
is that the solution of (2.1) is
where
and
are arbitrary mutually orthogonal vectors of equal length, lying on the plane S. Vector k should be normal to that plane,
and
.
In the above result, the sense of the
orientation and the direction of k were assumed to agree with
. Opposite orientation,
, that’s k and
compose left hand screw and
, will give solution
with
.
Summary:
For a plane S in three dimensions Maxwell equation (2.1) has two solutions
·
, with
,
, and the triple
is right hand screw oriented, that’s rotation of
to
by π/2 gives the movement of right hand screw in the direction of
.
·
, with
,
, and the triple
is left hand screw oriented, that’s rotation of
to
by π/2 gives movement of left hand screw in the direction of
or, equivalently, movement of right hand screw in the opposite direction,
.
·
and
, initial values of
and
, are arbitrary mutually orthogonal vectors of equal length, lying on the plane S. Vectors
are normal to that plane. The length of the “wave vectors”
is equal to the angular frequency
.
Maxwell Equation (2.1) is a linear one. Then any linear combination of
and
saving the structure of (2.2) will also be a solution.
Let’s write:
(2.3)
Then for arbitrary (real2) scalars
and
:
(2.4)
is solution of (2.1). The item in the second parenthesis is a weighted linear combination of two states with the same phase in the same plane but opposite sense of orientation. The states are strictly coupled, entangled if you prefer, because the bivector plane should be the same for both, does not matter what happens with that plane.
Arbitrary linear combination (2.4) can be rewritten as:
(2.5)
where
,
The triple of unit value basis orthonormal bivectors
is comprised of the
bivector, dual to the propagation direction vector;
is dual to an initial vector of the magnetic field;
is dual to the initial vector of the electric field. The expression (2.5) is a linear combination of two geometric algebra states, g-qubits.
Linear combination of the two equally weighted basic solutions of the Maxwell equation
and
,
with
reads:
(2.6)
where
and
. It can be written in
standard exponential form
.3
I will call such g-qubits spreons because they spread over the whole three-dimensional space for all values of time and instantly change under Clifford translations over the whole three-dimensional space for all values of time, along with the results of measurement of any observable.4
3. Measurement of Observables by Spreons
Measurement of any observable
(actually Hopf fibration) by a state
in the current formalism is, see, for example, [6] :
In the case of spreon (2.6):
,
,
,
and we get a
element spreading through the three-dimensional space for all values of the time parameter t. It is
element spreading through the three-dimensional space for all values of time parameter t:
(3.1)
Geometrically, that means that the measured observable is rotated by π/2 in the
plane, such that the
component becomes orthogonal to the plane
and remains unchanged. Two other components became orthogonal to
and
and continue rotating in
with angular velocity
. The factor
defines the dependency of those transformed values through all points of the three-dimensional space.
The current approach transcends common quantum computing schemes since the latter are principally based on qubit entanglement (whatever it is) and thus have tough problems of creating large sets of entangled qubits. In the current scheme, any test observable can be placed anywhere into the continuum of the
dependent values of the spreon state. The above formula (3.1) gives the result of measurements simultaneously at all points
.
Let us get back to the geometric algebra entanglements explained at the end of Sec.1. Assume we have two observables,
and
. Write them in exponential form:
where
,
,
,
,
where
,
,
,
,
The observable
can be written as:
,
and then its measurement by any g-qubit
reads:
that is, up to the factor
, the result of measurement of
, multiplied by the result of measurement of
.
4. Software Simulation of Analog Computing
The hardware creating superfields requires special implementation of a photonic/laser device that does not exist yet. Instead, we have a very convenient equivalent simulation scheme where the amount of simultaneously available space/time points of observable measured values is only restricted, for example, by the overall available Nvidia GPU number of threads.
Formula (3.1) gives the result of measuring any type of
element
by the sprefield state.
Let us consider pieces of code for possible implementation.
The CUDA code begins with allocating an initial array of (3.1) on the host at
. For example:
size_t memsize = 0;
struct cudaDeviceProp deviceProp;
checkCudaErrors(cudaGetDeviceProperties(&deviceProp, 0));
memsize = deviceProp.totalGlobalMem;
// We can never use all the memory so to keep things simple we aim to
// use around half the total memory
memsize /= 2;
int defaultDim = (int)floor(pow((memsize / (2.0 * sizeof(float))), 1.0 / 3.0));
// above is the size of cube for all considered points
float step = volumeSize.width / defaultDim;
float4* r = (float3*)malloc(memsize * sizeof(float4));
r[0].x = r[0].y = r[0].z = r[0].w = 0.0;
for (int k = 0; k < defaultDim; k++)
{
r[k].x = r[0].x + k * step;
r[k].y = r[0].y + k * step;
r[k].z = r[0].z + k * step;
}
float omega = 12560000.0; // an option from laser beams
float* factor = (float*)malloc(defaultDim * sizeof(float));
for (int k = 0; k < defaultDim; k++)
{
// below from formula (3.1)
factor[k] = 4. * (cosf(omega * r[k].x)) * (cosf(omega * r[k].x));
r[k].x = factor[k] * C3;
r[k].y = factor[k] * C2;
r[k].z = -factor[k] * C1;
r[k].w = factor[k] * C0;
}
Now, copy to the device memory:
Float4* d_r = NULL;
checkCudaErrors(cudaMalloc(&d_r, memsize));
checkCudaErrors(cudaMemcpy(d_r, r, memsize, cudaMemcpyHostToDevice));
Run the CUDA part of calculations and graphics output:
// create VBO
createVBO(&vbo, &cuda_vbo_resource, cudaGraphicsMapFlagsWriteDiscard);
// run the cuda part
runCuda(&cuda_vbo_resource);
// start rendering mainloop
glutMainLoop();
with specially written kernel function
vbo_kernel <<< grid, block >>> (pos, mesh_width, mesh_height, mesh_depth, time);
running inside the runCuda. The vbo_kernel executes in the parallel recalculation of all measured observable values sitting in the cube, changing with time and their graphics output.
5. Conclusions
The geometric algebra lift of conventional quantum mechanics qubits is the game-changing quantum leap forward potentially kicking from the quantum computing market big fishes (IBM, Microsoft, Google, dozens of smaller ones) investing billions in elaborating quantum computing devices. The approach brings into reality a kind of physical field spreading through the whole three-dimensional space and values of the time parameter. The fields can be modified instantly in all points of space and time values. All measured observable values are simultaneously available all together, not through looking one by one. In this way, the new type of quantum computer appeared to be a kind of analog computer keeping and instantly processing information by and on sets of objects possessing an infinite number of degrees of freedom. As a practical implementation, the multithread GPUs with the CUDA language functionality allow the creation of software simulating that kind of field processing numbers of space/time discrete points only restricted by the GPU threads capacity.
Specific problems that can get huge speedups with the suggested analog quantum computers are mainly related to the continuous media calculations and blockchain schemes [8] .
Thus, the innovative scheme of quantum computing with its availability in the near future, incomparable low cost value and scalable simulation software platform is becoming a tremendous challenge for companies like those mentioned above.
The corresponding anticipated objective is to create simulating software running on a GPU, preferably programmed in the Nvidia CUDA language. In that way, we will get a quantum computer as an analog computer simulating instantly parallelizable processing, and transformation of states identified by points on the three-dimensional spheres, that replace the Hilbert space formalism. The amount of simultaneously available space/time points of the observable measured values is then only restricted by the overall available GPU number of threads.
The simulating software functionality will follow the requirements of further developing quantum computing in all areas where its paradigm will give tremendous effectiveness, first of all in speed increasing by orders, from running years to seconds: material science, finances, drug development, cryptography, weather forecasting, to mention just some immediately seen ones.
NOTES
1It is universally possible due to the hedgehog theorem.
2Remember, in the current theory scalars are real ones. “Complex” scalars have no sense.
3Good to remember that the two basic solutions F+ and F− differ only by the sign of I3IS, which is caused by orientation of IS that in its turn defines if the triple
is right-hand screw or left-hand screw oriented.
4The two received solutions are similar in their form to the Majorana operators, though with a bivector instead of formal imaginary unit. Thus, not surprising that the following computational scheme represents more general software simulation which gets around the problem of physical implementation of Majorana fermions.