M. Navascues, A. Garcia-Saez, A. Acin, S. Pironio, M. B. Plenio
We show that the recent hierarchy of semidefinite programming approximations based on non-commutative polynomial optimization and reduced density matrix variational methods exhibits an interesting paradox when extended to the bosonic case: even though it can be proven that the hierarchy collapses after the first step, one finds numerically that higher order steps generate a sequence of lower bounds that converges to the optimal solution. We analyze this effect and compare it with similar behavior observed in implementations of semidefinite programming relaxations for classical polynomial minimization. We conclude that the method converges due to the rounding errors occurring during the execution of the numerical program, and show that convergence is lost as soon as computer precision is incremented. We support this conclusion by proving that for any element p of a Weyl algebra which is non-negative in the Schrodinger representation there exists another element p arbitrarily close to p that admits a sum of squares decomposition.
View original:
http://arxiv.org/abs/1203.3777
No comments:
Post a Comment