Non-Hermitian physics for optical manipulation uncovers inherent instability of large clusters

Intense light traps and binds small particles, offering unique control to the microscopic world. With incoming illumination and radiative losses, optical forces are inherently nonconservative, thus non-Hermitian. Contrary to conventional systems, the operator governing time evolution is real and asymmetric (i.e., non-Hermitian), which inevitably yield complex eigenvalues when driven beyond the exceptional points, where light pumps in energy that eventually “melts” the light-bound structures. Surprisingly, unstable complex eigenvalues are prevalent for clusters with ~10 or more particles, and in the many-particle limit, their presence is inevitable. As such, optical forces alone fail to bind a large cluster. Our conclusion does not contradict with the observation of large optically-bound cluster in a fluid, where the ambient damping can take away the excess energy and restore the stability. The non-Hermitian theory overturns the understanding of optical trapping and binding, and unveils the critical role played by non-Hermiticity and exceptional points, paving the way for large-scale manipulation.

Spectra, eigenstates and transport properties of a $\mathcal{PT}$-symmetric ring

We study, analytically and numerically, a simple $\mathcal{PT}$-symmetric tight-binding ring with an onsite energy $a$ at the gain and loss sites. We show that if $a\neq 0$, the system generically exhibits an unbroken PT -symmetric phase. We study the nature of the spectrum in terms of the singularities in the complex parameter space as well as the behavior of the eigenstates at large values of the gain and loss strength. We find that in addition to the usual exceptional points, there are “diabolical points”, and inverse exceptional points at which complex eigenvalues reconvert into real eigenvalues. We also study the transport through the system. We calculate the total flux from the source to the drain, and how it splits along the branches of the ring. We find that while usually the density flows from the source to the drain, for certain eigenstates a stationary “backflow” of density from the drain to the source along one of the branches can occur. We also identify two types of singular eigenstates, i.e. states that do not depend on the strength of the gain and loss, and classify them in terms of their transport properties.

Exceptional topological insulators

We introduce the exceptional topological insulator (ETI), a non-Hermitian topological state of matter that features exotic non-Hermitian surface states which can only exist within the three-dimensional topological bulk embedding. We show how this phase can evolve from a Weyl semimetal or Hermitian three-dimensional topological insulator close to criticality when quasiparticles acquire a finite lifetime. The ETI does not require any symmetry to be stabilized. It is characterized by a bulk energy point gap, and exhibits robust surface states that cover the bulk gap as a single sheet of complex eigenvalues or with a single exceptional point. The ETI can be induced universally in gapless solid-state systems, thereby setting a paradigm for non-Hermitian topological matter.

Generation of test matrices with specified eigenvalues using floating-point arithmetic

This paper concerns test matrices for numerical linear algebra using an error-free transformation of floating-point arithmetic. For specified eigenvalues given by a user, we propose methods of generating a matrix whose eigenvalues are exactly known based on, for example, Schur or Jordan normal form and a block diagonal form. It is also possible to produce a real matrix with specified complex eigenvalues. Such test matrices with exactly known eigenvalues are useful for numerical algorithms in checking the accuracy of computed results. In particular, exact errors of eigenvalues can be monitored. To generate test matrices, we first propose an error-free transformation for the product of three matrices YSX. We approximate S by ${S^{\prime }}$ S ′ to compute ${YS^{\prime }X}$ Y S ′ X without a rounding error. Next, the error-free transformation is applied to the generation of test matrices with exactly known eigenvalues. Note that the exactly known eigenvalues of the constructed matrix may differ from the anticipated given eigenvalues. Finally, numerical examples are introduced in checking the accuracy of numerical computations for symmetric and unsymmetric eigenvalue problems.

Concentrated matrix exponential distributions with real eigenvalues

Concentrated random variables are frequently used in representing deterministic delays in stochastic models. The squared coefficient of variation ( $\mathrm {SCV}$ ) of the most concentrated phase-type distribution of order $N$ is $1/N$ . To further reduce the $\mathrm {SCV}$ , concentrated matrix exponential (CME) distributions with complex eigenvalues were investigated recently. It was obtained that the $\mathrm {SCV}$ of an order $N$ CME distribution can be less than $n^{-2.1}$ for odd $N=2n+1$ orders, and the matrix exponential distribution, which exhibits such a low $\mathrm {SCV}$ has complex eigenvalues. In this paper, we consider CME distributions with real eigenvalues (CME-R). We present efficient numerical methods for identifying a CME-R distribution with smallest SCV for a given order $n$ . Our investigations show that the $\mathrm {SCV}$ of the most concentrated CME-R of order $N=2n+1$ is less than $n^{-1.85}$ . We also discuss how CME-R can be used for numerical inverse Laplace transformation, which is beneficial when the Laplace transform function is impossible to evaluate at complex points.

On eigenvalues of a matrix arising in energy-preserving/dissipative continuous-stage Runge-Kutta methods

In this short note, we define an s × s matrix Ks constructed from the Hilbert matrix H s = ( 1 i + j - 1 ) i , j = 1 s {H_s} = \left( {{1 \over {i + j - 1}}} \right)_{i,j = 1}^s and prove that it has at least one pair of complex eigenvalues when s ≥ 2. Ks is a matrix related to the AVF collocation method, which is an energy-preserving/dissipative numerical method for ordinary differential equations, and our result gives a matrix-theoretical proof that the method does not have large-grain parallelism when its order is larger than or equal to 4.

The Impact of Small Time Delays on the Onset of Oscillations and Synchrony in Brain Networks

The human brain constitutes one of the most advanced networks produced by nature, consisting of billions of neurons communicating with each other. However, this communication is not in real-time, with different communication or time-delays occurring between neurons in different brain areas. Here, we investigate the impacts of these delays by modeling large interacting neural circuits as neural-field systems which model the bulk activity of populations of neurons. By using a Master Stability Function analysis combined with numerical simulations, we find that delays (1) may actually stabilize brain dynamics by temporarily preventing the onset to oscillatory and pathologically synchronized dynamics and (2) may enhance or diminish synchronization depending on the underlying eigenvalue spectrum of the connectivity matrix. Real eigenvalues with large magnitudes result in increased synchronizability while complex eigenvalues with large magnitudes and positive real parts yield a decrease in synchronizability in the delay vs. instantaneously coupled case. This result applies to networks with fixed, constant delays, and was robust to networks with heterogeneous delays. In the case of real brain networks, where the eigenvalues are predominantly real, owing to the nearly symmetric nature of these weight matrices, biologically plausible, small delays, are likely to increase synchronization, rather than decreasing it.

Algebraic treatment of the Bateman Hamiltonian

We apply the algebraic method to the Bateman Hamiltonian and obtain its natural frequencies and ladder operators from the adjoint or regular matrix representation of that operator. Present analysis shows that the eigenfunctions compatible with the complex eigenvalues obtained earlier by other authors are not square integrable. In addition to this, the ladder operators annihilate an infinite number of such eigenfunctions.