The aim of this paper is to illustrate both analytically and numerically the interplay of two fundamentally distinct non-Hermitian mechanisms in the deep subwavelength regime. Considering a parity-time symmetric system of one-dimensional subwavelength resonators equipped with two kinds of non-Hermiticity – an imaginary gauge potential and on-site gain and loss – we prove that all but two eigenmodes of the system pass through exceptional points and decouple. By tuning the gain-to-loss ratio, the system changes from a phase with unbroken parity-time symmetry to a phase with broken parity-time symmetry. At the macroscopic level, this is observed as a transition from symmetrical eigenmodes to condensated eigenmodes at one edge of the structure. Mathematically, it arises from a topological state change. The results of this paper open the door to the justification of a variety of phenomena arising from the interplay between non-Hermitian reciprocal and nonreciprocal mechanisms not only in subwavelength wave physics but also in quantum mechanics, where the tight-binding model coupled with the nearest neighbour approximation can be analysed with the same tools as those developed here.

This paper investigates structural changes in the parameters of first-order autoregressive (AR) models by analyzing the edge eigenvalues of the precision matrices. Specifically, edge eigenvalues in the precision matrix are observed if and only if there is a structural change in the AR coefficients. We show that these edge eigenvalues correspond to the zeros of a determinantal equation. Additionally, we propose a consistent estimator for detecting outliers within the panel time series framework, supported by numerical experiments.

We determine the characteristic polynomials of the matrices $[q^{\,j-k}+t]_{1\le \,j,k\le n}$ and $[q^{\,j+k}+t]_{1\le \,j,k\le n}$ for any complex number $q\not =0,1$. As an application, for complex numbers $a,b,c$ with $b\not =0$ and $a^2\not =4b$, and the sequence $(w_m)_{m\in \mathbb Z}$ with $w_{m+1}=aw_m-bw_{m-1}$ for all $m\in \mathbb Z$, we determine the exact value of $\det [w_{\,j-k}+c\delta _{jk}]_{1\le \,j,k\le n}$.

Motivated by the work initiated by Chapman [‘Determinants of Legendre symbol matrices’, Acta Arith. 115 (2004), 231–244], we investigate some arithmetical properties of generalised Legendre matrices over finite fields. For example, letting $a_1,\ldots ,a_{(q-1)/2}$ be all the nonzero squares in the finite field $\mathbb {F}_q$ containing q elements with $2\nmid q$, we give the explicit value of the determinant $D_{(q-1)/2}=\det [(a_i+a_j)^{(q-3)/2}]_{1\le i,j\le (q-1)/2}$. In particular, if $q=p$ is a prime greater than $3$, then

$$ \begin{align*}\bigg(\frac{\det D_{(p-1)/2}}{p}\bigg)= \begin{cases} 1 & \mbox{if}\ p\equiv1\pmod4,\\ (-1)^{(h(-p)+1)/2} & \mbox{if}\ p\equiv 3\pmod4\ \text{and}\ p>3, \end{cases}\end{align*} $$

where $(\frac {\cdot }{p})$ is the Legendre symbol and $h(-p)$ is the class number of $\mathbb {Q}(\sqrt {-p})$.

Williamson’s theorem states that for any $2n \times 2n$ real positive definite matrix A, there exists a $2n \times 2n$ real symplectic matrix S such that $S^TAS=D \oplus D$, where D is an $n\times n$ diagonal matrix with positive diagonal entries known as the symplectic eigenvalues of A. Let H be any $2n \times 2n$ real symmetric matrix such that the perturbed matrix $A+H$ is also positive definite. In this paper, we show that any symplectic matrix $\tilde {S}$ diagonalizing $A+H$ in Williamson’s theorem is of the form $\tilde {S}=S Q+\mathcal {O}(\|H\|)$, where Q is a $2n \times 2n$ real symplectic as well as orthogonal matrix. Moreover, Q is in symplectic block diagonal form with the block sizes given by twice the multiplicities of the symplectic eigenvalues of A. Consequently, we show that $\tilde {S}$ and S can be chosen so that $\|\tilde {S}-S\|=\mathcal {O}(\|H\|)$. Our results hold even if A has repeated symplectic eigenvalues. This generalizes the stability result of symplectic matrices for non-repeated symplectic eigenvalues given by Idel, Gaona, and Wolf [Linear Algebra Appl., 525:45–58, 2017].

As an extension of Sylvester’s matrix, a tridiagonal matrix is investigated by determining both left and right eigenvectors. Orthogonality relations between left and right eigenvectors are derived. Two determinants of the matrices constructed by the left and right eigenvectors are evaluated in closed form.

We first establish a lower bound on the size and spectral radius of a graph G to guarantee that G contains a fractional perfect matching. Then, we determine an upper bound on the distance spectral radius of a graph G to ensure that G has a fractional perfect matching. Furthermore, we construct some extremal graphs to show all the bounds are best possible.

A $D_{\infty }$-topological Markov chain is a topological Markov chain provided with an action of the infinite dihedral group $D_{\infty }$. It is defined by two zero-one square matrices A and J satisfying $AJ=JA^{\textsf {T}}$ and $J^2=I$. A flip signature is obtained from symmetric bilinear forms with respect to J on the eventual kernel of A. We modify Williams’ decomposition theorem to prove the flip signature is a $D_{\infty }$-conjugacy invariant. We introduce natural $D_{\infty }$-actions on Ashley’s eight-by-eight and the full two-shift. The flip signatures show that Ashley’s eight-by-eight and the full two-shift equipped with the natural $D_{\infty }$-actions are not $D_{\infty }$-conjugate. We also discuss the notion of $D_{\infty }$-shift equivalence and the Lind zeta function.

Motivated by considerations of the quadratic orthogonal bisectional curvature, we address the question of when a weighted graph (with possibly negative weights) has nonnegative Dirichlet energy.

We show that the elements of the dual of the Euclidean distance matrix cone can be described via an inequality on a certain weighted sum of its eigenvalues.

A celebrated theorem of Douglas Lind states that a positive real number is equal to the spectral radius of some integral primitive matrix, if and only if, it is a Perron algebraic integer. Given a Perron number p, we prove that there is an integral irreducible matrix with spectral radius p, and with dimension bounded above in terms of the algebraic degree, the ratio of the first two largest Galois conjugates, and arithmetic information about the ring of integers of its number field. This arithmetic information can be taken to be either the discriminant or the minimal Hermite-like thickness. Equivalently, given a Perron number p, there is an irreducible shift of finite type with entropy $\log (p)$ defined as an edge shift on a graph whose number of vertices is bounded above in terms of the aforementioned data.

We prove an analogue of Alon’s spectral gap conjecture for random bipartite, biregular graphs. We use the Ihara–Bass formula to connect the non-backtracking spectrum to that of the adjacency matrix, employing the moment method to show there exists a spectral gap for the non-backtracking matrix. A by-product of our main theorem is that random rectangular zero-one matrices with fixed row and column sums are full rank with high probability. Finally, we illustrate applications to community detection, coding theory, and deterministic matrix completion.

We investigate a class of generalised stochastic complex matrices constructed from the class of all doubly stochastic matrices and a special class of circulant matrices. We determine the exact values of the structured singular values of all matrices in the class in terms of the constant row (column) sum.

We provide a description of the spectrum and compute the eigenvalues distribution of circulant Hankel matrices obtained as symmetrization of classical Toeplitz circulant matrices. Other types of circulant matrices such as simple and Cesàro circulant matrices are also considered.

Associated with every $2n\times 2n$ real positive definite matrix $A,$ there exist n positive numbers called the symplectic eigenvalues of $A,$ and a basis of $\mathbb {R}^{2n}$ called the symplectic eigenbasis of A corresponding to these numbers. In this paper, we discuss differentiability and analyticity of the symplectic eigenvalues and corresponding symplectic eigenbasis and compute their derivatives. We then derive an analogue of Lidskii’s theorem for symplectic eigenvalues as an application.

We study the $L^{q}$-spectrum of measures in the plane generated by certain nonlinear maps. In particular, we consider attractors of iterated function systems consisting of maps whose components are $C^{1+\alpha }$ and for which the Jacobian is a lower triangular matrix at every point subject to a natural domination condition on the entries. We calculate the $L^{q}$-spectrum of Bernoulli measures supported on such sets by using an appropriately defined analogue of the singular value function and an appropriate pressure function.

We propose a robust numerical algorithm for solving the nonlinear eigenvalue problem A(ƛ)x = 0. Our algorithm is based on the idea of finding the value of ƛ for which A(ƛ) is singular by computing the smallest eigenvalue or singular value of A(ƛ) viewed as a constant matrix. To further enhance computational efficiency, we introduce and use the concept of signed singular value. Our method is applicable when A(ƛ) is large and nonsymmetric and has strong nonlinearity. Numerical experiments on a nonlinear eigenvalue problem arising in the computation of scaling exponent in turbulent flow show robustness and effectiveness of our method.

The Toda equation and its variants are studied in the filed of integrable systems. One particularly generalized time discretisation of the Toda equation is known as the discrete hungry Toda (dhToda) equation, which has two main variants referred to as the dhTodaI equation and dhTodaII equation. The dhToda equations have both been shown to be applicable to the computation of eigenvalues of totally nonnegative (TN) matrices, which are matrices without negative minors. The dhTodaI equation has been investigated with respect to the properties of integrable systems, but the dhTodaII equation has not. Explicit solutions using determinants and matrix representations called Lax pairs are often considered as symbolic properties of discrete integrable systems. In this paper, we clarify the determinant solution and Lax pair of the dhTodaII equation by focusing on an infinite sequence. We show that the resulting determinant solution firmly covers the general solution to the dhTodaII equation, and provide an asymptotic analysis of the general solution as discrete-time variable goes to infinity.

Hexagonal grids are valuable in two-dimensional applications involving Laplacian. The methods and analysis are investigated in current work in both linear and nonlinear problems related to anisotropic Laplacian. Ordinary and compact hexagonal grid finite difference methods are developed by elementary arguments, and then analyzed by perturbation for standard Laplacian. In the anisotropic case, analysis is done through reduction to the standard one by using Fourier vectors of mixed types. These hexagonal seven-point methods, with established theoretic stabilities and accuracies, are numerically confirmed in linear and semi-linear anisotropic Poisson problems, and can be applied also in time-dependent problems and in many applications in two-dimensional irregular domains.

For symmetric eigenvalue problems, we constructed a three-term recurrence polynomial filter by means of Chebyshev polynomials. The new filtering technique does not need to solve linear systems and only needs matrix-vector products. It is a memory conserving filtering technique for its three-term recurrence relation. As an application, we use this filtering strategy to the Davidson method and propose the filtered-Davidson method. Through choosing suitable shifts, this method can gain cubic convergence rate locally. Theory and numerical experiments show the efficiency of the new filtering technique.