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.

Generalised centrohermitian and skew-centrohermitian matrices arise in a variety of applications in different fields. Based on the vibrating structure equation where M, D, G, K are given matrices with appropriate sizes and x is a column vector, we design a new vibrating structure mode. This mode can be discretised as the left and right inverse eigenvalue problem of a certain structured matrix. When the structured matrix is generalised centrohermitian, we discuss its left and right inverse eigenvalue problem with a submatrix constraint, and then get necessary and sufficient conditions such that the problem is solvable. A general representation of the solutions is presented, and an analytical expression for the solution of the optimal approximation problem in the Frobenius norm is obtained. Finally, the corresponding algorithm to compute the unique optimal approximate solution is presented, and we provide an illustrative numerical example.

The hyperbolic quadratic eigenvalue problem (HQEP) was shown to admit Courant–Fischer type min–max principles in 1955 by Duffin and Cauchy type interlacing inequalities in 2010 by Veselić. It can be regarded as the closest analog (among all kinds of quadratic eigenvalue problem) to the standard Hermitian eigenvalue problem (among all kinds of standard eigenvalue problem). In this paper, we conduct a systematic study on the HQEP both theoretically and numerically. On the theoretical front, we generalize Wielandt–Lidskii type min–max principles and, as a special case, Fan type trace min/max principles and establish Weyl type and Wielandt–Lidskii–Mirsky type perturbation results when an HQEP is perturbed to another HQEP. On the numerical front, we justify the natural generalization of the Rayleigh–Ritz procedure with existing principles and our new optimization principles, and, as consequences of these principles, we extend various current optimization approaches—steepest descent/ascent and nonlinear conjugate gradient type methods for the Hermitian eigenvalue problem—to calculate a few extreme eigenvalues (of both positive and negative type). A detailed convergence analysis is given for the steepest descent/ascent methods. The analysis reveals the intrinsic quantities that control convergence rates and consequently yields ways of constructing effective preconditioners. Numerical examples are presented to demonstrate the proposed theory and algorithms.