This article focuses on the discrete double-curl operator arising in the Maxwell equation that models three dimensional photonic crystals with face centered cubic lattice. The discrete double-curl operator is the degenerate coefficient matrix of the generalized eigenvalue problems (GEVP) due to the Maxwell equation. We derive an eigendecomposition of the degenerate coefficient matrix and explore an explicit form of orthogonal basis for the range and null spaces of this matrix. To solve the GEVP, we apply these theoretical results to project the GEVP to a standard eigenvalue problem (SEVP), which involves only the eigenspace associated with the nonzero eigenvalues of the GEVP and therefore the zero eigenvalues are excluded and will not degrade the computational efficiency. This projected SEVP can be solved efficiently by the inverse Lanczos method. The linear systems within the inverse Lanczos method are well-conditioned and can be solved efficiently by the conjugate gradient method without using a preconditioner. We also demonstrate how two forms of matrix-vector multiplications, which are the most costly part of the inverse Lanczos method, can be computed by fast Fourier transformation due to the eigendecomposition to significantly reduce the computation cost. Integrating all of these findings and techniques, we obtain a fast eigenvalue solver. The solver has been implemented by MATLAB and successfully solves each of a set of $5.184$ million dimension eigenvalue problems within $50$ to $104$ minutes on a workstation with two Intel Quad-C ore Xeon X5687 3.6 GHz CPUs.

We provide a refinement of the Poincaré inequality on the torus $\mathbb{T}^{d}$ : there exists a set $\mathcal{B} \subset \mathbb{T} ^{d}$ of directions such that for every $\alpha \in \mathcal{B}$ there is a $c_{\alpha } > 0$ with $$\begin{aligned} \|\nabla f\|_{L^{2}(\mathbb{T}^{d})}^{d-1} \| \langle \nabla f, \alpha \rangle \|_{L^{2}(\mathbb{T}^{d})} \geq c_{\alpha }\|f\| _{L^{2}(\mathbb{T}^{d})}^{d} \quad \mbox{for all}~f\in H^{1}\bigl( \mathbb{T}^{d}\bigr)~ \mbox{with mean 0.} \end{aligned}$$ The derivative $\langle \nabla f, \alpha \rangle $ does not detect any oscillation in directions orthogonal to $\alpha $ , however, for certain $\alpha $ the geodesic flow in direction $\alpha $ is sufficiently mixing to compensate for that defect. On the two-dimensional torus $\mathbb{T}^{2}$ the inequality holds for $\alpha = (1, \sqrt{2})$ but is not true for $\alpha = (1,e)$ . Similar results should hold at a great level of generality on very general domains.

We introduce a non-linear injective transformation τ from the set of non-vanishing normalized Hausdorff moment sequences to the set of normalized Stieltjes moment sequences by the formula$T$[($a$_{$n$})_{$n$=1}^{∞}]_{$n$}= 1/$a$_{1}...$a$_{$n$}. Special cases of this transformation have appeared in various papers on exponential functionals of Lévy processes, partly motivated by mathematical finance. We give several examples of moment sequences arising from the transformation and provide the corresponding measures, some of which are related to$q$-series.

Consider an age-dependent, single-species branching process defined by a progeny number distribution, and a lifetime distribution associated with each independent particle. In this paper, we focus on the associated inverse problem where one wishes to formally solve for the progeny number distribution or the lifetime distribution that defines the Bellman-Harris branching process. We derive results for the existence and uniqueness (the identifiability) of these two distributions given one of two types of information: the extinction time probability of the entire process (extinction time distribution), or the distribution of the total number of particles at one fixed time. We demonstrate that perfect knowledge of the distribution of extinction times allows us to formally determine either the progeny number distribution or the lifetime distribution. Furthermore, we show that these constructions are unique. We then consider “data” consisting of a perfectly known total number distribution given at one specific time. For a process with known progeny number distribution and exponentially distributed lifetimes, we show that the rate parameter is identifiable. For general lifetime distributions, we also show that the progeny distribution is globally unique. Our results are presented through four theorems, each describing the constructions in the four distinct cases.

No abstract uploaded!