Search

We consider large non-Hermitian random matrices $X$ with complex, independent, identically distributed centred entries and show that the linear statistics of their eigenvalues are asymptotically Gaussian for test functions having $2+\epsilon$ derivatives. Previously this result was known only for a few special cases; either the test functions were required to be analytic [Rider, Silverstein 2006], or the distribution of the matrix elements needed to be Gaussian [Rider, Virag 2007], or at least match the Gaussian up to the first four moments [Tao, Vu 2016; Kopel 2015]. We find the exact dependence of the limiting variance on the fourth cumulant that was not known before. The proof relies on two novel ingredients: (i) a local law for a product of two resolvents of the Hermitisation of $X$ with different spectral parameters and (ii) a coupling of several weakly dependent Dyson Brownian Motions. These methods are also the key inputs for our analogous results on the linear eigenvalue statistics of real matrices $X$ that are presented in the companion paper [Cipolloni, Erdős, Schroder 2019].

In the block version of the classical two-sided Jacobi method for the Hermitian eigenvalue problem, the off-diagonal elements of iterated matrix A ( k ) converge to zero. However, this fact alone does not necessarily guarantee that A ( k ) converges to a fixed diagonal matrix. The same is true for the matrix of accumulated unitary transformations Q ( k ) . We prove that under certain assumptions A ( k ) indeed converges to a fixed diagonal matrix, whose diagonal elements are the eigenvalues of the input matrix A. Next it is shown that for a simple eigenvalue the corresponding column of Q ( k ) converges to the corresponding eigenvector. For a multiple eigenvalue or a cluster of eigenvalues, we prove that the orthogonal projectors constructed from the corresponding columns of Q ( k ) converge to the orthogonal projector onto the eigenspace corresponding to those eigenvalues. Moreover, the appropriate convergence bounds are obtained for all discussed cases. Convergence results are also valid for the parallel block-Jacobi method with dynamic ordering. The developed theory is illustrated by numerical example.

In this article, we propose the Fast Algorithms for Maxwell’s Equations (FAME) package for solving Maxwell’s equations for modeling three-dimensional photonic crystals. FAME combines the null-space free method with fast Fourier transform (FFT)-based matrix-vector multiplications to solve the generalized eigenvalue problems (GEPs) arising from Yee’s discretization. The GEPs are transformed into a null-space free standard eigenvalue problem with a Hermitian positive-definite coefficient matrix. The computation times for FFT-based matrix-vector multiplications with matrices of dimension 7 million are only 0.33 and 3.6 × 10 − 3 seconds using MATLAB with an Intel Xeon CPU and CUDA C++ programming with a single NVIDIA Tesla P100 GPU, respectively. Such multiplications significantly reduce the computational costs of the conjugate gradient method for solving linear systems. We successfully use FAME on a single P100 GPU to solve a set of GEPs with matrices of dimension more than 19 million, in 127 to 191 seconds per problem. These results demonstrate the potential of our proposed package to enable large-scale numerical simulations for novel physical discoveries and engineering applications of photonic crystals.

In this paper, we study T-eigenvalues of third-order tensors. Definitions of the T-eigenvalues and Hermitian tensors are proposed. We present a commutative tensor family. We prove some T-eigenvalue inequalities for Hermitian tensors, including extensions of Weyl's theorem and Cauchy's interlacing theorem from the matrix case to the tensor case. Finally, we introduce the stability of the T-eigenvalues and study the Lyapunov equation for tensors.

For Hermitian and non-Hermitian Hamiltonian matrices H, we present the Schr¨odinger equation for qudit (spin-j system, N-level atom) with the state vector |ψ〉 in a new form of the linear eigenvalue equation for the matrix ℋ = (H ⊗ 1N) and the probability eigenvector |p〉 identified with quantum states in the probability representation of quantum mechanics. We discuss the possibility to experimentally detect the difference between the system states described by the solutions, corresponding to the Schrodinger equation with Hermitian and non-Hermitian Hamiltonians, by measuring the probabilities of artificial spin-1/2 projections m = ± 1/2, sets of which are identified with qudit states. We show that different symmetries of systems, including 𝒫𝒯 -symmetry and broken 𝒫𝒯 -symmetry, are determined by a set of N complex eigenvalues of the Hamiltonian matrix H.

In this work, we study the adiabatic phases and shortcut for a pseudo-Hermitian interacting spin-1/2 composite system in which one system is driven by a non-Hermitian magnetic field. The biorthonormal eigenbasis and the condition for real eigenvalues have been derived and discussed. By studying the Berry phases and the transitionless quantum driving (TQD) Hamiltonian, the roles of correlation between the subsystems and the non-Hermicity are discussed. The non-Hermitian dynamics and the Berry phases of the subsystem, as well as their influences on the TQD Hamiltonian are studied in both no correlation and large correlation cases. These results mean that we can modify the correlation between the subsystems to control the adiabatic shortcuts of the subsystems. These schemes can provide a feasible way for studying the adiabatic shortcuts of the composite system and the non-Hermitian system.

We generalize von Neumann's well-known trace inequality, as well as related eigenvalue inequalities for hermitian matrices, to Schatten-class operators between complex Hilbert spaces of infinite dimension. To this end, we exploit some recent results on the $C$-numerical range of Schatten-class operators. For the readers' convenience, we sketched the proof of these results in the Appendix., 16 pages

Elastography is a ill-posed inverse problem that aims at recovering the Lame and density of the domain of interest from finite number of observations, we consider as model the Helmoltz equation. We present an implementation for solving the Helmholtz inverse problem in two dimensions via an optimization problem based on Bayesian approach. In addition, the accuracy of the method is also investigated with respect to the amount of information taken from the generalized Hermitian Eigenvalue problem and by comparing the maximum a posterior estimate to the true parameter distribution in simulated experiments.

The Ritz upper bound to eigenvalues of Hermitian operators is essential for many applications in science. It is a staple of quantum chemistry and physics computations. The lower bound devised by Temple in 1928 [G. Temple, Proc. R. Soc. A Math. Phys. Eng. Sci. 119, 276–293 (1928)] is not, since it converges too slowly. The need for a good lower-bound theorem and algorithm cannot be overstated, since an upper bound alone is not sufficient for determining differences between eigenvalues such as tunneling splittings and spectral features. In this paper, after 90 y, we derive a generalization and improvement of Temple’s lower bound. Numerical examples based on implementation of the Lanczos tridiagonalization are provided for nontrivial lattice model Hamiltonians, exemplifying convergence over a range of 13 orders of magnitude. This lower bound is typically at least one order of magnitude better than Temple’s result. Its rate of convergence is comparable to that of the Ritz upper bound. It is not limited to ground states. These results complement Ritz’s upper bound and may turn the computation of lower bounds into a staple of eigenvalue and spectral problems in physics and chemistry.

We study certain pairs of subspaces V and W of C^n we call supports that consist of eigenspaces of the eigenvalues ±‖M‖ of a minimal hermitian matrix M(‖M‖ ≤‖M+D‖ for all real diagonals D). For any pair of orthogonal subspaces we define a non negative invariant δ called the adequacy to measure how close they are to form a support and to detect one. This function δ is the minimum of another map F defined in a product of spheres of hermitian matrices. We study the gradient, Hessian and critical points of F in order to approximate δ. These results allow us to prove that the set of supports has interior points in the space of flag manifolds. Fil: Mendoza, Alberto. Universidad Simón Bolívar; Venezuela Fil: Recht, Lázaro. Universidad Simón Bolívar; Venezuela. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Saavedra 15. Instituto Argentino de Matemática Alberto Calderón; Argentina Fil: Varela, Alejandro. Universidad Nacional de General Sarmiento. Instituto de Ciencias; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Saavedra 15. Instituto Argentino de Matemática Alberto Calderón; Argentina

This paper considers the problem of estimating the number of sinusoidal signals in additive sub-Gaussian noise. We develop a novel non-asymptotic analysis of the eigenvalues of a Hermitian Toeplitz data matrix with finite noisy measurements. With this analysis, we propose a new threshold-based estimation algorithm which is guaranteed to correctly recover the true source number with high probability under certain conditions. The analysis is applicable to any independent sub-Gaussian additive noise. In particular, exact knowledge of the likelihood function is not needed and the noise terms may not be identically distributed. The theoretical claims are demonstrated by extensive numerical experiments.

We demonstrate that crystal defects can act as a probe of intrinsic non-Hermitian topology. In particular, in point-gapped systems with periodic boundary conditions, a pair of dislocations may induce a non-Hermitian skin effect, where an extensive number of Hamiltonian eigenstates localize at only one of the two dislocations. An example of such a phase are two-dimensional systems exhibiting weak non-Hermitian topology, which are adiabatically related to a decoupled stack of Hatano-Nelson chains. Moreover, we show that strong two-dimensional point-gap topology may also result in a dislocation response, even when there is no skin effect present with open boundary conditions. For both cases, we directly relate their bulk topology to a stable dislocation non-Hermitian skin effect. Finally, and in stark contrast to the Hermitian case, we find that gapless non-Hermitian systems hosting bulk exceptional points also give rise to a well-localized dislocation response., 6 pages, 4 figures, supplement included, accepted manuscript

The geometry of Hamiltonian's eigenstates is encoded in the quantum geometric tensor (QGT). It contains both the Berry curvature, central to the description of topological matter and the quantum metric. So far the full QGT has been measured only in Hermitian systems, where the role of the quantum metric is mostly shown to determine corrections to physical effects. On the contrary, in non-Hermitian systems, and in particular near exceptional points, the quantum metric is expected to diverge and to often play a dominant role, for example on the enhanced sensing and on wave packet dynamics. In this work, we report the first experimental measurement of the quantum metric in a non-Hermitian system. The specific platform under study is an organic microcavity with exciton-polariton eigenstates, which demonstrate exceptional points. We measure the quantum metric's divergence and we determine the scaling exponent $n=-1.01\pm0.08$, which is in agreement with theoretical predictions for the second-order exceptional points.

Time evolution operators of a strongly ionizing medium are calculated by a time-dependent unitary transformation (TDUT) method. The TDUT method has been employed in a quantum mechanical system composed of discrete states. This method is especially helpful for solving molecular rotational dynamics in quasi-adiabatic regimes because the strict unitary nature of the propagation operator allows us to set the temporal step size to large, a tight limitation on the temporal step size (δt&lt, &lt, 1) can be circumvented by the strict unitary nature. On the other hand, in a strongly ionizing system where the Hamiltonian is not Hermitian, the same approach cannot be directly applied because it is demanding to define a set of field-dressed eigenstates. In this study, the TDUT method was applied to the ionizing regime using the Kramers–Henneberger frame, in which the strong-field-dressed discrete eigenstates are given by the field-free discrete eigenstates in a moving frame. Although the present work verifies the method for a one-dimensional atom as a prototype, the method can be applied to three-dimensional atoms, and molecules exposed to strong laser fields.

Efficiently processing basic linear algebra subroutines is of great importance for a wide range of computational problems. In this paper, we consider techniques to implement matrix functions on a quantum computer. We embed given matrices into 3 times larger Hermitian matrices and assume as input a given set of unitary operators generated by the embedding matrices. With the matrix embedding formula, we give Trotter-based quantum subroutines for elementary matrix operations include addition, multiplication, Kronecker sum, tensor product, Hadamard product, and arbitrary real eigenvalue single-matrix functions. We then discuss the composed matrix functions in terms of the estimation of scalar quantities such as inner products, traces, determinants, and Schatten p-norms with bounded errors. We thus provide a framework for compiling instructions for linear algebraic computations into gate sequences on actual quantum computers. The framework for calculating the matrix functions is more efficient than the best classical counterpart for a set of matrices.

We introduce a powerful analytic method to study the statistics of the number $\mathcal{N}_{\textbf{A}}(\gamma)$ of eigenvalues inside any contour $\gamma \in \mathbb{C}$ for infinitely large non-Hermitian random matrices ${\textbf A}$. Our generic approach can be applied to different random matrix ensembles, even when the analytic expression for the joint distribution of eigenvalues is not known. We illustrate the method on the adjacency matrices of weighted random graphs with asymmetric couplings, for which standard random-matrix tools are inapplicable. The main outcome is an effective theory that determines the cumulant generating function of $\mathcal{N}_{\textbf{A}}$ via a path integral along $\gamma$, with the path probability distribution following from the solution of a self-consistent equation. We derive the expressions for the mean and the variance of $\mathcal{N}_{\textbf{A}}$ as well as for the rate function governing rare fluctuations of ${\mathcal{N}}_{\textbf{A}}{(\gamma)}$. All theoretical results are compared with direct diagonalization of finite random matrices, exhibiting an excellent agreement., Comment: 6 pages, 2 figures. SI as ancillary file

A common wisdom in quantum mechanics is that the Hamiltonian has to be Hermitian in order to ensure a real eigenvalue spectrum. Yet, parity–time (PT)-symmetric Hamiltonians are sufficient for real eigenvalues and therefore constitute a complex extension of quantum mechanics beyond the constraints of Hermiticity. However, as only single-particle or classical wave physics has been exploited so far, an experimental demonstration of the true quantum nature of PT symmetry has been elusive. In our work, we demonstrate two-particle quantum interference in a PT-symmetric system. We employ integrated photonic waveguides to reveal that the quantum dynamics of indistinguishable photons shows strongly counterintuitive features. To substantiate our experimental data, we analytically solve the quantum master equation using Lie algebra methods. The ideas and results presented here pave the way for non-local PT-symmetric quantum mechanics as a novel building block for future quantum devices. Parity–time symmetry in second quantization is demonstrated in an integrated non-Hermitian coupled waveguide structure. A counterintuitive shift of the position of the Hong–Ou–Mandel dip is observed in integrated lossy waveguide structures.

We prove a sharp Hardy-type inequality for the Dirac operator. We exploit this inequality to obtain spectral properties of the Dirac operator perturbed with Hermitian matrix-valued potentials $${\mathbf {V}}$$ of Coulomb type: we characterise its eigenvalues in terms of the Birman–Schwinger principle and we bound its discrete spectrum from below, showing that the ground-state energy is reached if and only if $${\mathbf {V}}$$ verifies some rigidity conditions. In the particular case of an electrostatic potential, these imply that $${\mathbf {V}}$$ is the Coulomb potential.

The standard Yee's scheme for the Maxwell eigenvalue problem places the discrete electric field variable at the midpoints of the edges of the grid cells. It performs well when the permittivity is a scalar field. However, when the permittivity is a Hermitian full tensor field it would generate un-physical complex eigenvalues or frequencies. In this paper, we propose a finite element method which can be interpreted as a modified Yee's scheme to overcome this difficulty. This interpretation enables us to create a fast FFT eigensolver that can compute very effectively the band structure of the anisotropic photonic crystal with SC and FCC lattices. Furthermore, we overcome the usual large null space associated with the Maxwell eigenvalue problem by deriving a null-space free discrete eigenvalue problem which involves a crucial Hermitian positive definite linear system to be solved in each of the iteration steps. It is demonstrated that the CG method without preconditioning converges in 37 iterations even when the dimension of a matrix is as large as 5 , 184 , 000 .

We compute the low lying eigenvalues of the Hermitian Dirac operator in lattice QCD with N f = 2 + 1 + 1 twisted mass fermions . We discuss whether these eigenvalues are in the ϵ -regime or the p -regime of Wilson chiral perturbation theory ( χ PT) for twisted mass fermions. Reaching the deep ϵ -regime is practically unfeasible with presently typical simulation parameters, but still the few lowest eigenvalues of the employed ensemble evince some characteristic ϵ -regime features. With this conclusion in mind, we develop a fitting strategy to extract two low energy constants from analytical ϵ -regime predictions at a fixed index. Thus, we obtain results for the chiral condensate and the low energy constant W 8 . We also discuss how to improve both the theoretical calculation and the lattice computation.

In classical Whitham modulation theory, the transition of the dispersionless Whitham equations from hyperbolic to elliptic is associated with a pair of nonzero purely imaginary eigenvalues coalescing and becoming a complex quartet, suggesting that a Krein signature is operational. However, there is no natural symplectic structure. Instead, we find that the operational signature is the “sign characteristic” of real eigenvalues of Hermitian matrix pencils. Its role in classical Whitham single‐phase theory is elaborated for illustration. However, the main setting where the sign characteristic becomes important is in multiphase modulation. It is shown that a necessary condition for two coalescing characteristics to become unstable (the generalization of the hyperbolic to elliptic transition) is that the characteristics have opposite sign characteristic. For example the theory is applied to multiphase modulation of the two‐phase traveling wave solutions of coupled nonlinear Schrodinger equation.

Let [math] and [math] be two [math] Hermitian matrices. Assume that the eigenvalues [math] of [math] are known, as well as the eigenvalues [math] of [math] . What can be said about the eigenvalues of the sum [math] ? This is Horn’s problem. We revisit this question from a probabilistic viewpoint. The set of Hermitian matrices with spectrum [math] is an orbit [math] for the natural action of the unitary group [math] on the space of [math] Hermitian matrices. Assume that the random Hermitian matrix [math] is uniformly distributed on the orbit [math] and, independently, the random Hermitian matrix [math] is uniformly distributed on [math] . We establish a formula for the joint distribution of the eigenvalues of the sum [math] . The proof involves orbital measures with their Fourier transforms, and Heckman’s measures.

Quantum systems governed by non-Hermitian Hamiltonians with PT symmetry are special in having real energy eigenvalues bounded below and unitary time evolution. We argue that PT symmetry may also be important and present at the level of Hermitian quantum field theories because of the process of renormalisation. In some quantum field theories renormalisation leads to PT -symmetric effective Lagrangians. We show how PT symmetry may allow interpretations that evade ghosts and instabilities present in an interpretation of the theory within a Hermitian framework. From the study of examples PT -symmetric interpretation is naturally built into a path integral formulation of quantum field theory; there is no requirement to calculate explicitly the PT norm that occurs in Hamiltonian quantum theory. We discuss examples where PT -symmetric field theories emerge from Hermitian field theories due to effects of renormalisation. We also consider the effects of renormalisation on field theories that are non-Hermitian but PT -symmetric from the start.

The complex-valued quantum mechanics considers quantum motion on the complex plane instead of on the real axis, and studies the variations of a particle complex position, momentum and energy along a complex trajectory. On the basis of quantum Hamilton-Jacobi formalism in the complex space, we point out that having complex-valued motion is a universal property of quantum systems, because every quantum system is actually accompanied with an intrinsic complex Hamiltonian originating from the equation. It is revealed that the conventional real-valued quantum mechanics is a special case of the complex-valued quantum mechanics in that the eigen-structures of real and complex quantum systems, such as their eigenvalues, eigenfunctions and eigen-trajectories, are invariant under linear complex mapping. In other words, there is indeed no distinction between Hermitian systems, PT-symmetric systems, and non PT-symmetric systems when viewed from a complex domain. Their eigen-structures can be made coincident through linear transformation of complex coordinates.

Parity-time ($\mathcal{PT}$) symmetric systems are classical, gain-loss systems whose dynamics are governed by non-Hermitian Hamiltonians with exceptional-point (EP) degeneracies. The eigenvalues of a $\mathcal{PT}$-symmetric Hamiltonian change from real to complex conjugates at a critical value of gain-loss strength that is called the $\mathcal{PT}$ breaking threshold. Here, we obtain the $\mathcal{PT}$-threshold for a one-dimensional, finite Kitaev chain -- a prototype for a p-wave superconductor -- in the presence of a single pair of gain and loss potentials as a function of the superconducting order parameter, on-site potential, and the distance between the gain and loss sites. In addition to a robust, non-local threshold, we find a rich phase diagram for the threshold that can be qualitatively understood in terms of the band-structure of the Hermitian Kitaev mo del. In particular, for an even chain with zero on-site potential, we find a re-entrant $\mathcal{PT}$-symmetric phase bounded by second-order EP contours. Our numerical results are supplemented by analytical calculations for small system sizes., 8 pages, 6 figures

We revisit the topic of two-state quantum systems using the Clifford Algebra in three dimensions $$Cl_3$$ . In this description, both the quantum states and Hermitian operators are written as elements of $$Cl_3$$ . By writing the quantum states as elements of the minimal left ideals of this algebra, we compute the energy eigenvalues and eigenvectors for the Hamiltonian of an arbitrary two-state system. The geometric interpretation of the Hermitian operators enables us to introduce an algebraic method to diagonalize these operators in $$Cl_3$$ . We then use this approach to revisit the problem of a spin-1/2 particle interacting with an external arbitrary constant magnetic field, obtaining the same results as in the conventional theory. However, Clifford algebra reveals the underlying geometry of these systems, which reduces to the Larmor precession in an arbitrary plane of $$Cl_3$$ .

Bulk-boundary correspondence, connecting the bulk topology and the edge states, is an essential principle of the topological phases. However, the bulk-boundary correspondence is broken down in general non-Hermitian systems. In this paper, we construct one-dimensional non-Hermitian Su-Schrieffer-Heeger model with periodic driving that exhibits non-Hermitian skin effect: all the eigenstates are localized at the boundary of the systems, whether the bulk states or the zero and the $\pi$ modes. To capture the topological properties, the non-Bloch winding numbers are defined by the non-Bloch periodized evolution operators based on the generalized Brillouin zone. Furthermore, the non-Hermitian bulk-boundary correspondence is established: the non-Bloch winding numbers ($W_{0,\pi}$) characterize the edge states with quasienergies $\epsilon=0, \pi$. In our non-Hermitian system, a novel phenomenon can emerge that the robust edge states can appear even when the Floquet bands are topological trivial with zero non-Bloch band invariant, which is defined in terms of the non-Bloch effective Hamiltonian. We also show that the relation between the non-Bloch winding numbers ($W_{0,\pi}$) and the non-Bloch band invariant ($\mathcal{W}$): $\mathcal{W}= W_{0}- W_{\pi}$., Comment: 8 pages, 4 figures

Much has been learned about universal properties of the eigenstate entanglement entropy for one-dimensional lattice models, which is described by a Hermitian Hamiltonian. While very less of it has been understood for non-Hermitian systems. In the present work we study a non-Hermitian, non-interacting model of fermions which is invariant under combined $PT$ transformation. Our models show a phase transition from $PT$ unbroken phase to broken phase as we tune the hermiticity breaking parameter. Entanglement entropy of such systems can be defined in two different ways, depending on whether we consider only right (or equivalently only left) eigenstates or a combination of both left and right eigenstates which form a complete set of bi-orthonormal eigenstates. We demonstrate that the entanglement entropy of the ground state and also of the typical excited states show some unique features in both of these phases of the system. Most strikingly, entanglement entropy obtained taking a combination of both left and right eigenstates shows an exponential divergence with system size at the transition point. While in the $PT$-unbroken phase, the entanglement entropy obtained from only the right (or equivalently left) eigenstates shows identical behavior as of an equivalent Hermitian system which is connected to the non-Hermitian system by a similarity transformation., 9 pages, 7 figures

We consider large non-Hermitian real or complex random matrices X with independent, identically distributed centred entries. We prove that their local eigenvalue statistics near the spectral edge, the unit circle, coincide with those of the Ginibre ensemble, i.e. when the matrix elements of X are Gaussian. This result is the non-Hermitian counterpart of the universality of the Tracy-Widom distribution at the spectral edges of the Wigner ensemble., Probability Theory and Related Fields, 179, ISSN:0178-8051, ISSN:1432-2064

We study quantum tomography from a continuous measurement record obtained by measuring expectation values of a set of Hermitian operators obtained from unitary evolution of an initial observable. For this purpose, we consider the application of a random unitary, diagonal in a fixed basis at each time step and quantify the information gain in tomography using Fisher information of the measurement record and the Shannon entropy associated with the eigenvalues of covariance matrix of the estimation. Surprisingly, very high fidelity of reconstruction is obtained using random unitaries diagonal in a fixed basis even though the measurement record is not informationally complete. We then compare this with the information generated and fidelities obtained by application of a different Haar random unitary at each time step. We give an upper bound on the maximal information that can be obtained in tomography and show that a covariance matrix taken from the Wishart-Laguerre ensemble of random matrices and the associated Marchenko-Pastur distribution saturates this bound. We find that physically, this corresponds to an application of a different Haar random unitary at each time step. We show that repeated application of random diagonal unitaries gives a covariance matrix in tomographic estimation that corresponds to a new ensemble of random matrices. We analytically and numerically estimate eigenvalues of this ensemble and show the information gain to be bounded from below by the Porter-Thomas distribution., 12 pages, 9 figures

We point out that in non-Hermitian systems, Jordan decomposition should replace eigendecomposition so that all ``good approximate eigenstates'' of the system are identified. These states can be resonantly excited. As a concrete example, we study the location and field distribution of zero-energy corner states in a non-Hermitian quadrupole insulator (QI) and split the parameter space into three distinct regimes according to properties of the corner states: near-Hermitian QI, intermediate phase, and trivial insulator. In the newly discovered intermediate phase, the Hamiltonian becomes defective, and the counterintuitive and delicate response of the system to external drives can be perfectly explained by the Jordan decomposition.

Skin effect, experimentally discovered in one dimension, describes the physical phenomenon that on an open chain, an extensive number of eigenstates of a non-Hermitian hamiltonian are localized at the end(s) of the chain. Here in two and higher dimensions, we establish a theorem that the skin effect exists, if and only if periodic-boundary spectrum of the hamiltonian covers a finite area on the complex plane. This theorem establishes the universality of the effect, because the above condition is satisfied in almost every generic non-Hermitian hamiltonian, and, unlike in one dimension, is compatible with all spatial symmetries. We propose two new types of skin effect in two and higher dimensions: the corner-skin effect where all eigenstates are localized at one corner of the system, and the geometry-dependent-skin effect where skin modes disappear for systems of a particular shape, but appear on generic polygons. An immediate corollary of our theorem is that any non-Hermitian system having exceptional points (lines) in two (three) dimensions exhibits skin effect, making this phenomenon accessible to experiments in photonic crystals, Weyl semimetals, and Kondo insulators., Comment: 21 pages, 10 figures

Distant boundaries in linear non-Hermitian lattices can dramatically change energy eigenvalues and corresponding eigenstates in a nonlocal way. This effect is known as non-Hermitian skin effect (NHSE). Combining non-Hermitian skin effect with nonlinear effects can give rise to a host of novel phenomenas, which may be used for nonlinear structure designs. Here we study nonlinear non-Hermitian skin effect and explore nonlocal and substantial effects of edges on stationary nonlinear solutions. We show that fractal and continuum bands arise in a long lattice governed by a nonreciprocal discrete nonlinear Schrodinger equation. We show that stationary solutions are localized at the edge in the continuum band. We consider a non-Hermitian Ablowitz-Ladik model and show that nonlinear exceptional point disappears if the lattice is infinitely long., Comment: to appear in pla

We consider a realization of fractional supersymmetric of quantum mechanics of order r, where the Hamiltonian and supercharges involve reflection operators. It is shown that the Hamiltonian has r-fold degenerate spectrum and the eigenvalues of hermitian supercharges are zeros of the associated Hermite polynomials of Askey and Wimp. Also it is shown that the associated eigenfunctions involve lacunary Hermite polynomials.

An important theorem in Gaussian quantum information tells us that we can diagonalise the covariance matrix of any Gaussian state via a symplectic transformation. Whilst the diagonal form is easy to find, the process for finding the diagonalising symplectic can be more difficult, and a common, existing method requires taking matrix powers, which can be demanding analytically. Inspired by a recently presented technique for finding the eigenvectors of a Hermitian matrix from certain submatrix eigenvalues, we derive a similar method for finding the diagonalising symplectic from certain submatrix determinants, which could prove useful in Gaussian quantum information., Comment: 10 pages, supplementary files available at https://github.com/softquanta/symplectic_decomposition

The J\'{a}nossy density for a determinantal point process is the probability density that an interval $I$ contains exactly $p$ points except for those at $k$ designated loci. The J\'{a}nossy density associated with an integrable kernel $\mathbf{K}\doteq (\varphi(x)\psi(y)-\psi(x)\varphi(y))/(x-y)$ is shown to be expressed as a Fredholm determinant $\mathrm{Det}(\mathbb{I}-\tilde{\mathbf{K}}|_I)$ of a transformed kernel $\tilde{\mathbf{K}}\doteq (\tilde{\varphi}(x)\tilde{\psi}(y)-\tilde{\psi}(x)\tilde{\varphi}(y))/(x-y)$. We observe that $\tilde{\mathbf{K}}$ satisfies Tracy and Widom's criteria if $\mathbf{K}$ does, because of the structure that the map $(\varphi, \psi)\mapsto (\tilde{\varphi}, \tilde{\psi})$ is a meromorphic $\mathrm{SL}(2,\mathbb{R})$ gauge transformation between covariantly constant sections. This observation enables application of the Tracy--Widom method to J\'{a}nossy densities, expressed in terms of a solution to a system of differential equations in the endpoints of the interval. Our approach does not explicitly refer to isomonodromic systems associated with Painlev\'{e} equations employed in the preceding works. As illustrative examples we compute J\'{a}nossy densities with $k=1, p=0$ for Airy and Bessel kernels, related to the joint distributions of the two largest eigenvalues of random Hermitian matrices and of the two smallest singular values of random complex matrices., Comment: 18 pages, 8 figs, 2 Mathematica nb's attached. (v2) Additional remarks on SL(r) extension etc, improved agreement with quadrature approx. Version to appear in PTEP

We present a new math-physics modeling approach, called canonical quantization with numerical mode-decomposition, for capturing the physics of how incoming photons interact with finite-sized dispersive media, which is not describable by the previous Fano-diagonalization methods. The main procedure is to (1) study a system where electromagnetic (EM) fields are coupled to non-uniformly distributed Lorentz oscillators in Hamiltonian mechanics, (2) derive a generalized Hermitian eigenvalue problem for conjugate pairs in coordinate space, (3) apply computational electromagnetics methods to find a countably/finite set of time-harmonic eigenmodes that diagonalizes the Hamiltonian, and (4) perform the subsequent canonical quantization with mode-decomposition. Moreover, we provide several numerical simulations that capture the physics of full quantum effects, impossible by classical Maxwell's equations, such as non-local dispersion cancellation of an entangled photon pair and Hong-Ou-Mandel (HOM) effect in a dispersive beam splitter., Comment: submitted to Physical Review A (under review)

Non-Hermitian skin effect, namely that the eigenvalues and eigenstates of a non-Hermitian tight-binding Hamiltonian have significant differences under open or periodic boundary conditions, is a remarkable phenomenon of non-Hermitian systems. Inspired by the presence of the non-Hermitian skin effect, we study the evolution of wave-packets in non-Hermitian systems, which can be determined using the single-particle Green's function. Surprisingly, we find that in the thermodynamical limit, the Green's function does not depend on boundary conditions, despite the presence of skin effect. We proffer a general proof for this statement in arbitrary dimension with finite hopping range, with an explicit illustration in the non-Hermitian Su-Schrieffer-Heeger model. We also explore its applications in non-interacting open quantum systems described by the master equation, where we demonstrate that the evolution of the density matrix is independent of the boundary condition., Comment: 5 pages, 3 figures

The non-Hermitian skin effect, i.e., eigenstate condensation at the edges in lattices with open boundaries, is an exotic manifestation of non-Hermitian systems. In Bloch theory, an effective non-Hermitian Hamiltonian is generally used to describe dissipation, which, however, is not norm preserving and neglects quantum jumps. Here it is shown that in a self-consistent description of the dissipative dynamics in a one-band lattice, based on the stochastic Schrödinger equation or Lindblad master equation with a collective jump operator, the skin effect and its dynamical features are washed out. Nevertheless, both short- and long-time relaxation dynamics provide a hidden signature of the skin effect found in the semiclassical limit. In particular, relaxation toward a maximally mixed state with the largest von Neumann entropy in a lattice with open boundaries is a manifestation of the semiclassical skin effect., The author acknowledges the Spanish State Research Agency through the Severo Ochoa and María de Maeztu Program for Centers and Units of Excellence in R&D (Grant No. MDM-2017-0711).

Standard approach to dynamical random matrix models relies on the description of trajectories of eigenvalues. Using the analogy from optics, based on the duality between the Fermat principle(trajectories) and the Huygens principle (wavefronts), we formulate the Hamilton-Jacobi dynamics for large random matrix models. The resulting equations describe a broad class of random matrix models in a unified way, including normal (Hermitian or unitary) as well as strictly non-normal dynamics. HJ formalism applied to Brownian bridge dynamics allows one for calculations of the asymptotics of the Harish-Chandra-Itzykson-Zuber integrals., Comment: 5 + 9 pages, published version

The spectral statistics of Hermitian random Toeplitz matrices with independent and identically distributed elements are investigated numerically. It is found that eigenvalue statistics of complex Toeplitz matrices are surprisingly well approximated by the semi-Poisson distribution belonging to intermediate-type statistics observed in certain pseudointegrable billiards. The origin of intermediate behavior could be attributed to the fact that Fourier transformed random Toeplitz matrices have the same slow decay outside the main diagonal as critical random matrix ensembles. The statistical properties of the full spectrum of real random Toeplitz matrices are close to the Poisson distribution, but each of their constituent subspectra is again well described by the semi-Poisson distribution. The findings indicate that intermediate statistics in general and the semi-Poisson distribution in particular are more universal than considered before.

The energy band attraction (EBA) caused by the nonorthogonal eigenvectors is a unique phenomenon in the non-Hermitian (NH) system. However, restricted by the required tight-binding approximation and meticulously engineered complex potentials, such an effect has never been experimentally demonstrated before. Here by a suitable design of all-dielectric Mie resonators in a parallel-plate transmission line, we for the first time verify the photonic analog of the EBA effects both theoretically and experimentally. The evolution of the EBA effect in a two-level NH system from gapped bands to gapless bands to flat bands is observed by precisely tuning the loss of the Mie resonators. The transmission spectra can be theoretically connected to the eigenvalues and eigenvectors of the NH Hamiltonian. Furthermore we extend our methods to a graphenelike two-dimensional NH system. Our works show a metamaterial approach toward NH topological photonics and foster a deeper understanding of band theory in open systems.

We propose a scheme for simulating a momentum dependent non-Hermitian Hamiltonian term by adding loss in the hopping process in a one-dimensional array of optical cavities, which leads to a one-band topological nontrivial non-Hermitian model. We find that there is a kind of transition from all real energy spectrum to all imaginary one under open boundary conditions: the eigenvalues are all real or imaginary dependent on which amplitude of the Hermitian or non-Hermitian part is dominant even though the eigenvalues of the bulk system are complex under the periodic boundary condition. In addition, for such a non-Hermitian Hamiltonian, the left eigenvectors are the mirror symmetry of the right eigenvectors. Finally, by using the input-output formula, we show that the distribution of the boundary states can be directly observed by detecting the optical transmission in both cases.

In non-Hermitian (NH) quantum mechanics, Hamiltonians are studied whose eigenvalues are not necessarily real since the condition of hermiticity is not imposed. Certain symmetries of NH operators can ensure that some or all of the eigenvalues are real and thus suitable for the description of physical systems whose energies are always real. While the mathematics of NH quantum mechanics is well developed, applications of the theory to real quantum systems are scarce, and no closed system is known whose Hamiltonian is NH. Here, we consider the elementary textbook example of a NH Hamiltonian matrix, and we show how it naturally emerges as a simplifying concept in the modeling of molecular electronic devices. We analyze the consequences of non-Hermiticity and exceptional points in the spectrum of NH operators for the molecular conductance and the spectral density of simple models for molecules on surfaces.

The direct tracking of a moving radio transmitter using spatially distributed receivers is considered based on delay and Doppler-shift induced signal waveforms. One key step for direct tracking is to approximate the generalized likelihood function corresponding to the largest eigenvalue of a Hermitian matrix. Although stochastic approximations can be set arbitrarily close to the optimal solution, these methods generally have high computational cost and tend to be inefficient when applied to high dimensional problems. In this paper, we develop a Laplace approximation approach for direct tracking using Taylor series expansions. The proposed method belongs to the deterministic approximation domain and therefore it is computational efficient. The simulations validate the performance of the proposed method.

The relevance in Physics of non-Hermitian operators with real eigenvalues is being widely recognized not only in quantum mechanics but also in other areas, such as quantum optics, quantum fluid dynamics and quantum field theory. %stochastic processesand so on. In this note, a quantum system described by a non-Hermitian Hamiltonian, which is constituted by two types of interacting bosons, is investigated. The real eigenvalues of the Hamiltonian are explicitly determined, as well as complete biorthogonal sets of eigenfunctions of the Hamiltonian and its adjoint. The diagonal representation of $H$ is obtained using pseudo-bosonic operators., 16 pages, 1 figure

Schur-like forms are developed for matrices that have a symmetry structure with respect to an indefinite inner product induced by a Hermitian and unitary Gram matrix. It is characterized under which conditions these forms can be computed by structure-preserving unitary transformations. The main results combines and generalizes the two well-known results from the literature that on the one hand any normal matrix can be unitarily diagonalized and on the other hand a Hamiltonian matrix can be transformed to Hamiltonian Schur form via a unitary similarity transformation if and only if its purely imaginary eigenvalues satisfy certain conditions that involve the sign characteristic of the matrix under consideration.

In the hard edge scaling limit of the Jacobi unitary ensemble generated by the weight $x^{\alpha}(1-x)^{\beta},~x\in[0,1],~\alpha,\beta>0$, the probability that all eigenvalues of Hermitian matrices from this ensemble lie in the interval $[t,1]$ is given by the Fredholm determinant of the Bessel kernel. We derive the constant in the asymptotics of this Bessel-kernel determinant. A specialization of the results gives the constant in the asymptotics of the probability that the interval $(-a,a),a>0,$ is free of eigenvalues in the Jacobi unitary ensemble with the symmetric weight $(1-x^2)^{\beta}, x\in[-1,1]$., Comment: 19 pages

Unlike their fermionic counterparts, the dynamics of Hermitian quadratic bosonic Hamiltonians are governed by a generally non-Hermitian Bogoliubov-de Gennes effective Hamiltonian. This underlying non-Hermiticity gives rise to a dynamically stable regime, whereby all observables undergo bounded evolution in time, and a dynamically unstable one, whereby evolution is unbounded for at least some observables. We show that stability-to-instability transitions may be classified in terms of a suitably generalized $\mathcal{P}\mathcal{T}$ symmetry, which can be broken when diagonalizability is lost at exceptional points in parameter space, but also when degenerate real eigenvalues split off the real axis while the system remains diagonalizable. By leveraging tools from Krein stability theory in indefinite inner-product spaces, we introduce an indicator of stability phase transitions, which naturally extends the notion of phase rigidity from non-Hermitian quantum mechanics to the bosonic setting. As a paradigmatic example, we fully characterize the stability phase diagram of a bosonic analogue to the Kitaev-Majorana chain under a wide class of boundary conditions. In particular, we establish a connection between phase-dependent transport properties and the onset of instability, and argue that stable regions in parameter space become of measure zero in the thermodynamic limit. Our analysis also reveals that boundary conditions that support Majorana zero modes in the fermionic Kitaev chain are precisely the same that support stability in the bosonic chain., 30+7 pages, 11 figures, uses iopart.cls

Absence of level repulsion between extended states in random non-Hermitian systems is demonstrated. As a result, the general Wigner-Dyson distributions of level spacing of diffusive metals in the usual Hermitian systems is replaced by the Poisson distribution for quasiparticle level spacing of non-Hermitian disordered metals in the thermodynamic limit of infinite system size. This is a very surprising result because Poisson statistics is universally true for the Anderson insulators where energy eigenstates do not overlap with each other so that energy levels are independent from each other.For disordered metals where different eigenstates overlap with each other, one should expect different levels trying to stay away from each other so that the Poisson distribution should not apply there. Our results show that the larger non-Hermitian energy (dissipation) can invalidate level repulsion principle that holds dearly in quantum mechanics. Thus, our theory provides a unified picture for recent discovery of so called "level attraction" in various systems. It provides also a theoretical basis for manipulating energy levels., 10 pages, 8 figures