A popular and classical method for finding the best rank one approximation of a real tensor is the higher order power method (HOPM). It is known in the literature that the iterative sequence generated by HOPM converges globally, while the convergence rate can be superlinear, linear or sublinear. In this paper, we examine the convergence rate of HOPM in solving the best rank one approximation problem of real tensors. We first show that the iterative sequence of HOPM always converges globally and provide an explicit eventual sublinear convergence rate. The sublinear convergence rate estimate is in terms of the dimension and the order of the underlying tensor space. Then, we examine the concept of nondegenerate singular vector tuples and show that, if the sequence of HOPM converges to a nondegenerate singular vector tuple, then the global convergence rate is R-linear. We show that, for almost all tensors (in the sense of Lebesgue measure), all the singular vector tuples are nondegenerate, and so, the HOPM “typically” exhibits global R-linear convergence rate. Moreover, without any regularity assumption, we establish that the sequence generated by HOPM always converges globally and R-linearly for orthogonally decomposable tensors with order at least 3. We achieved this by showing that each nonzero singular vector tuple of an orthogonally decomposable tensor with order at least 3 is nondegenerate.

In this paper we derive new fractional error bounds for polynomial systems with exponents explicitly determined by the dimension of the underlying space and the number/degree of the involved polynomials. Our major result extends the existing error bounds from the system involving only a single polynomial to a general polynomial system and do not require any regularity assumptions. In this way we resolve, in particular, some open questions posed in the literature. The developed techniques are largely based on variational analysis and generalized differentiation, which allow us to establish, e.g., a nonsmooth extension of the seminal Łojasiewicz’s gradient inequality to maxima of polynomials with explicitly determined exponents. Our major applications concern quantitative Hölderian stability of solution maps for parameterized polynomial optimization problems and nonlinear complementarity systems with polynomial data as well as high-order semismooth properties of the eigenvalues of symmetric tensors.

No abstract uploaded!

No abstract uploaded!

A familiar anomaly affects SU(2) gauge theory in four dimensions: a theory with an odd number of fermion multiplets in the spin 1/2 representation of the gauge group, and more generally in representations of spin 2r+1/2, is inconsistent. We describe here a more subtle anomaly that can affect SU(2) gauge theory in four dimensions under the condition that fermions transform with half-integer spin under SU(2) and bosons with integer spin. Such a theory, formulated in a way that requires no choice of spin structure, and with an odd number of fermion multiplets in representations of spin 4r+3/2, is inconsistent. The theory is consistent if one picks a spin or spin_c structure. Under Higgsing to U(1), the new SU(2) anomaly reduces to a known anomaly of "all-fermion electrodynamics." Like that theory, an SU(2) theory with an odd number of fermion multiplets in representations of spin 4r+3/2 can provide a boundary state for a five-dimensional gapped theory whose partition function on a closed five-manifold Y is (−1)^{∫_Y w2w3}. All statements have analogs with SU(2) replaced by Sp(2N). There is also an analog in five dimensions.