Prerpints.org logo
Preprint
Article

Entanglement Distillation Optimization using Fuzzy Relations for Quantum State Tomography

Altmetrics

Downloads

112

Views

42

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

01 June 2023

Posted:

02 June 2023

You are already at the latest version

Alerts
Abstract
Practical entanglement distillation is a critical component in quantum information theory. Entanglement distillation is often utilized for designing quantum computer networks and quantum repeaters. The practical entanglement distillation problem is formulated as a bilevel optimization problem. A fuzzy formulation is introduced to estimate the quantum state (density matrix) from pseudo-likelihood functions (i.e., quantum state tomography). A scale-independent relationship between fuzzy relations in terms of the pseudo-likelihood functions is obtained. The entanglement distillation optimization problem was solved using the combined coupled map lattice and dual annealing approach. Comparative analysis of the results is then conducted against a standard dual annealing algorithmic implementation.
Keywords: 
Subject: Computer Science and Mathematics  -   Computer Networks and Communications

1. Introduction

Entanglement distillation is a critical aspect of quantum information systems. The key idea of entanglement distillation is to preserve or restore the quality of diluted entanglement states of quantum information transmitted over large distances. Decoherence effects during transmission causes the dilution of entanglement states. Past theoretical and experimental research works have focused on the studies of quantum distillation frameworks [1,2,3]. Recently, pairs of single photons (entangled in multiple degrees of freedom) were used to experimentally determine the domain of distillable states and their relative fidelity [4]. In Ecker et al., (2021) [4] comparative studies were also done on various distillation schemes to gain a deeper understanding in terms of resilient quantum network design. A proof-of-concept experiment was recently conducted to study the application of filtering protocols (in atomic ensembles) for constructing quantum repeater nodes [5]. The experiment was conducted in a crystal (rare-earth-ion-doped). In that setting, the entanglement states were prepared. The relationship between bit thread, entanglement distillation, and entanglement purification (in the holographic framework) was recently studied [6]. In Lin et al., (2021) [6], a bit thread interpretation of the one-shot entanglement distillation tensor network was provided. It was shown that the holographic entanglement purification process could be viewed as a special case of a type of surface growth scheme. The objective of the study in Lin et al., (2021) [6] was to develop an accurate framework for describing physical entanglement structures. Another interesting work is seen in the theoretical investigation of entanglement distillability presented in Qian et al., (2021) [7]. In that work, the authors studied the mentioned subject with regards to the undistillability conjecture of specific Werner states - where they gained deeper understanding of the distillability problem. Recent developments of entanglement distillation are seen in Kondra et al., (2021) [8], He et al., (2021) [9], Gour and Scandolo, (2021) [10], Yan et al., (2022) [11], Shchukin and van Loock, (2022) [12] and Riera-Sàbat et al., (2021) [13].
In the past few years, much research works in quantum information systems have been directed towards optimization. This can be seen in the work of Gyongyosi and Imre (2020) [14]. In that work, the optimization of the read-out procedure of local unitaries of a high-retrieval efficiency quantum memory was done. The authors also studied the retrieval efficiency of quantum memory and the signal-to-noise ratio. A similar work is seen in Gyongyosi et al., (2020) [15]; where an approach for obtaining the optimal quantum state and computational path evaluation for gate-model quantum computing devices was proposed. Efficient algorithms for attaining generation time and transmission fidelity of entangled pairs between the end nodes of quantum chains is seen in Brand et al., (2020) [16].
Quantum state tomography (QST) is an effective method to reconstruct or estimate quantum states from measurements of identical quantum states. This directly impacts the development of reliable quantum resources as well as quantum devices for quantum information processing. A recent interesting study on information gain in QST is seen in the work of Sahu et al., (2022) [17]. In that work, the authors discuss the effects of chaotic dynamics on the information transfer during QST. Their investigations also uncovered operational insights into the mechanisms of fidelity gain during actual quantum information tomography protocols. In Schmale et al., (2022) [18], the authors improve the accuracy of observables measurement from tomographic data. This was done by using a QST scheme that approximates a probability distribution over a measurement (informationally complete) within a variational manifold characterized by a convolutional neural network. In that work, the authors regenerated high- classical fidelity states which performed more efficiently as compared to standard techniques (e.g., maximum likelihood estimation). Another recent approach in QST is observed in the work of Farooq et al., (2022) [19]. In that work, the authors developed a formulation for pure quantum state reconstruction via eigenvalue decomposition. In Farooq et al., (2022) [19], the authors demonstrated that the proposed approach is robust against depolarizing noise (high-strength white noise) - where the quantum states were reconstructed accurately similar to the noiseless case. A computational-focused strategy for QST is seen in the work Ahmad et al., (2022) [20]. In that research, the authors proposed a simultaneous perturbation stochastic approximation algorithm with high-speed convergence. The central idea in that work was to develop a computational technique for estimating quantum states under limited computational conditions using a Barzilai–Borwein two-point step size gradient method. In Choi (2022) [21], a low-complexity effective QST scheme was proposed. The proposed scheme requires the measurement of only three observables for systems of any scale. The key concept is coupling of the system and the ‘pointer’ of single qubit. In that work, the authors also developed a similarly scalable efficient maximum likelihood computational framework for approximating states from incomplete statistical data.
Complex nonlinear dynamical systems have been successfully modeled using chaotic maps such as coupled map lattices (CML). For instance, in Lu et al., (2022) [22] an urban rail transit system under cascading failure was modeled with respect to network vulnerability using CMLs. In that work, the authors quantified the relation between passenger flow and the perturbation index. In addition to investigations of the anti-risk resistance capability of stations, the dependency of cascading failures on location and type of stations were identified. In Stenzinger and Tragtenberg, (2022) [23], CMLs were for a biomedical-focused application. In that work, the authors utilized spatiotemporal chaos in the CML to simulate cardiac reentry. Among their key research findings were clinical manifestations of certain types of tachycardia in the electrocardiogram as well as a novel type of dynamical pattern (with wavefronts comprised of harmonized bursts and cardiac plateaus). In Wang and Liu, (2021) [24], a one-dimensional, two-parameter with a wide-range system mixed coupled map lattice model was proposed for image encryption. The authors of that work performed simulations to establish the effectiveness of the proposed encryption algorithm for grayscale and color images. In addition, they also carried out security tests to ensure that the proposed method could hold against conventional security attacks. CMLs have also been recently employed in the field of wireless communication. For instance, in Xie et al., (2022) [25], a CML model was developed to assess the vulnerability of unmanned aerial vehicle (UAV) network. The authors in that work showed that precision interference on critical UAV nodes could inflict significant damage to the entire UAV network. In addition, the authors in Xie et al., (2022) [25] also discovered that network vulnerability increases with the intensity of external interference.
This work reformulates the practical entanglement distillation problem as a bilevel optimization model. This work extends the practical entanglement distillation problem to account for approximating the quantum state (density matrix) using QST using fuzzy relations. The bilevel practical entanglement distillation optimization problem is then solved using the combined CML and the dual annealing approach. The central idea is to leverage on the complex behavior of the CML to enhance the optimization capabilities of the dual annealing algorithm. This paper is organized as follows: the second section describes the model formulation for the bilevel practical entanglement distillation problem. In this section, the novel QST formulation using fuzzy relations is presented. The third section discusses the CML and its hybridization with the dual annealing algorithm. The fourth section presents analysis on the results generated by the numerical experiments. This paper concludes with some final remarks and recommendations for future research works.

2. Entanglement Distillation with Fuzzy Relations for QST

In this work, a bipartite entanglement distillation model is considered – where the central idea is to convert a state, ρ A B (density matrix form) into a state which is close to a maximally entangled state utilizing only local operations and classical communication. This communication takes place between two nodes of a communication network A and B. This can be represented mathematically as follows:
F =   Φ d η A ^ B ^ Φ d   such   that   | Φ d = 1 d ∑ i = 0 d − 1 | i A ^ | i B ^
where F ∈ ( 0,1 ) is the fidelity – i.e., closeness of the converted state to the maximally entangled state. A and B are the input registers while A ^ and B ^ are the output registers. d is the dimension of the quantum state and η A ^ B ^ is the converted state ( ρ A B → η A ^ B ^ ). | Φ d is the maximally entangled state across output registers A ^ and B ^ . Therefore, as an example if the dimension, d = A = B = 2 , then ρ A B would be as follows:
ρ A B = 1 − p | 01 01 | + p | Φ 2 Φ 2 |
In contrast to theoretical entanglement distillation, practical entanglement distillation frameworks allow for the possibility of failure [26]. Therefore, the fidelity parameter would only be relevant to the analysis if the entanglement distillation is a success. This then creates a multilevel scenario with the probability of distillation success, P δ ∈ ( 0,1 ) is cascaded by the fidelity parameter, F. Since the local quantum memory utilized to store the quantum entanglement is imperfect/non-ideal, this entanglement cannot be preserved for an arbitrary amount of time. Hence, the probability of distillation success, P ( δ ) would control the rate at which high-fidelity entanglement between the different nodes in the network is aimed at.
Practical entanglement distillation uses schemes involving the application of local operation and measurement on A and B registers. This is then followed by a measurement outcome (single exchange) using classical communication to determine distillation success/failure. In this work, the practical quantum distillation scheme presented in Rozpędek et al., (2018) [27] is considered. The central idea is to search for the following optimal parameters: (i) output dimension, d, (ii) input state, ρ A B as well as (iii) quantum channels (i.e., quantum operations). To optimize quantum operations, Choi isomorphism could be utilized - where a one-to-one correspondence between quantum channels and quantum states is established (with certain properties) (Jiang et al., 2013). This way, the isomorphism would carry over all the information from the original channel to the Choi state. Following the problem formulation in Rozpędek et al., (2018) [27], the bilevel optimization formulation of entangle distillation is follows:
M a x i m i z e →   F = A B P δ   T r | Φ d Φ d | A ^ B ^ ⊗   ρ A ′ B ′ T C ^ 1 , A ^ A ′ ⊗ C ^ 1 , B ^ B ′
subject to,
M a x i m i z e   →   P δ = A B   T r ρ A ′ B ′ T C ^ 1 , A ′ ⊗ C ^ 1 , B ′  
such that,
ρ A ′ B ′ T   ≥ 0
C ^ 1 , A ^ A ′ ≥ 0   ,   C ^ 1 , B ^ B ′ ≥ 0
C ^ 1 , A ′ ≤ I A ′ A   ,   C ^ 1 , B ′ ≤ I B ′ B ,   A = B = d ≥ 0   : d ∈ N
where A ′ and B ′ are Choi state equivalent for output registers A and B. Similarly, C ^ 1 , A ^ A ′   ,   C ^ 1 , B ^ B ′ ,   C ^ 1 , A ′ and C ^ 1 , B ′ are matrices depicting Choi states which correspond to quantum channels. The symbol ⊗   represents the Kronecker product and the dimensions of the identity matrices, I A ′ and I B ′ depend on the dimensions of the registers A ′ and B ′ .
In designing practical quantum information processing systems, quantum state tomography (QST) is employed to approximate a certain state of the density matrix. This is done by using measurement results of repeated state preparation. Unfortunately, QST becomes very challenging as the system scales up in terms of size - i.e., dimensions, d . Current methods for QST formulation include linear inversion, maximum likelihood, Bayesian, hybrid Bayesian-Monte Carlo and neural network approaches [28,29]. In this work, a fuzzy approach is developed and employed for QST. In Recasens, (2022) [30] it was demonstrated that symmetric positive semi-definite matrices could be constructed from the application of fuzzy relations on certain sets (fuzzy subsets). Consider sets: X and Y with both having cardinality of n and μ are fuzzy subsets of X and Y with the additive generator of the t -norm being: t N ( x ) = ( 1 − x ) N . Then the following matrix is constructed:
ρ ( x , y ) = 1 − ( 1 − μ ( x ) ) N − ( 1 − μ ( y ) ) N 1 / N   where   N = n + 1 l o g 2 e
In the case of estimating the density matrix, ρ A ′ B ′ in equations (2) and (3), ρ A ′ B ′ ∼ ρ is considered as symmetric and positive semi-definite with the additional normalization condition: T r   ( ρ A ′ B ′ ) = 1 . In this view, the μ ( x ) and μ ( y ) are fuzzy sets that act as pseudo-likelihood functions obtained from measurement data (empirical). Since x and y are quantum states | Φ d Φ d | for d ≥ 0 , it can be stated that x = y . The measurements of these quantum states in the form of fuzzy memberships may not be equal: μ 1 ( x = y ) ≠ μ 2 ( x = y ) and the normalization condition: ∑ A ′ B ′ ρ A ′ B ′ = 1 holds. Using the simplified notations for the pseudolikelihood measurements μ 1 ( x ) ~ μ 1 and μ 2 ( x ) ~ μ 2 , Equation (4) then becomes: ρ ( x , y ) = 1 − ( 1 − μ 1 ) N − ( 1 − μ 2 ) N 1 / N . An interesting result is obtained from the symmetric semi-definite density matrix construction when the differential form of equation (4) is obtained:
∂ ρ ∂ μ 2 = − 1 − μ 2 N − 1 1 − μ 1 N − 1 ∂ ρ ∂ μ 1   w h e r e   N = n + 1 l o g 2 e   ( 5 )
Equation (5) is then obtained by taking the partial differential for the density matrix in equation (4) and assuming μ and ρ to be continuous parameters. Similarly, the analysis is extended to seek for an analogue to approximate the factor in equation (5) using the binomial formula and Maclaurin series:
− 1 − μ 2 N − 1 1 − μ 1 N − 1 = − ∑ k ≥ 0 N − 1 k − μ 2 k ∑ k ≥ 0 N − 1 k − μ 1 k =   =   − ∑ k ≥ 0 − μ 2 k / k ! ∑ k ≥ 0 − μ 1 k / k !   ≈ − e − μ 2 e − μ 1  
Therefore, a generalized partial differential equation (PDE) independent of scale (i.e., cardinality or dimension) is obtained:
∂ ρ ∂ μ 2 + e   ( μ 1 − μ 2 ) ∂ ρ ∂ μ 1 = 0  
The measurements μ 1 and μ 2 are treated as continuous parameters and the solution to the scale-free exact solution to the PDE in equation (7) is obtained:
ρ μ 1 , μ 2 = c 0 e − μ 1 − c 0 e − μ 2
where c 0 is a constant. Considering μ 1 and μ 2   to be subsequent measurements, pseudo-likelihood measurements is conjectured to have the Markov property such that:
μ 2 =   P μ 1
where P is an n × n transition or stochastic matrix. With this conjecture, the density matrix, ρ ( x ) is approximated using data from the current pseudo-likelihood measurement, μ 2 and the previous pseudo-likelihood function, μ 1 . The current pseudo-likelihood measurement, μ 2 could be generated using the stochastic matrix, P as in equation (9). QST is highly efficient when estimating density matrices, ρ for small quantum systems. However, as the system scales-up and ρ increases in dimensions, effective QST becomes a challenging feat [28]. The fuzzy formulation for pseudo-likelihood relations presented in equations (4)-(9) aim to enable effective estimation of density matrices, ρ at larger scales. Treating density matrices, ρ as symmetric and positive semi-definite are critical requirements for this analysis.

3. Chaotic Dual Annealing Optimization

In recent times, optimization problems have grown in terms of complexity – i.e., nonlinearity, nonconvexity, multilevel and uncertainty in parameters. This directly results in stagnation in the solution method upon implementation – where the algorithmic technique gets trapped in a local optimum. This causes the task of optimizing such problems to become exceedingly difficult. An effective approach to tackle such issues is the utilization of stochastic optimization approaches. One such approach is simulated annealing (SA) [31]. SA is based on the idea of the physical concept of annealing where a solid is repeatedly heated and cooled until reaches a configurational state of minimal energy. In this work we use the dual annealing optimization approach - which is a coupled technique inspired by the classical simulated annealing as well as the fast simulated annealing approaches [32]. The dual annealing optimization approach used in this work implements a distorted Cauchy-Lorentz visiting distribution, g ∆ x t with q (the distribution’s shape controlling parameter):
g ∆ x t = k T t − A 1 1 + q − 1 A 2 A 3 ,   A 1 = D 3 − q , A 2 = ∆ x ( t ) 2 T ( t ) 2 3 − q ,   A 3 = 1 q − 1 + D − 1 2  
where k and D are constants, t is artificial time, Δ x ( t ) is the trial jump distance and T ( t ) is the artificial temperature. The candidate solution’s acceptance probability is set by defining an acceptance parameter, q ′ < 1 is determined as follows:
P = m i n 1 , f ( q ′ )   where   f ( q )   =   1 − ( 1 − q ′ ) β Δ E 1 1 − q ′
The decrements of the artificial temperature, T ( t ) is given by the following relation:
T t = T 1 2 q − 1 − 1 1 + t q − 1 − 1
In this work, the dual annealing technique is coupled with coupled map lattices (CML) to enhance its optimization capability by simulating chaotic non-linear dynamics. CMLs have applications in diverse fields of research and implementation: system identification, process modeling, design optimization, internet of things (IoT) and cryptography [24,33,34,35]. For a specific discreet-time dynamical map, CMLs consist of discreet interacting elements. These elements make up a multidimensional dynamical system. Hence, the CML model could be designed to simulate a dynamical system - where its macroscopic field variables are defined on a lattice. Independent processes then decompose the dynamical system. In this sense, each process represents a procedure in the CML. Nonlinear transformations are generated using each state variable (i.e., Lattice points) and coupling terms, which are executed iteratively. Since CMLs commonly generate chaotic behavior, a logistic map is introduced in this work to induce such dynamics by setting the constant, θ > 3.57 . The logistic map is given as follows:
x i + 1 = θ x i 1 − x i  
where i ∈ [ 1 , m ] is the iteration count. The recursive map for CML is as follows:
x i + 1 = ε r x i 1 − x i j + ( 1 − ε ) r x i 1 − x i j − 1
where s is the index for the vertices on the lattice and the coupling parameter, ε ∈ [ 0,1 ] . The second term in CML recursive map is for the neighboring lattice points. In this work a unidirectional coupling approach is taken with the convex map. The proposed techniques that combine the CML and the dual annealing optimization approaches is given in the algorithm below:
Algorithm: Combined Coupled Map Lattice and Dual Annealing Optimization Approach
START
1. Initialize dual annealing parameters: maximum iterations ,   T 1 , q , q ′ .
2. Initialize CML parameters: Logistic map parameter, θ number of lattice vertices (N), coupling parameter ( ε ).
3. Solve P ( δ ) using dual annealing algorithm to find optimal quantum state, ρ A ′ B ′ T and Choi states: C ^ 1 , A ′ and C ^ 1 , B ′
4.Initiate and run CML simulation.
5. IF Iteration=1:
Randomly initialize lattice values, x ( i = 1 )
6. IF Iteration>1 and fitness function = 1:
Back-substitute the improved computed lattice values
7. IF Iteration>1 and fitness function = 0:
Back-substitute the lattice values from the previous iteration ( i − 1 )
(so it does not lead the CML algorithm towards non-optimal regions)
8. Compute new lattice values ( x ) using CML (with logistic map)
9. Check fitness of solution candidates (binary decision function)
10. IF Candidate solution is improved:
10.1. Fitness function is Satisfied:
Fitness = 1, HALT PROGRAM
10.2. Fitness function is NOT Satisfied:
Fitness = 1, Repeat Step 4
11. IF Candidate solution is NOT improved:
11.1. Fitness is Satisfied:
Fitness function = 0, HALT PROGRAM
11.2. Fitness is NOT Satisfied:
Fitness function = 0, Repeat Step 5
12. Determine statistical moments on CML simulation data – mean and variance.
13. Using statistical moments on a Gaussian distribution, simulate random values for quantum state ρ A ′ B ′ T .
14. Using standard PRNG, simulate random values for Choi states: C ^ 1 , A ^ A ′ and C ^ 1 , B ^ B ′ .
15. Solve for F in the upper-level problem.
16. Re-initialize Stackelberg game framework until fitness function cannot be further improved.
END
The bilevel optimization problem presented in Equation (3) is solved within a Stackelberg game-theoretic framework [34]. The fidelity objective function, F is the upper level/leader while the probability of success, P ( δ ) lower level/follower. The leader’s strategy is to optimize the objective of the sub-problem ( P ( δ ) . This influences the follower’s strategy – optimizing the objective (F). The numerical framework iteratively solves every level of the optimization problem as a Stackelberg game reaches the optimal solution. In this work, the entanglement distillation problem was solved using: Combined Coupled Map Lattice and Dual Annealing technique (CML-DA) technique and dual annealing method using pseudo-random number generators (PRNG-DA). The sub-problem problem is solved using DA-PRNG approach by searching for the optimal dimension, d that maximizes, P ( δ ) . The quantum state, ρ A ′ B ′ T and the Choi states: C ^ 1 , A ′ and C ^ 1 , B ′ are generated using the PRNG. Consequently, using the obtained dimension, d, the density state ρ A ′ B ′   T and P δ , the upper-level problem is solved by searching for the best Choi states: C ^ 1 , A ^ A ′ and C ^ 1 , B ^ B ′ using the PRNG. As for the proposed CML-DA technique, the Stackelberg framework is employed, and each level of the problem is solved iteratively. In the CML-DA approach, the quantum state of the qubit before conversion ( ρ A ′ B ′   T ) is estimated using the CML and PRNG techniques. The parameter settings employed for the numerical experiments in this work for both the CML and DA techniques are given in Table 1 and Table 2:

4. Computational Analysis

The Stackelberg game-theoretic framework was utilized for the entanglement distillation problem with a QST using fuzzy relations to estimate the quantum state. The pseudolikelihood measurements were simulated and the optimization problem was then solved using: the combined CML and dual annealing technique (CML-DA) and the conventional dual annealing approach with pseudo-random number generators (PRNG-DA). The numerical experiments were performed using the Python programming language on Google Collaboratory platform on a cloud with Python 3 Google Compute Engine (RAM 12.68 GB and Disk space:107.72 GB). Each technique was executed with a total of 40 times - where each time the technique was run 3 times and the best solution was taken for each execution. Therefore, both techniques were individually executed a total of 120 times. The computational results obtained using both techniques were measured using the weighted hypervolume indicator for the bilevel problem:
w H V I = w 1 x * − x + w 2 x o * − x o
where the optimal solution candidate is ( x * , x 0 * ) and the nadir point is ( x , x o ). The weights w 1 and w 2 enables the relative importance of the contribution of the upper-level problem (1) and lower-level problem (2). In these experiments the weights: w 1 = 0.7 and w 1 = 0.3 . The nadir point is for the upper-level problem (or fidelity objective) (1) and lower-level subproblem (or probability of success) (2) is ( x = 1   a n d   x o = 1 ). The larger the value of the wHVI metric, the better the optimization performance. The rated individual solutions obtained using the PRNG-DA approach is given in Table 3 and the optimal quantum state of the best solution (pure state density matrix) is shown in Figure 1:
The individual solutions obtained using the CML-DA technique is ranked and shown in Table 4. The generated optimal quantum state (pure density matrix) for the individual best solution is given in Figure 2:
The best individual solutions reached by the PRNG-DA and CML-DA techniques have a quantum state, ρ A ′ B ′ T with the dimensions: d = 4 and d = 8 respectively. The difference in terms of the level of optimality reached by the individual best solutions measured using the wHVI was about 7.722 % . The scatter plot in Figure 3 shows the spread of the individual solutions generated using the PRNG-DA and CML-DA:
In Figure 3, it can be seen the shaded area highlights the optimal region in the objective space – where the optimality of the fidelity objective is given more importance as compared to the cascaded objective function: probability of success. Thus, more solutions generated by the CML-DA technique fall into the optimal region as compared to the PRNG-DA approach. This is reflected in the value of the overall optimality of the solutions measured using the weighted HVI. The overall optimality of the CML-DA and PRNG-DA approaches are 27.109 and 25.718 respectively. The CML-DA outperforms the PRNG-DA in terms of overall optimization performance by approximately 5.267%. The key factor influencing the techniques in this study is the CML component. The chaotic dynamics exhibited by the CML enables wide-range dynamical behavior as compared the standard PRNG. This boosts the performance of the CML-DA approach enabling it to perform a more thorough search of the objective landscape as compared to standard techniques (PRNG-based approaches). The CML is thus able to adapt better to the multilevel structure of the optimization problem and avoiding scenarios that may cause stagnation. In other words, the CML manages to escape local optima traps in the objective space. The programs developed using the mentioned computational techniques did not encounter convergent issues and performed smoothly during execution. The execution time taken for overall solution generation using the CML-DA technique was 156.804 seconds while the PRNG-DA approach took about 27 seconds. Thus, the trade-off experienced by the CML-DA in terms of optimization efficiency is the execution time as compared to the PRNG-DA approach. The additional algorithmic complexity embedded in the CML chaotic simulator comes with an additional computational cost. Nevertheless, both techniques showed stable performance when solving the entanglement distillation problem with the fuzzy relations for QST.

5. Conclusions & Recommendations

In this work, the entanglement distillation optimization was formulated in a bilevel structure. Using fuzzy relations, the quantum state of the system (density matrix, ρ A ′ B ′ T ) represented as a symmetric positive semidefinite matrix was estimated using pseudo-likelihood measurements. A scale-independent partial differential equation (PDE) was derived (see equation (7)). The solution to the PDE was employed for the QST of the quantum state. This multilevel optimization formulation was then solved using the CML-DA and PRNG-DA approaches within a Stackelberg game theoretic-framework. Comparative analysis showed that the CML-DA technique proved more efficient as compared to the PRNG-DA approach. Due to the CML-DAs capacity for wide-range chaotic behavior, it was able to navigate the objective space efficiently as compared to conventional techniques such as the PRNG-DA.
Future works could be directed on testing the performance of other novel meta-heuristics or evolutionary algorithms for practical entanglement distillation. In addition, research efforts could also be focused on reformulating the practical entanglement distillation problem using pseudo-likelihood measurements; to include environmental factors influencing the density matrix estimation in QST. This could be done by considering the quantum Liouville equation or the Gorini–Kossakowski–Sudarshan–Lindblad equation [36,37].

Conflicts of Interest

On behalf of all authors, the corresponding author states that there is no conflict of interest.

References

  1. Kalb, N., Reiserer, A.A., Humphreys, P.C., Bakermans, J.J., Kamerling, S.J., Nickerson, N.H., Benjamin, S.C., Twitchen, D.J., Markham, M. and Hanson, R., 2017. Entanglement distillation between solid-state quantum network nodes. Science, 356(6341), pp.928-932. [CrossRef]
  2. Li, M., Fei, S. and Li-Jost, X., 2011. Bell inequality, separability and entanglement distillation. Chinese Science Bulletin, 56(10), pp.945-954. [CrossRef]
  3. Ruan, L., Dai, W. and Win, M.Z., 2018. Adaptive recurrence quantum entanglement distillation for two-Kraus-operator channels. Physical Review A, 97(5), p.052332. [CrossRef]
  4. Ecker, S., Sohr, P., Bulla, L., Huber, M., Bohmann, M. and Ursin, R., 2021. Experimental single-copy entanglement distillation. Physical Review Letters, 127(4), p.040506. [CrossRef]
  5. Liu, C., Tu, T., Li, P.Y., Liu, X., Zhu, X.Y., Zhou, Z.Q., Li, C.F. and Guo, G.C., 2022. Towards entanglement distillation between atomic ensembles using high-fidelity spin operations. Communications Physics, 5(1), pp.1-9. [CrossRef]
  6. Lin, Y.Y., Sun, J.R. and Sun, Y., 2021. Bit thread, entanglement distillation, and entanglement of purification. Physical Review D, 103(12), p.126002.
  7. Qian, L., Chen, L., Chu, D. and Shen, Y., 2021. A matrix inequality for entanglement distillation problem. Linear Algebra and its Applications, 616, pp.139-177. [CrossRef]
  8. Kondra, T.V., Datta, C. and Streltsov, A., 2021. Catalytic transformations of pure entangled states. Physical Review Letters, 127(15), p.150503. [CrossRef]
  9. He, M., Malaney, R. and Burnett, B.A., 2021. Noiseless linear amplifiers for multimode states. Physical Review A, 103(1), p.012414. [CrossRef]
  10. Gour, G. and Scandolo, C.M., 2021. Entanglement of a bipartite channel. Physical Review A, 103(6), p.062422. [CrossRef]
  11. Yan, P.S., Zhou, L., Zhong, W. and Sheng, Y.B., 2022. Measurement-based logical qubit entanglement purification. Physical Review A, 105(6), p.062418. [CrossRef]
  12. Shchukin, E. and van Loock, P., 2022. Optimal entanglement swapping in quantum repeaters. Physical Review Letters, 128(15), p.150502.
  13. Riera-Sàbat, F., Sekatski, P., Pirker, A. and Dür, W., 2021. Entanglement-assisted entanglement purification. Physical Review Letters, 127(4), p.040502. [CrossRef]
  14. Gyongyosi, L. and Imre, S., 2020. Optimizing high-efficiency quantum memory with quantum machine learning for near-term quantum devices. Scientific reports, 10(1), pp.1-24. [CrossRef]
  15. Gyongyosi, L., 2020. Quantum state optimization and computational pathway evaluation for gate-model quantum computers. Scientific reports, 10(1), pp.1-12. [CrossRef]
  16. Brand, S., Coopmans, T. and Elkouss, D., 2020. Efficient computation of the waiting time and fidelity in quantum repeater chains. IEEE Journal on Selected Areas in Communications, 38(3), pp.619-639. [CrossRef]
  17. Sahu, A., Sreeram, P.G. and Madhok, V., 2022. Effect of chaos on information gain in quantum tomography. Physical Review E, 106(2), p.024209. [CrossRef]
  18. Schmale, T., Reh, M. and Gärttner, M., 2022. Efficient quantum state tomography with convolutional neural networks. npj Quantum Information, 8(1), pp.1-8. [CrossRef]
  19. Farooq, A., Khalid, U., ur Rehman, J. and Shin, H., 2022. Robust Quantum State Tomography Method for Quantum Sensing. Sensors, 22(7), p.2669. [CrossRef]
  20. Ahmad, S.T., Farooq, A. and Shin, H., 2022. Self-guided quantum state tomography for limited resources. Scientific Reports, 12(1), pp.1-8. [CrossRef]
  21. Choi, M.S., 2022. Single-qubit reaped quantum state tomography. Scientific reports, 12(1), pp.1-9. [CrossRef]
  22. Lu, Q.C., Zhang, L., Xu, P.C., Cui, X. and Li, J., 2022. Modeling network vulnerability of urban rail transit under cascading failures: A Coupled Map Lattices approach. Reliability Engineering & System Safety, 221, p.108320. [CrossRef]
  23. Stenzinger, R.V. and Tragtenberg, M.H.R., 2022. Cardiac reentry modeled by spatiotemporal chaos in a coupled map lattice. The European Physical Journal Special Topics, pp.1-12. [CrossRef]
  24. Wang, X. and Liu, P., 2021. A new full chaos coupled mapping lattice and its application in privacy image encryption. IEEE Transactions on Circuits and Systems I: Regular Papers, 69(3), pp.1291-1301. [CrossRef]
  25. Xie, X., Li, J., Huang, Z., Yang, Q. and Kwak, K.S., 2022. Coupled-Map-Lattices-Based Vulnerability Assessment of UAV Network in Interference Scenarios. Wireless Communications and Mobile Computing, 2022. [CrossRef]
  26. Xu, F., Ma, X., Zhang, Q., Lo, H.K. and Pan, J.W., 2020. Secure quantum key distribution with realistic devices. Reviews of Modern Physics, 92(2), p.025002. [CrossRef]
  27. Rozpędek, F., Schiet, T., Elkouss, D., Doherty, A.C. and Wehner, S., 2018. Optimizing practical entanglement distillation. Physical Review A, 97(6), p.062333. [CrossRef]
  28. Lukens, J.M., Law, K.J., Jasra, A. and Lougovski, P., 2020. A practical and efficient approach for Bayesian quantum state estimation. New Journal of Physics, 22(6), p.063038. [CrossRef]
  29. Quek, Y., Fort, S. and Ng, H.K., 2021. Adaptive quantum state tomography with neural networks. npj Quantum Information, 7(1), pp.1-7. [CrossRef]
  30. Recasens, J., 2022. On the relationship between positive semi-definite matrices and t-norms. Fuzzy Sets and Systems, 446, pp.26-37. [CrossRef]
  31. Guilmeau, T., Chouzenoux, E. and Elvira, V., 2021, July. Simulated annealing: A review and a new scheme. In 2021 IEEE Statistical Signal Processing Workshop (SSP) (pp. 101-105). IEEE.
  32. Xiang, Y., Gubian, S. and Martin, F., 2017. Generalized simulated annealing. Computational Optimization in Engineering-Paradigms and Applications, pp.25-46.
  33. Zia, U., McCartney, M., Scotney, B., Martinez, J. and Sajjad, A., 2022. A novel pseudo-random number generator for IoT based on a coupled map lattice system using the generalised symmetric map. SN Applied Sciences, 4(2), pp.1-17. [CrossRef]
  34. Ganesan, T., Vasant, P. and Litvinchev, I., 2022. Chaotic simulator for bilevel optimization of virtual machine placements in cloud computing. Journal of the Operations Research Society of China, 10(4), pp.703-723. [CrossRef]
  35. Wang, X., Wang, X., Teng, L. and Jiang, D., 2022. A novel meaningful image encryption algorithm based on newly-designed coupled map lattice and adaptive embedding. Optik, 270, p.170073.
  36. Kim, H.W. and Rhee, Y.M., 2020. Two-oscillator mapping modification of the Poisson bracket mapping equation formulation of the quantum–classical Liouville equation. The Journal of Chemical Physics, 153(21), p.214103. [CrossRef]
  37. Manzano, D., 2020. A short introduction to the Lindblad master equation. Aip Advances, 10(2), p.025106. [CrossRef]
Preprints 75430 g001
Figure 1. The quantum state, ρ A ′ B ′ T for the best individual solution generated using the PRNG-DA approach.
Figure 1. The quantum state, ρ A ′ B ′ T for the best individual solution generated using the PRNG-DA approach.
Preprints 75430 g001
Preprints 75430 g002
Figure 2. The quantum state, ρ A ′ B ′ T for the best individual solution generated using the CML-DA approach.
Figure 2. The quantum state, ρ A ′ B ′ T for the best individual solution generated using the CML-DA approach.
Preprints 75430 g002
Preprints 75430 g003
Figure 3. Spread objective values for individual solutions generated using the PRNG-DA and CML-DA.
Figure 3. Spread objective values for individual solutions generated using the PRNG-DA and CML-DA.
Preprints 75430 g003
Table 1. Parameter settings for the DA technique.
Table 1. Parameter settings for the DA technique.
Parameter Value
Initial Temperature, T ( 1 ) 15
Restart temperature ratio 0.0002
Shape controlling parameter, q 2.62
Acceptance parameter, q ′ -5
Maximum iterations 100
Table 2. Parameter settings for the CML.
Table 2. Parameter settings for the CML.
Parameter Value
Logistic map parameter, θ 3.7
Maximum iteration 100
Lattice vertex count (N) 100
Coupling parameter ( ε ) 0.5
Parametric interval ( Δ r ) 0.05
Table 3. Rated individual solutions obtained using the PRNG-DA approach.
Table 3. Rated individual solutions obtained using the PRNG-DA approach.
Parameters Best Median Worst
d 4 4 4
P ( δ ) 0.9897 0.9878 0.9839
F 0.8873 0.5201 0.1218
Iterations 341 397 557
  ρ A ′ B ′ T [0.10023951, 0.18421193, 0.02052628, 0.69502228] [0.1290744, 0.46015852, 0.36462488, 0.0461422] [0.07538315, 0.81589332, 0.07227743, 0.0364461]
wHVI 0.918 0.6604 0.3805
Table 4. Rated individual solutions obtained using the CML-DA approach.
Table 4. Rated individual solutions obtained using the CML-DA approach.
Parameters Best Median Worst
d 9 9 4
P ( δ ) 0.9745 0.9813 0.9699
F 0.9991 0.5105 0.0435
Iterations 315 495 443
  ρ A ′ B ′ T [5.91809983e-03, 2.26849719e-04, 2.75340776e-02, 1.54353405e-01, 5.56746335e-02, 4.33851377e-01, 4.06979776e-03, 3.17812492e-01, 5.59267748e-04] [0.00190849, 0.05472652, 0.03569338, 0.23633382, 0.34003655, 0.20481441, 0.00174142, 0.08249384, 0.04225155]] [0.33783083, 0.09176041, 0.19728488, 0.37312388]
wHVI 0.9917 0.6517 0.3214
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated