1. Introduction
Quantum tomography, as a significant framework in quantum characterization, verification, and validation (QCVV), is indispensable for manufacturing and developing quantum computers [
1,
2]. The term “tomography” originates from the Greek words “tomos”, meaning “description” or “writing”. In the context of quantum tomography as depicted in
Figure 1, the term is used to describe a process that involves reconstructing a quantum state by obtaining information about its various slices or projections. The results of tomography are generally considered as the foundation of quantum device calibration [
3,
4,
5,
6,
7], quantum error mitigation (QEM) [
8], and quantum error correction (QEC) [
9,
10,
11], providing detailed information about qubits, quantum gates, as well as the quantum noise. In this framework, an experimenter prepares a set of experiments consisting of quantum states, circuits, and measurements. Then, a set of data is collected by executing the prepared experiments. Results of interest are derived by performing estimation algorithms.
Based on the common skeleton but different targets, the main techniques are quantum state tomography (QST) [
6,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22], quantum process tomography (QPT) [
23,
24,
25,
26,
27,
28,
29,
30,
31], and gate set tomography (GST) [
5,
32,
33,
34,
35], process tensor tomography (PTT) [
36,
37], and instrument set tomography (IST) [
38]. The effectiveness of these methods in characterizing the target components within specified assumptions has been proven. Furthermore, some methods have been used in developing real quantum hardware [
4,
39,
40,
41].
Tomography tools typically operate under the assumption that components other than the target are knowable. This implies that the experimenter possesses complete knowledge of non-target components or can acquire their information implicitly through learning algorithms. Relaxing knowable assumptions indicates the progress of comprehensiveness and practicality at a cost of complexity in terms of the number of experiments required and the post-processing.
The tomography first emerged by proposing QST [
23,
24,
25,
26,
27,
28,
29,
30,
31] with the assumption that quantum gates and measurements are known to the experimenter. Therefore, one can construct a complete basis of the Hilbert–Schmidt space corresponding to the target quantum state. The black-box quantum state is mathematically reconstructed utilizing measurement probabilities and the knowledge of quantum gates and measurements. Analogously, the QPT [
23,
24,
25,
26] follows this approach, where the target is a black box quantum process, and the knowable assumptions are made to non-target quantum gates, quantum states, and measurements. Utilizing non-target quantum gates to compose symmetry state preparation and measurement (SPAM) circuits and combining them with the known quantum states and measurements, complete bases are constructed before and after the black-box quantum process, respectively. Consequently, the QPT tries to recover a completely positive (CP) trace-preserving (TP) map between the input and output quantum states to represent the target process.
Developing from the QPT, GST was proposed based on the fact that all components are noisy. While applying tomography to quantum gates, QPT becomes inaccurate in the presence of SPAM errors. GST addressed this issue by independently modeling quantum gates, initial states, and measurements as tomography targets. Therefore, the GST explicitly separates the SPAM error into intrinsic and extrinsic errors. Intrinsic errors are inherent in the state preparation and measurement process, while extrinsic errors are introduced by performing quantum gates in SPAM circuits, facilitating QEC and QEM. Furthermore, GST exhibits gauge freedom up to a similarity transformation by a gauge matrix. This implies that an initialization error cannot be easily distinguished from a faulty measurement. As a result, GST provides a systematic method to simultaneously characterize quantum gates, initial states, and measurements. However, the scaling with the system size is polynomially worse than QPT due to the requirement to characterize multiple gates at once.
QPT and GST do not account for the dynamic influence of the environment on quantum processes. Nevertheless, compelling evidence, as highlighted in recent research [
37], suggests the existence of multi-time system–environment (SE) correlations. Neglecting these SE correlations may introduce systematic errors in the tomography process. For instance, using QPT or GST to perform tomography on a black-box quantum process with SE correlations in the state preparation may cause the theoretical violation of CP constraints [
37]. This can disrupt tomography under the Markovian model and degrade the effectiveness of QEC codes.
PTT [
36,
37] initially aims to operationally characterize SE correlations. Based on the quantum stochastic process [
42], a multi-time-step non-Markovian system can be modeled by a principal quantum system [
43], an environment with the same dimensions as the system, as SE unitaries act on both the system and environment dimensions between adjacent time steps, and instruments that the experimenter applies to the system at time steps. In NISQ settings, those instruments, quantum gates, and measurements are the only components accessible to the experimenter. The PTT exploits the process tensor to represent all inaccessible parts consisting of the initial SE state and SE unitaries. By combining informationally complete instruments that span the space of CP and trace non-increasing (TNI) quantum operations, the process tensor is recovered. However, the PTT requires an exponential number of experiments with respect to the Markovian order, which limits large-scale implementations.
Growing out of PTT, the IST [
38] was proposed to perform GST in non-Markovian situations. It tackles the inconsistency in characterizing SE correlations and instruments. For example, two inconsistent process tensors may be generated by two sets of imperfect instruments. Both SE correlations and instruments are targets. The IST successfully estimates the instrument and the SE correlations that are consistent with the measurement probabilities. However, it still suffers from the exponential number of experiments and the high complexity of post-processing.
In this review, we provide a fundamental introduction to quantum tomography techniques, including QST, QPT, GST, PTT, and IST. We briefly introduce the basics and direction of the development of these techniques for readers to facilitate systematic comprehension. This review contains many abbreviations. To facilitate prompt reference for readers, we have compiled the key abbreviations in
Table 1, aiming to enhance accessibility and comprehension.
This review is structured as follows. In
Section 2, we briefly introduce mathematical representations for quantum tomography. Then, the basic frameworks for Markovian and non-Markovian quantum tomography are introduced in
Section 3 and
Section 4, respectively. In
Section 5, we introduce the maximum likelihood estimation method as an indispensable statistical tool in the standard toolkit. Finally, further developments to enhance the performance of tomography are discussed in
Section 6.
2. Preliminaries
Prior to delving into quantum tomography techniques, it is essential to establish mathematical representations for physical entities. This section introduces the mathematical expressions for quantum states, gates, and measurements. Subsequently, a calligraphic notation will be employed to depict the physical entities associated with quantum operations, such as quantum gates (
) and measurements (
). The quantum state is written in the Roman alphabet, such as
. We primarily utilize the Pauli transfer matrix (PTM) formalism. Then, the Choi–Jamiolkowski isomorphism (CJI) representation is additionally introduced [
37,
44].
In the PTM representation, a quantum state
on a
d-dimensional Hilbert space
is represented by a superoperator
as a vector in the
-dimensional Hilbert–Schmidt space. We specify the basis of Hilbert–Schmidt space to be Pauli basis without loss of generality. Hence, the elements of
represent the coefficients of the corresponding normalized Pauli matrix [
5],
where
is the
k-th normalized Pauli matrix. Furthermore, a quantum measurement, which is described as an operator
M, can be represented by a set of row vectors [
5], i.e.,
where
. Let
denote
, the probability of performing measurement
M on
, be represented as [
5]
A quantum process
can be conveniently represented by a PTM defined as
This definition describes the mapping between input and output quantum states. Notably, the output state
is represented by the left product of input state
that
Hence, the experimental probability result of the experiment, which involves a quantum state
, quantum process
, and a measurement
, can be represented as
The PTM constructs a mapping from the input state to the output state. It is easy to check and constrain the trace characteristics by examining the first row of PTM. For example, the trace-preserving (TP) constraints indicate the first row of PTM to be
[
5]. However, it is difficult to intuitively represent the CP property.
Therefore, we also introduce the CJI representations. For a quantum process
, the CJI matrix is defined as [
37,
44]
where the output state can be determined by
Since the CJI matrix represents a quantum state without normalization, it is also referred to as the Choi state. The CP constraints for are equivalent to the positive semi-definiteness of the .
Additionally, PTM representation
A can be easily transformed into Choi–Jamiolkowski isomorphism (CJI) [
37,
44] representation
, and vice versa, that
3. Basic Markovian Tomography Framework
Considering a Markovian quantum experiment conducted by an experimenter with specified components consisting of quantum state
, circuit
, and measurement
, the probability of the experiment, whose samples are the only accessible outputs to the experimenter, is given by [
5]
where the PTM representation is
Current quantum devices generally provide limited initial quantum states and quantum measurements. Therefore, the specified quantum state and measurement are referred to as the state preparation and measurement (SPAM) circuits that act on the initial states and native measurements. Then, the probability can be rewritten in the following form [
5]
where
and
represent the state preparation circuit and measurement circuit, respectively. In the following, the quantum state and measurement are referred to as the SPAM circuits acting on the initial state and native measurement, respectively. Moreover, the specified quantum states and measurements with index in experiments are implemented by specifying the SPAM circuits without loss of generality [
5]
The basic Markovian quantum tomography framework can be summarized in the following steps:
- (1)
Experiment preparation. Prepare a set of experiments, where each experiment consists of the quantum state, circuit, and measurement.
- (2)
Data collection. Execute the prepared experiments and record the measurement samples.
- (3)
Tomography reconstruction. Reconstruct the tomographic target by performing a post-process algorithm based on the collected data.
In general, quantum tomography techniques assume that components except for tomographic targets are perfectly implemented as the knowledge to the experimenter, or implicitly learnable by learning algorithms. Experimenters usually prepare experiments that form a tomographically complete basis for the target to obtain sufficient information.
In this section, based on the linear inversion, we will introduce the basic methods of QST, QPT, and GST, i.e., LQST, LQPT, and LGST, respectively.
3.1. Basic QST
The task of QST is to find the explicit representation of an unknown quantum state
[
6,
13,
14,
15,
16,
17,
18]. The construction of a tomographically complete basis for the state is constructed by combining quantum circuits and measurements, which can be solely accomplished by measurements with identity circuits.
Let the set of measurement
be tomographically complete so that
are linear independent. The probability of sampling a desired output value of
is
which indicates the decomposition of the unknown quantum state
. Since the measurements are known to the experimenter, the dual set of
can be determined as
by performing the linear inversion, where
.
Then, the unknown quantum state can be reconstructed by
where
is the estimated probability,
is the number of the desired output value of
in samples, and
is the total sampling times.
The LQST requires at least experiments to collect complete information on the unknown quantum state, and the linear inverse process on a matrix. The reconstructed state is represented by a -dimensional real vector.
3.2. Basic QPT
The task of QPT is to find the explicit representation of an unknown quantum process
, where
represents the space of bounded linear operator on the operand space [
23,
24,
25,
27,
28,
29,
30,
31]. The process
can be treated as a map from the space of the input state to the space of the output state. Therefore, to construct a tomographically complete basis for the unknown quantum process, QPT requires a set of quantum states
and a set of quantum measurements
that simultaneously span
, i.e., both
and
are linear independent.
Given
and
, the probability is
which indicates the coefficient of the map from
to
. Based on the known knowledge of quantum states and measurements, the dual sets of
and
can be determined by
and
, respectively, by performing the linear inversion. Then, the quantum process can be reconstructed by
where
is the estimated probability.
It can be observed that the construction of a tomographically complete basis requires experiments and processes the linear inversion on the matrix twice. The reconstructed quantum process is saved as a -dimensional PTM.
3.3. Basic GST
QST and QPT typically assume the perfect implementation of non-target components, serving as knowledge for the experimenter. However, nothing is absolutely perfect. While using QPT to characterize a set of quantum gates that are available to a quantum device, errors introduced by SPAM circuits impact the self-consistency of tomography, which motivate the development of GST [
5,
32,
33,
34,
35]. Compared to QST and QPT, the goal of GST is to self-consistently and completely characterize a set of
quantum gates
, an initial quantum state
, and a native quantum measurement
. Based on the probability defined in (
11), a quantum gate set is defined as [
5]
where
is a two-outcome measurement and
is the
null gate, which does nothing for no time. Based on the given gate set
, SPAM gates are defined as
where
is composed of gates in the set
, that is,
where
is the length of the
i-th SPAM gate. Then, we can also obtain
tomographically complete input quantum states
and
measurement operators
. Therefore, we can obtain
probabilities,
Likewise, we can rewrite the above equality in the matrix form,
where
Experimentally measuring the values
, we can obtain matrices
, whose elements satisfy
Furthermore, we define
and
, which can be experimentally estimated with and satisfy
where
is measurable.
Then, the unknown gate set can be recovered by performing convex optimization technology. Based on high-quality quantum devices, it is known a priori that the measured gates will differ from an ideal (target) set of gates by some very small error. Therefore, we have a target set
and define
Then, we can find the estimation of matrix
B via the following convex optimization problem,
where
and
. Lastly, we can recover the unknown gate set with
4. Basic Non-Markovian Tomography Framework
No system is isolated. Based on the quantum stochastic process theory [
42,
45], an open quantum system subjected to
k-time-step interventions by quantum operations (instruments), such as quantum processes and measurements, can be modeled as shown in
Figure 2. It utilizes a principal quantum system and an environment to represent the entire SE dimensions. At time step
t, the experimenter applies an instrument
from the
t-available set
to the system. The instrument transforms the system state and outputs a value
as the observation of the experimenter. Without loss of generality, the output
is treated as the index of the instrument applied at time step
t. Notably, there exists an SE unitary between adjacent time steps to represent the non-Markovian SE evolution. Therefore, the probability of obtaining a sequence of output value
is represented by
where we simplify the
and
into
and
, respectively, without confusing them.
Non-negligible non-Markovian SE correlations of open quantum systems theoretically impact the tomography techniques designed under Markovian settings. A simple example to demonstrate the impact is the theoretical violation of CP constraints given in Example 1 in [
33,
37] (c.f. Example A1 in
Appendix A). The requirement of characterization of non-Markovian quantum correlations motivates non-Markovian quantum tomography techniques, such as process tensor tomography (PTT) and instrument set tomography (IST). Non-Markovian quantum tomography also exploits the three-step skeleton as shown in
Section 3, where each experiment in the first step of
Experiment Preparation consists of instruments.
The target in the non-Markovian tomography should be further exclaimed. Based on the criterion of accessibility, the components are classified into accessible and inaccessible parts. Specifically, the only accessible components to the experimenter are the instruments acting on the principal quantum system, while the initial SE state and SE unitaries remain inaccessible. Quantum stochastic process theory [
42,
45] utilizes the process tensor to represent all inaccessible parts as a map from the product of
k spaces of CP and trace non-increasing (TNI) bounded operators on quantum states to the space of output quantum states. In this review, we principally adopt the definitions of process tensor
that
with physical properties:
- (P1)
Linearity. , for any .
- (P2)
Complete positivity. , where is the identity process on the ancilla, for any instruments act on the system and the ancilla.
- (P3)
Containment. is contained in , where .
Recently, the most acceptable mathematical representation of process tensor is the CJI matrix
where ★ represents the link product defined in [
37],
is the Choi state of SE unitary
. This formulation indicates that the SE correlations of an open quantum system with
k time step instrument interventions can be represented by a Choi state consisting of
subsystems labeled by
, respectively. In this way, the output state can be determined by
where
, and
represent the Choi state of
. It is easy to verify the linearity and complete positivity of
. Furthermore, the process tensor of a sub-time span can be derived as
for any
because of the causality, where
, and partial trace with label
means partial trace other dimensions except for •.
In this section, based on the linear inversion, we will introduce the basic methods of PTT and IST, i.e., LPTT and LIST, respectively.
4.1. Basic PTT
Linear inverse process tensor tomography (LPTT) provides a common and quick method for PTT [
36,
42]. The process tensor can be treated as a map from the space of
k-time-step instruments to the probability. The tomography of process tensor requires the informationally complete (IC) available sets at all time steps, which means that the instruments in each available set span the space of the bounded linear operator on the quantum state. Furthermore, we assume that instruments are perfectly implemented as the knowledge to the experimenter.
The experimenter prepares a set of experiments labeled by that consist of all combinations of instruments. By performing the experiments, the experimenter collects estimated probabilities .
Let
be the dual set of
t-available set
that
. The process tensor can be reconstructed by
Given a sequence of output
, it can be verified that
Since the IC available set at each time step consists of instruments, LPTT requires experiments to characterize the entire process tensor. However, the linear inverse procedures are independently applied on each time step to determine the dual set of the available set. This indicates k independent linear inverse procedures on the complex matrix. Nevertheless, it still requires the reconstruction of the matrix of the process tensor.
4.2. Basic IST
The emergence of instrument set tomography (IST) results from non-avoidable imperfect implementations of instruments while performing PTT [
38]. The differences between the knowledge and the practical implementations of instruments may lead to the systematic inconsistency of tomographic results. A simple example to show the inconsistency is that PTT generates two different process tensors utilizing two sets of imperfectly implemented IC instruments, respectively. Moreover, it is not always reasonable to require the quantum device to be characterized as this can formulate an IC available set at each time step. Recent superconductive quantum devices can hardly provide quantum gates that span the space of unitary when a time step is specified to be a time slot for applying a quantum gate [
33,
38]. Tomographic results should be responsible for all available quantum instruments that the quantum device can provide, which results in the requirement of the self-consistency of the IST.
Equation (
38) indicates that, given the initial SE state, SE unitaries, and instruments, the probability can be determined. By exploiting the process tensor to represent all inaccessible components, the probability can be determined by specifying the instruments and the process tensor. Therefore, the full and reduced instrument set can be defined as
respectively, where
represents available sets,
is the set of SE unitaries, and
represents the initial SE state.
LIST consists of two sub-procedures. First, LIST fixes the linear relationship of available sets by measurement probabilities and knowledge of instruments. For time step
t, the probability described in (
39) can be reformed as
where
and
represent the output values before and after time step
t in a
k-time step experiment, respectively. This indicates the decomposition of
on the inaccessible non-orthogonal basis
By adjusting
and
, the LIST tries to make the basis
complete. Utilizing a unique integer
to denote
as the index of the basis matrix, the LIST connects the measured probabilities and the basis as
where
and
represent the vectorization of the
and
, respectively. If
is invertible, we can obtain the instruments
where
and
.
Assuming that the quantum instruments are well implemented and approximately close to the ideal instruments, the inaccessible basis matrix can be optimized by
where
represents the knowledge of instruments to the experimenter. Then, the PTMs of instruments are recovered by the devectorization of determined
in
.
Then, the LIST performs the process tensor reconstruction to recover the process tensor. Utilized instruments are the maximum linear independent set of instruments in the corresponding available sets. Details of process tensor reconstruction can be referred to in
Section 4.1.
Note that the requirement of basis matrix optimization implies gauge freedom, which is also present in the GST. The optimal basis matrix is non-unique and results in the non-uniqueness of tomographic results. This implies that correlations introduced before the time step cannot be distinguished by correlations introduced after.
5. Maximum Likelihood Estimation-Based Methods
Maximum likelihood estimation (MLE)-based methods [
5,
32,
37,
38,
46,
47] grew out of the drawbacks of basic linear inverse methods. It should be highlighted that the precision of linear inverse methods highly relies on the accuracy of estimated probability generated by measurement sampling, which is of the order
, where
is the number of samples. The probability estimation of finite measurement sampling introduces the sampling error, leading to a biased estimated probability. Subsequently, the result of the tomography target may not be physical, such as the CP violation of the estimated quantum state in QST. Experimenters may be interested in utilizing more experiment data to enhance the performance of the tomography. However, linear inverse methods lack the capability of extracting additional information from overcomplete data.
Addressing these issues, statistical methods are proposed based on the MLE, which enables the exploitation of overcomplete data to reduce the error of estimation. These methods compose the measurement data into a likelihood function with respect to the estimated target which requires maximization. The parameterization of the tomographic target also varies from the method. Typically, quantum states, processes, and measurements can be directly parameterized by -dimensional real column vectors, -dimensional real matrices, and -dimensional real row vectors, respectively, as their PTM representations.
The parameterization of SE correlations in non-Markovian situations is more complicated. While representing SE correlations by the process tensor, it can be modeled by the Cholesky decomposition of the Choi state which naturally satisfies the CP constraints. The PTM representation of the process tensor is also available as described in
Section 4. Furthermore, the SE correlations can also be represented by the initial SE state and SE unitaries, where the SE state can be represented by a
-dimensional real column vector and each SE unitary can be modeled by a
-dimensional real vector
whose
i-th element represents the rotation angle on the
-th Pauli matrix, i.e., the estimated SE unitary is
where
.
Although several kinds of likelihood functions are proposed, we apply the most common likelihood as described in
Appendix B, where maximizing
is equivalent to minimizing
with respect to the desired estimator
as described in their basic framework. Remarkably, constraints are introduced in MLE methods to guarantee that the result is physical. The details of constraints for tomography components are listed below:
State: A quantum state
is constrained to be completely positive, meaning its density matrix must be positive semi-definite with a unit trace,
This guarantees that the probabilities of the system are positive with summation 1. Note that TP constraints can be efficiently represented by
Process: A quantum process is constrained to be CPTP. The CP constraints require the Choi state of
to be completely positive,
which indicates that the process remains the CP property of the quantum state. Meanwhile, entries of the first row of its PTM are limited as
to guarantee the TP property, which indicates that the summation of probabilities of the quantum state is 1.
Measurement: A quantum measurement
and its complementary are constrained by a complete positive that
representing the positive outcome of probabilities.
Process tensor: The process tensor
should be constrained by CP and causality. The CP constraints require the Choi state of the process tensor to be positive semi-definite,
Causality constraints result from the fact that future events do not influence past statistics, as a generalization of trace preservation, which in the Choi state representation is
By superoperators of PTM representation, the causality constraints can be conveniently represented by
Then, an optimization algorithm is performed to optimize the optimization problem
The workflow of MLE methods is demonstrated in
Figure 3.
MLE methods fill the blank of constrained statistical methods within the standard quantum tomography toolkits and obtained wild acceptance from researchers. They provide physical implementable results by introducing constraints. The capability of utilizing overcomplete data promotes the accuracy of the tomography by extracting additional information compared with basic linear inverse methods. However, MLE methods are quite computationally resource-consuming. In the classic memory aspect, all MLE methods require saving exponentially large data with respect to the system dimension d to perform the optimization algorithm. In the time aspect, MLE methods require more experiments than the linear inverse methods which is exponentially large to d. Moreover, solving constrained optimization problems with an exponentially large number of parameters with respect to d is extremely time-consuming. Therefore, the recent applications of MLE methods are limited to small-scale quantum systems.
6. Further Developments
The advancement of quantum tomography performance necessitates further development built upon standard methods, including linear inverse and maximum likelihood estimation (MLE) techniques. First, the effectiveness of frequency-based standard methods requires further consideration. The estimated probabilities are biased with finite sampling, which results in the possibility of the violation of physical implementability as it may occur in linear inverse methods. It seems that this problem has been solved by MLE methods. However, it can be verified that the MLE outputs the same result as linear inverse methods when using the same set of measurement data without constraints. The introduction of constraints forces the output to satisfy the physical requirement, which means that the result lies on the boundary of feasible solutions, as shown in
Figure 4. This indicates that no reasonable error bar can be applied in this estimation.
Bayesian tomographic methods are proposed for this systematic problem [
14,
26,
48,
49]. These methods propose meaningful formations for the estimation error of the target. The tomographic targets are assumed to obey a probability distribution. After optimization, the final target is set as the expectation of the distribution, which naturally has a reasonable error bar corresponding to the variance. Nevertheless, Bayesian methods are quite time-consuming due to the construction of distribution. Efficient Bayesian estimation for quantum tomography needs further research.
Then, the efficiency indicated by the accuracy gained from the unit number of experiments and the unit time of the post-process algorithm should be further improved. As previously introduced, the linear inverse provides relatively quick methods for tomography, but the accuracy is bounded by the reciprocal of the number of samples without the capability of utilizing additional data to improve the accuracy. Moreover, sampling errors may cause the physical violations of results. The MLE methods have the capability to utilize additional data with physical constraints. However, performing optimization algorithms to solve constrained optimization problems with an exponentially large number of parameters is time-consuming, which limits these methods to small-scale systems. Recently, researchers have proposed several methods to enhance the efficiency of quantum tomography.
Following the trend of the development of machine learning, a quantum tomography scheme utilizing the neural network has been proposed [
50,
51,
52]. Typically, training and testing datasets consist of measurement results and prior tomographic results. Then, the neural network explicitly learns the map from measurement data to the tomographic target, and implicitly learns the knowledge of non-target components. To a certain extent, the absence of knowledge of non-target components mitigates the impact of imperfect implementations. Remarkably, when the training has been finished, neural-network-based methods have obtained higher accuracy with the same number of probability data inputs compared with the linear inverse method. It should be noted that neural-network-based methods may not generate physically implementable results due to the absence of constraints. However, physical constraints can be introduced into the cost function while performing neural network training, which results in time-consuming neural network training.
Another approach to enhance efficiency involves leveraging the prior symmetrical structure of the tomographic target. In other words, incorporating certain assumptions about the tomographic target can prove beneficial. One well-known assumption is employed in compressed tomography [
13,
53,
54], which posits that the tomographic target exhibits a low-rank characteristic. This methodology has evolved from matrix completion techniques, where the complete target matrix is reconstructed from incomplete elements using random sampling and a nuclear norm minimization algorithm. For instance, consider an unknown quantum state
with
. By employing this approach, an estimate
can be derived with only
measurement settings. Additionally, other prior symmetry structural assumptions, such as matrix product states [
55], permutationally invariant states [
56], and similar concepts, also significantly contribute to enhancing efficiency.