CC BY
1
Mathematical modeling of processes in quantum computer elements based on methods of quantum theory to improve their efficiency
A.A. Biryukov1, M.A. Shleenkov2 1 Samara State Transport University, 443066, Samara, Russia, Svoboda str. 2V;
2 Samara National Research University, 443086, Samara, Russia, Moskovskoe shosse 34
Abstract
The paper studies entangled states of two qubits interacting with each other and with an electromagnetic field. The state of the qubits is determined by a statistical density matrix. The degree of entanglement of the state is characterized by the Peres-Gorodeckii (PG) parameter. The statistical density matrix and its evolution are determined in the energy representation within the framework of the path integral formalism. The obtained equations determine the dependence of the PG parameter on the parameters of qubit dipole-dipole interaction and the acting electromagnetic field. The results of numerical calculations are presented in graphs for the PG parameter. It is shown that it is possible to choose parameters corresponding to qubit states with a high degree of entanglement (0.99).
Keywords: qubits, quantum entanglement, Peres-Gorodeckii parameter, path integration.
Citation: Biryukov AA, Shleenkov MA. Mathematical modeling of processes in quantum computer elements based on methods of quantum theory to improve their efficiency. Computer Optics 2024; 48(6): 889-893. DOI: 10.18287/2412-6179-CO-1489.
Introduction
Richard Feynman in 1982 proposed constructing computers from quantum mechanical elements as a means of improving computer and communication device performance [1]. Since then, active efforts have been made to create computers and other information devices based on quantum elements. The main components of a quantum computer are systems consisting of two identical qubits. Entangled states of two qubits are of fundamental importance in the operation of a quantum computer and in quantum information science (cryptography [2], teleportation [3]). However, in real constructions, qubits in an entangled state can only be maintained for a limited time. Therefore, active theoretical and experimental research is currently underway to find conditions that can prolong the time of quantum entanglement of qubits.
Quantum entanglement is important for understanding quantum theory and the prospects of its applications. Fundamental principles and primary methods of creating and detecting entangled states in various quantum systems can be found, for example, in the review [4]. Several dimensionless parameters have been introduced for quantitative measurement of quantum entanglement, such as the PG parameter [5], negativity [6], and state entropy [7].
In reference [8], the entangled states of cold ions confined in a linear trap and interacting with laser beams are discussed. It is shown that the decay of the entangled state in this system is negligible, and quantum register measurement (quantum reading) can be performed with high efficiency. Reference [9] considers the stabilization of the entangled state of two-level atoms interacting with a classical laser field with the consideration of spontaneous emission. The exact solution obtained shows that the optimal entanglement that can be achieved in the
Lamb-Dicke limit is 0.43. Reference [10] considers the decay of quantum entanglement of two-level systems, each interacting with a thermal reservoir at finite temperature. For a wide class of initially entangled states, the authors show that the system always disentangles in finite time.
References [11, 12] theoretically investigate the effect of dipole-dipole interaction on quantum entanglement and demonstrate the potential use of dipole-dipole interaction as a tool for controlling the degree of quantum entanglement. Reference [13] investigates the influence of chaotic behavior of coherent state parameters on the dynamics of entanglement in a system of two-level atoms and photons in a cavity. Reference [14] shows that the presence of a high-intensity laser field leads to a high degree of entanglement of two qubits. The paper [15] explores within the framework standard Tavis-Cummings model, entanglement of qubits under the influence of a non-stationary single-mode field and a Kerr-like medium. The work [16] examines the entanglement between two superconducting qubits in which there is a dipole - dipole interaction and interaction with the thermal field. The authors of article [17] analyze the dynamics of entangled states of qubits interacting with the external environment using the method of correlation functions. The paper [18] describes original experiments in which two molecules were presented in an entangled quantum state. A review of research on entangled qubit states allows us to conclude that the influence of an external electromagnetic field and dipole-dipole interaction between qubits is a promising tool for controlling the degree of quantum entanglement.
The article [19] notes the extraordinary progress in the modern development of quantum computing technologies, which indicates the prospect of widespread use of quantum computers.
In this work, we construct a mathematical model describing two interacting qubits with an external electromagnetic field. The model allows us to determine the PG entanglement parameter and its dependence on the qubit interaction parameters among themselves and with the external electromagnetic field. The density matrix and its evolution are determined in the energy representation within the framework of the path integral formalism. This formalism allows interactions of any intensity to be considered.
Mathematical model of two qubits interacting with each other and an electromagnetic field
We will construct a mathematical model of a system in which two qubits interact with an external electromagnetic field and with each other, and determine the degree of entanglement of the qubit states as a function of the interaction parameters. By varying these parameters, we can alter the degree of entanglement of the qubits.
The Hamiltonian of this model is given by:
H Ml = H q + F(x),
F(x) = VQF (x) + VQQ,
(1)
where Hq - Hamiltonian of two non-interacting qubits; VQf(x) - interaction operator between qubits and the electromagnetic field; Vqq - dipole-dipole interaction operator of qubits with each other.
We describe the studied system by the statistical operator p in the interaction representation [20]:
p(t ) = U d (t )p(0)UD (t),
(2)
where
U d (t ) = fexp
. t
- I jVD, (x)d x
VD (x) = exp
IHQ x
n
F(x)exp
- ! HQ X n
where T is time-ordering operator.
The explicit form of the statistical density matrix will be calculated in the energy representation [21]. To achieve this, we will take the eigenstates of the two-qubit Hamiltonian as the basis with the following properties:
Hq | n, m) = Enm | n, m), i
% | n, m)<m, n |= 1,
n,m=0
<m',n | n,m) = Sn„Sm„,
(3)
where En, m En + Em is the energy of the noninteracting qubit system; |n, m) = |n)®|m) are the eigenstates of the two-qubit Hamiltonian, they describe the state of the system; each eigenvector |n), |m) takes two values, i.e., n = 0,1, m = 0,1. In this basis, equation (3) takes the form of equation (2).
i
i
Pmfmf (t) = % IL < mfnf | UD(t)| nm )P,
"in ,"in =0min ,min=0
minninninmin
(0)<rninriin | UD (t) | n m f ),
(4)
where
Pm/n/n>m> (t) = <mf , nf | P(t) | nf, mf ), Pm n n m (0) = <mn, nn | p(0) | rim, mn ).
in in in in
The kernel of the evolution operator is computed based on the summation over all possible trajectories in the energy space:
<mf, nf | Ud (t)| nin, mn ) =
i i K
= % % n<m*,nk | Ud(tk,4-i)| nk-i,mk-i),
n0,--,nK-i=0m).--'mK-i=0 k=i
(5)
with tk> tk-i and n0 = nin, m0 = miri at the initial time fe = 0 and nK=n/, mK = m/ at the final time tN = t > 0.
The computation of the matrix elements of the evolution operator <n/, m/1 ^D(t)|n,n, m„) based on equation (5) will be more accurate the smaller the time intervals tk - tk-i we consider.
Based on equations (4), (5), we will establish recursive relations for the matrix elements of the statistical operator at different time points:
(tk) =
mknknkmk i
= % % <mk, nk | Ud (tk - tk-i) | nk-i, mk-i )* (6)
nk-i"-i=0mk-i.mk-i=0
(tk -i)< mk-i, nk-i | U D (tk - tk-i) | "m ),
XP,
mk-i"k -i"k-imk-i
where the statistical operator is normalized at each time moment, i.e., Spp(t) = 1.
Equations (4 - 6) represent an algorithm for calculating the statistical density ma"tnx Pn/m/n'/m'/ (t) at any time moment t.
For small (tk - tk-i) the evolution operator UD(tk - tk-i) according to (2) can be represented as
^D (tk, tk-i) =
= i -1 exp
7 H Q x h
V(x)exp
- 7 HQ x h
At,
(7)
where At = tk - tk-i, x = 1/2(tk + tk-i).
Based on (3), (7), we can express the kernel of the evolution operator (8) as:
<mk, nk ,| Ud (tk, tk-i) | "k-i, mk-i) ~ S„, n . Sm m . +
"k"kmk ,mk-i
<mk,"k | V>(x) | «k-1, mk-i) exp (ra"k«k-i +ramkmk.
-i
(8)
where ftra„k„k-1 = E„k - E„k-1
i.e., ra«k"k-i, ®mkmk-
Irn = E - E
mkmk-i mk mk-i
are the transition frequencies of the qubits between their states.
It is clear that the adjoint matrix of the evolution operator is also constructed in the same way
<mk_i«k_i | UD (tk - tk_i) | nkmk ) .
Note that for physical models, constructing the evolution operator (9) in explicit form for a small time interval (tk- t—) is much simpler than constructing the evolution operator {mk, nk | Uo(tk, tk-i)\nk-i, mk-1) for a finite time interval (tk - tk-1).
The explicit form of the statistical matrix is constructed based on the recurrence relation (6). For this purpose, we specify an explicit form of the statistical matrix at the initial moment of time Pm,nn,nn'inm',n(0), a transition matrix (9) over a small time interval (tk - tk-1 = t) and in accordance with equation (6) we obtain the statistical operator Pm,nn,nn'inm',n(t) at the moment of time t (with the normalization condition Spp(t) = 1). Applying recurrent relation (6) N times, we will construct an explicit form of the statistical operator pmffm'/X), for a point in time
N
t = X(t* - tk-1), t0=0, tN = t .
k=1
This method is convenient for numerical calculations. Calculation accuracy increases with shorter time intervals
(tk - tk-1).
To quantitatively assess the quantum entanglement of the two qubits, we use the PG parameter [3, 4] with the measure e:
s(t ) = £(| X, (t ^ (t)),
(9)
where Xi - are the eigenvalues of the following equation:
det((mf,nf | p(/) | nf',mf')pT -XI) = 0, (10)
where I is the identity matrix; the index pT denotes matrix partial transformation. The entanglement is maximal when e = 1 and minimal when e = 0.
The explicit form of interaction potentials Vqf(x), Vqq is constructed based on the physical model. The proposed system of equations (6), (8), (9), (10) allows one to numerically generate plots depicting the dependence of e on the qubit interaction parameters and parameters of the external electromagnetic field.
Numerical modeling of entangled states of two qubits
Based on the proposed mathematical model, we describe the entangled states of two identical qubits in a specific physical model.
Consider a system of two qubits, in which the interaction operators are represented by the following expressions.
The operator Vqq of qubit dipole-dipole interaction:
Vqq = gdxd2 ,
(11)
where dx, d2 are the operators of the dipole moments of the first and second qubits, respectively, and g is the interaction constant.
The operator VQF (x) of qubit interaction with a single-mode electromagnetic field with amplitude E0 and frequency Q:
VQF (x) = djEo cos (Qt + ) + <i2 E0cos (Qt + ^2), (12)
where ^1, ^2 - are the phases of oscillations of the electromagnetic field acting on the first and second qubits, respectively.
The matrix elements of the evolution operator (8) considering (1), (11), (12) take the form
<mk, nk | UD (tk, tk-1) | nk-1, mk= 8„k ^ ^ ,mk-1 + +j—g<mk,nk | dd21 nk-1,mk-)exp[i (co„k„k_j + romkmk-1 )x]^At-
+j(-E0
{mk, nk | d^ nk-1, mk-1) cos (Qr + ^ )exp [i (ra„k„k_j +®mkmk_x ) + {mk,nk | d21 nk-1,m^1)cos(Qx + ^)exp[i(„^ +ramkmk-1 )x]Ut.
(13)
We also construct the expression for
{nk-1mk-11 UD (tk - tk-1) | n'kMk ) .
Let's number our basis vectors in the following order:
|nk -1, mk-1 )^|0,0),|0,1),|1,0),|1,1),
(
{mk, nk | Ud (At) | nk-1, mk-1) =
Q R At -i-e
Q R At -i-e
(14)
Q R At
{mk,nk |^{0,0|,{1,01,{0,11,{1,11.
(15)
In this basis representation, we express the matrix of the evolution operator with elements (13) explicitly for the case when Q = ra:
-iQQ At
Q R At -i-e
Q R At
-iQQ At
-i ^ e-«2
0
. QR At + -i-e
. Qr At + -i-e
(16)
where At = tk-4-i, Qr = dEo/ h is the Rabi oscillation frequency, Qq=gd2 / h is oscillation frequency for qubit interaction, d is absolute value of qubit dipole moments (nk |d|nk-l) = dl(1 - S„k, nk-i), (mk \d2\mk-l) = d2(1 - 8m*, mk-l), absolute values of dipoles are equal to each other d1= d2 = d. We took into account ranm > 0 for n > m, ranm < 0 for n < m and ranm = - ranm. Terms with frequencies Q + ra, 2raoi, 2raio are neglected when we construct (16) in rotating wave approximation.
The corresponding structure has the adjoint matrix. In our modeling, we assume = ^2 = 0.
Dimensionless times, t
Fig. 1. The dependence of the PG parameter on the dimensionless time parameter
The fig. 1 represents the results of numerical modeling of the dependence of the PG parameter e of our model on the dimensionless time parameter x = Qo-t and interaction parameters Qr and Qq. Where Qo is constant that has the dimension of frequency and must be specified when considering specific experimental conditions.
The first numerical experiment was conducted with the qubit dipole-dipole interaction enabled (Qq = 0.05-Qo) and in the absence of an external electromagnetic field (Qr = 0.00-Qo). The graph depicting the dependence of the PG parameter on time is shown as a dashed line. The results indicate continuous oscillations of qubit entanglement in this case, as the PG parameter e oscillates from the minimum value of 0 to the maximum value of 1.
The second numerical experiment was conducted with the qubit dipole-dipole interaction enabled (Qq = 0.05-Qo) and in the presence of an external electromagnetic field (Qr = 7.50-Qo). The results indicate oscillations of qubit entanglement, which stabilize at a state of high entanglement when the PG parameter takes a value of around 0.99. Note that an external field with a smaller amplitude changes the value of the entanglement parameter, but does not establish its stable value.
In the third experiment, the external electromagnetic field (Qr = 7.50-Qo) is turned on at the moment x = 0 and turned off the moment x = 100. Turning off the field leads to a decrease parameter e and the destruction of the entagled state. A graph of this behavior e over time is presented with a dash-dotted line.
Conclusion
The proposed mathematical model allows for the study of the entanglement parameter e of a two-qubit system for different interaction potentials.
It is shown that the interaction between qubits leads to the entanglement of their states. Additional interaction with an external field leads to the stabilization of the entangled state.
The proposed mathematical model allows for the exploration of entangled states of qubits under the influence of external fields of any structure.
References
[1] Feynman RP. Simulating physics with computers. Int J Theor Phys 1982; 21: 467-488. DOI: 10.1007/BF02650179.
[2] Klobus W, Grudka A, Baumgartner A, Tomaszewski D, Schonenberger C, Martinek J. Entanglement witnessing and quantum cryptography with nonideal ferromagnetic detectors. Phys Rev B 2014; 89: 125404. DOI: 10.1103/PhysRevB.89.125404.
[3] Torres JM, Bernad JZ, Alber G. Quantum teleportation and entanglement swapping of matter qubits with coherent multiphoton states. Phys Rev A 2014; 90: 012304. DOI: 10.1103/PhysRevA.90.012304.
[4] Horodecki R, Horodecki P, Horodecki M, Horodecki K. Quantum entanglement. Rev Mod Phys 2009; 81: 865-942. DOI: 10.1103/RevModPhys.81.865.
[5] Peres A. Separability criterion for density matrices. Phys Rev Lett 1996; 77: 1413-1415. DOI: 10.1103/PhysRevLett.77.1413.
[6] Vidal G, Werner RF. Computable measure of entanglement. Phys Rev A 2002; 65: 032314. DOI: 10.1103/PhysRevA.65.032314.
[7] Bashkirov EK, Sochkova EYu. Entanglement in two-atom model with degenerate Raman transitions. Journal of Samara State Technical University, Ser Physical and Mathematical Sciences 2011; 2: 135-141. DOI: 10.14498/vsgtu934.
[8] Cirac JI, Zoller P. Quantum computations with cold trapped ions. Phys Rev Lett 1995; 74: 4091-4094. DOI: 10.1103/PhysRevLett.74.4091.
[9] Cakir O, Klyachko AA, Shumovsky AS. Steady-state entanglement of two atoms created by classical driving field. Phys Rev A 2005; 71: 034303. DOI: 10.1103/PhysRevA.71.034303.
[10] Al-Qasimi A, James DFV. Sudden death of entanglement at finite temperature. Phys Rev A 2008; 77: 012117. DOI: 10.1103/PhysRevA.77.012117.
[11] Bashkirov EK, Mastyugin MS. Entanglement of two flux qubits interacting with thermal field. Proc SPIE 2014; 9031: 903110. DOI: 10.1117/12.2051291.
[12] Bashkirov EK, Stupatskaya MP. The entanglement of two dipole-dipole coupled atoms induced by nondegenerate two-mode thermal noise. Laser Phys 209; 19: 525-530. DOI: 10.1134/S1054660X09030281.
[13] Gorokhov AV. Chaos and entanglement in atomic systems interacting with photons. Proc SPIE 2013; 8699: 86990Y. DOI: 10.1117/12.2019110.
[14] Soltani M, Ezatabadipour H, Jalali J, Darabi P, Azizi E, Rashedi G. Control of entanglement between two dissipative non-interacting qubits in a common heat bath by a laser field. Eur Phys J D 2013; 67: 256-264. DOI: 10.1140/epjd/e2013-40350-8.
[15] Abdel-Khalek S, Berrada K, Khalil EM, Obada A-SF, Reda E, Eleuch H. Quantumness measures for a system of two qubits interacting with a field in the presence of the time-dependent interaction and kerr medium. Entropy 2021; 23: 635. DOI: 10.3390/e23050635.
[16] Bashkirov EK. Entanglement of two dipole-coupled qubits induced by a detuned thermal field. J Phys: Conf Ser 2018; 1096: 012063. DOI: 10.188/1742-6596/1096/1/012063.
[17] Bratus E, Pastur L. The dynamics of quantum correlations of two qubits in a common environment. J Math Phys Anal Geom 2020; 16(3): 228-262. DOI: 10.15407/mag16.03.228.
[18] Holland CM, Lu Y, Cheuk LW. On-demand entanglement of molecules in a reconfigurable optical tweezer array. Science 2023; 382(6675): 1143-1147. DOI: 10.1126/science.adf4272.
[19] Ayoade O, Rivas P, Orduz J. Artificial intelligence computing at the quantum level. Data 2022; 7(3): 28. DOI: 10.3390/data7030028.
[20] Dirac PAM. Principles of quantum mechanics. 4th ed. Oxford University Press; 1982. ISBN: 978-0198512080.
[21] Biryukov AA, Degtyareva YV, Shleenkov MA. Calculating the probabilities of quantum transitions in atoms and molecules numerically through functional integration. Bulletin of the Russian Academy of Sciences: Physics 2018; 82(12): 1565-1569. DOI: 10.3103/S1062873818120055.
Authors' information
Alexander Aleksandrovich Biryukov (born in 1942) graduated from Saratov State University in 1966, completed his postgraduate studies at the same university in 1969, and defended his dissertation in 1970 in the field of Theoretical Physics to obtain the degree of candidate of Physical and Mathematical sciences. From 1969 to 2016, he held various positions at Samara State University, including Associate Professor, Professor, Dean of Physics faculty, and Head of the General and Theoretical Physics department. From 2016 to 2020, he worked as a professor at Samara National Research University. Since 2020, he has been serving as a professor at Samara State University of Transport. He has authored of more than 200 scientific and methodological publications covering quantum theory, quantum optics, stochastic processes theory, thermodynamics, and physics teaching methods. E-mail: [email protected] ORCID ID: 0000-0003-3955-1726.
Mark Aleksandrovich Shleenkov (born in 1988) pursued his studies at Samara State University (now known as Samara National Research University) and obtained his Ph.D. degree in Laser Physics in 2013. Currently, he holds the position of Associate Professor at the General and Theoretical Physics department at Samara National Research University. E-mail: [email protected] ORCID ID: 0000-0002-0930-4348.
Code o/State Categories Scienti/ic and Technical Information (in Russian- GRNTI): 29.05.15, 50.09.31 Received December 26, 2023. The final version - March 17, 2024.
