Cover Image

This paper presents a framework for incorporating arbitrary implicit multistep schemes into the lattice Boltzmann method. While the temporal discretization of the lattice Boltzmann equation is usually derived using a second-order trapezoidal rule, it appears natural to augment the time discretization by using multistep methods. The effect of incorporating multistep methods into the lattice Boltzmann method is studied in terms of accuracy and stability. Numerical tests for the third-order accurate Adams-Moulton method and the second-order backward differentiation formula show that the temporal order of the method can be increased when the stability properties of multistep methods are considered in accordance with the second Dahlquist barrier.     

The degenerate C. Neumann system I: symmetry reduction and convexity

The C. Neumann system describes a particle on the sphere S ⁿ under the influence of a potential that is a quadratic form. We study the case that the quadratic form has ℓ +1 distinct eigenvalues with multiplicity. Each group of m _σ equal eigenvalues gives rise to an O(m _σ )-symmetry in configuration space. The combined symmetry group G is a direct product of ℓ + 1 such factors, and its cotangent lift has an Ad*-equivariant momentum mapping. Regular reduction leads to the Rosochatius system on S ^ℓ , which has the same form as the Neumann system albeit for an additional effective potential.     

Generating Hyperbolic Singularities in Semitoric Systems Via Hopf Bifurcations

Let \((M,\Omega )\) be a connected symplectic 4-manifold and let \(F=(J,H) :M\rightarrow \mathbb {R}^2\) be a completely integrable system on M with only non-degenerate singularities. Assume that F does not have singularities with hyperbolic blocks and that \(p_1,\ldots ,p_n\) are the focus–focus singularities of F. For each subset \(S=\{i_1,\ldots ,i_j\}\), we will show how to modify F locally around any \(p_i, i \in S\), in order to create a new integrable system \(\widetilde{F}=(J, \widetilde{H}) :M \rightarrow \mathbb {R}^2\) such that its classical spectrum \(\widetilde{F}(M)\) contains j smooth curves of singular values corresponding to non-degenerate transversally hyperbolic singularities of \(\widetilde{F}\). Moreover the focus–focus singularities of \(\widetilde{F}\) are precisely \(p_i\), \(i \in \{1,\ldots ,n\} \setminus S\). The proof is based on Eliasson’s linearization theorem for non-degenerate singularities, and properties of the Hamiltonian Hopf bifurcation.     

Boundary and Lens Rigidity of Lorentzian Surfaces

Let be a Lorentzian metric on the plane that agrees with the standard metric outside a compact set and so that there are no conjugate points along any time-like geodesic of . Then and are isometric. Further, if and are two dimensional compact time oriented Lorentzian manifolds with space--like boundaries and so that all time-like geodesics of maximize the distances between their points and and are ``boundary isometric', then there is a conformal diffeomorphism between and and they have the same areas. Similar results hold in higher dimensions under an extra assumption on the volumes of the manifolds.

     

Reaction matrix for Dirichlet billiards with attached waveguides

The reaction matrix of a cavity with attached waveguides connects scattering properties to properties of a corresponding closed billiard for which the waveguides are cut off by straight walls. On the one hand, this matrix is directly related to the S-matrix, on the other hand it can be expressed by a spectral sum over all eigenfunctions of the closed system. However, in the physically relevant situation where these eigenfunctions vanish on the impenetrable boundaries of the closed billiard, the spectral sum for the reaction matrix, as it was used before, fails to converge and does not reliably reproduce the scattering properties. We derive here a convergent representation of the reaction matrix in terms of eigenmodes satisfying Dirichlet boundary conditions and demonstrate its validity in the rectangular and the Sinai billiards.     

Suprema of Chaos Processes and the Restricted Isometry Property

We present a new bound for suprema of a special type of chaos process indexed by a set of matrices, which is based on a chaining method. As applications we show significantly improved estimates for the restricted isometry constants of partial random circulant matrices and time‐frequency structured random matrices. In both cases the required condition on the number m of rows in terms of the sparsity s and the vector length n is m ? s log² s log² n. © 2014 Wiley Periodicals, Inc.     

Low rank matrix recovery from rank one measurements

We study the recovery of Hermitian low rank matrices $X\in {\mathbb{C}}^{n\times n}$ from undersampled measurements via nuclear norm minimization. We consider the particular scenario where the measurements are Frobenius inner products with random rank-one matrices of the form $a_j{a}_j^{?}$ for some measurement vectors $a_1,\dots ,a_m$, i.e., the measurements are given by $b_j=\mathrm{tr}(X{a}_j{a}_j^{?})$. The case where the matrix $X=x{x}^{?}$ to be recovered is of rank one reduces to the problem of phaseless estimation (from measurements $b_j=|〈x,a_j〉{|}^2$) via the PhaseLift approach, which has been introduced recently. We derive bounds for the number m of measurements that guarantee successful uniform recovery of Hermitian rank r matrices, either for the vectors $a_j$, $j=1,\dots ,m$, being chosen independently at random according to a standard Gaussian distribution, or $a_j$ being sampled independently from an (approximate) complex projective t-design with $t=4$. In the Gaussian case, we require $m\ge Crn$ measurements, while in the case of 4-designs we need $m\ge \mathit{Cr}n\log ?(n)$. Our results are uniform in the sense that one random choice of the measurement vectors $a_j$ guarantees recovery of all rank r-matrices simultaneously with high probability. Moreover, we prove robustness of recovery under perturbation of the measurements by noise. The result for approximate 4-designs generalizes and improves a recent bound on phase retrieval due to Gross, Krahmer and Kueng. In addition, it has applications in quantum state tomography. Our proofs employ the so-called bowling scheme which is based on recent ideas by Mendelson and Koltchinskii.