We consider the problem of quantum multiparameter estimation with experimental constraints and formulate the solution in terms of a convex optimization. Specifically, we outline an efficient method to identify the optimal strategy for estimating multiple unknown parameters of a quantum process and apply this method to a realistic example. The example is two electron-spin qubits coupled through the dipole and exchange interactions with unknown coupling parameters---explicitly, the position vector relating the two qubits and the magnitude of the exchange interaction are unknown. This coupling Hamiltonian generates a unitary evolution which, when combined with arbitrary single-qubit operations, produces a universal set of quantum gates. However, the unknown parameters must be known precisely to generate high-fidelity gates. We use the Cram\'er-Rao bound on the variance of a point estimator to construct the optimal series of experiments to estimate these free parameters and present a complete analysis of the optimal experimental configuration. Our method of transforming the constrained optimal parameter estimation problem into a convex optimization is powerful and widely applicable to other systems.
We consider the problem of quantum multiparameter estimation with experimental constraints and formulate the solution in terms of a convex optimization. Specifically, we outline an efficient method to identify the optimal strategy for estimating multiple unknown parameters of a quantum process and apply this method to a realistic example. The example is two electron-spin qubits coupled through the dipole and exchange interactions with unknown coupling parameters---explicitly, the position vector relating the two qubits and the magnitude of the exchange interaction are unknown. This coupling Hamiltonian generates a unitary evolution which, when combined with arbitrary single-qubit operations, produces a universal set of quantum gates. However, the unknown parameters must be known precisely to generate high-fidelity gates. We use the Cram'er-Rao bound on the variance of a point estimator to construct the optimal series of experiments to estimate these free parameters and present a complete analysis of the optimal experimental configuration. Our method of transforming the constrained optimal parameter estimation problem into a convex optimization is powerful and widely applicable to other systems.
Estimation of physical parameters encoded in a Hamiltonian is a central task in quantum sensing and learning. While the ultimate precision limit for estimating a single parameter coupled to a single generator is well established, the corresponding bound for estimating a function of multiple parameters-each coupled to distinct and possibly non-commuting generators-remains unknown in general. Here, we derive the ultimate quantum limit and present an estimation protocol for any function of parameters in a general Hamiltonian that attains this bound. We show that, although the task is fundamentally a multiparameter problem, our tight bound reduces to an optimized single-parameter quantum Cramér-Rao bound, even for arbitrary generator sets. Our result unifies and extends previous works, providing a general framework for optimal function estimation in quantum systems.
Estimating multiple unknown parameters simultaneously is essential for practical quantum sensing. However, it faces a fundamental challenge: the optimal strategy for estimating one parameter is often incompatible with that for another, making it impossible to simultaneously achieve the ultimate precision limits for all parameters. Here we develop a general and efficient computational framework that jointly optimizes probe states, control operations, and measurements across different strategy families, including parallel, sequential, and those with indefinite causal order. Our approach provides exact semidefinite-program formulations for several precision bounds, including the Holevo, Nagaoka-Hayashi, and quantum Cramér-Rao bounds. We demonstrate the capabilities of the framework in multiparameter magnetometry and frequency estimation, identifying optimal protocols within each class and revealing a strict hierarchy among the achievable performances of different classes in the multiparameter regime. The framework also directly incorporates resource constraints, such as energy budgets, enabling systematic investigation of experimentally realistic sensing scenarios. Furthermore, we develop a finite-memory optimization method for sequential strategies with restricted ancillary-memory dimension. By decomposing the protocol into initial probe preparation and intermediate control operations, this method provides a practical route to designing resource-constrained sequential sensing schemes. Our work establishes a versatile computational tool for determining fundamental precision limits and designing optimal quantum-sensing protocols in complex multiparameter settings.
A number of problems in quantum state and system identification are addressed. Specifically, it is shown that the maximum likelihood estimation (MLE) approach, already known to apply to quantum state tomography, is also applicable to quantum process tomography (estimating the Kraus operator sum representation (OSR)), Hamiltonian parameter estimation, and the related problems of state and process (OSR) distribution estimation. Except for Hamiltonian parameter estimation, the other MLE problems are formally of the same type of convex optimization problem and therefore can be solved very efficiently to within any desired accuracy. Associated with each of these estimation problems, and the focus of the paper, is an optimal experiment design (OED) problem invoked by the Cramer-Rao Inequality: find the number of experiments to be performed in a particular system configuration to maximize estimation accuracy; a configuration being any number of combinations of sample times, hardware settings, prepared initial states, etc. We show that in all of the estimation problems, including Hamiltonian parameter estimation, the optimal experiment design can be obtained by solving a convex optimization problem. Software to solve the MLE and OED convex optimization problems is available upon request from the first author.
In this technical note, we study the problem of parameter estimation for quantum channels, which are (completely positive) maps acting on quantum states. We focus on an important subclass of channels called Pauli channels which are characterized by a certain set of vectors (directions) and a number of scalar parameters. For the case where the directions are known, a special parametrization turns the parameter estimation problem into a convex optimization problem. For the case of unknown directions we give a simple algorithm to estimate the directions for qubit Pauli channels. These results assume that the identification experiment configuration is given, namely, a set of quantum states are given on which the channel acts, and a set of positive measurement (operator)s are fixed from which information is gathered. In the second part of the technical note we consider the problem of determining the experiment design, namely, determining states and measurements for optimal parameter identification of the channels. We formalize this problem as a maximization problem for Fisher information and, assuming known channel directions, prove that this problem is convex. We also prove that the optimal states to be used in experiments are pure and the optimal measurements are extremal. For qubit Pauli channels we prove that both the optimal pure input states and projective measurements should be directed towards the channel directions. We illustrate the results of the technical note with numerical examples and simulations.
We investigate multiparameter quantum estimation in a Raman-coupled two-qubit system at thermal equilibrium. Analytical expressions for the quantum Fisher information matrix are derived to characterize the simultaneous estimation of the temperature and Raman coupling strength. The corresponding quantum Cramér--Rao bounds are obtained and compared with those of individual estimation strategies. Our results reveal optimal operating regimes determined by the interplay between thermal fluctuations and coherent interactions. In particular, quantum thermometry exhibits a well-defined optimal temperature window, whereas the estimation of the Raman coupling strength is significantly enhanced in the low-temperature and weak-coupling regime. We further show that simultaneous estimation can outperform independent estimation within appropriate parameter regions, highlighting the advantages of multiparameter quantum metrology. These results provide analytical insights into the ultimate precision limits of Raman-coupled two-qubit systems and identify promising operating regimes for quantum sensing and quantum thermometry.
Multiparameter estimation remains a fundamental challenge in quantum estimation theory because the optimal measurements associated with different parameters are generally incompatible. In this work, we develop a Bayesian framework for the simultaneous estimation of multiple parameters. Our approach builds on Personick's method, in which the estimation problem is reduced to a set of Lyapunov equations whose solutions define optimal observables for the individual parameters. Since these observables generally do not commute, we construct a convex combination of the solutions, parameterized by a set of variational parameters. The spectral decomposition of the resulting operator defines a parametrized projection-valued measure. The resulting projective measurement determines the likelihood function, from which the minimum mean-square error estimators and the corresponding Bayesian mean-square errors are obtained, following standard Bayesian procedures, as functions of the variational parameters. To determine their optimal values, we formulate a minimax optimization problem based on the normalized Bayesian mean-square errors, introducing a min-max normalization procedure that enables a meaningful comparison of estimation errors associated with different parameters. This optimization yields a variationally optimized projective measurement for simultaneous Bayesian estimation. Two qubit examples, involving phase estimation and the estimation of parameters defining convex combinations of unitary operations, demonstrate the construction of optimized projective measurements and the corresponding minimum mean-square error estimators within the proposed framework for finite data sets.