Quantum image representations based on density matrices in open quantum systems

So far, research on quantum image representation has gone through more than 20 years. During this time, the quantum image representation models used have almost all been based on state vectors. However, in practical problems, the environment and the principal quantum system cannot be separated, and isolated quantum systems do not exist in principle. This case is often referred to as an open quantum system. In open quantum systems, many problems involve density matrices, such as the calculation of Von Neumann entropy, the quantization of coherence, and the operator-sum representations of quantum operations. Therefore, the existing quantum image representation models are only suitable for closed quantum systems. To this end, the paper proposes three models that can not only represent quantum images in an open quantum system but also decompose the evolution process of quantum images utilizing operator-sum decomposition. These three models are the representation model of quantum gray-scale images, the tensor product representation model of quantum color images, and the representation model of quantum color images based on mixed states in the Bloch sphere, respectively. All these image representation models have strong correlations among them and are very diﬀerent from their classical analogues. Between them, the biggest diﬀerence is that the paper employs density matrices, inspired by incoherent-coherent states, to represent quantum images rather than classical state vectors. By means of one of the representation models proposed in the paper, we ﬁnally demonstrate the evolution process of the quantum image going through the amplitude damping channel.


Introduction
In recent decades, the combination of quantum computing [1] and digital image processing [2] has been extensively studied.Quantum image processing (QIP) [3], an area focused on extending conventional image processing tasks and operations to the quantum computing framework, is the new sub-area that has emerged in that regard.The most fundamental problem in this field is the problem of quantum image representation (PQIR) [4][5][6][7], the development process of which can be divided into two stages.The first stage starts with the Qubit Lattice representation of quantum images proposed by Benegas-Andraka and Boshi [8] in 2003.This was closely followed by Entangled Image [9] and Real Ket [10] representation models.In the second stage, the flexible representation of quantum images (FRQI) [11] was proposed in 2011, followed by the novel enhanced quantum representation of digital images (NEQR) [12] was proposed in 2013, and they are generally considered pioneers of this phase.
However, these representations and applications of quantum images mentioned above are all in closed quantum systems, which require the principal system not to interact with the external environment.In the classical system-environment division, the environmental perturbations to the system can be arbitrarily reduced until their effects are negligible through rational design and division.But the quantum system-environment division is different from that of the classical system-environment, everything that happens in the principal system is always in extensive entanglement with the external environment and even if the energy and matter transfer on the boundary is small enough to be negligible, this entanglement still has a huge impact.Thus, the environment and the principal system can never be separated, and isolated quantum systems cease to exist in principle.
In real physical systems, no matter how weak the coupling that can prevent the principal system from being isolated, the evolution of an open quantum system [13] is eventually plagued by non-unitary features like decoherence and dissipation.The so-called quantum decoherence [14] here is a purely quantum-mechanical effect whereby the system loses its ability to exhibit coherent behavior by getting entangled with the ambient degrees of freedom.Decoherence stands as a serious obstacle common to all applications relying on the capability of maintaining and exploiting quantum coherence.According to the DiVincenzo criteria [15], we require relatively long decoherence times to ensure coherence is preserved during gate operations.Therefore, when the quantum systems we are interested in are no longer isolated or closed but become open quantum systems, it is necessary to quantify the coherence of quantum images in the evolution process.This allows us to better understand and control the effects of these non-unitary features (which can be considered as noise in quantum information processing) on the system's evolution, and take appropriate measures to reduce their negative impact.Hence, it is of significant importance to find better ways to represent quantum images for the future development of quantum image processing.
In the framework of quantum Internet [16], the study of methods and techniques for open-system quantum communication is crucial for building efficient quantum communication networks.Distributed gate models [17] can be used to construct complex and powerful quantum computing networks, where multiple gate model quantum computers can operate and exchange information in a distributed manner.In quantum image transmission, especially under the open system, the effects of noise and distortion need to be considered.At this time the capacity of the quantum noise channel [18] is an important consideration, i.e., the maximum amount of quantum information that can be transmitted in a channel affected by quantum noise.Therefore, when designing and implementing a distributed gate model quantum computing setup, one needs to consider how to ensure the integrity and accuracy of quantum image transmission.The interplay of these aspects promotes the development of the fields of quantum computing and quantum communication, helps to build more reliable quantum communication systems, and promotes the application of quantum technology in various fields.
Inspired by the above problem, the paper considers the PQIR in open quantum systems and gives three specific representation models that can be used to quantify coherence [19][20][21] and analyze the dynamic conditions of coherence completely unaffected by noise in open quantum systems.Since the principal system is entangled with the environment in an open quantum system, it is usually described by a mixed state rather than a pure state of the wave function.On this occasion, the following advantages of the state density matrix are demonstrated: • The density matrix can accurately and completely characterize all measurement properties of open quantum systems.
• Some rigorous frameworks have been proposed to quantify coherence and determine intuitive and computable measures of coherence.
• It can be used to calculate Von Neumann entropy, entropy exchange, and operatorsum representations of quantum operations.Therefore, the paper utilizes the density matrix to store the gray value of each pixel and completes the preparation process by performing rotation operations on the surface of the Bloch sphere.
The rest of the paper is organized as follows: -Sect. 2 presents the preliminary knowledge, -Sect.3 gives three novel models for representing gray-scale images and RGB color images in open quantum systems and gives the process of quantum image preparation, -Sect.4 introduces the environments and quantum noise channels and gives an example of a quantum image coupling with the environment.
Finally, we draw conclusions and outline possible future research tracks.

Related work
To date, various quantum image-representation model have been proposed to store and process image information, in quantum computing, such as qubit lattice, entangled image, real ket representation, flexible representation of quantum images (FRQI), novel enhanced quantum representation of digital images (NEQR).
In FRQI, a normalized superposition is used to store used to store position and grayscale information for all the pixels in an image.The main drawback of this model is that FRQI uses only one qubit to store the gray-scale information of each pixel in the image, which cannot process complex color operations.To improve FRQI, the following proposed NEQR model uses two entangled qubit sequences to store the gray-scale and position information and stores the whole image in the superposition of the two qubit sequences.However, since NEQR can only represent square images of size 2 n × 2 n , it cannot be applied to more rectangular images.To solve this problem, NEQR was improved into INEQR [22] to represent quantum images of size 2 n 1 × 2 n 2 Furthermore, since in image representation and processing, we usually treat quantum images as two-dimensional data, which are typically represented by floating-point numbers.Therefore, in 2020, Zhang [23] et al. improved the representation of quantum images by increasing the pixel size from positive integers to floating-point representations.This enhancement allows for a finer and more flexible representation of quantum images of arbitrary size, enabling the representation of arbitrary sized 2-D data and accommodating element values of arbitrary precision.Additionally, this method not only handles twodimensional data of arbitrary size but also facilitates the processing of three-dimensional data through dimensionality reduction.Furthermore, fewer qubits are needed after the dimensionality reduction process compared to the direct preparation of three-dimensional data.This proposal integrates information about a H × W 2-D data into a quantum state having it formula in (2.1): here, |YX is the location information, |S YX F can store a floating-point number of 2-D data.Therefore, these latest models are also improvements on the FRQI as well as the NEQR.Subsequently, numerous research results in the field of quantum computing have emerged, including the QIC algorithm [24][25][26], floating point quantum algorithm [27,28], image scrambling [29], quantum encryption [30], quantum neural networks [31,32] and, so on.
However, all the above representations are based on closed quantum systems, inspired by [33], when a quantum image transmits quantum information through a noisy channel (i.e., a channel that is not isolated from the environment), assuming that the quantum image is the primary system and undergoes a dynamic evolution process, which can be expressed as the transmission of a quantum image through a noisy quantum channel, the evolution process from initial state to final state is described using density matrices.At the same time, the fact of mathematically describing the evolution of the quantum image in a noisy channel need to use the operators-sum representations of the noisy channel, which inevitably concern the density matrix.Also the use of density matrix description paves the way for the quantization of coherence as well as frozen coherence in future researches.

Preliminaries
In this section, we introduce the basic concepts and some knowledge about the Bloch sphere and present the definition of the incoherent-coherent state [34][35][36].

Bloch sphere
Comparing the 0 and 1 of classical bits, the two possible states of a single qubit are |0 and |1 .In three-dimensional Bloch Sphere, a qubit can be written in the form as: where | cos θ 2 | 2 and |e iφ sin θ 2 | 2 respectively represent the probability of becoming |0 and |1 , also they satisfy the normalization condition: , and e iφ is the relative phase factor.In the Bloch sphere, a point P can be determined by two angles θ and φ, i.e., the state vector corresponds to a point on the surface of the Bloch sphere, as it is shown in Fig. 1.
Before we describe the specific use of density matrices in QIP, it is important to introduce the related knowledge of density matrices.To do so, let us first consider the density matrix Next, let Bloch vectors υ = (υ 1 , υ 2 , υ 3 ), where υ 1 = sin θ cos φ, υ 2 = sin θ sin φ, υ 3 = cos θ , then (3.3) can be rewritten as follows: (3.3) Now, since the amplitude of the Bloch vector is unity for a pure state, i.e., υ 2 1 + υ 2 2 + υ 2 3 = 1, this is perfectly consistent with the definition of a pure state and we may conclude that the points on the surface of the Bloch sphere depict pure states.Reference [1] also shows that mixed states can be characterized by the condition υ 2 1 + υ 2 2 + υ 2 3 < 1, thus the single-qubit mixed states can be represented by the points inside the Bloch sphere.

Incoherence-coherence state
Quantum coherence represents a basic feature of quantum systems that is not present in the classical world.Coherence is regarded as a precious resource that cannot be generated or increased under a restricted class of operations known as incoherent operations.Unlike other resources, quantum coherence is a basis dependent quantity, a quantum state ρ is said to be incoherent in a given reference basis {|i }, if the state is diagonal in this basis, i.e., if ρ = i p i |i i| with some probabilities p i .In Chitambar's research assisted coherent distillation task [34], this task arises naturally in bipartite systems where both parties work together to generate the maximal possible coherence on one of the subsystems.For a bipartite system, the reference basis is assumed to be a tensor product of local bases.So in the bipartite system, if a quantum state is called an incoherent-coherent state, then such a state has the following form: Despite the fundamental importance of quantum coherence, only very recently have relevant first steps been introduced to a rigorous framework for the quantification of coherence and identified some intuitive and computable measures of coherence [19].In fact, for the incoherent-coherent state, since one postulate for coherent measurements requires that C(ρ), the coherence of ρ, does not increase on average under selective incoherent operations, then for which we choose the incoherent operation for the local projective measurements {|i i|} on the first subsystem.Another postulate for coherent measurements requires that C(ρ) is a convex function of density matrices, which implies Therefore, the coherence of incoherent-coherent and coherent-incoherent states is the average of the coherent parts, which proves that the coherence of the incoherent-coherent states can be quantized.

Image representation in open quantum systems
All real-world quantum systems always interact with their surroundings to a greater or lesser extent, so that quantum decoherence occurs, such the quantum systems are often called open quantum systems.Representation of quantum images in open quantum systems becomes a crucial step in exploring quantum computers with noise for practical applications.But the quantum image representation models used so far, such as FRQI and NEQR given in (4.1) and (4.2) respectively, are all represented in the form of state vectors as follows: Therefore, the desired density matrices corresponding to (4.1) and (4.2) can be obtained, respectively, by It is intuitively obvious that both ρ FRQI and ρ NEQR are so complex that their storage and preparation would require a large number of resources, which question demands us to find new tools to represent quantum images in open quantum systems.

Representation of quantum gray-scale images
In this subsection, combining the encoding method given by FRQI with the incoherentcoherent state, a novel representation for gray-scale images based on density matrices is proposed.Furthermore, the time complexity of the preparation process is also discussed.
The representative expression of a quantum image in the size of 2 n × 2 n can be written as follows: where |Y X| encodes the information about position, and ρ 1 (θ YX ) is a density matrix related to θ YX , which encodes the gray-scale value corresponding to pixel (Y , X) and is specifically expressed as follows: It is obvious that formula (4.4) is equivalent to formula (3.3) with φ at 0, i.e., ρ 1 (θ YX ) = ρ(θ YX , 0).This interesting consequence inspires us to represent the intensity in the grayscale image by employing the quantum states in the Bloch sphere with φ at 0. As shown in Fig. 2, ρ 1 (θ YX ) = ρ(θ YX , 0) is the quantum state in the dashed part of the Bloch sphere.Therefore, only the dashed part of the Bloch sphere is used to encode the gray scale.In quantum computing, computers are usually initialized in well-prepared states.Therefore, the preparation process of transforming the quantum computer from the initialization state to the desired quantum state is necessary.The differences between the evolution of the state vector and that of the density matrix during the preparation process can be specifically described as There are two ways to prepare the desired quantum states in (4.5).One is based on the relevance of the evolution process between the state vector and the density matrix in (4.7).As a result, in the preparation process, the newly proposed gray-scale image representation, which is encoded based on the density matrix, can replace the state vector in FRQI.Also the same evolutionary operator R y is adopted in both preparation processes, so it can be concluded that the time complexity of both methods is the same.
The other method is that we take advantage of the fact that the dashed part on the Bloch sphere shown in Fig. 2 represents the intensity of the quantum gray-scale image.So, we can apply the rotation operators to realize the preparation of pixel gray-scale information.Recall that Any quantum state on the Bloch sphere can be realized by the rotation operators R z (φ) and R x (θ ) acting on the |0 , i.e., When discussing the preparation process of gray-scale images in this method, it is worth noting that in equation (4.8), when φ = 0, the effect of R z (φ + π/2)R x (θ ) acting on |0 0| is Figure 4 The quantum gray-scale image circuit of Fig. 3 and its equivalent form equivalent to that of R y (θ ) acting on |0 0|, i.e., Therefore, the quantum gray-scale image circuit of Fig. 3 and its equivalent form can be shown in Fig. 4. In addition, the position information setting in the circuit will be described in detail in Sect.4.4, Fig. 9.
Recalling the facts that the first preparation method employs the rotation operator R y (θ ), and the second one utilizes the composite of the rotation operator R z (φ)R x (θ ), it follows from the above equation that the two preparation methods are equivalent, Therefore, the time complexity of the preparation process of (4.5) is similar to that of FRQI.
In contrast, when using the FRQI model or the NEQR model to represent a quantum image in an open quantum system, 4n qubits are required to generate the density matrix which stores the position information as far as the position information is concerned, consuming a large amount of resources.However, only 2n qubits are needed to store the position information when the density matrix represents the quantum image, saving 50% of qubits.Moreover, as for the computational efficiency, which refers to the complexity of the preparation process, i.e., the number of quantum universal gates, by the way the density matrix shown in Eq. (4.7) evolves, we can learn that the time complexity of the preparation process in Eq. (4.5) is similar to that of FRQI.As for the use of the latest models simply changing the storage of gray values on the basis of FRQI as well as NEQR, the preparation process is similar.Therefore, this approach effectively saves precious qubits and has a unique advantage in dealing with the dynamical evolution of open quantum systems.

Tensor product representation of quantum color images
The RGB color representation model is one of the most commonly used multi-channel color models, which produces new color information by combining three channels of color information: red, green and blue.Using the RGB color model, a color image can be decomposed into three gray-scale channel images.As mentioned above, the FRQI model has widespread appeal in quantum color images.Inspired by the representation technique in Sect.4.1, i.e., the gray value at each pixel can be encoded by a pure state ρ(θ , φ) with φ = 0 on the surface of the Bloch sphere.A question arose naturally: Whether the encoding technique for gray-scale images can be extended to RGB color images?
The way to deal with this problem is as follows: The gray value information at each pixel is still encoded with ρ(θ , φ), but the difference is that φ encodes the information of red, green and blue.Specifically speaking, ρ(θ , 0) encodes color red; ρ(θ , π/2) encodes color green; and ρ(θ , π) encodes color blue.In all these cases, the angle θ is still used to encode the brightness of the color.
The representation of a 2 n × 2 n RGB color image can be written as follows: As can be seen from the above formula, |Y X| encodes the information about the position and encodes the RGB color information corresponding to the position |Y X| in quantum images, and at this time ) and ρ(θ 3 YX , π) are used to encode the information of red, green and blue, respectively.According to (3.3), a simple calculation leads to the following conclusion: ) 12) The preparation of gray-scale information can be implemented by applying rotation operators that act on a single qubit on the Bloch sphere according to Sect.4.1, and the quantum circuits for (4.10), (4.11) and (4.12) are shown in Fig. 5.The method above is a slight modification of that shown in Sect.4.1.More specifically, we use the three curves on the Bloch sphere to represent the states of the R, G, and B channels, and the preparation process for each curve is the same.When φ = 0, it is observed that the representation is the

Mixed state representation of quantum color images
This subsection provides another representation of quantum color images in open systems.The color of each pixel in a quantum color image is determined by the three components R, G, and B, and each component has a value range of 0 to 255.First, the R, G, and B components of the quantum color images are normalized, and then υ 1 , υ 2 and υ 3 are applied to encode R, G, and B, respectively.The most intuitive purpose of normalizing R, G, and B components is to reduce the dependence of the color represented by the RGB value on the brightness of the image.Therefore, the information about (R, G, B) of the pixel position (Y , X) in an arbitrary quantum image can be written as Thus, the corresponding (R, G, B) information for the pixel position (Y , X) in a quantum color image can be encoded by The second representative expression of a quantum color image in the size of 2 n × 2 n can be written as follows: The density matrix ρ YX in (4.14) encodes the information of color corresponding to the pixel position |Y X|, where , I is the identical matrix, σ = (σ 1 , σ 2 , σ 3 ) and σ 1 , σ 2 , σ 3 are the Pauli matrices, respectively.
It is worth noting that ρ YX is the density matrix of a certain mixed state ensemble, and so the point corresponding to ρ YX will be in the interior or on the surface of the Bloch sphere.) is its complexity in preparing the mixed state ρ YX , which limits its practical applicability and confines it to theoretical research.In the paper [37], we learned that, in the general situation, it is possible to purify the mixed state formed by a quantum system Q in interaction with its environment by introducing a reference system R. Specifically, the system R is a dynamically isolated system and has a zero internal Hamiltonian, so that the system R and the quantum system Q can be formed into a joint system RQ to reach a purely entangled state |ψ RQ .Then, the pure state ensemble of quantum system Q can be obtained through the purification process.
And in the paper, the preparation method of the mixed state is explored by utilizing the relation equation between the mixed state and the pure state on the Bloch sphere.Specifically, to realize the preparation of ρ YX , a commonly used approach is employed, where the point corresponding to ρ YX within the Bloch sphere is mapped onto the surface of the Bloch sphere.This mapping allows us to represent ρ YX using a set of pure states.The precise understanding of the relationship between mixed states and pure states is further clarified in Lemma 4.1, which is presented below.

Lemma 4.1 A mixed state characterized by any point C in the Bloch sphere can be decomposed into two pure states characterized by two points A and B on the Bloch sphere, respectively, where C is on the line connecting A and B.
Proof From the geometric relationship in Fig. 6, we can have

Recalling (3.3) we know
In addition, This finishes the proof.
Problem 2 Although Lemma 4.1 gives the answer to Question 1, another question arises naturally, the two pure states, representing the mixed state, corresponding to A and B are obviously not unique, which means that the code based on such pure state decomposition is not uniquely decodable.
Next, the method to solve Problem 2 is given, which is divided into the following two steps.
Step 1: The RGB values of the quantum color image are translated into the corresponding points inside the Bloch sphere.For the (R, G, B) information corresponding to the given pixel position (Y , X), the vector − → υ YX = (υ 1 YX , υ 2 YX , υ 3 YX ) can be obtained using equation (4.13).Then, it is possible to calculate (r YX , θ YX , φ YX ) for the corresponding points inside the Bloch sphere.
Step 2: At fixed point A, the mixed state B(υ 1 , υ 2 , υ 3 ) inside the Bloch sphere can be uniquely decomposed into pure states.
In order to obtain such a result, we first give a fixed point A(0, -1, 0) on the Bloch sphere (see Fig. 7).Then, the half-line − → AB intersects the surface of the Bloch sphere at the point C(υ 1 , υ 2 , υ 3 ), where C is the point related to B that we are looking for.According to the following simultaneous equations to find the intersection, the coordinates of point C, of the half-line − → AB with the Bloch sphere we have Since C is on the surface of the Bloch sphere, through geometric relationships we have 2υ 1 (υ 2 + 1) . (4.17) In order to establish a unique decodable relation between B and C, the parameter λ is also introduced, where λ = |AB|/|AC|.Here, |AB| and |AC| represent the length of line segments AB and AC respectively.So far, the mixed state B is jointly determined by pure states A, C and the parameter λ.The whole mapping procedure of a unique decodable code representing a mixed state by two pure states can be seen in Fig. 7.
To better understand the encoding method presented above, two examples are given next.
Example 2 Calculate the coordinates of point C on the surface of the Bloch sphere and the parameter λ, which jointly determine (R, G, B) information as (255, 0, 0).
It is easy to know that the coordinates B(υ 1 , υ 2 , υ 3 ) of the mixed state are ( √ 3/3, 0, 0).So we have Finally, based on the above discussion, an alternative representation in the size of 2 n × 2 n RGB color image is as shown in (4.18): The new model proposed here is called the mixed state representation of quantum color images, abbreviated as the MSR-QCI model.Supposing that the accuracy of the parameter λ YX is 2 -q , i.e., λ YX ∈ [0, 2 -q -1], then the parameter λ YX corresponding to the pixel (Y , X) is encoded as follows: In (4.18), |Y X| is used to represent position information of pixel, ρ YX refers to the representation of the mapping of the mixed states within the Bloch sphere to the pure states of the corresponding sphere, while ρ YX ⊗ |λ YX λ YX | encodes the color information.

Preparation process of MSR-QCI model
In order to obtain the desired quantum image on the quantum computer, the preparation process of the MSR-QCI model is necessary and is described in this subsection.From the representation of the MSR-QCI model, n + q + 1 qubits are needed to construct the quantum image model for a 2 n × 2 n RGB color image.The first step is to prepare n + q + 1 qubits and to set all of them to |0 0|.The initial state can be expressed by Recalling the discussion in Sect.4.3, the desired pure state C(θ 1 , φ 1 ) and parameter |λ can be calculated in advance.The preparation process for the MSR-QCI model is now divided into three steps.
Step 1: In this step, single qubit gates I and H are used to prepare position information, which converts the initial state ρ 0 to the intermediate state ρ 1 , where The whole quantum operation in step 1 can be expressed by U 1 : (4.20) interprets the transformation from the initial state ρ 0 to the intermediate state ρ 1 .After this step, the position information for all the pixels is stored in the MSR-QCI model. (4.20) Step 2: To prepare the quantum color image, it remains to set the information of color for every pixel.In effect, this step can also apply the rotation of qubit on the surface of the Bloch sphere to realize the preparation of pixel color information.Given pixel position (Y , X), define the quantum controlled rotation operation U YX as follows: where Acting YX (θ 1 , φ 1 ) on the initial state |0 0| gives us: Then U YX acting on ρ 1 can obtain the following result: Therefore, the whole rotation operator U 2 can be defined as Through the quantum operator U 2 , the intermediate state ρ 1 is transformed into ρ 2 .
Step 3: The final step is to prepare the parameter |λ YX at each pixel position (Y , X), which is expressed as a fixed-point decimal with an accuracy of 2 -q and is stored in the quantum circuits using q qubits.In effect, we also need to divide the last step into 2 2n suboperations to store the parameter information of each pixel.Fixed pixel position (Y , X), the corresponding quantum sub-operation YX is defined as follows: where ϒ YX is a quantum operator as shown in (4.23), which is the value-setting operator for pixel position (Y , X): From(4.23), if Otherwise, it is a quantum identity gate that will do nothing on the quantum state.Therefore, the quantum transformation of ϒ YX to set parameter |λ YX for the pixel is as follows: Therefore, the intermediate state ρ 2 is transformed as in (4.24): It can be seen from the above formula that every operator YX only sets the parameter |λ YX of its corresponding pixel position.Therefore, the whole work U 3 of step 3 consisting of all the sub-operators is as follows: After the three steps described above, the entire preparation is done.(4.26) illustrates the representation of a quantum color image of size 2 × 2 in the MSR model, and Fig. 8 shows the detailed quantum circuits for this 2 × 2 image MSR model.
In quantum computation, computers are usually initialized in well-prepared states.When we discussed the time complexity of quantum image preparation, a complex operation needs to be broken down into simple gates.Thus, the controlled operation in which the control qubit is the position qubit and the target qubit is the color qubit in Fig. 8 has its equivalent form as shown in Fig. 9.
The total number of single-qubit and double-qubit operations used to prepare the MSR-QCI mode is Similarly, the time complexity of the preparation for the tensor product representation of quantum color images is also obtained.From the model (4.9), n + 3 qubits are needed to construct the quantum image model for a 2 n × 2 n image.To begin with, for the quantum operation of U 1 = I ⊗3 ⊗ H ⊗2n , its time complexity is obviously O(3 + 2n).

Environments and quantum noise channels
A natural way to describe the dynamics of open quantum systems is to think of them as interactions between the system of interest (here called the principal system) and the environment, which together form a closed quantum system.The dynamics of such a closed quantum system is described by a unitary transformation, and we can think of this unitary transformation as a box.For us, we are not concerned with the internal mechanism of the box, which can be implemented by quantum circuits or some Hamiltonian systems, or even anything else.In other words, suppose we have a system in state ρ, and this state is sent into a box which is coupled to an environment.Then, the final state ε(ρ) of this system may not be related by a unitary transformation to the initial state ρ.
As illustrated in Fig. 11, we now assume that the input state of the system-environment is a product state ρ ⊗ ρ env .After the U transformation of the box, the system no longer interacts with the environment, because we perform a partial trace over the environment to obtain the reduced state of the system alone.
This approach is easy to relate to the real world, but it has the disadvantage of being mathematically inconvenient, whereas operator-sum representations for understanding quantum operations, although exactly similar to it, provide a powerful mathematical representation of quantum operations.
If the representation of equation (5.1) is based on the operator on the Hilbert space of the principal system alone, let |e k be a set of orthonormal basis of the (finite dimensional) state space of the environment, and then let ρ env = |e 0 e 0 | be the initial state of the environment.Suppose that the environment starts from a pure state, even though the environment is started in a mixed state, we can use it as an intermediate step in the calculation by introducing an additional external system that purifies the environment and does not make any difference to the dynamics experienced by the principal system.Therefore equation (5.1) can be rewritten as: where E k = e k |U|e 0 is an operator on the state space of the principal system.The above equation is called the operator-sum representation of ε, and the operator E k is called the operator element of the quantum operation ε.
We can draw a physical explanation for the operator-sum representation, assuming that after the unitary U transformation, the measurement of the environment is performed on bases |e k .Applying the principle of implied measurement, it can be seen that such a measurement only affects the state of the environment without changing the state of the principal system.According to this idea, we can deal with noise in the quantum noise channels without affecting the state of the principal system.We regard unnecessary interactions in open quantum systems as noise in quantum information processing, and the processed quantum images as the principal system.In order to build a useful quantum information processing system, we need to understand and control these noise systems without affecting the principal system.
For instance, the amplitude damping channel, which is a quantum operation that describes the energy dissipation of a system to the environment.It is typically an atomic process and can be thought of as the spontaneous emission of a photon into the environment by the decay of an excited state of a two-level atom in the presence of an electromagnetic field.This channel can be described with the map: When there is no excitation present, the system and the environment remain unaltered, while when an excitation is present in the system, it can either remain itself with probability (1p) or decay to |0 prin |1 env with probability p, producing an excitation in the environment.So the operator elements of amplitude damping operator ε are where γ = sin 2 θ can be considered as the probability of losing a photon.If we transmit the quantum gray-scale image in the amplitude damping channel, the quantum image will be in the principal system.After coupling with the environment, the partial trace for the environment can be obtained as follows: Figure 12 Frozen coherence pre-processing process and specifically speaking, If we compare the above result with the original density matrix of the qubit given in (4.6), then we can easily observe that the effect of the coupling with the environment is the introduction of (1γ ) to the density matrix.Meanwhile, the diagonal terms in the density matrix will evolve faster than the off-diagonal terms, which implies that decoherence occurs after the quantum image model is coupled to the noise channel.This coupling will result in the loss of information from the principal system into the environment.Therefore, the model (5.In an open quantum system, a quantum image may face the phenomenon of "decoherence" during transmission due to interference from the external environment, resulting in the quantum state become incoherence state.However, by pre-processing the quantum image, we can make its gray scale information satisfy the condition of frozen coherence, which greatly protects the coherence during transmission / process and maintains the clarity of the image at the same time.This process can be depicted as shown in Fig. 12.
This performance enhancement in coherence protection will not only enable more secure and reliable quantum communication / process, leading to higher resolution and clearer images, but will also advance the development of some application fields such as quantum cryptography and so on.
The quantum image representation models proposed in the paper are also applicable to various other noise channels, such as bit-flip, bit-phase flip, phase-flip channels, and the phase damping channel.Their action on the quantum image model is described below in terms of a parameter q ∈ [0, 1] which encodes the strength of the noise.The bit-flip, bit-phase flip, and phase-flip channels can be represented in operator form by with k = 1, 2 and 3, respectively, and σ j being the j-th Pauli matrix.Finally, the phase damping channel can be characterized by

Conclusion
The focus of this paper is to address the problem of quantum image representation encountered in the field of open quantum systems.In the paper, the shortcomings of these existing quantum image models are first identified by analyzing the quantum systems that exist in reality.Since the quantum system inevitably interacts with other quantum systems around it and becomes an open quantum system, the previous image representation methods are only applicable to closed quantum systems.To overcome this problem, we introduce a density matrix to store the gray values of all the pixels in the image instead of the probability amplitude or the information of the qubit sequence.Inspired by incoherentcoherent states, we propose three models for quantum image representation.When encountering the preparation of mixed states, unlike the previous approach of introducing a reference system to make the mixed system a purely entangled state, we deal with the preparation of mixed states through the mapping relation between mixed and pure states on the Bloch sphere, and utilize the rotation operator acting on the Bloch sphere in the preparation process.These models can better capture the coherence and entanglement between pixels in an image, leading to more accurate image representations.However, our proposed model has some limitations, such as the geometrical relations required during the preparation of the MSR-QCI model, which is so complex that the preparation process is relatively difficult.In addition, we investigate the transmission process of quantum images in noisy channels.This is of great significance to the field of quantum communication and image transmission.The results can help us understand how quantum images are affected by noise and provide corresponding error correction and freezing methods.Through these studies, we give a quantum image representation with practical significance in open quantum systems and provide a theoretical basis and practical method for quantum image transmission.With comparative analysis in our research, we demonstrate the advantages of the new approach over traditional quantum image representation methods, this provides an important background for future related work.Specifically, the use of density matrices representing the quantum images can advance the quantization of quantum images in the coherence domain, enabling the features of coherence changes to be captured in a mathematically rigorous manner.In quantum communication or evolution, we want the transmitted or evolved quantum information to maintain a certain degree of coherence, i.e., not to be distorted by the interference of the external environment that causes the information to be distorted.And frozen coherence helps to obtain a relatively long decoherence time to ensure that the coherence is maintained during gate operation.By studying the frozen coherence, one can design more stable and reliable quantum communication systems and promote the development of quantum technology in practical applications.
The quantification of coherence generally uses three functions: l 1 -norm coherence, relative entropy, and quantum fidelity [38].Similar to the frozen coherence process that we will discuss in the future, the general conditions for frozen coherence based on quantum fidelity could be explored based on the change of quantum fidelity coherence distance before and after quantum image evolution.These different methods of coherence quantization provide us with a more comprehensive view, allowing us to better understand the properties of quantum images and providing new ideas and methods for future researches.

Figure 3 A 2 ×
Figure 3 A 2 × 2 example image and its representative expression in an open quantum system

Figure 5
Figure 5The quantum circuits preparing RGB color information

Figure 6 Problem 1
Figure 6 Relationship between point C and points A and B on the Bloch sphere

Figure 7
Figure 7 Unique decodable code representing a mixed state by two pure states

Figure 9
Figure 9 Equivalent form of controlled operation in quantum circuits

Figure 11
Figure 11 An open quantum system model 3) provides the foundation for the quantification and protection of coherence in open quantum systems.
Here, ρ B i are arbitrary quantum states on quantum subsystem B, and the states |i A belong to the local incoherent basis of quantum subsystem A. Similarly, it is called a coherentincoherent state if it can be written as ρ = i p i ρ A i ⊗ |i i| B .Incoherent-coherent and coherent-incoherent states are incoherent in one subsystem.