Search

In this paper, we develop an $n$ dimensional volumetric stretch energy ($n$-VSE) functional for the volume-/mass-preserving parameterization of the $n$-manifolds topologically equivalent to $n$-ball. The $n$-VSE has a lower bound and equal to it if and only if the map is volume-/mass-preserving. This motivates us to minimize the $n$-VSE to achieve the ideal volume-/mass-preserving parameterization. In the discrete case, we also guarantee the relation between the lower bound and the volume-/mass-preservation, and propose the spherical and ball volume-/mass-preserving parameterization algorithms. The numerical experiments indicate the accuracy and robustness of the proposed algorithms. The modified algorithms are applied to the manifold registration and deformation, showing the versatility of $n$-VSE.

In this paper, we study the Helmholtz transmission eigenvalue problem for inhomogeneous anisotropic media with the index of refraction $n(x)\equiv 1$ in two and three dimension. Starting with a nonlinear fourth order formulation established by Cakoni, Colton and Haddar [2009], by introducing some auxiliary variables, we present an equivalent mixed formulation for this problem, followed up with the finite element discretization. Using the proposed scheme, we rigorously show that the optimal convergence rate for the transmission eigenvalues both on convex and nonconvex domains can be expected. Moreover, by this scheme, we will obtain a sparse generalized eigenvalue problem whose size is so demanding even with a coarse mesh that its smallest few real eigenvalues fail to be solved by the shift and invert method. We partially overcome this critical issue by deflating the almost all of the $\infty$ eigenvalue of huge multiplicity, resulting in a drastic reduction of the matrix size without deteriorating the sparsity. Extensive numerical examples are reported to demonstrate the effectiveness and efficiency of the proposed scheme.

This paper focuses on studying the bifurcation analysis of the eigenstructure of the $\gamma$-parameterized generalized eigenvalue problem ($\gamma$-GEP) arising in three-dimensional (3D) source-free Maxwell's equations with Pasteur media, where $\gamma$ is the magnetoelectric chirality parameter. For the weakly coupled case, namely, $\gamma < \gamma_{*} \equiv$ critical value, the $\gamma$-GEP is positive definite, which has been well-studied by Chern et.\ al, 2015. For the strongly coupled case, namely, $\gamma > \gamma_{*}$, the $\gamma$-GEP is no longer positive definite, introducing a totally different and complicated structure. For the critical strongly coupled case, numerical computations for electromagnetic fields have been presented by Huang et.\ al, 2018. In this paper, we build several theoretical results on the eigenstructure behavior of the $\gamma$-GEPs. We prove that the $\gamma$-GEP is regular for any $\gamma > 0$, and the $\gamma$-GEP has $2 \times 2$ Jordan blocks of infinite eigenvalues at the critical value $\gamma_{*}$. Then, we show that the $2 \times 2$ Jordan block will split into a complex conjugate eigenvalue pair that rapidly goes down and up and then collides at some real point near the origin. Next, it will bifurcate into two real eigenvalues, with one moving toward the left and the other to the right along the real axis as $\gamma$ increases. A newly formed state whose energy is smaller than the ground state can be created as $\gamma$ is larger than the critical value. This stunning feature of the physical phenomenon would be very helpful in practical applications. Therefore, the purpose of this paper is to clarify the corresponding theoretical eigenstructure of 3D Maxwell's equations with Pasteur media., Comment: 26 pages, 5 figures

We study the non-existence, existence and multiplicity of positive solutions to the following nonlinear Kirchhoff equation:% \begin{equation*} \left\{ \begin{array}{l} -M\left( \int_{\mathbb{R}^{3}}\left\vert \nabla u\right\vert ^{2}dx\right) \Delta u+\mu V\left( x\right) u=Q(x)\left\vert u\right\vert ^{p-2}u+\lambda f\left( x\right) u\text{ in }\mathbb{R}^{N}, \\ u\in H^{1}\left( \mathbb{R}^{N}\right) ,% \end{array}% \right. \end{equation*}% where $N\geq 3,20\right) ,$ the potential $V$ is a nonnegative function in $\mathbb{R}% ^{N}$ and the weight function $Q\in L^{\infty }\left( \mathbb{R}^{N}\right) $ with changes sign in $\overline{\Omega }:=\left\{ V=0\right\} .$ We mainly prove the existence of at least two positive solutions in the cases that $% \left( i\right) $ $20$ sufficiently large, where $\lambda _{1}\left( f_{\Omega }\right) $ is the first eigenvalue of $-\Delta $ in $% H_{0}^{1}\left( \Omega \right) $ with weight function $f_{\Omega }:=f|_{% \overline{\Omega }},$ whose corresponding positive principal eigenfunction is denoted by $\phi _{1}.$ Furthermore, we also investigated the non-existence and existence of positive solutions if $a,\lambda $ belongs to different intervals., Comment: 34 pages

Calculation of band structure of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalue. In this paper, we try to address this problem in such a broad context that band structure of three dimensional isotropic photonic crystals with all 14 Bravais lattices can be efficiently computed in a unified framework. We uncover the delicate machinery behind several key results of our work and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments in our framework. In this work particular effort is made on reformulating the Bloch boundary condition for all 14 Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by systematic use of commutativity among them, which has been overlooked previously, and reducing eigen-decomposition of double-curl operator to the canonical form of a 3x3 complex skew-symmetric matrix under unitary congruence. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm., Comment: 37 pages, 20 figures

To solve the Bethe-Salpeter eigenvalue problem with distinct sizes, two efficient methods, called {\Gamma}QR algorithm and {\Gamma}-Lanczos algorithm, are proposed in this paper. Both algorithms preserve the special structure of the initial matrix $H=\begin{bmatrix}A & B-\overline{B} & -\overline{A}\end{bmatrix}$, resulting the computed eigenvalues and the associated eigenvectors still hold the properties similar to those of $H$. Theorems are given to demonstrate the validity of the proposed two algorithms in theory. Numerical results are presented to illustrate the superiorities of our methods., Comment: 31pages, 4figures

In this paper, we study intermittent behaviors of coupled piecewise-expanding map lattices with two nodes and a weak coupling. We show that the successive phase transition between ordered and disordered phases occurs for almost every orbit. That is, we prove $\liminf_{n\rightarrow \infty}| x_1(n)-x_2(n)|=0$ and $\limsup_{n\rightarrow \infty}| x_1(n)-x_2(n)|\ge c_0>0$, where $x_1(n), x_2(n)$ correspond to the coordinates of two nodes at the iterative step $n$. We also prove the same conclusion for weakly coupled tent-map lattices with any multi-nodes., Comment: 33 pages, 4 figures