Pair Of Eigenvalues Research Articles

Exactly solvable Stuart-Landau models in arbitrary dimensions.

We use Clifford's geometric algebra to extend the Stuart-Landau system to dimensions D>2 and give an exact solution of the oscillator equationsin the general case. At the supercritical Hopf bifurcation marked by a transition from stable fixed-point dynamics to oscillatory motion, the Jacobian matrix evaluated at the fixed point has N=⌊D/2⌋ pairs of complex conjugate eigenvalues which cross the imaginary axis simultaneously. For odd D there is an additional purely real eigenvalue that does the same. Oscillatory dynamics is asymptotically confined to a hypersphere S^{D-1} and is characterised by extreme multistability, namely the coexistence of an infinite number of limiting orbits, each of which has the geometry of a torus T^{N} on which the motion is either periodic or quasiperiodic. We also comment on similar Clifford extensions of other limit cycle oscillator systems and their generalizations.

Medium induced mixing, spatial modulations, and critical modes in QCD

The mixing between the chiral condensate and the density in hot and dense quantum chromodynamics (QCD) matter is familiar. We show that the mixing relevant for the ground state is considerably more extensive, and in particular also involves gluonic degrees of freedom. As a result, the Hessian of the QCD effective action is non-Hermitian, but retains a symmetry under combined charge and complex conjugation. This can lead to complex-conjugate pairs of eigenvalues of this Hessian, signaling regimes with spatially modulated correlations. Furthermore, based on the analytic structure of the quark determinant at a chiral critical point, we demonstrate that the corresponding massless critical mode is composed of the chiral condensate, the density and the Polyakov loops. Due to an avoided crossing, the critical mode turns out to be disconnected from the chiral condensate in vacuum. We present general arguments for all these features and illustrate them through explicit model calculations. Published by the American Physical Society 2024

The Extended Diamond Difference - Constant Nodal Method with Decoupled Cell Iteration Scheme in Two-Dimensional Discrete Ordinate Transport Problems

Analytical solutions for the Boltzmann transport equation are typically limited to specific scenarios, prompting the research transport community to extensively explore numerical approaches with special attention to the discrete ordinates (SN) formulation. The spectral nodal methods are widely recognized for their capability to preserve the analytical solution of the SN mathematical model. However, these formulations exhibit notable complexity with high-order polynomial approximations for the transverse leakage terms, motivating researchers to pursue simpler hybrid approaches by partially preserving the analytical solution. Both spectral nodal and hybrid spectral nodal methods introduce nonstandard auxiliary equations not amenable to standard iterative schemes due to a highly coupling in angular directions. Following these ideas, we present a finite element – spectral nodal hybrid method referred to as Extended Diamond Difference-Constant Nodal (EDD-CN) method for two-dimensional fixed-source discrete ordinate transport problems with constant nodal approximation for the leakage terms. The proposed method seeks to preserve eigenfunctions associated with the dominant pair of eigenvalues in the transverse integrated SN equations. In addition, we propose a novel iterative scheme called Decoupled Cell Iteration (DCI) scheme which decouples the angular directions in the nonstandard auxiliary equations and follows conventional iteration transport sweeps for solving SN problems.

Asymptotics of Solutions of Autonomous Singularly Perturbed Equations when the Stability of the Equilibrium Position Changes at Several Points

Considers a system consisting of 2n equations with fast variables and one equation with a slow variable. The first approximation matrix of the system of fast variables has 2n pairwise complex conjugate eigenvalues. The stability of the equilibrium position, a system of fast variables, is influenced by all eigenvalues. Early studies considered cases where the stability of the equilibrium position is affected by only one pair of complex conjugate eigenvalues. The problem of delaying the solution of a system of fast variables near an unstable equilibrium position has been solved. The delay time for the decision is determined.

On Spectral Radius and Energy of a Graph with Self-Loops

The spectral radius of a square matrix is the maximum among absolute values of its eigenvalues. Suppose a square matrix is nonnegative; then, by Perron–Frobenius theory, it will be one among its eigenvalues. In this paper, Perron–Frobenius theory for adjacency matrix of graph with self-loops AGS will be explored. Specifically, it discusses the nontrivial existence of Perron–Frobenius eigenvalue and eigenvector pair in the matrix AGS−σnI, where σ denotes the number of self-loops. Also, Koolen–Moulton type bound for the energy of graph GS is explored. In addition, the existence of a graph with self-loops for every odd energy is proved.

Computation of Normal Form and Unfolding of Codimension-3 Zero-Hopf–Hopf Bifurcation

The computation of the normal form as well as its unfolding is a key step to understand the topological structure of a bifurcation. Though a lot of results have been obtained, it still remains unsolved for higher co-dimensional bifurcations. The main purpose of this paper is devoted to the computation of a codimension-3 zero-Hopf–Hopf bifurcation, at which a zero as well as two pairs of pure imaginary eigenvalues can be found from the matrix evaluated at the equilibrium point. Different distributions of eigenvalues are considered, which may behave in a non-semisimple form for 1:1 internal resonance. Based on the combination of center manifold and normal form theory, all the coefficients of normal forms and nonlinear transformations are derived explicitly in terms of parameters of the original vector field, which are obtained via a recursive procedure. Accordingly, a user friendly computer program using a symbolic computation language Maple is developed to compute the coefficients up to an arbitrary order for a specific vector field with zero-Hopf–Hopf bifurcation. Furthermore, universal unfolding parameters are derived in terms of the perturbation of physical parameters, which can be employed to investigate the local behaviors in the neighborhood of the bifurcation point. Here, we emphasize that though different norm forms based on different choices may exist, their topological structures are the same, corresponding to qualitatively equivalent dynamics.

Gaussian-process-regression-based method for the localization of exceptional points in complex resonance spectra

Resonances in open quantum systems depending on at least two controllable parameters can show the phenomenon of exceptional points (EPs), where not only the eigenvalues but also the eigenvectors of two or more resonances coalesce. Their exact localization in the parameter space is challenging, in particular in systems, where the computation of the quantum spectra and resonances is numerically very expensive. We introduce an efficient machine learning algorithm to find EPs based on Gaussian process regression (GPR). The GPR-model is trained with an initial set of eigenvalue pairs belonging to an EP and used for a first estimation of the EP position via a numerically cheap root search. The estimate is then improved iteratively by adding selected exact eigenvalue pairs as training points to the GPR-model. The GPR-based method is developed and tested on a simple low-dimensional matrix model and then applied to a challenging real physical system, viz., the localization of EPs in the resonance spectra of excitons in cuprous oxide in external electric and magnetic fields. The precise computation of EPs, by taking into account the complete valence band structure and central-cell corrections of the crystal, can be the basis for the experimental observation of EPs in this system.

Equilibria and bifurcations in contact dynamics

We provide a systematic study of equilibria of contact vector fields and the bifurcations that occur generically in 1-parameter families, and express the conclusions in terms of the Hamiltonian functions that generate the vector fields. Equilibria occur at points where the zero-level set of the Hamiltonian function is either singular or is tangent to the contact structure. The eigenvalues at an equilibrium have an interesting structure: there is always one particular real eigenvalue of any equilibrium, related to the contact structure, that we call the principal coefficient, while the other eigenvalues arise in quadruplets, similar to the symplectic case except they are translated by a real number equal to half the principal coefficient. There are two types of codimension 1 equilibria, named Type I, arising where the zero-set of the Hamiltonian is singular, and Type II where it is not, but there is a degeneracy related again to the principal coefficient and the contact of the zero-level set of the Hamiltonian with the contact structure. Both give rise generically to saddle-node bifurcations. Some special features include: (i) for Type II singularities, Hopf bifurcations cannot occur in dimension 3, but they may in dimension 5 or more; (ii) for Type I singularities, a fold-Hopf bifurcation can occur with codimension 1 in any dimension, and (iii) again for Type I, and in dimension at least 5, a fold-multi-Hopf bifurcation (where several pairs of eigenvalues pass through the imaginary axis simultaneously together with one through the origin) may also occur with codimension 1.

MULTIPERIODIC OSCILLATIONS IN SECOND-ORDER LYAPUNOV'S SYSTEMS WITH A DIFFERENTIATION OPERATOR IN THE DIRECTION OF A QUASI-PERIODIC VECTOR FIELD

Abstract. The work is devoted to studying oscillatory solutions of Hamiltonian systems of two equations with a differentiation operator in the directions of a vector field with a quasi-periodically varying velocity, in which the frequencies satisfy the condition of strong incommensurability. The system is quasi–linear, has a pair of conjugate pure imaginary eigenvalues, admits the first integral, and is real analytical in time variables in the neighbourhood of the real axis of the complex plane, and space variables on the real plane. It can be called the Lyapunov system. Each equation of a vector field is represented by an independent quasi-periodic scalar function. According to one study by J. Moser the integrals of such equations, under the condition of strong incommensurability, are described by the sum of linear and quasi-periodic functions. To investigate the existence of oscillatory solutions in the system under consideration, Lyapunov's method is extended to this system, considering the result of J. Moser. The differentiation operator with a quasi-periodic property is considered for the first time, and periodic or constant coefficients are often considered in connection with various issues of the theory of multi-frequency oscillations. Thus, the result of the conducted research is an analogue of the result on Lyapunov's periodic systems in the case of quasi-periodic systems with a special differentiation operator.

Eigenvalues bifurcating from the continuum in two-dimensional potentials generating non-Hermitian gauge fields

It has been recently shown that complex two-dimensional (2D) potentials Vɛ(x,y)=V(y+iɛη(x)) can be used to emulate non-Hermitian matrix gauge fields in optical waveguides. Here x and y are the transverse coordinates, V(y) and η(x) are real functions, ɛ>0 is a small parameter, and i is the imaginary unit. The real potential V(y) is required to have at least two discrete eigenvalues in the corresponding 1D Schrödinger operator. When both transverse directions are taken into account, these eigenvalues become thresholds embedded in the continuous spectrum of the 2D operator. Small nonzero ɛ corresponds to a non-Hermitian perturbation which can result in a bifurcation of each threshold into an eigenvalue. Accurate analysis of these eigenvalues is important for understanding the behavior and stability of optical waves propagating in the artificial non-Hermitian gauge potential. Bifurcations of complex eigenvalues out of the continuum is the main object of the present study. Using recent mathematical results from the rigorous analysis of elliptic operators, we obtain simple asymptotic expansions in ɛ that describe the behavior of bifurcating eigenvalues. The lowest threshold can bifurcate into a single eigenvalue, while every other threshold can bifurcate into a pair of complex eigenvalues. These bifurcations can be controlled by the Fourier transform of function η(x) evaluated at certain isolated points of the reciprocal space. When the bifurcation does not occur, the continuous spectrum of 2D operator contains a quasi-bound-state which is characterized by a strongly localized central peak coupled to small-amplitude but nondecaying tails. The analysis is applied to the case examples of parabolic and double-well potentials V(y). In the latter case, the bifurcation of complex eigenvalues can be dampened if the two wells are widely separated.

Theoretical and experimental characterization of non-Markovian anti-parity-time systems

Open systems with anti-parity-time ({{{{{{{mathcal{APT}}}}}}}}) or {{{{{{{mathcal{PT}}}}}}}} symmetry exhibit a rich phenomenology absent in their Hermitian counterparts. To date all model systems and their diverse realizations across classical and quantum platforms have been local in time, i.e., Markovian. Here we propose a non-Markovian system with anti-{{{{{{{mathcal{PT}}}}}}}}-symmetry where a single time-delay encodes the retention of memory, and experimentally demonstrate its consequences with two time-delay coupled semiconductor lasers. A transcendental characteristic equation with infinitely many eigenvalue pairs sets our model apart. We show that a sequence of amplifying-to-decaying dominant mode transitions is induced by the time delay in our minimal model. The signatures of these transitions quantitatively match results obtained from four, coupled, nonlinear rate equations for laser dynamics, and are experimentally observed as constant-width sideband oscillations in the laser intensity profiles. Our work introduces a paradigmatic non-Hermitian system with memory, paves the way for its realization in classical systems, and may apply to time-delayed feedback-control for quantum systems.

4 × 4 Irreducible sign pattern matrices that require four distinct eigenvalues

A sign pattern matrix is a matrix whose entries are from the set {+,−,0}. For a sign pattern matrix A, the qualitative class of A, denoted Q(A), is the set of all real matrices whose entries have signs given by the corresponding entries of A. An n×n sign pattern matrix A requires all distinct eigenvalues if every real matrix in Q(A) has n distinct eigenvalues. Li and Harris (2002) [13] characterized the 2×2 and 3×3 irreducible sign pattern matrices that require all distinct eigenvalues, and established some useful general results on n×n sign patterns that require all distinct eigenvalues. In this paper, we characterize 4×4 irreducible sign patterns that require four distinct eigenvalues. This is done by characterizing 4×4 irreducible sign patterns that require four distinct real eigenvalues, that require four distinct nonreal real eigenvalues, or that require two distinct real eigenvalues and a pair of conjugate nonreal eigenvalues. The last case turns out to be much more involved. Some interesting open problems are presented. Three important tools that are used in the paper are the following: the discriminant of a polynomial; the fact that if a square sign pattern matrix A requires all distinct eigenvalues then A requires a fixed number of real eigenvalues; and the known result that if A is a “k-cycle” sign pattern then for each B∈Q(A), the k nonzero eigenvalues of B are evenly distributed on a circle in the complex plane centered at the origin.

Synchronization in networked systems with large parameter heterogeneity

Systems that synchronize in nature are intrinsically different from one another, with possibly large differences from system to system. While a vast part of the literature has investigated the emergence of network synchronization for the case of small parametric mismatches, we consider the general case that parameter mismatches may be large. We present a unified stability analysis that predicts why the range of stability of the synchronous solution either increases or decreases with parameter heterogeneity for a given network. We introduce a parametric approach, based on the definition of a curvature contribution function, which allows us to estimate the effect of mismatches on the stability of the synchronous solution in terms of contributions of pairs of eigenvalues of the Laplacian. For cases in which synchronization occurs in a bounded interval of a parameter, we study the effects of parameter heterogeneity on both transitions (asynchronous to synchronous and synchronous to asynchronous.).

Benjamin–Feir Instability of Stokes Waves in Finite Depth

Whitham and Benjamin predicted in 1967 that small-amplitude periodic traveling Stokes waves of the 2d-gravity water waves equations are linearly unstable with respect to long-wave perturbations, if the depth {mathtt h} is larger than a critical threshold texttt{h}_{scriptscriptstyle {textsc {WB}}}approx 1.363 . In this paper, we completely describe, for any finite value of mathtt h >0 , the four eigenvalues close to zero of the linearized equations at the Stokes wave, as the Floquet exponent mu is turned on. We prove, in particular, the existence of a unique depth texttt{h}_{scriptscriptstyle {textsc {WB}}}, which coincides with the one predicted by Whitham and Benjamin, such that, for any 0< mathtt h < texttt{h}_{scriptscriptstyle {textsc {WB}}}, the eigenvalues close to zero are purely imaginary and, for any mathtt h > texttt{h}_{scriptscriptstyle {textsc {WB}}}, a pair of non-purely imaginary eigenvalues depicts a closed figure “8”, parameterized by the Floquet exponent. As {mathtt h} rightarrow texttt{h}_{scriptscriptstyle {textsc {WB}}}^{, +} the “8” collapses to the origin of the complex plane. The complete bifurcation diagram of the spectrum is not deduced as in deep water, since the limits texttt{h}rightarrow +infty (deep water) and mu rightarrow 0 (long waves) do not commute. In finite depth, the four eigenvalues have all the same size mathcal {O}(mu ), unlike in deep water, and the analysis of their splitting is much more delicate, requiring, as a new ingredient, a non-perturbative step of block-diagonalization. Along the whole proof, the explicit dependence of the matrix entries with respect to the depth texttt{h} is carefully tracked.

Bound states in Bose-Einstein condensates with radially-periodic spin-orbit coupling

We consider Bose-Einstein condensate (BEC) subject to the action of spin-orbit-coupling (SOC) periodically modulated in the radial direction. In contrast to the commonly known principle that periodic potentials do not create bound states, the binary BEC maintains multiple localized modes in the linear limit, with their chemical potential falling into spectral gaps of the (numerically found) radial band structure induced by the spatial modulation of the SOC. In the presence of the repulsive nonlinearity, the SOC modulation supports fundamental gap solitons of the semi-vortex types, as well as higher-order vortex gap solitons. The localization degree and stability of the gap solitons strongly depend on the location of their chemical potential in the gap. Stability intervals for vortex gap solitons in a broad range of the intrinsic vorticity, from −2 to 3, are identified. Thus, the analysis reveals the previously unexplored mechanism of linear and nonlinear localization provided by the spatially periodic modulation of SOC, which may be extended to other settings, such as those for optical beams and polaritons. Unlike the commonly known quartets of eigenvalues for small perturbations, in the present system the instability is accounted for by shifted complex eigenvalue pairs.

Constructible reality condition of pseudo entropy via pseudo-Hermiticity

As a generalization of entanglement entropy, pseudo entropy is not always real. The real-valued pseudo entropy has promising applications in holography and quantum phase transition. We apply the notion of pseudo-Hermiticity to formulate the reality condition of pseudo entropy. We find the general form of the transition matrix for which the eigenvalues of the reduced transition matrix possess real or complex pairs of eigenvalues. Further, we find a class of transition matrices for which the pseudo (Rényi) entropies are non-negative. Some known examples which give real pseudo entropy in quantum field theories can be explained in our framework. Our results offer a novel method to generate the transition matrix with real pseudo entropy. Finally, we show the reality condition for pseudo entropy is related to the Tomita-Takesaki modular theory for quantum field theory.

Generalized Mercier stability criterion for stellarators

The Mercier criterion is a well-known stability criterion for tokamaks. It is derived from a 2 × 2 matrix eigenvalue problem arising from the expansion of resonant solutions about a singular surface where m−nq=0, with m and n being the poloidal and toroidal mode numbers, respectively, and q being the safety factor. The stability criterion is that the eigenvalues must be real, otherwise, the solution oscillates, violating the Newcomb crossing criterion. Because of the non-axisymmetry of stellarators, different toroidal as well as poloidal harmonics couple to each other. It follows that each singular surface can have multiple resonant harmonics, with multiplicity M≥1. The corresponding matrix eigenvalue problem involves a 2M×2M matrix, resulting in M pairs of positive and negative eigenvalues. The generalized stability criterion is that all eigenvalues must be real. While the original Mercier criterion can be expressed in terms of quadratures of equilibrium quantities over the singular surface, which can be evaluated anywhere, the generalized Mercier criterion can only be evaluated on rational q surfaces with a given set of resonant harmonics.

The Hamiltonian formalism of the inverse problem

A big-isotropic structure is a generalization of the notion of Dirac structure, due to Vaisman. We discuss the inverse problem of deciding if a vector field is Hamiltonian having a big-isotropic structure as underlying geometry. In [1] we have considered this question for the special case of replicator equations. Here we generalize that approach to any vector field that can be written in the form X(u)=(Bη)(u), where B=B(u) is a matrix and η=η(u) is a vector. For a linear system we show that, if the representing matrix of the system has at least one pair of positive-negative non-zero eigenvalues, or a zero eigenvalue with at least one 3-dimensional Jordan block associated to it, then the linear system has a Hamiltonian description with respect to a big-isotropic structure. As a byproduct, we find a class of linear systems with zero eigenvalue that are Hamiltonian with respect to a big-isotropic structure but not a Dirac structure. Moreover, we prove that every linear Hamiltonian system, having a big-isotropic structure as underlying geometry, is completely integrable in the sense of Zung [12]. For linear systems, in the both Hamiltonian formulation and complete integrability cases, explicit descriptions of the geometric structures, the Hamiltonian functions, the first integrals and the commuting flows are provided.

The 0:1 resonance bifurcation associated with the supercritical Hamiltonian pitchfork bifurcation

We consider the non-semisimple 0:1 resonance (i.e., the unperturbed equilibrium has two purely imaginary eigenvalues ( and ) and a non-semisimple double-zero one) Hamiltonian bifurcation with one distinguished parameter, which corresponds to the supercritical Hamiltonian pitchfork bifurcation. Based on BCKV singularity theory established by Broer et al. (Z. Angew. Math. Phys. 44:389-432, 1993), this bifurcation essentially triggered by the reversible universal unfolding with respect to BCKV-restricted morphisms of the planar non-semisimple singularity (the is regarded as distinguished parameter with respect to the external parameter λ). We first give the plane bifurcation diagram of the integrable Hamiltonian on each level of integral in detail, which is related to the usual supercritical Hamiltonian pitchfork bifurcation. Then, we use the -symmetry generated by the additional pair of imaginary eigenvalues to reconstruct the above plane bifurcation phenomenon caused by the zero eigenvalue pair into the case with two degrees of freedom. Finally, we prove the persistence of typical bifurcation scenarios (e.g., 2-dimensional invariant tori and the symmetric homoclinic orbit) under the small Hamiltonian perturbations, as proposed by Broer et al. (Z Angew Math Phys 44: 389-432, 1993). An example system (the coupled Duffing oscillator) with strong linear coupling and weak local nonlinearity is given for this bifurcation.

Inseparable Gershgorin discs and the existence of conjugate complex eigenvalues of real matrices

We investigate the converse of the known fact that if the Gershgorin discs of a real n-by-n matrix may be separated by positive diagonal similarity, then the eigenvalues are real. In the 2-by-2 case, with appropriate signs for the off-diagonal entries, we find that the converse is correct, which raises several questions. First, in the 3-by-3 case, the converse is not generally correct, but, empirically, it is frequently true. Then, in the n-by-n case, n ≥ 3 , we find that if all the 2-by-2 principal submatrices have inseparable discs (‘strongly inseparable discs’), the full matrix must have a nontrivial pair of conjugate complex eigenvalues (i.e. cannot have all real eigenvalues). This hypothesis cannot generally be weakened.