Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Consider the following optimization problem: Given a quantum state $\sigma$, a constant $b$ and a Hermitian operator $A$, find

$\underset{\rho} \max F(\rho,\sigma)$

subject to $\text{Tr}(\rho A)>b$,

where $F(\rho,\sigma)=||\sqrt{\rho} \sqrt{\sigma}||_1$ is the fidelity function.

Does anyone know of a fast algorithm for this problem? Even better, an available package that implements such an algorithm?

share|cite|improve this question
up vote 4 down vote accepted

I have luckily found an answer to this question:

The fidelity is concave in both its arguments i.e.

$F(\lambda \rho_1+(1-\lambda)\rho_2,\lambda \sigma_1+(1-\lambda)\sigma_2 )\geq \lambda F(\rho_1,\rho_2)+(1-\lambda)F(\sigma_1,\sigma_2)$.

This means that this maximization problem can be dealt with efficiently with the very well developed tools of convex optimization. The algorithm can be implemented, in particular, with the 'cvx' package for MatLab available in One has only to be careful and write the fidelity function as

$F(\rho,\sigma)=\text{Tr}(\sqrt{\sqrt{\rho}\sigma \sqrt{\rho}})$.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.