The performance of quantum gate operations is experimentally determined by how correctly operational parameters can be determined and set, and how stable these parameters can be maintained. In addition, gates acting on different sets of qubits require unique sets of control parameters. Thus, an efficient multidimensional parameter estimation procedure is crucial to calibrate even medium-sized quantum processors. Here, we develop and characterize an efficient calibration protocol to automatically estimate and adjust experimental parameters of the widely used two-qubit Mølmer-Sørensen entangling gate operation in a trapped-ion quantum information processor. The protocol exploits Bayesian parameter estimation methods that include a stopping criterion based on a desired gate infidelity. We experimentally demonstrate a tune-up procedure that leads to a residual median gate infidelity due to miscalibration of $1.3(1)\times {10}^{-3}$, requiring $1200\pm 500$ experimental cycles, while completing the entire gate calibration procedure in less than one minute, which provides a significant speedup over commonly used manual tune-up routines. This approach is applicable to other quantum information processor architectures with known or sufficiently characterized theoretical models.

Realisation of experiments even on small and medium-scale quantum computers requires an optimisation of several parameters to achieve high-fidelity operations. As the size of the quantum register increases, the characterisation of quantum states becomes more difficult since the number of parameters to be measured grows as well and finding efficient observables in order to estimate the parameters of the model becomes a crucial task. Here we propose a method relying on application of Bayesian inference that can be used to determine systematic, unknown phase shifts of multi-qubit states. This method offers important advantages as compared to Ramsey-type protocols. First, application of Bayesian inference allows the selection of an adaptive basis for the measurements which yields the optimal amount of information about the phase shifts of the state. Secondly, this method can process the outcomes of different observables at the same time. This leads to a substantial decrease in the resources needed for the estimation of phases, speeding up the state characterisation and optimisation in experimental implementations. The proposed Bayesian inference method can be applied in various physical platforms that are currently used as quantum processors.

Combinatorial optimization problems are ubiquitous in industry. In addition to finding a solution with minimum cost, problems of high relevance involve a number of constraints that the solution must satisfy. Variational quantum algorithms (VQAs) have emerged as promising candidates for solving these problems in the noisy intermediate-scale quantum stage. However, the constraints are often complex enough to make their efficient mapping to quantum hardware difficult or even infeasible. An alternative standard approach is to transform the optimization problem to include these constraints as penalty terms, but this method involves additional hyperparameters and does not ensure that the constraints are satisfied due to the existence of local minima. In this paper, we introduce a new method for solving combinatorial optimization problems with challenging constraints using VQAs. We propose the multi-objective variational constrained optimizer (MOVCO) to classically update the variational parameters by a multiobjective optimization performed by a genetic algorithm. This optimization allows the algorithm to progressively sample only states within the in-constraints space, while optimizing the energy of these states. We test our proposal on a real-world problem with great relevance in finance: the cash handling problem. We introduce a novel mathematical formulation for this problem, and compare the performance of MOVCO versus a penalty based optimization. Our empirical results show a significant improvement in terms of the cost of the achieved solutions, but especially in the avoidance of local minima that do not satisfy any of the mandatory constraints.

A major challenge for the realization of useful universal quantum computers is achieving high fidelity two-qubit entangling gate operations. However, calibration errors can affect the quantum gate operations and limit their fidelity. To reduce such errors it is desirable to have an analytical understanding and quantitative predictions of the effects that miscalibrations of gate parameters have on the gate performance. In this work, we study a systematic perturbative expansion in miscalibrated parameters of the Mølmer-Sørensen entangling gate, which is widely used in trapped-ion quantum processors. Our analytical treatment particularly focuses on systematic center-line detuning miscalibrations. Via a unitary Magnus expansion, we compute the gate evolution operator, which allows us to obtain relevant key properties such as relative phases, electronic populations, quantum state purity and fidelities. These quantities, subsequently, are used to assess the performance of the gate using the fidelity of entangled states as performance metric. We verify the predictions from our model by benchmarking them against measurements in a trapped-ion quantum processor. The method and the results presented here can help design and calibrate high-fidelity gate operations of large-scale quantum computers.

Population annealing is a variant of the simulated annealing algorithm that improves the quality of the thermalization process in systems with rough free-energy landscapes by introducing a resampling process. We consider the diluted Sherrington-Kirkpatrick Ising model using population annealing to study its efficiency in finding solutions to combinatorial optimization problems. From this study, we find an easy-hard-easy transition in the model hardness as the problem instances become more diluted and associate this behavior to the clusterization and connectivity of the underlying Erdős-Rényi graphs. We calculate the efficiency of obtaining minimum energy configurations and find that population annealing outperforms simulated annealing for the cases close to this hardness peak while reaching similar efficiencies in the easy limits. Finally, it is known that population annealing can be used to define an adaptive inverse temperature annealing schedule. We compare this adaptive method to a linear schedule and find that the adaptive method achieves improved efficiencies while being robust against final temperature miscalibrations.

We introduce SA-FDR, a novel algorithm for $\ell_0$-norm feature selection that considers this task as a combinatorial optimisation problem and solves it by using simulated annealing to perform a global search over the space of feature subsets. The optimisation is guided by the Fisher discriminant ratio, which we use as a computationally efficient proxy for model quality in classification tasks. Our experiments, conducted on datasets with up to hundreds of thousands of samples and hundreds of features, demonstrate that SA-FDR consistently selects more compact feature subsets while achieving a high predictive accuracy. This ability to recover informative yet minimal sets of features stems from its capacity to capture inter-feature dependencies often missed by greedy optimisation approaches. As a result, SA-FDR provides a flexible and effective solution for designing interpretable models in high-dimensional settings, particularly when model sparsity, interpretability, and performance are crucial.