Source author record

Joshua Feinberg

Joshua Feinberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Pseudo-Hermitian Random Matrix Models: General Formalism

Pseudo-hermitian matrices are matrices hermitian with respect to an indefinite metric. They can be thought of as the truncation of pseudo-hermitian operators, defined over some Krein space, together with the associated metric, to a finite dimensional subspace. As such, they can be used, in the usual spirit of random matrix theory, to model chaotic or disordered PT-symmetric quantum systems, or their gain-loss-balanced classical analogs, in the phase of broken PT-symmetry. The eigenvalues of pseudo-hermitian matrices are either real, or come in complex-conjugate pairs. In this paper we introduce a family of pseudo-hermitian random matrix models, depending parametrically on their metric. We apply the diagrammatic method to obtain its averaged resolvent and density of eigenvalues as explicit functions of the metric, in the limit of large matrix size $N$. As a concrete example, which is essentially an ensemble of elements of the non-compact unitary Lie algebra, we choose a particularly simple set of metrics, and compute the resulting resolvent and density of eigenvalues in closed form. The spectrum consists of a finite fraction of complex eigenvalues, which occupy uniformly two two-dimensional blobs, symmetric with respect to the real axis, as well as the complimentary fraction of real eigenvalues, condensed in a finite segment, with a known non-uniform density. The numbers of complex and real eigenvalues depend on the signature of the metric, that is, the numbers of its positive and negative eigenvalues. We have also carried thorough numerical analysis of the model for these particular metrics. Our numerical results converge rapidly towards the asymptotic analytical large-$N$ expressions.

preprint2021arXiv

Dynamics of Disordered Mechanical Systems with Large Connectivity, Free Probability Theory, and Quasi-Hermitian Random Matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. We calculate the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. In addition to matrix model analysis, we study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

preprint2021arXiv

Which Metrics Are Consistent with a Given Pseudo-Hermitian Matrix?

Given a diagonalizable $N\times N$ matrix $H$, whose non-degenerate spectrum consists of $p$ pairs of complex conjugate eigenvalues and additional $N-2p$ real eigenvalues, we determine all metrics $M$, of all possible signatures, with respect to which $H$ is pseudo-hermitian. In particular, we show that any compatible $M$ must have $p$ pairs of opposite eigenvalues in its spectrum so that $p$ is the minimal number of both positive and negative eigenvalues of $M$. We provide explicit parametrization of the space of all admissible metrics and show that it is topologically a $p$-dimensional torus tensored with an appropriate power of the group $Z_2$.

preprint2016arXiv

Chebyshev-polynomial expansion of the localization length of Hermitian and non-Hermitian random chains

We carry Chebyshev-polynomial expansion of the inverse localization length of Hermitian and non-Hermitian random chains as function of energy. For Hermitian models, the expansion produces numerically this energy-dependent function in one run of the algorithm. This is in strong contrast to the standard transfer-matrix method, which produces the inverse localization length for a fixed energy in each run. For non-Hermitian models, as in the transfer-matrix method, our algorithm computes the inverse localization length for a fixed (complex) energy. We also find a formula of the Chebyshev-polynomial expansion of the density of states of non-Hermitian models. As explained in more detail in the Introduction, our algorithm for non-Hermitian models may be the only available efficient algorithm for finding the density of states of models with interactions.

preprint2009arXiv

Probabilistic Interpretation of Resonant States

We provide probabilistic interpretation of resonant states. This we do by showing that the integral of the modulus square of resonance wave functions (i.e., the conventional norm) over a properly expanding spatial domain is independent of time, and therefore leads to probability conservation. This is in contrast with the conventional employment of a bi-orthogonal basis that precludes probabilistic interpretation, since wave functions of resonant states diverge exponentially in space. On the other hand, resonant states decay exponentially in time, because momentum leaks out of the central scattering area. This momentum leakage is also the reason for the spatial exponential divergence of resonant state. It is by combining the opposite temporal and spatial behaviors of resonant states that we arrive at our probabilistic interpretation of these states. The physical need to normalize resonant wave functions over an expanding spatial domain arises because particles leak out of the region which contains the potential range and escape to infinity, and one has to include them in the total count of particle number.

preprint2009arXiv

Statistics of Resonances in One Dimensional Continuous Systems

We study the average density of resonances (DOR) of a disordered one-dimensional continuous open system. The disordered system is semi-infinite, with white-noise random potential, and it is coupled to the external world by a semi-infinite continuous perfect lead. Our main result is an integral representation for the DOR which involves the probability density function of the logarithmic derivative of the wave function at the contact point.

preprint2007arXiv

Does the complex deformation of the Riemann equation exhibit shocks?

The Riemann equation $u_t+uu_x=0$, which describes a one-dimensional accelerationless perfect fluid, possesses solutions that typically develop shocks in a finite time. This equation is $\cP\cT$ symmetric. A one-parameter $\cP\cT$-invariant complex deformation of this equation, $u_t-iu(iu_x)^ε= 0$ ($ε$ real), is solved exactly using the method of characteristic strips, and it is shown that for real initial conditions, shocks cannot develop unless $ε$ is an odd integer.

preprint2006arXiv

Non-Hermitean Random Matrix Theory: Summation of Planar Diagrams, the "Single-Ring" Theorem and the Disk-Annulus Phase Transition

I review aspects of work done in collaboration with A. Zee and R. Scalettar \cite{fz1,fz2,fsz} on complex non-hermitean random matrices. I open by explaining why the bag of tools used regularly in analyzing hermitean random matrices cannot be applied directly to analyze non-hermitean matrices, and then introduce the Method of Hermitization, which solves this problem. Then, for rotationally invariant ensembles, I derive a master equation for the average density of eigenvalues in the complex plane, in the limit of infinitely large matrices. This is achieved by resumming all the planar diagrams which appear in the perturbative expansion of the hermitized Green function. Remarkably, this resummation can be carried {\em explicitly} for any rotationally invariant ensemble. I prove that in the limit of infinitely large matrices, the shape of the eigenvalue distribution is either a disk or an annulus. This is the celebrated ``Single-Ring'' Theorem. Which of these shapes is realized is determined by the parameters (coupling constants) which determine the ensemble. By varying these parameters a phase transition may occur between the two possible shapes. I briefly discuss the universal features of this transition. As the analysis of this problem relies heavily on summation of planar Feynman diagrams, I take special effort at presenting a pedagogical exposition of the diagrammatic method, which some readers may find useful.

preprint2003arXiv

All about the Static Fermion Bags in the Gross-Neveu Model

We review in detail the construction of {\em all} stable static fermion bags in the 1+1 dimensional Gross-Neveu model with $N$ flavors of Dirac fermions, in the large $N$ limit. In addition to the well known kink and topologically trivial solitons (which correspond, respectively, to the spinor and antisymmetric tensor representations of O(2N)), there are also threshold bound states of a kink and a topologically trivial soliton: the heavier topological solitons (HTS). The mass of any of these newly discovered HTS's is the sum of masses of its solitonic constituents, and it corresponds to the tensor product of their O(2N) representations. Thus, it is marginally stable (at least in the large $N$ limit). Furthermore, its mass is independent of the distance between the centers of its constituents, which serves as a flat collective coordinate, or a modulus. There are no additional stable static solitons in the Gross-Neveu model. We provide detailed derivation of the profiles, masses and fermion number contents of these static solitons. For pedagogical clarity, and in order for this paper to be self-contained, we also included detailed appendices on supersymmetric quantum mechanics and on reflectionless potentials in one spatial dimension, which are intimately related with the theory of static fermion bags. In particular, we present a novel simple explicit formula for the diagonal resolvent of a reflectionless Schrödinger operator with an arbitrary number of bound states. In additional appendices we summarize the relevant group representation theoretic facts, and also provide a simple calculation of the mass of the kinks.

preprint2001arXiv

"Single Ring Theorem" and the Disk-Annulus Phase Transition

Recently, an analytic method was developed to study in the large $N$ limit non-hermitean random matrices that are drawn from a large class of circularly symmetric non-Gaussian probability distributions, thus extending the existing Gaussian non-hermitean literature. One obtains an explicit algebraic equation for the integrated density of eigenvalues from which the Green's function and averaged density of eigenvalues could be calculated in a simple manner. Thus, that formalism may be thought of as the non-hermitean analog of the method due to Brézin, Itzykson, Parisi and Zuber for analyzing hermitean non-Gaussian random matrices. A somewhat surprising result is the so called "Single Ring" theorem, namely, that the domain of the eigenvalue distribution in the complex plane is either a disk or an annulus. In this paper we extend previous results and provide simple new explicit expressions for the radii of the eigenvalue distiobution and for the value of the eigenvalue density at the edges of the eigenvalue distribution of the non-hermitean matrix in terms of moments of the eigenvalue distribution of the associated hermitean matrix. We then present several numerical verifications of the previously obtained analytic results for the quartic ensemble and its phase transition from a disk shaped eigenvalue distribution to an annular distribution. Finally, we demonstrate numerically the "Single Ring" theorem for the sextic potential, namely, the potential of lowest degree for which the "Single Ring" theorem has non-trivial consequences.

preprint2001arXiv

Generalized Supersymmetric Quantum Mechanics and Reflectionless Fermion Bags in 1+1 Dimensions

We study static fermion bags in the 1+1 dimensional Gross-Neveu and Nambu-Jona-Lasinio models. It has been known, from the work of Dashen, Hasslacher and Neveu (DHN), followed by Shei's work, in the 1970's, that the self-consistent static fermion bags in these models are reflectionless. The works of DHN and of Shei were based on inverse scattering theory. Several years ago, we offered an alternative argument to establish the reflectionless nature of these fermion bags, which was based on analysis of the spatial asymptotic behavior of the resolvent of the Dirac operator in the background of a static bag, subjected to the appropriate boundary conditions. We also calculated the masses of fermion bags based on the resolvent and the Gelfand-Dikii identity. Based on arguments taken from a certain generalized one dimensional supersymmetric quantum mechanics, which underlies the spectral theory of these Dirac operators, we now realize that our analysis of the asymptotic behavior of the resolvent was incomplete. We offer here a critique of our asymptotic argument.

preprint1997arXiv

Self-Isospectral Periodic Potentials and Supersymmetric Quantum Mechanics

We discuss supersymmetric quantum mechanical models with periodic potentials. The important new feature is that it is possible for both isospectral potentials to support zero modes, in contrast to the standard nonperiodic case where either one or neither (but not both) of the isospectral pair has a zero mode. Thus it is possible to have supersymmetry unbroken and yet also have a vanishing Witten index. We present some explicit exactly soluble examples for which the isospectral potentials have identical band spectra, and which are ``self-isospectral'' in the sense that the potentials have identical shape, but are translated by one half period relative to one another.

preprint1994arXiv

On Kinks and Bound States in the Gross-Neveu Model

We investigate static space dependent $\sigx=\lag\barψψ\rag$ saddle point configurations in the two dimensional Gross-Neveu model in the large N limit. We solve the saddle point condition for $\sigx$ explicitly by employing supersymmetric quantum mechanics and using simple properties of the diagonal resolvent of one dimensional Schrödinger operators rather than inverse scattering techniques. The resulting solutions in the sector of unbroken supersymmetry are the Callan-Coleman-Gross-Zee kink configurations. We thus provide a direct and clean construction of these kinks. In the sector of broken supersymmetry we derive the DHN saddle point configurations. Our method of finding such non-trivial static configurations may be applied also in other two dimensional field theories.

preprint1993arXiv

Dissipation of the String Embedding Dimension in the Singlet Sector of a Matrix Model at Large $α^\prime$

The one dimensional Fermi gas of matrix eigenvalues of the Marinari-Parisi model at positive values of the cosmological constant is generalised.The number of matrix eigenvalues (i.e. gas particles) is varied while keeping the effective potential fixed. This model exhibits a transition from a phase whose continuum behaviour is that of $c=1$ conformal matter coupled to gravity to the well known pure gravity phase of the original model.The former phase is character- ised by an extremely large Regge slope $α^\prime$ which scales as $β^{2/5}$ causing the scaling regions of the two phases to overlap. In this way a continuous flow from one phase to the other is made poss- ible. This phase transition occurs in the singlet sector of the matrix model. The density of states and the two puncture correlator at non zero momenta are calculated on the sphere and are found to behave very differ- ently in the two phases,a fact which demonstrates the phase transition. We comment on a possible relation between this transition and large α^\prime$ semiclassical expansions in the continuum.